refactor(dftu): prune redundant #include in module_dftu

Remove unused #include directives across 12 files in
source/source_lcao/module_dftu/. Each removal was verified by
symbol-usage grep against the file body (excluding the include lines
themselves) and by compiling the dftu + hamilt_lcao CMake targets.

Removed includes per file:
- dftu_lcao.cpp: parameter.h, constants.h, global_function.h,
  inverse_matrix.h, memory_recorder.h, magnetism.h, charge.h and
  unused <cstdint>/<cmath>/<cstdio>/<cstring>/<fstream>/<iomanip>/
  <iostream>/<sstream>; add focused matrix.h (for
  ModuleBase::matrix::operator() via occ_mat), tool_quit.h, tool_title.h
- dftu_folding.cpp: parameter.h, hcontainer.h, hcontainer_funcs.h
- dftu_force.cpp: parameter.h, constants.h, inverse_matrix.h,
  elecstate_lcao.h, magnetism.h, charge.h, <cmath>, <fstream>,
  <iomanip>, <iostream>, <sstream>, <stdio.h>, <string.h>
- dftu_hamilt.cpp: parameter.h
- dftu_lcao_op.cpp: hcontainer_funcs.h, <unordered_set>
- dftu_lcao_op_legacy.cpp: dftu_lcao.h (already provided by its own .h)
- dftu_occup.cpp: parameter.h
- dftu_tools.cpp: timer.h (no direct ModuleBase::timer use), parameter.h
- dftu_yukawa.cpp: parameter.h, <complex>, <fstream>, <iomanip>,
  <iostream>, <sstream>, <cstdio>, <cstring>
- dftu_lcao.h: charge_mixing.h (still provided transitively via
  dftu_base.h), force_stress_arrays.h (not directly used here)
- dftu_lcao_op.h: density_matrix.h (not directly used here)
- dftu_lcao_op_legacy.h: timer.h (not directly used here)

Note: scalapack_connector.h was kept in dftu_force.cpp / dftu_hamilt.cpp
/ dftu_occup.cpp -- these files call ScalapackConnector::gemm.

Verification: OMP_NUM_THREADS=1 make -C build_std_para dftu hamilt_lcao
=> EXIT=0, no errors, no warnings. Full abacus build deferred to user;
external modules that include dftu_lcao.h and use Charge_Mixing
(chgmixing.cpp, esolver_ks.cpp/h) should still resolve via
dftu_lcao.h -> dftu_base.h -> charge_mixing.h.

Author: abacus_fixer

source/source_lcao/module_dftu/dftu_folding.cpp

@@ -2,11 +2,8 @@
 #include "dftu_folding.h"
 #include "dftu_lcao.h"
 #include "source_base/timer.h"
-#include "source_io/module_parameter/parameter.h"
 #include "source_cell/module_neighbor/sltk_grid_driver.h"
 #include "source_lcao/hamilt_lcao.h"
