exchange '/bignxyz' in recip2real and real2recip

Qianruipku · Qianruipku · commit b426c1fa8dae · 2022-01-09T20:06:11.000+08:00
diff --git a/source/module_pw/pw_transform.cpp b/source/module_pw/pw_transform.cpp
@@ -35,7 +35,7 @@ void PW_Basis:: real2recip(std::complex<double> * in, std::complex<double> * out
 
     for(int ig = 0 ; ig < this->npw ; ++ig)
     {
-        out[ig] = this->ft.aux1[this->ig2isz[ig]];
+        out[ig] = this->ft.aux1[this->ig2isz[ig]] / double(this->bignxyz);
     }
     return;
     ModuleBase::timer::tick("PW_Basis", "real2recip");
@@ -72,7 +72,7 @@ void PW_Basis:: real2recip(double * in, std::complex<double> * out)
 
     for(int ig = 0 ; ig < this->npw ; ++ig)
     {
-        out[ig] = this->ft.aux1[this->ig2isz[ig]];
+        out[ig] = this->ft.aux1[this->ig2isz[ig]] / double(this->bignxyz);
     }
     ModuleBase::timer::tick("PW_Basis", "real2recip_gamma_only");
     return;
@@ -105,7 +105,7 @@ void PW_Basis:: recip2real(std::complex<double> * in, std::complex<double> * out
     
     for(int ir = 0 ; ir < this->nrxx ; ++ir)
     {
-        out[ir] = this->ft.aux1[ir] / double(this->bignxyz);
+        out[ir] = this->ft.aux1[ir];
     }
     ModuleBase::timer::tick("PW_Basis", "recip2real");
 
@@ -144,7 +144,7 @@ void PW_Basis:: recip2real(std::complex<double> * in, double * out)
     {
         for(int ipy = 0 ; ipy < npy ; ++ipy)
         {
-            out[ix*npy + ipy] = this->ft.r_rspace[ix*npy*2 + ipy] / double(this->bignxyz);
+            out[ix*npy + ipy] = this->ft.r_rspace[ix*npy*2 + ipy];
         }
     }
     ModuleBase::timer::tick("PW_Basis", "recip2real_gamma_only");
@@ -167,7 +167,7 @@ void PW_Basis:: real2recip(std::complex<float> * in, std::complex<float> * out)
 
     for(int ig = 0 ; ig < this->npw ; ++ig)
     {
-        out[ig] = this->ft.auxf1[this->ig2isz[ig]];
+        out[ig] = this->ft.auxf1[this->ig2isz[ig]] / float(this->bignxyz);
     }
     return;
 }
@@ -192,7 +192,7 @@ void PW_Basis:: real2recip(float * in, std::complex<float> * out)
 
     for(int ig = 0 ; ig < this->npw ; ++ig)
     {
-        out[ig] = this->ft.auxf1[this->ig2isz[ig]];
+        out[ig] = this->ft.auxf1[this->ig2isz[ig]] / float(this->bignxyz);
     }
     return;
 }
@@ -214,7 +214,7 @@ void PW_Basis:: recip2real(std::complex<float> * in, std::complex<float> * out)
     
     for(int ir = 0 ; ir < this->nrxx ; ++ir)
     {
-        out[ir] = this->ft.auxf1[ir] / float(this->bignxyz);
+        out[ir] = this->ft.auxf1[ir];
     }
 
     return;
@@ -240,7 +240,7 @@ void PW_Basis:: recip2real(std::complex<float> * in, float * out)
     {
         for(int ipy = 0 ; ipy < npy ; ++ipy)
         {
-            out[ix*npy + ipy] = this->ft.rf_rspace[ix*npy*2 + ipy] / float(this->bignxyz);
+            out[ix*npy + ipy] = this->ft.rf_rspace[ix*npy*2 + ipy];
         }
     }
     return;
