Control{
    cpumax = 25 min
    condition = fixed_charge 
}
accuracy{
    ksampling{
        method = file 
    }
    xctype = ggapbe 
    cutoff_wf = 25.0 Rydberg
    cutoff_cd = 225.0 Rydberg
    num_bands = 10 
    force_convergence{
        max_force = 0.1e-3 Hartree/Bohr
    }
    scf_convergence{
        delta_total_energy = 1e-8 Hartree
    }
    initial_charge_density = file 
    hybrid_functional{
        sw_hybrid_functional = on
        sw_screened_exchange = on
        functional_type = hse06 
        sw_change_axis = on 
        alpha = 0.25 
        omega = 0.106 
        sw_rspace = on 
    }
    initial_wavefunctions = file
    ek_convergence{
    !    num_max_iteration = 500 
        sw_eval_eig_diff = on 
        delta_eigenvalue = 1.e-5 
        succession = 2 
        num_extra_bands = 5 
    }
}
structure{
    element_list{
        #tag element atomicnumber mass zeta
         Ga 31 127096.476 0.3
         N 7 25532.4958 0.0
    }
    atom_list{
        coordinate_system = internal 
        atoms{
            #units angstrom
            #tag element rx ry rz
             N 0.33333333337736604 0.6666666667436242 0.002132533348762255
             N 0.6666666666655773 0.33333333327561526 0.502132533323024
             Ga 0.33333333337736604 0.6666666667436242 0.37924698407950147
             Ga 0.6666666666655773 0.33333333327561526 0.8792469840537632
        }
    }
    unit_cell{
        #units angstrom
        a_vector = 3.175994091150693 0.0 0.0 
        b_vector = -1.5879970456018053 2.750491565226408 0.0 
        c_vector = 0.0 0.0 5.139990392653881 
    }
    symmetry{
        tspace{
            lattice_system = hexagonal 
        }
        method = automatic 
    }
    unit_cell_type = primitive 
    magnetic_state = ferro 
}
!structure_evolution{
!    method = cg 
!}
wavefunction_solver{
    solvers{
        #tag sol till_n submat
         lmMSD 1 on
         RMM3 -1 on
    }
    rmm{
        edelta_change_to_rmm = 1.0e-3 
    }
    submat{
        before_renewal = on 
    }
}
charge_mixing{
    mixing_methods{
        #tag no method rmxs rmxe istr nbmix update
         1 pulay 0.10 0.10 3 15 renew
    }
}
postprocessing{
   sw_band_symmetry_analysis = on
}