44#include " source_base/parallel_common.h"
55#include " source_base/timer.h"
66#include " source_base/tool_check.h"
7- #include " source_io/module_parameter/parameter.h"
87
98#include < algorithm>
109#include < cassert>
@@ -58,6 +57,7 @@ void LCAO_Orbitals::init(
5857 const double & lcao_rmax_in,
5958 const bool & deepks_setorb,
6059 const int & out_mat_r,
60+ const bool & out_element_info,
6161 const bool & force_flag,
6262 const int & my_rank
6363)
@@ -86,7 +86,7 @@ void LCAO_Orbitals::init(
8686#ifdef __MPI
8787 bcast_files (ntype, my_rank);
8888#endif
89- Read_Orbitals (ofs_in, ntype, lmax, deepks_setorb, out_mat_r, force_flag, my_rank);
89+ Read_Orbitals (ofs_in, ntype, lmax, deepks_setorb, out_mat_r, out_element_info, force_flag, my_rank);
9090 return ;
9191}
9292
@@ -156,6 +156,7 @@ void LCAO_Orbitals::Read_Orbitals(std::ofstream& ofs_in,
156156 const int & lmax_in,
157157 const bool & deepks_setorb,
158158 const int & out_mat_r,
159+ const bool & out_element_info,
159160 const bool & force_flag, // mohan add 2021-05-07
160161 const int & my_rank) // mohan add 2021-04-26
161162{
@@ -247,7 +248,7 @@ void LCAO_Orbitals::Read_Orbitals(std::ofstream& ofs_in,
247248 this ->Phi = new Numerical_Orbital[ntype];
248249 for (int it = 0 ; it < ntype; it++)
249250 {
250- this ->Read_PAO (ofs_in, it, force_flag, my_rank);
251+ this ->Read_PAO (ofs_in, it, out_element_info, force_flag, my_rank);
251252 // caoyu add 2021-05-24 to reconstruct atom_arrange::set_sr_NL
252253 this ->rcutmax_Phi = std::max (this ->rcutmax_Phi , this ->Phi [it].getRcut ());
253254 }
@@ -263,7 +264,7 @@ void LCAO_Orbitals::Read_Orbitals(std::ofstream& ofs_in,
263264
264265 delete[] this ->Alpha ;
265266 this ->Alpha = new Numerical_Orbital[1 ]; // not related to atom type -- remain to be discussed
266- this ->Read_Descriptor (ofs_in, force_flag, my_rank);
267+ this ->Read_Descriptor (ofs_in, out_element_info, force_flag, my_rank);
267268 }
268269
269270 ModuleBase::timer::end (" LCAO_Orbitals" , " Read_Orbitals" );
@@ -279,6 +280,7 @@ void LCAO_Orbitals::Read_Orbitals(std::ofstream& ofs_in,
279280// -------------------------------------------------------
280281void LCAO_Orbitals::Read_PAO (std::ofstream& ofs_in,
281282 const int & it,
283+ const bool & out_element_info,
282284 const bool & force_flag, // mohan add 2021-05-07
283285 const int & my_rank) // mohan add 2021-04-26
284286{
@@ -311,7 +313,7 @@ void LCAO_Orbitals::Read_PAO(std::ofstream& ofs_in,
311313 int lmaxt = 0 ;
312314 int nchimaxt = 0 ;
313315
314- this ->read_orb_file (ofs_in, in_ao, it, lmaxt, nchimaxt, this ->Phi , force_flag, my_rank);
316+ this ->read_orb_file (ofs_in, in_ao, it, lmaxt, nchimaxt, this ->Phi , out_element_info, force_flag, my_rank);
315317
316318 // lmax and nchimax for all types
317319 this ->lmax = std::max (this ->lmax , lmaxt);
@@ -323,6 +325,7 @@ void LCAO_Orbitals::Read_PAO(std::ofstream& ofs_in,
323325
324326// caoyu add 2021-3-16
325327void LCAO_Orbitals::Read_Descriptor (std::ofstream& ofs_in,
328+ const bool & out_element_info,
326329 const bool & force_flag, // mohan add 2021-05-07
327330 const int & my_rank) // read descriptor basis
328331{
@@ -355,7 +358,7 @@ void LCAO_Orbitals::Read_Descriptor(std::ofstream& ofs_in,
355358 this ->lmax_d = 0 ;
356359 this ->nchimax_d = 0 ;
357360
358- this ->read_orb_file (ofs_in, in_de, 0 , this ->lmax_d , this ->nchimax_d , this ->Alpha , force_flag, my_rank);
361+ this ->read_orb_file (ofs_in, in_de, 0 , this ->lmax_d , this ->nchimax_d , this ->Alpha , out_element_info, force_flag, my_rank);
359362
360363 in_de.close ();
361364
@@ -368,6 +371,7 @@ void LCAO_Orbitals::read_orb_file(std::ofstream& ofs_in, // GlobalV::ofs_running
368371 int & lmax,
369372 int & nchimax,
370373 Numerical_Orbital* ao,
374+ const bool & out_element_info,
371375 const bool & force_flag,
372376 const int & my_rank)
373377{
@@ -588,7 +592,7 @@ void LCAO_Orbitals::read_orb_file(std::ofstream& ofs_in, // GlobalV::ofs_running
588592 this ->kmesh ,
589593 this ->dk ,
590594 this ->dr_uniform ,
591- PARAM . inp . out_element_info ,
595+ out_element_info,
592596 true ,
593597 force_flag); // delta k mesh in reciprocal space
594598
0 commit comments