diff --git a/source/module_pw/pw_transform_k.cpp b/source/module_pw/pw_transform_k.cpp
@@ -31,7 +31,7 @@ void PW_Basis_K:: real2recip(std::complex<double> * in, std::complex<double> * o
 
     for(int igl = 0 ; igl < this->npwk[ik] ; ++igl)
     {
-        out[igl] = this->ft.aux1[this->igl2isz_k[igl+ik*this->npwk_max]];
+        out[igl] = this->ft.aux1[this->igl2isz_k[igl+ik*this->npwk_max]] / double(this->bignxyz);
     }
     return;
     ModuleBase::timer::tick("PW_Basis_K", "real2recip");
@@ -68,7 +68,7 @@ void PW_Basis_K:: real2recip(double * in, std::complex<double> * out, int ik)
 
     for(int igl = 0 ; igl < this->npwk[ik] ; ++igl)
     {
-        out[igl] = this->ft.aux1[this->igl2isz_k[igl+ik*this->npwk_max]];
+        out[igl] = this->ft.aux1[this->igl2isz_k[igl+ik*this->npwk_max]] / double(this->bignxyz);
     }
     ModuleBase::timer::tick("PW_Basis_K", "real2recip_gamma_only");
     return;
@@ -97,7 +97,7 @@ void PW_Basis_K:: recip2real(std::complex<double> * in, std::complex<double> * o
     
     for(int ir = 0 ; ir < this->nrxx ; ++ir)
     {
-        out[ir] = this->ft.aux1[ir] / double(this->bignxyz);
+        out[ir] = this->ft.aux1[ir];
     }
     ModuleBase::timer::tick("PW_Basis_K", "recip2real");
 
@@ -136,7 +136,7 @@ void PW_Basis_K:: recip2real(std::complex<double> * in, double * out, int ik)
     {
         for(int ipy = 0 ; ipy < npy ; ++ipy)
         {
-            out[ix*npy + ipy] = this->ft.r_rspace[ix*npy*2 + ipy] / double(this->bignxyz);
+            out[ix*npy + ipy] = this->ft.r_rspace[ix*npy*2 + ipy];
         }
     }
     ModuleBase::timer::tick("PW_Basis_K", "recip2real_gamma_only");
@@ -159,7 +159,7 @@ void PW_Basis_K:: real2recip(std::complex<float> * in, std::complex<float> * out
 
     for(int igl = 0 ; igl < this->npwk[ik] ; ++igl)
     {
-        out[igl] = this->ft.auxf1[this->igl2isz_k[igl+ik*this->npwk_max]];
+        out[igl] = this->ft.auxf1[this->igl2isz_k[igl+ik*this->npwk_max]] / float(this->bignxyz);
     }
     return;
 }
@@ -184,7 +184,7 @@ void PW_Basis_K:: real2recip(float * in, std::complex<float> * out, int ik)
 
     for(int igl = 0 ; igl < this->npwk[ik] ; ++igl)
     {
-        out[igl] = this->ft.auxf1[this->igl2isz_k[igl+ik*this->npwk_max]];
+        out[igl] = this->ft.auxf1[this->igl2isz_k[igl+ik*this->npwk_max]] / float(this->bignxyz);
     }
     return;
 }
@@ -206,7 +206,7 @@ void PW_Basis_K:: recip2real(std::complex<float> * in, std::complex<float> * out
     
     for(int ir = 0 ; ir < this->nrxx ; ++ir)
     {
-        out[ir] = this->ft.auxf1[ir] / float(this->bignxyz);
+        out[ir] = this->ft.auxf1[ir];
     }
 
     return;
@@ -232,7 +232,7 @@ void PW_Basis_K:: recip2real(std::complex<float> * in, float * out, int ik)
     {
         for(int ipy = 0 ; ipy < npy ; ++ipy)
         {
-            out[ix*npy + ipy] = this->ft.rf_rspace[ix*npy*2 + ipy] / float(this->bignxyz);
+            out[ix*npy + ipy] = this->ft.rf_rspace[ix*npy*2 + ipy];
         }
     }
     return;
diff --git a/source/module_pw/test/test1-2.cpp b/source/module_pw/test/test1-2.cpp
@@ -87,7 +87,6 @@ TEST_F(PWTEST,test1_2)
                 ModuleBase::Vector3<double> real_r(ix, iy, iz);
                 double phase_im = -delta_g * real_r;
                 complex<double> phase(0,ModuleBase::TWO_PI * phase_im);
-                tmp[ixy * nz + iz] /= nxyz;
                 tmp[ixy * nz + iz] *= exp(phase);
             }
         }
diff --git a/source/module_pw/test/test1-2f.cpp b/source/module_pw/test/test1-2f.cpp
@@ -87,7 +87,6 @@ TEST_F(PWTEST,test1_2f)
                 ModuleBase::Vector3<double> real_r(ix, iy, iz);
                 double phase_im = -delta_g * real_r;
                 complex<double> phase(0,ModuleBase::TWO_PI * phase_im);
-                tmp[ixy * nz + iz] /= nxyz;
                 tmp[ixy * nz + iz] *= exp(phase);
             }
         }
