fix process

Author: A-006

source/module_basis/module_pw/module_fft/fft_bundle.cpp

@@ -57,10 +57,7 @@ void FFT_Bundle::initfft(int nx_in,
         }
 #endif
     }
-    if (this->precision == "double" || this->precision == "mixing")
-    {
         double_flag = true;
-    }
 #if defined(__DSP)
     if (device == "dsp")
     {

source/module_basis/module_pw/pw_basis.cpp

@@ -15,7 +15,11 @@ PW_Basis::PW_Basis()
 
 PW_Basis::PW_Basis(std::string device_, std::string precision_) : device(std::move(device_)), precision(std::move(precision_)) {
     classname="PW_Basis";
-    this->fft_bundle.setfft("cpu",this->precision);
+    #if (not defined(__ENABLE_FLOAT_FFTW) and (defined(__CUDA)))
+        this->fft_bundle.setfft("cpu","double");
+    #else
+        this->fft_bundle.setfft("cpu",this->precision);
+    #endif
     this->double_data_ = (this->precision == "double") || (this->precision == "mixing");
     this->float_data_ = (this->precision == "single") || (this->precision == "mixing");
 }

source/module_basis/module_pw/pw_transform.cpp

@@ -210,7 +210,7 @@ void PW_Basis::recip2real(const std::complex<FPTYPE>* in, FPTYPE* out, const boo
 #endif
     for (int i = 0; i < this->nst * this->nz; ++i)
     {
-        fft_bundle.get_auxg_data<FPTYPE>()[i] = std::complex<double>(0, 0);
+        fft_bundle.get_auxg_data<FPTYPE>()[i] = std::complex<FPTYPE>(0, 0);
     }
 
 #ifdef _OPENMP

source/module_hamilt_pw/hamilt_pwdft/structure_factor.cpp

@@ -68,8 +68,7 @@ void Structure_Factor::setup_structure_factor(const UnitCell* Ucell, const Paral
 //	std::ofstream ofs( outstr.c_str() ) ;
     bool usebspline;
     if(nbspline > 0) {   usebspline = true;
-    } else {    usebspline = false;
-}
+    } else {    usebspline = false;}
 
     if(usebspline)
     {
@@ -147,6 +146,7 @@ void Structure_Factor::setup_structure_factor(const UnitCell* Ucell, const Paral
             inat++;
         }
     }
+    
     if (device == "gpu") {
         if (PARAM.globalv.has_float_data) {
             resmem_cd_op()(this->c_eigts1, Ucell->nat * (2 * rho_basis->nx + 1));
@@ -164,7 +164,9 @@ void Structure_Factor::setup_structure_factor(const UnitCell* Ucell, const Paral
         syncmem_z2z_h2d_op()(this->z_eigts3, this->eigts3.c, Ucell->nat * (2 * rho_basis->nz + 1));
     }
     else {
+        std::cout<<"the test is "<<std::endl;
         if (PARAM.globalv.has_float_data) {
+            std::cout<<"here set the float data"<<std::endl;
             resmem_ch_op()(this->c_eigts1, Ucell->nat * (2 * rho_basis->nx + 1));
             resmem_ch_op()(this->c_eigts2, Ucell->nat * (2 * rho_basis->ny + 1));
             resmem_ch_op()(this->c_eigts3, Ucell->nat * (2 * rho_basis->nz + 1));

source/module_io/read_set_globalv.cpp

@@ -72,7 +72,7 @@ void ReadInput::set_globalv(const Input_para& inp, System_para& sys)
     bool float_cond = false;
 #endif
     sys.has_double_data = (inp.precision == "double") || (inp.precision == "mixing") || float_cond;
-    sys.has_float_data = (inp.precision == "float") || (inp.precision == "mixing") || float_cond;
+    sys.has_float_data = (inp.precision == "single") || (inp.precision == "mixing") || float_cond;
 }
 
 /// @note Here para.inp has not been synchronized of all ranks.