@@ -239,7 +239,7 @@ void ESolver_KS_PW<T, Device>::before_scf(const int istep)
239239 {
240240 std::stringstream ss;
241241 ss << PARAM.globalv .global_out_dir << " SPIN" << is + 1 << " _CHG_INI.cube" ;
242- ModuleIO::write_cube (GlobalC::Pgrid,
242+ ModuleIO::write_grid (GlobalC::Pgrid,
243243 this ->pelec ->charge ->rho [is],
244244 is,
245245 PARAM.inp .nspin ,
@@ -257,7 +257,7 @@ void ESolver_KS_PW<T, Device>::before_scf(const int istep)
257257 {
258258 std::stringstream ss;
259259 ss << PARAM.globalv .global_out_dir << " SPIN" << is + 1 << " _POT_INI.cube" ;
260- ModuleIO::write_cube (GlobalC::Pgrid,
260+ ModuleIO::write_grid (GlobalC::Pgrid,
261261 this ->pelec ->pot ->get_effective_v (is),
262262 is,
263263 PARAM.inp .nspin ,
@@ -462,7 +462,7 @@ void ESolver_KS_PW<T, Device>::iter_finish(int& iter)
462462 data = this ->pelec ->charge ->rho_save [is];
463463 }
464464 std::string fn = PARAM.globalv .global_out_dir + " /tmp_SPIN" + std::to_string (is + 1 ) + " _CHG.cube" ;
465- ModuleIO::write_cube (GlobalC::Pgrid,
465+ ModuleIO::write_grid (GlobalC::Pgrid,
466466 data,
467467 is,
468468 PARAM.inp .nspin ,
@@ -475,7 +475,7 @@ void ESolver_KS_PW<T, Device>::iter_finish(int& iter)
475475 if (XC_Functional::get_func_type () == 3 || XC_Functional::get_func_type () == 5 )
476476 {
477477 fn = PARAM.globalv .global_out_dir + " /tmp_SPIN" + std::to_string (is + 1 ) + " _TAU.cube" ;
478- ModuleIO::write_cube (GlobalC::Pgrid,
478+ ModuleIO::write_grid (GlobalC::Pgrid,
479479 this ->pelec ->charge ->kin_r_save [is],
480480 is,
481481 PARAM.inp .nspin ,
0 commit comments