diff --git a/source/module_pw/test/test1-3.cpp b/source/module_pw/test/test1-3.cpp
@@ -89,7 +89,6 @@ TEST_F(PWTEST,test1_3)
                 ModuleBase::Vector3<double> real_r(ix, iy, iz);
                 double phase_im = -delta_g * real_r;
                 complex<double> phase(0,ModuleBase::TWO_PI * phase_im);
-                tmp[ixy * nz + iz] /= nxyz;
                 tmp[ixy * nz + iz] *= exp(phase);
             }
         }
diff --git a/source/module_pw/test/test1-3f.cpp b/source/module_pw/test/test1-3f.cpp
@@ -89,7 +89,6 @@ TEST_F(PWTEST,test1_3f)
                 ModuleBase::Vector3<double> real_r(ix, iy, iz);
                 double phase_im = -delta_g * real_r;
                 complex<double> phase(0,ModuleBase::TWO_PI * phase_im);
-                tmp[ixy * nz + iz] /= nxyz;
                 tmp[ixy * nz + iz] *= exp(phase);
             }
         }
@@ -112,7 +111,7 @@ TEST_F(PWTEST,test1_3f)
     {
         for(int iz = 0 ; iz < nplane ; ++iz)
         {
-            EXPECT_NEAR(tmp[ixy * nz + startiz + iz].real(),rhor[ixy*nplane+iz],1e-6);
+            EXPECT_NEAR(tmp[ixy * nz + startiz + iz].real(),rhor[ixy*nplane+iz],1e-4);
         }
     }
     
diff --git a/source/module_pw/test/test1-4.cpp b/source/module_pw/test/test1-4.cpp
@@ -98,7 +98,6 @@ TEST_F(PWTEST,test1_4)
                     ModuleBase::Vector3<double> real_r(ix, iy, iz);
                     double phase_im = -delta_g * real_r;
                     complex<double> phase(0,ModuleBase::TWO_PI * phase_im);
-                    tmp[ixy * nz + iz] /= nxyz;
                     tmp[ixy * nz + iz] *= exp(phase);
                 }
             }
diff --git a/source/module_pw/test/test1-4f.cpp b/source/module_pw/test/test1-4f.cpp
@@ -98,7 +98,6 @@ TEST_F(PWTEST,test1_4f)
                     ModuleBase::Vector3<double> real_r(ix, iy, iz);
                     double phase_im = -delta_g * real_r;
                     complex<double> phase(0,ModuleBase::TWO_PI * phase_im);
-                    tmp[ixy * nz + iz] /= nxyz;
                     tmp[ixy * nz + iz] *= exp(phase);
                 }
             }
@@ -119,8 +118,8 @@ TEST_F(PWTEST,test1_4f)
         {
             for(int iz = 0 ; iz < nplane ; ++iz)
             {
-                EXPECT_NEAR(tmp[ixy * nz + startiz + iz].real(),rhor[ixy*nplane+iz].real(),1e-6);
-                EXPECT_NEAR(tmp[ixy * nz + startiz + iz].imag(),rhor[ixy*nplane+iz].imag(),1e-6);
+                EXPECT_NEAR(tmp[ixy * nz + startiz + iz].real(),rhor[ixy*nplane+iz].real(),1e-4);
+                EXPECT_NEAR(tmp[ixy * nz + startiz + iz].imag(),rhor[ixy*nplane+iz].imag(),1e-4);
             }
         }
 
diff --git a/source/module_pw/test/test2-2.cpp b/source/module_pw/test/test2-2.cpp
@@ -88,7 +88,6 @@ TEST_F(PWTEST,test2_2)
                 ModuleBase::Vector3<double> real_r(ix, iy, iz);
                 double phase_im = -delta_g * real_r;
                 complex<double> phase(0,ModuleBase::TWO_PI * phase_im);
-                tmp[ixy * nz + iz] /= nxyz;
                 tmp[ixy * nz + iz] *= exp(phase);
             }
         }
diff --git a/source/module_pw/test/test2-3.cpp b/source/module_pw/test/test2-3.cpp
@@ -89,7 +89,6 @@ TEST_F(PWTEST,test2_3)
                 ModuleBase::Vector3<double> real_r(ix, iy, iz);
                 double phase_im = -delta_g * real_r;
                 complex<double> phase(0,ModuleBase::TWO_PI * phase_im);