-#include "source_hamilt/module_hcontainer/hcontainer.h"
-#include "source_hamilt/module_hcontainer/hcontainer_funcs.h"
 
 namespace DFTU_LCAO {
 

source/source_lcao/module_dftu/dftu_force.cpp

@@ -1,27 +1,13 @@
-#include "source_io/module_parameter/parameter.h"
-
 #ifdef __LCAO
 #include "dftu_force.h"
 #include "dftu_folding.h"
 #include "dftu_lcao.h"
-#include "source_base/constants.h"
 #include "source_base/global_function.h"
-#include "source_base/inverse_matrix.h"
 #include "source_base/module_external/scalapack_connector.h"
 #include "source_base/parallel_reduce.h"
 #include "source_base/timer.h"
-#include "source_estate/elecstate_lcao.h"
-#include "source_cell/magnetism.h"
-#include "source_estate/module_charge/charge.h"
 
-#include <cmath>
 #include <complex>
-#include <fstream>
-#include <iomanip>
-#include <iostream>
-#include <sstream>
-#include <stdio.h>
-#include <string.h>
 #include <string>
 
 

source/source_lcao/module_dftu/dftu_hamilt.cpp

@@ -1,6 +1,5 @@
 #include "dftu_lcao.h"
 #include "source_base/module_external/scalapack_connector.h"
-#include "source_io/module_parameter/parameter.h"
 #include "source_base/timer.h"
 
 

source/source_lcao/module_dftu/dftu_lcao.cpp

@@ -1,23 +1,11 @@
 #include "dftu_lcao.h"
 
-#include "source_io/module_parameter/parameter.h"
-#include "source_base/constants.h"
-#include "source_base/global_function.h"
-#include "source_base/inverse_matrix.h"
-#include "source_base/memory_recorder.h"
+#include "source_base/matrix.h"  // occ_mat uses ModuleBase::matrix::operator()
+#include "source_base/tool_quit.h"
+#include "source_base/tool_title.h"
 #include "source_base/timer.h"
-#include "source_cell/magnetism.h"
-#include "source_estate/module_charge/charge.h"
 
-#include <cstdint>
-#include <cmath>
 #include <complex>
-#include <cstdio>
-#include <cstring>
-#include <fstream>
-#include <iomanip>
-#include <iostream>
-#include <sstream>
 #include <vector>
 
  // mohan add 2025-11-06

source/source_lcao/module_dftu/dftu_lcao.h

@@ -4,14 +4,12 @@
 #include "source_cell/klist.h"
 #include "source_cell/unitcell.h"
 #include "source_basis/module_ao/parallel_orbitals.h"
-#include "source_estate/module_charge/charge_mixing.h"
 #include "source_pw/module_pwdft/dftu_base.h"
 #ifdef __LCAO
 #include "source_basis/module_ao/orb_read.h"
 #include "source_hamilt/hamilt.h"
 #include "source_hamilt/module_hcontainer/hcontainer.h"
 #include "source_estate/module_dm/density_matrix.h"
-#include "source_lcao/force_stress_arrays.h" // mohan add 2024-06-15
 #endif
 
 #include <string>

source/source_lcao/module_dftu/dftu_lcao_op.cpp

@@ -4,11 +4,7 @@
 #include "source_base/tool_title.h"
 #include "source_cell/module_neighbor/sltk_grid_driver.h"
 #include "source_lcao/module_operator_lcao/operator_lcao.h"
-#include "source_hamilt/module_hcontainer/hcontainer_funcs.h"
 #include "source_io/module_parameter/parameter.h"
-#ifdef _OPENMP
-#include <unordered_set>
-#endif
 #include "source_base/parallel_reduce.h"
 
 template <typename TK, typename TR>

source/source_lcao/module_dftu/dftu_lcao_op.h

@@ -4,7 +4,6 @@
 #include "source_basis/module_nao/two_center_integrator.h"
 #include "source_cell/module_neighbor/sltk_grid_driver.h"
 #include "source_cell/unitcell.h"
-#include "source_estate/module_dm/density_matrix.h"
 #include "source_lcao/module_operator_lcao/operator_lcao.h"
 #include "source_lcao/module_dftu/dftu_lcao.h"
 #include "source_hamilt/module_hcontainer/hcontainer.h"

source/source_lcao/module_dftu/dftu_lcao_op_legacy.cpp

@@ -1,7 +1,6 @@
 #include "dftu_lcao_op_legacy.h"
 #include "source_base/timer.h"
 #include "source_base/tool_title.h"
-#include "source_lcao/module_dftu/dftu_lcao.h"
 
 namespace hamilt
 {

source/source_lcao/module_dftu/dftu_lcao_op_legacy.h

@@ -1,7 +1,6 @@
 #ifndef OPDFTULCAO_H
 #define OPDFTULCAO_H
 
-#include "source_base/timer.h"
 #include "source_lcao/module_operator_lcao/operator_lcao.h"
 #include "source_lcao/module_dftu/dftu_lcao.h" // mohan add 20251107
 

source/source_lcao/module_dftu/dftu_occup.cpp

@@ -1,11 +1,10 @@
 #include "dftu_lcao.h"
 #include "dftu_folding.h"
 #include "source_base/timer.h"
-#include "source_io/module_parameter/parameter.h"
+#include "source_base/module_external/scalapack_connector.h"
 #ifdef __LCAO
 #include "source_lcao/hamilt_lcao.h"
 #endif
-#include "source_base/module_external/scalapack_connector.h"
 
 // copy_occ_mat(), zero_occ_mat(), mix_occ_mat(), set_occ_mat(ucell),
 // get_occ_mat_flat(), set_occ_mat_flat()