type transformation

lyb9812 · lyb9812 · commit a917e0b58ea0 · 2025-10-12T08:41:16.000+08:00
diff --git a/source/source_pw/module_ofdft/evolve_ofdft.cpp b/source/source_pw/module_ofdft/evolve_ofdft.cpp
@@ -185,7 +185,7 @@ void Evolve_OFDFT::cal_CD_potential(std::vector<std::complex<double>> psi_,
 
         for (int ir = 0; ir < pw_rho->nrxx; ++ir)
         {
-            rpot(0, ir) -= mCD_para*std::real(2.0*rCDPotential[ir]*std::pow(ModuleBase::PI,3) / (2*std::pow(kF_r[ir],2)));
+            rpot(0, ir) -= mCD_para*2.0*std::real(rCDPotential[ir])*std::pow(ModuleBase::PI,3) / (2.0*std::pow(std::real(kF_r[ir]),2));
         }
         delete[] recipCurrent_x;
         delete[] recipCurrent_y;

Original file line number	Diff line number	Diff line change
`@@ -185,7 +185,7 @@ void Evolve_OFDFT::cal_CD_potential(std::vector<std::complex<double>> psi_,`
`185`	`185`
`186`	`186`	`for (int ir = 0; ir < pw_rho->nrxx; ++ir)`
`187`	`187`	`{`
`188`		`- rpot(0, ir) -= mCD_parastd::real(2.0rCDPotential[ir]std::pow(ModuleBase::PI,3) / (2std::pow(kF_r[ir],2)));`
	`188`	`+ rpot(0, ir) -= mCD_para2.0std::real(rCDPotential[ir])std::pow(ModuleBase::PI,3) / (2.0std::pow(std::real(kF_r[ir]),2));`
`189`	`189`	`}`
`190`	`190`	`delete[] recipCurrent_x;`
`191`	`191`	`delete[] recipCurrent_y;`