-                tmp[ixy * nz + iz] /= nxyz;
                 tmp[ixy * nz + iz] *= exp(phase);
             }
         }

Original file line number	Diff line number	Diff line change
`@@ -35,7 +35,7 @@ void PW_Basis:: real2recip(std::complex<double> * in, std::complex<double> * out`
`35`	`35`
`36`	`36`	`for(int ig = 0 ; ig < this->npw ; ++ig)`
`37`	`37`	`{`
`38`		`- out[ig] = this->ft.aux1[this->ig2isz[ig]];`
	`38`	`+ out[ig] = this->ft.aux1[this->ig2isz[ig]] / double(this->bignxyz);`
`39`	`39`	`}`
`40`	`40`	`return;`
`41`	`41`	`ModuleBase::timer::tick("PW_Basis", "real2recip");`
`@@ -72,7 +72,7 @@ void PW_Basis:: real2recip(double * in, std::complex<double> * out)`
`72`	`72`
`73`	`73`	`for(int ig = 0 ; ig < this->npw ; ++ig)`
`74`	`74`	`{`
`75`		`- out[ig] = this->ft.aux1[this->ig2isz[ig]];`
	`75`	`+ out[ig] = this->ft.aux1[this->ig2isz[ig]] / double(this->bignxyz);`
`76`	`76`	`}`
`77`	`77`	`ModuleBase::timer::tick("PW_Basis", "real2recip_gamma_only");`
`78`	`78`	`return;`
`@@ -105,7 +105,7 @@ void PW_Basis:: recip2real(std::complex<double> * in, std::complex<double> * out`
`105`	`105`
`106`	`106`	`for(int ir = 0 ; ir < this->nrxx ; ++ir)`
`107`	`107`	`{`
`108`		`- out[ir] = this->ft.aux1[ir] / double(this->bignxyz);`
	`108`	`+ out[ir] = this->ft.aux1[ir];`
`109`	`109`	`}`
`110`	`110`	`ModuleBase::timer::tick("PW_Basis", "recip2real");`
`111`	`111`
`@@ -144,7 +144,7 @@ void PW_Basis:: recip2real(std::complex<double> * in, double * out)`
`144`	`144`	`{`
`145`	`145`	`for(int ipy = 0 ; ipy < npy ; ++ipy)`
`146`	`146`	`{`
`147`		`- out[ixnpy + ipy] = this->ft.r_rspace[ixnpy*2 + ipy] / double(this->bignxyz);`
	`147`	`+ out[ixnpy + ipy] = this->ft.r_rspace[ixnpy*2 + ipy];`
`148`	`148`	`}`
`149`	`149`	`}`
`150`	`150`	`ModuleBase::timer::tick("PW_Basis", "recip2real_gamma_only");`
`@@ -167,7 +167,7 @@ void PW_Basis:: real2recip(std::complex<float> * in, std::complex<float> * out)`
`167`	`167`
`168`	`168`	`for(int ig = 0 ; ig < this->npw ; ++ig)`
`169`	`169`	`{`
`170`		`- out[ig] = this->ft.auxf1[this->ig2isz[ig]];`
	`170`	`+ out[ig] = this->ft.auxf1[this->ig2isz[ig]] / float(this->bignxyz);`
`171`	`171`	`}`
`172`	`172`	`return;`
`173`	`173`	`}`
`@@ -192,7 +192,7 @@ void PW_Basis:: real2recip(float * in, std::complex<float> * out)`
`192`	`192`
`193`	`193`	`for(int ig = 0 ; ig < this->npw ; ++ig)`
`194`	`194`	`{`
`195`		`- out[ig] = this->ft.auxf1[this->ig2isz[ig]];`
	`195`	`+ out[ig] = this->ft.auxf1[this->ig2isz[ig]] / float(this->bignxyz);`
`196`	`196`	`}`
`197`	`197`	`return;`
`198`	`198`	`}`
`@@ -214,7 +214,7 @@ void PW_Basis:: recip2real(std::complex<float> * in, std::complex<float> * out)`
`214`	`214`
`215`	`215`	`for(int ir = 0 ; ir < this->nrxx ; ++ir)`
`216`	`216`	`{`
`217`		`- out[ir] = this->ft.auxf1[ir] / float(this->bignxyz);`
	`217`	`+ out[ir] = this->ft.auxf1[ir];`
