diff --git a/source/source_base/kernels/math_kernel_op.cpp b/source/source_base/kernels/math_kernel_op.cpp
index 9f57b51cb8..dad790e463 100644
--- a/source/source_base/kernels/math_kernel_op.cpp
+++ b/source/source_base/kernels/math_kernel_op.cpp
@@ -1,4 +1,5 @@
 #include "source_base/kernels/math_kernel_op.h"
+#include "source_base/module_external/blas_connector.h"
 
 #include <iomanip>
 #include <iostream>
diff --git a/source/source_base/kernels/math_kernel_op.h b/source/source_base/kernels/math_kernel_op.h
index cb35663f7d..9476d32fdc 100644
--- a/source/source_base/kernels/math_kernel_op.h
+++ b/source/source_base/kernels/math_kernel_op.h
@@ -5,7 +5,6 @@
 
 #include "source_base/macros.h"
 
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/parallel_reduce.h"
 
 #include "source_base/module_device/memory_op.h"
diff --git a/source/source_base/kernels/math_kernel_op_vec.cpp b/source/source_base/kernels/math_kernel_op_vec.cpp
index 6d35a956e1..8353b82660 100644
--- a/source/source_base/kernels/math_kernel_op_vec.cpp
+++ b/source/source_base/kernels/math_kernel_op_vec.cpp
@@ -1,4 +1,6 @@
 #include "source_base/kernels/math_kernel_op.h"
