fix bug of #include

Author: PeizeLin

source/source_base/kernels/cuda/math_kernel_op_vec.cu

@@ -2,6 +2,7 @@
 
 #include <base/macros/macros.h>
 #include <thrust/complex.h>
+#include "source_base/parallel_reduce.h"
 
 template <>
 struct GetTypeReal<thrust::complex<float>> {

source/source_base/kernels/rocm/math_kernel_op_vec.hip.cu

@@ -2,6 +2,7 @@
 
 #include <base/macros/macros.h>
 #include <thrust/complex.h>
+#include "source_base/parallel_reduce.h"
 template <>
 struct GetTypeReal<thrust::complex<float>> {
     using type = float; /**< The return type specialization for std::complex<double>. */

source/source_lcao/FORCE_STRESS.cpp

@@ -24,6 +24,8 @@
 #include "source_lcao/module_operator_lcao/overlap.h"
 #include "source_lcao/pulay_fs.h"
 
+#include "source_base/parallel_reduce.h"
+
 #ifdef __MPI
 #include "source_base/parallel_reduce.h"
 #endif

source/source_lcao/module_deltaspin/cal_mw.cpp

@@ -11,6 +11,8 @@
 #include "source_estate/elecstate_lcao.h"
 #include "source_lcao/hamilt_lcao.h"
 #include "source_lcao/module_operator_lcao/dspin_lcao.h"
+#endif
+
 #include "source_base/parallel_reduce.h"
 
 template <>
@@ -51,8 +53,6 @@ void spinconstrain::SpinConstrain<std::complex<double>>::cal_mi_lcao(const int&
     ModuleBase::timer::end("spinconstrain::SpinConstrain", "cal_mi_lcao");
 }
 
-#endif
-
 template <>
 void spinconstrain::SpinConstrain<std::complex<double>>::cal_mi_pw()
 {