`218`	`218`	`}`
`219`	`219`
`220`	`220`	`return;`
`@@ -240,7 +240,7 @@ void PW_Basis:: recip2real(std::complex<float> * in, float * out)`
`240`	`240`	`{`
`241`	`241`	`for(int ipy = 0 ; ipy < npy ; ++ipy)`
`242`	`242`	`{`
`243`		`- out[ixnpy + ipy] = this->ft.rf_rspace[ixnpy*2 + ipy] / float(this->bignxyz);`
	`243`	`+ out[ixnpy + ipy] = this->ft.rf_rspace[ixnpy*2 + ipy];`
`244`	`244`	`}`
`245`	`245`	`}`
`246`	`246`	`return;`
Original file line number	Diff line number	Diff line change
`@@ -31,7 +31,7 @@ void PW_Basis_K:: real2recip(std::complex<double> * in, std::complex<double> * o`
`31`	`31`
`32`	`32`	`for(int igl = 0 ; igl < this->npwk[ik] ; ++igl)`
`33`	`33`	`{`
`34`		`- out[igl] = this->ft.aux1[this->igl2isz_k[igl+ik*this->npwk_max]];`
	`34`	`+ out[igl] = this->ft.aux1[this->igl2isz_k[igl+ik*this->npwk_max]] / double(this->bignxyz);`
`35`	`35`	`}`
`36`	`36`	`return;`
`37`	`37`	`ModuleBase::timer::tick("PW_Basis_K", "real2recip");`
`@@ -68,7 +68,7 @@ void PW_Basis_K:: real2recip(double * in, std::complex<double> * out, int ik)`
`68`	`68`
`69`	`69`	`for(int igl = 0 ; igl < this->npwk[ik] ; ++igl)`
`70`	`70`	`{`
`71`		`- out[igl] = this->ft.aux1[this->igl2isz_k[igl+ik*this->npwk_max]];`
	`71`	`+ out[igl] = this->ft.aux1[this->igl2isz_k[igl+ik*this->npwk_max]] / double(this->bignxyz);`
`72`	`72`	`}`
`73`	`73`	`ModuleBase::timer::tick("PW_Basis_K", "real2recip_gamma_only");`
`74`	`74`	`return;`
`@@ -97,7 +97,7 @@ void PW_Basis_K:: recip2real(std::complex<double> * in, std::complex<double> * o`
`97`	`97`
`98`	`98`	`for(int ir = 0 ; ir < this->nrxx ; ++ir)`
`99`	`99`	`{`
`100`		`- out[ir] = this->ft.aux1[ir] / double(this->bignxyz);`
	`100`	`+ out[ir] = this->ft.aux1[ir];`
`101`	`101`	`}`
`102`	`102`	`ModuleBase::timer::tick("PW_Basis_K", "recip2real");`
`103`	`103`
`@@ -136,7 +136,7 @@ void PW_Basis_K:: recip2real(std::complex<double> * in, double * out, int ik)`
`136`	`136`	`{`
`137`	`137`	`for(int ipy = 0 ; ipy < npy ; ++ipy)`
`138`	`138`	`{`
`139`		`- out[ixnpy + ipy] = this->ft.r_rspace[ixnpy*2 + ipy] / double(this->bignxyz);`
	`139`	`+ out[ixnpy + ipy] = this->ft.r_rspace[ixnpy*2 + ipy];`
`140`	`140`	`}`
`141`	`141`	`}`
`142`	`142`	`ModuleBase::timer::tick("PW_Basis_K", "recip2real_gamma_only");`
`@@ -159,7 +159,7 @@ void PW_Basis_K:: real2recip(std::complex<float> * in, std::complex<float> * out`
`159`	`159`
`160`	`160`	`for(int igl = 0 ; igl < this->npwk[ik] ; ++igl)`
`161`	`161`	`{`
`162`		`- out[igl] = this->ft.auxf1[this->igl2isz_k[igl+ik*this->npwk_max]];`
	`162`	`+ out[igl] = this->ft.auxf1[this->igl2isz_k[igl+ik*this->npwk_max]] / float(this->bignxyz);`
`163`	`163`	`}`
`164`	`164`	`return;`
`165`	`165`	`}`
`@@ -184,7 +184,7 @@ void PW_Basis_K:: real2recip(float * in, std::complex<float> * out, int ik)`
`184`	`184`
`185`	`185`	`for(int igl = 0 ; igl < this->npwk[ik] ; ++igl)`
`186`	`186`	`{`
`187`		`- out[igl] = this->ft.auxf1[this->igl2isz_k[igl+ik*this->npwk_max]];`
	`187`	`+ out[igl] = this->ft.auxf1[this->igl2isz_k[igl+ik*this->npwk_max]] / float(this->bignxyz);`