+#include "source_base/module_external/blas_connector.h"
+
 
 namespace ModuleBase
 {
diff --git a/source/source_base/module_grid/batch.cpp b/source/source_base/module_grid/batch.cpp
index b11c22e4ce..718d9277d3 100644
--- a/source/source_base/module_grid/batch.cpp
+++ b/source/source_base/module_grid/batch.cpp
@@ -4,7 +4,6 @@
 #include <cassert>
 #include <iterator>
 
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/module_external/lapack_connector.h"
 
 namespace {
diff --git a/source/source_hsolver/diago_bpcg.cpp b/source/source_hsolver/diago_bpcg.cpp
index 54d76d6d75..45c9b8b24d 100644
--- a/source/source_hsolver/diago_bpcg.cpp
+++ b/source/source_hsolver/diago_bpcg.cpp
@@ -1,7 +1,6 @@
 #include "source_hsolver/diago_bpcg.h"
 
 #include "diago_iter_assist.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/global_function.h"
 #include "source_base/kernels/math_kernel_op.h"
 #include "source_base/parallel_comm.h" // different MPI worlds
diff --git a/source/source_hsolver/diago_iter_assist.cpp b/source/source_hsolver/diago_iter_assist.cpp
index 49bc307324..5ce64062f2 100644
--- a/source/source_hsolver/diago_iter_assist.cpp
+++ b/source/source_hsolver/diago_iter_assist.cpp
@@ -1,6 +1,5 @@
 #include "diago_iter_assist.h"
 #include "source_io/module_parameter/parameter.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/complexmatrix.h"
 #include "source_base/constants.h"
 #include "source_base/global_variable.h"
diff --git a/source/source_lcao/module_gint/gint_gamma_vl.cpp b/source/source_lcao/module_gint/gint_gamma_vl.cpp
index 2546bb8cbe..161a3e7083 100644
--- a/source/source_lcao/module_gint/gint_gamma_vl.cpp
+++ b/source/source_lcao/module_gint/gint_gamma_vl.cpp
@@ -4,7 +4,6 @@
 #include "gint_gamma.h"
 #include "gint_tools.h"
 #include "grid_technique.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/memory.h"
 #include "source_base/timer.h"
 #include "source_basis/module_ao/ORB_read.h"
diff --git a/source/source_lcao/module_gint/gint_rho_old.cpp b/source/source_lcao/module_gint/gint_rho_old.cpp
index 00e26cbd21..b3027d6b12 100644
--- a/source/source_lcao/module_gint/gint_rho_old.cpp
+++ b/source/source_lcao/module_gint/gint_rho_old.cpp
@@ -1,7 +1,6 @@
 #include "gint_k.h"
 #include "gint_tools.h"
 #include "grid_technique.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/global_function.h"
 #include "source_base/global_variable.h"
 #include "source_base/timer.h"
diff --git a/source/source_lcao/module_gint/gint_vl_old.cpp b/source/source_lcao/module_gint/gint_vl_old.cpp
index a0793a6ce3..9ebc341d7f 100644
--- a/source/source_lcao/module_gint/gint_vl_old.cpp
+++ b/source/source_lcao/module_gint/gint_vl_old.cpp
@@ -5,7 +5,6 @@
 #include "grid_technique.h"
 #include "source_base/ylm.h"
 #include "source_pw/module_pwdft/global.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/timer.h"
 #include "source_base/array_pool.h"
 #include "source_base/vector3.h"
diff --git a/source/source_lcao/module_gint/temp_gint/gint_vl.cpp b/source/source_lcao/module_gint/temp_gint/gint_vl.cpp
index f4470d2873..3cdd3c4549 100644
--- a/source/source_lcao/module_gint/temp_gint/gint_vl.cpp
+++ b/source/source_lcao/module_gint/temp_gint/gint_vl.cpp
@@ -1,4 +1,3 @@
-#include "source_base/module_external/blas_connector.h"
 #include "gint_common.h"
 #include "gint_vl.h"
 #include "phi_operator.h"
diff --git a/source/source_lcao/module_gint/temp_gint/gint_vl_metagga.cpp b/source/source_lcao/module_gint/temp_gint/gint_vl_metagga.cpp
index 6ba02e9bcd..56dd2ff6b0 100644
--- a/source/source_lcao/module_gint/temp_gint/gint_vl_metagga.cpp
+++ b/source/source_lcao/module_gint/temp_gint/gint_vl_metagga.cpp
@@ -1,4 +1,3 @@
-#include "source_base/module_external/blas_connector.h"
 #include "gint_common.h"
 #include "gint_vl_metagga.h"
 #include "phi_operator.h"
diff --git a/source/source_lcao/module_gint/temp_gint/gint_vl_metagga_nspin4.cpp b/source/source_lcao/module_gint/temp_gint/gint_vl_metagga_nspin4.cpp
index 3066195ed4..d5880600b0 100644
--- a/source/source_lcao/module_gint/temp_gint/gint_vl_metagga_nspin4.cpp
+++ b/source/source_lcao/module_gint/temp_gint/gint_vl_metagga_nspin4.cpp
@@ -1,6 +1,5 @@
 #include "source_base/global_function.h"
 #include "source_lcao/module_hcontainer/hcontainer.h"
-#include "source_base/module_external/blas_connector.h"
 #include "gint_common.h"
 #include "gint_vl_metagga_nspin4.h"
 #include "phi_operator.h"
diff --git a/source/source_lcao/module_gint/temp_gint/gint_vl_nspin4.cpp b/source/source_lcao/module_gint/temp_gint/gint_vl_nspin4.cpp
index 18f1230cc7..56aca5edeb 100644
--- a/source/source_lcao/module_gint/temp_gint/gint_vl_nspin4.cpp
+++ b/source/source_lcao/module_gint/temp_gint/gint_vl_nspin4.cpp
@@ -1,6 +1,5 @@
 #include "source_base/global_function.h"
 #include "source_lcao/module_hcontainer/hcontainer.h"
-#include "source_base/module_external/blas_connector.h"
 #include "gint_common.h"
 #include "gint_vl_nspin4.h"
 #include "phi_operator.h"
diff --git a/source/source_lcao/module_lr/operator_casida/operator_lr_hxc.cpp b/source/source_lcao/module_lr/operator_casida/operator_lr_hxc.cpp
index 58c30ad8ae..56d2e4fda7 100644
--- a/source/source_lcao/module_lr/operator_casida/operator_lr_hxc.cpp
+++ b/source/source_lcao/module_lr/operator_casida/operator_lr_hxc.cpp
@@ -1,7 +1,6 @@
 #include "operator_lr_hxc.h"
 #include <vector>
 #include "source_io/module_parameter/parameter.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/timer.h"
 #include "source_lcao/module_lr/utils/lr_util.h"
 #include "source_lcao/module_lr/utils/lr_util_hcontainer.h"
diff --git a/source/source_lcao/module_operator_lcao/dspin_lcao.cpp b/source/source_lcao/module_operator_lcao/dspin_lcao.cpp
index 07f2217c10..d4ebcaaf7e 100644
--- a/source/source_lcao/module_operator_lcao/dspin_lcao.cpp
+++ b/source/source_lcao/module_operator_lcao/dspin_lcao.cpp
@@ -1,6 +1,5 @@
 #include "dspin_lcao.h"
 #include "source_lcao/module_deltaspin/spin_constrain.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/timer.h"
 #include "source_base/memory.h"
 #include "source_base/tool_title.h"
diff --git a/source/source_lcao/module_ri/module_exx_symmetry/symmetry_rotation_R.hpp b/source/source_lcao/module_ri/module_exx_symmetry/symmetry_rotation_R.hpp
index 2fe48a61fb..0aef2a2ba6 100644
--- a/source/source_lcao/module_ri/module_exx_symmetry/symmetry_rotation_R.hpp
+++ b/source/source_lcao/module_ri/module_exx_symmetry/symmetry_rotation_R.hpp
@@ -1,6 +1,5 @@
 #include "symmetry_rotation.h"
 #include "source_lcao/module_ri/RI_Util.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/timer.h"
 #include <array>
 #include <RI/global/Global_Func-2.h>
diff --git a/source/source_pw/module_pwdft/hamilt_pw.cpp b/source/source_pw/module_pwdft/hamilt_pw.cpp
index 57854d256c..86782ed665 100644
--- a/source/source_pw/module_pwdft/hamilt_pw.cpp
+++ b/source/source_pw/module_pwdft/hamilt_pw.cpp
@@ -1,7 +1,6 @@
 #include "hamilt_pw.h"
 
 #include "source_io/module_parameter/parameter.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/global_function.h"
 #include "source_base/global_variable.h"
 #include "source_pw/module_pwdft/global.h"
diff --git a/source/source_pw/module_pwdft/onsite_projector.cpp b/source/source_pw/module_pwdft/onsite_projector.cpp
index 49df42e08b..4c03ce9b74 100644
--- a/source/source_pw/module_pwdft/onsite_projector.cpp
+++ b/source/source_pw/module_pwdft/onsite_projector.cpp
@@ -7,7 +7,6 @@
 #include "source_pw/module_pwdft/onsite_projector.h"
 
 #include "source_base/projgen.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/kernels/math_kernel_op.h"
 #ifdef __MPI
 #include "source_base/parallel_reduce.h"
diff --git a/source/source_pw/module_pwdft/operator_pw/nonlocal_pw.cpp b/source/source_pw/module_pwdft/operator_pw/nonlocal_pw.cpp
index eb6c2e839c..2a87f3c1a7 100644
--- a/source/source_pw/module_pwdft/operator_pw/nonlocal_pw.cpp
+++ b/source/source_pw/module_pwdft/operator_pw/nonlocal_pw.cpp
@@ -1,7 +1,6 @@
 #include "nonlocal_pw.h"
 
 #include "source_io/module_parameter/parameter.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/timer.h"
 #include "source_base/parallel_reduce.h"
 #include "source_base/tool_quit.h"
diff --git a/source/source_pw/module_pwdft/operator_pw/onsite_proj_pw.cpp b/source/source_pw/module_pwdft/operator_pw/onsite_proj_pw.cpp
index 2ec2e6f09f..56b6d3756a 100644
--- a/source/source_pw/module_pwdft/operator_pw/onsite_proj_pw.cpp
+++ b/source/source_pw/module_pwdft/operator_pw/onsite_proj_pw.cpp
@@ -1,6 +1,5 @@
 #include "onsite_proj_pw.h"
 
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/timer.h"
 #include "source_base/parallel_reduce.h"
 #include "source_base/tool_quit.h"
diff --git a/source/source_pw/module_pwdft/operator_pw/op_exx_pw.h b/source/source_pw/module_pwdft/operator_pw/op_exx_pw.h
index 72b0a209ae..ee9501df66 100644
--- a/source/source_pw/module_pwdft/operator_pw/op_exx_pw.h
+++ b/source/source_pw/module_pwdft/operator_pw/op_exx_pw.h
@@ -2,7 +2,6 @@
 #define OPEXXPW_H
 
 #include "operator_pw.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/kernels/math_kernel_op.h"
 #include "source_base/macros.h"
 #include "source_base/matrix.h"
diff --git a/source/source_pw/module_stodft/sto_che.cpp b/source/source_pw/module_stodft/sto_che.cpp
index 858d5cf820..f59cbd08fc 100644
--- a/source/source_pw/module_stodft/sto_che.cpp
+++ b/source/source_pw/module_stodft/sto_che.cpp
@@ -1,5 +1,4 @@
 #include "sto_che.h"
-#include "source_base/module_external/blas_connector.h"
 #include "source_base/module_device/device.h"
 #include "source_base/kernels/math_kernel_op.h"
 #include "source_base/module_container/ATen/kernels/blas.h"