@@ -10,7 +10,7 @@ void write_elf(
1010 const int & nbz,
1111#endif
1212 const std::string& out_dir,
13- const int & istep ,
13+ const int & istep_in ,
1414 const int & nspin,
1515 const double * const * rho,
1616 const double * const * tau,
@@ -94,14 +94,14 @@ void write_elf(
9494
9595 if (nspin == 1 )
9696 {
97- std::string fn = out_dir + " /ELF .cube" ;
97+ std::string fn = out_dir + " /elf .cube" ;
9898
9999 int is = -1 ;
100100 ModuleIO::write_vdata_palgrid (pgrid,
101101 elf[0 ].data (),
102102 is,
103103 nspin,
104- istep ,
104+ istep_in ,
105105 fn,
106106 ef_tmp,
107107 ucell_,
@@ -112,14 +112,17 @@ void write_elf(
112112 {
113113 for (int is = 0 ; is < nspin; ++is)
114114 {
115- std::string fn_temp = out_dir + " /ELF_SPIN" + std::to_string (is + 1 ) + " .cube" ;
115+ std::string fn_temp = out_dir + " /elf" ;
116+
117+ fn_temp += fn_temp + std::to_string (is + 1 ) + " .cube" ;
118+
116119 int ispin = is + 1 ;
117120
118121 ModuleIO::write_vdata_palgrid (pgrid,
119122 elf[is].data (),
120123 ispin,
121124 nspin,
122- istep ,
125+ istep_in ,
123126 fn_temp,
124127 ef_tmp,
125128 ucell_,
@@ -133,19 +136,19 @@ void write_elf(
133136 elf_tot[ir] = (tau[0 ][ir] + tau[1 ][ir] - tau_vw[0 ][ir] - tau_vw[1 ][ir]) / (tau_TF[0 ][ir] + tau_TF[1 ][ir]);
134137 elf_tot[ir] = 1 . / (1 . + elf_tot[ir] * elf_tot[ir]);
135138 }
136- std::string fn = out_dir + " /ELF .cube" ;
139+ std::string fn = out_dir + " /elf .cube" ;
137140
138141 int is = -1 ;
139142 ModuleIO::write_vdata_palgrid (pgrid,
140143 elf_tot.data (),
141144 is,
142145 nspin,
143- istep ,
146+ istep_in ,
144147 fn,
145148 ef_tmp,
146149 ucell_,
147150 precision,
148151 out_fermi);
149152 }
150153}
151- }
154+ }
0 commit comments