`188`	`188`	`}`
`189`	`189`	`return;`
`190`	`190`	`}`
`@@ -206,7 +206,7 @@ void PW_Basis_K:: recip2real(std::complex<float> * in, std::complex<float> * out`
`206`	`206`
`207`	`207`	`for(int ir = 0 ; ir < this->nrxx ; ++ir)`
`208`	`208`	`{`
`209`		`- out[ir] = this->ft.auxf1[ir] / float(this->bignxyz);`
	`209`	`+ out[ir] = this->ft.auxf1[ir];`
`210`	`210`	`}`
`211`	`211`
`212`	`212`	`return;`
`@@ -232,7 +232,7 @@ void PW_Basis_K:: recip2real(std::complex<float> * in, float * out, int ik)`
`232`	`232`	`{`
`233`	`233`	`for(int ipy = 0 ; ipy < npy ; ++ipy)`
`234`	`234`	`{`
`235`		`- out[ixnpy + ipy] = this->ft.rf_rspace[ixnpy*2 + ipy] / float(this->bignxyz);`
	`235`	`+ out[ixnpy + ipy] = this->ft.rf_rspace[ixnpy*2 + ipy];`
`236`	`236`	`}`
`237`	`237`	`}`
`238`	`238`	`return;`
Original file line number	Diff line number	Diff line change
`@@ -87,7 +87,6 @@ TEST_F(PWTEST,test1_2)`
`87`	`87`	`ModuleBase::Vector3<double> real_r(ix, iy, iz);`
`88`	`88`	`double phase_im = -delta_g * real_r;`
`89`	`89`	`complex<double> phase(0,ModuleBase::TWO_PI * phase_im);`
`90`		`- tmp[ixy * nz + iz] /= nxyz;`
`91`	`90`	`tmp[ixy * nz + iz] *= exp(phase);`
`92`	`91`	`}`
`93`	`92`	`}`
Original file line number	Diff line number	Diff line change
`@@ -89,7 +89,6 @@ TEST_F(PWTEST,test1_3)`
`89`	`89`	`ModuleBase::Vector3<double> real_r(ix, iy, iz);`
`90`	`90`	`double phase_im = -delta_g * real_r;`
`91`	`91`	`complex<double> phase(0,ModuleBase::TWO_PI * phase_im);`
`92`		`- tmp[ixy * nz + iz] /= nxyz;`
`93`	`92`	`tmp[ixy * nz + iz] *= exp(phase);`
`94`	`93`	`}`
`95`	`94`	`}`
Original file line number	Diff line number	Diff line change
`@@ -98,7 +98,6 @@ TEST_F(PWTEST,test1_4)`
`98`	`98`	`ModuleBase::Vector3<double> real_r(ix, iy, iz);`
`99`	`99`	`double phase_im = -delta_g * real_r;`
`100`	`100`	`complex<double> phase(0,ModuleBase::TWO_PI * phase_im);`
`101`		`- tmp[ixy * nz + iz] /= nxyz;`
`102`	`101`	`tmp[ixy * nz + iz] *= exp(phase);`
`103`	`102`	`}`
`104`	`103`	`}`
Original file line number	Diff line number	Diff line change
`@@ -88,7 +88,6 @@ TEST_F(PWTEST,test2_2)`
`88`	`88`	`ModuleBase::Vector3<double> real_r(ix, iy, iz);`
`89`	`89`	`double phase_im = -delta_g * real_r;`
`90`	`90`	`complex<double> phase(0,ModuleBase::TWO_PI * phase_im);`
`91`		`- tmp[ixy * nz + iz] /= nxyz;`
`92`	`91`	`tmp[ixy * nz + iz] *= exp(phase);`
`93`	`92`	`}`
`94`	`93`	`}`
Original file line number	Diff line number	Diff line change
`@@ -89,7 +89,6 @@ TEST_F(PWTEST,test2_3)`
`89`	`89`	`ModuleBase::Vector3<double> real_r(ix, iy, iz);`
`90`	`90`	`double phase_im = -delta_g * real_r;`
`91`	`91`	`complex<double> phase(0,ModuleBase::TWO_PI * phase_im);`
`92`		`- tmp[ixy * nz + iz] /= nxyz;`
`93`	`92`	`tmp[ixy * nz + iz] *= exp(phase);`
`94`	`93`	`}`
`95`	`94`	`}`