Refactor: remove dead cross module includes

source/source_cell/read_atoms.cpp

@@ -8,16 +8,11 @@
 #include "source_io/module_parameter/parameter.h"
 #include "print_cell.h"
 #include "read_stru.h"
-#include "source_estate/read_orb.h"
 #include "source_base/timer.h"
 #include "source_base/constants.h"
 #include "source_base/formatter.h"
 #include "source_base/mathzone.h"
 
-#ifdef __LCAO
-#include "source_basis/module_ao/ORB_read.h" // to use 'ORB' -- mohan 2021-01-30
-#endif
-
 bool unitcell::read_atom_positions(UnitCell& ucell,
                          std::ifstream &ifpos,
                          std::ofstream &ofs_running,

source/source_esolver/esolver.cpp

@@ -19,7 +19,6 @@
 #include "esolver_lj.h"
 #include "esolver_of.h"
 #include "esolver_of_tddft.h"
-#include "source_io/module_parameter/md_parameter.h"
 
 #include <stdexcept>
 

source/source_esolver/esolver_dm2rho.cpp

@@ -9,7 +9,6 @@
 #include "source_lcao/module_operator_lcao/operator_lcao.h"
 #include "source_io/module_output/cube_io.h"
 #include "source_io/module_ml/io_npz.h"
-#include "source_io/module_output/print_info.h"
 #include "source_lcao/rho_tau_lcao.h" // mohan add 2025-10-24
 
 namespace ModuleESolver

source/source_esolver/esolver_double_xc.cpp

@@ -15,7 +15,6 @@
 #include "source_lcao/hamilt_lcao.h"
 #include "source_hsolver/hsolver_lcao.h"
 #include "source_io/module_parameter/parameter.h"
-#include "source_io/module_hs/write_HS.h" // use ModuleIO::write_hsk()
 #include "source_lcao/setup_deepks.h" // use deepks, mohan add 2025-10-10
 
 namespace ModuleESolver

source/source_esolver/esolver_fp.cpp

@@ -6,9 +6,6 @@
 #include "source_hamilt/module_ewald/H_Ewald_pw.h"
 #include "source_hamilt/module_vdw/vdw.h"
 #include "source_io/module_output/cif_io.h"
-#include "source_io/module_output/cube_io.h" // use write_vdata_palgrid
-#include "source_io/module_json/init_info.h"
-#include "source_io/module_json/output_info.h"
 #include "source_io/module_output/output_log.h"
 #include "source_io/module_output/print_info.h"
 #include "source_io/module_chgpot/rhog_io.h"

source/source_esolver/esolver_ks.cpp

@@ -2,7 +2,6 @@
 #include "source_base/timer_wrapper.h"
 
 // for jason output information
-#include "source_io/module_json/init_info.h"
 #include "source_io/module_json/output_info.h"
 
 #include "source_estate/update_pot.h" // mohan add 20251016
@@ -14,7 +13,6 @@
 #include "source_hamilt/module_xc/xc_functional.h"
 #include "source_io/module_output/output_log.h" // use write_head
 #include "source_estate/elecstate_print.h" // print_etot
-#include "source_io/module_output/print_info.h" // print_parameters
 #include "source_lcao/module_dftu/dftu.h" // mohan add 2025-11-07
 
 namespace ModuleESolver

source/source_esolver/esolver_ks_lcao_tddft.cpp

@@ -4,8 +4,6 @@
 //----------------IO-----------------
 #include "source_base/global_variable.h"
 #include "source_io/module_ctrl/ctrl_output_td.h"
-#include "source_io/module_current/td_current_io.h"
-#include "source_io/module_dipole/dipole_io.h"
 #include "source_io/module_output/output_log.h"
 #include "source_io/module_wf/read_wfc_nao.h"
 //------LCAO HSolver ElecState-------

source/source_esolver/esolver_ks_lcaopw.cpp

@@ -2,15 +2,13 @@
 
 #include "source_pw/module_pwdft/elecond.h"
 #include "source_io/module_parameter/input_conv.h"
-#include "source_io/module_output/output_log.h"
 
 #include <iostream>
 
 //--------------temporary----------------------------
 #include "source_estate/module_charge/symmetry_rho.h"
 #include "source_estate/occupy.h"
 #include "source_hamilt/module_ewald/H_Ewald_pw.h"
-#include "source_io/module_output/print_info.h"
 //-----force-------------------
 #include "source_pw/module_pwdft/forces.h"
 //-----stress------------------
@@ -23,11 +21,6 @@
 #include "source_hsolver/hsolver_lcaopw.h"
 #include "source_hsolver/kernels/hegvd_op.h"
 #include "source_base/kernels/math_kernel_op.h"
-#include "source_io/module_unk/berryphase.h"
-#include "source_io/module_bessel/numerical_basis.h"
-#include "source_io/module_bessel/numerical_descriptor.h"
-#include "source_io/module_wannier/to_wannier90_pw.h"
-#include "source_io/module_chgpot/write_elecstat_pot.h"
 #include "source_io/module_parameter/parameter.h"
 #include "source_hamilt/module_xc/xc_functional.h"
 

source/source_esolver/esolver_of.cpp

@@ -1,14 +1,10 @@
 #include "esolver_of.h"
 
 #include "source_io/module_parameter/parameter.h"
-#include "source_io/module_output/cube_io.h"
-#include "source_io/module_output/output_log.h"
-#include "source_io/module_chgpot/write_elecstat_pot.h"
 //-----------temporary-------------------------
 #include "source_base/global_function.h"
 #include "source_estate/module_charge/symmetry_rho.h"
 #include "source_hamilt/module_ewald/H_Ewald_pw.h"
-#include "source_io/module_output/print_info.h"
 #include "source_estate/cal_ux.h"
 #include "source_pw/module_pwdft/forces.h"
 #include "source_pw/module_ofdft/of_stress_pw.h"

source/source_esolver/esolver_of_tddft.cpp

@@ -1,14 +1,10 @@
 #include "esolver_of_tddft.h"
 
 #include "source_io/module_parameter/parameter.h"
-#include "source_io/module_output/cube_io.h"
-#include "source_io/module_output/output_log.h"
-#include "source_io/module_chgpot/write_elecstat_pot.h"
 //-----------temporary-------------------------
 #include "source_base/global_function.h"
 #include "source_estate/module_charge/symmetry_rho.h"
 #include "source_hamilt/module_ewald/H_Ewald_pw.h"
-#include "source_io/module_output/print_info.h"
 #include "source_estate/cal_ux.h"
 //-----force-------------------
 #include "source_pw/module_pwdft/forces.h"

source/source_esolver/lcao_others.cpp

@@ -10,20 +10,16 @@
 #include "source_cell/module_neighbor/sltk_grid_driver.h"
 #include "source_estate/elecstate_lcao.h"
 #include "source_estate/module_dm/cal_dm_psi.h"
-#include "source_io/module_unk/berryphase.h"
 #include "source_io/module_chgpot/get_pchg_lcao.h"
 #include "source_io/module_wf/get_wf_lcao.h"
 #include "source_io/module_parameter/parameter.h"
-#include "source_io/module_wf/read_wfc_nao.h"
 #include "source_io/module_hs/write_HS_R.h"
-#include "source_io/module_chgpot/write_elecstat_pot.h"
 #include "source_lcao/LCAO_domain.h"
 #include "source_lcao/module_deltaspin/spin_constrain.h"
 #include "source_lcao/module_operator_lcao/op_exx_lcao.h"
 #include "source_lcao/module_operator_lcao/operator_lcao.h"
 
 #ifdef __EXX
-#include "source_io/module_restart/restart_exx_csr.h"
 #endif
 
 // mohan add 2025-03-06

source/source_esolver/pw_others.cpp

@@ -1,6 +1,5 @@
 #include "esolver_ks_pw.h"
 #include "source_base/module_device/device.h"
-#include "source_io/module_bessel/numerical_basis.h"
 #include "source_io/module_bessel/numerical_descriptor.h"
 
 #include "source_base/formatter.h"

source/source_estate/elecstate_energy_terms.cpp

@@ -3,8 +3,6 @@
 #include "source_estate/module_pot/H_Hartree_pw.h"
 #include "source_estate/module_pot/efield.h"
 #include "source_estate/module_pot/gatefield.h"
-#include "source_lcao/module_deepks/LCAO_deepks.h"
-#include "source_lcao/module_deltaspin/spin_constrain.h"
 #include "source_lcao/module_dftu/dftu.h" // mohan add 2025-11-06
 
 namespace elecstate

source/source_estate/elecstate_print.cpp

@@ -6,7 +6,6 @@
 #include "source_estate/module_pot/efield.h"
 #include "source_estate/module_pot/gatefield.h"
 #include "source_hamilt/module_xc/xc_functional.h"
-#include "source_lcao/module_deepks/LCAO_deepks.h"
 #include "source_io/module_parameter/parameter.h"
 #include "occupy.h"
 namespace elecstate

source/source_io/module_parameter/input_conv.cpp

@@ -42,7 +42,6 @@
 #include "source_estate/module_pot/gatefield.h"
 #include "source_hsolver/hsolver_lcao.h"
 #include "source_hsolver/hsolver_pw.h"
-#include "source_md/md_func.h"
 #include "source_relax/bfgs_basic.h"
 #include "source_relax/ions_move_cg.h"
 

source/source_lcao/FORCE_gamma.cpp

@@ -13,7 +13,6 @@
 #include "source_estate/elecstate_lcao.h"
 #include "source_lcao/LCAO_domain.h"
 #include "source_lcao/pulay_fs.h"
-#include "source_io/module_hs/write_HS.h"
 
 template <>
 void Force_LCAO<double>::allocate(const UnitCell& ucell,

source/source_lcao/FORCE_k.cpp

@@ -10,7 +10,6 @@
 #include "source_estate/module_dm/cal_dm_psi.h"
 #include "source_lcao/LCAO_domain.h"
 #include "source_lcao/pulay_fs.h"
-#include "source_io/module_hs/write_HS.h"
 #include "source_io/module_parameter/parameter.h"
 
 #include <map>

source/source_lcao/module_rdmft/rdmft_tools.cpp

@@ -10,7 +10,6 @@
 #include "source_estate/module_pot/H_Hartree_pw.h"
 #include "source_estate/module_pot/pot_local.h"
 #include "source_estate/module_pot/pot_xc.h"
-#include "source_pw/module_pwdft/structure_factor.h"
 #include "source_lcao/module_gint/gint_interface.h"
 #include "source_io/module_parameter/parameter.h"
 

source/source_lcao/wavefunc_in_pw.cpp

@@ -5,7 +5,6 @@
 #include "source_base/math_sphbes.h"
 #include "source_base/math_polyint.h"
 #include "source_base/math_ylmreal.h"
-#include "source_pw/module_pwdft/soc.h"
 
 void Wavefunc_in_pw::make_table_q(
 	const UnitCell &ucell,

source/source_pw/module_pwdft/forces_cc.cpp

@@ -2,7 +2,6 @@
 #include "stress_func.h"
 #include "source_base/parallel_reduce.h"
 #include "source_io/module_parameter/parameter.h"
-#include "source_io/module_output/output_log.h"
 // new
 #include "source_base/complexmatrix.h"
 #include "source_base/libm/libm.h"

source/source_pw/module_pwdft/forces_scc.cpp

@@ -1,6 +1,5 @@
 #include "forces.h"
 #include "source_base/parallel_reduce.h"
-#include "source_io/module_output/output_log.h"
 #include "stress_func.h"
 // new
 #include "source_base/complexmatrix.h"