delete -march-native in Makefile.gnu or else we can not use valgrind

Author: Qianruipku

source/module_pw/pw_gatherscatter.cpp

@@ -23,7 +23,7 @@ void PW_Basis:: gatherp_scatters(complex<double> *in, complex<double> *out)
         return;
     }
 #ifdef __MPI
-    std::complex<double> * tmp;
+    std::complex<double> * tmp = NULL;
     if(this->poolrank == 0) tmp = new std::complex<double> [this->nz * this->nstot];
 
     //gather planes of different processors
@@ -40,7 +40,7 @@ void PW_Basis:: gatherp_scatters(complex<double> *in, complex<double> *out)
     MPI_Scatterv(tmp, this->nstnz_per, this->startnsz_per,mpicomplex,out,
                     this->nstnz,mpicomplex,0, POOL_WORLD); 
 
-    if(this->poolrank == 0) delete[] tmp;
+    if(tmp!=NULL) delete[] tmp;
 #endif
     return;
 }
@@ -107,7 +107,7 @@ void PW_Basis:: gathers_scatterp(complex<double> *in, complex<double> *out)
     }
 #ifdef __MPI
     if(this->poolnproc == 1) return;
-    std::complex<double> * tmp;
+    std::complex<double> * tmp = NULL;
     if(this->poolrank == 0) tmp = new std::complex<double> [this->nz * this->nstot];
 
     //scatter sticks to different processors
@@ -121,7 +121,7 @@ void PW_Basis:: gathers_scatterp(complex<double> *in, complex<double> *out)
         MPI_Scatterv(&tmp[istot*this->nz], this->numz,this->startz, mpicomplex, &out[ixy*this->nplane], 
                     this->nplane,mpicomplex,0,POOL_WORLD);
     }
-    if(this->poolrank == 0) delete[] tmp;
+    if(tmp!=NULL) delete[] tmp;
 #endif
     return;
 }

source/module_pw/unittest/Makefile.gnu

@@ -69,7 +69,7 @@ FFTW_LIB         = -L${FFTW_LIB_DIR} -lfftw3 -Wl,-rpath=${FFTW_LIB_DIR}
 # CUDA_LIB			= -L${CUDA_LIB_DIR} -lcufft -lcublas -lcudart
 
 LIBS = ${FFTW_LIB} ${CUDA_LIB}
-OPTS = -I${FFTW_INCLUDE_DIR} ${HONG} -Ofast -march=native -std=c++11 -Wall -g 
+OPTS = -I${FFTW_INCLUDE_DIR} ${HONG} -Ofast -std=c++11 -Wall -g 
 #==========================
 # MAKING OPTIONS
 #==========================

source/module_pw/unittest/test2.cpp

@@ -55,7 +55,7 @@ int main(int argc,char **argv)
     GT = latvec.Inverse();
 	G  = GT.Transpose();
 	GGT = G * GT;
-    complex<double> *tmp;
+    complex<double> *tmp = NULL;
     if(myrank == 0)
     {
         tmp = new complex<double> [nx*ny*nz];
@@ -157,6 +157,6 @@ int main(int argc,char **argv)
     MPI_Barrier(MPI_COMM_WORLD);     
     delete [] rhog;
     delete [] rhor;
-    if(myrank == 0) delete [] tmp; 
+    if(tmp!=NULL) delete []tmp; 
     return 0;
 }

source/module_pw/unittest/test3.cpp

@@ -55,7 +55,7 @@ int main(int argc,char **argv)
     GT = latvec.Inverse();
 	G  = GT.Transpose();
 	GGT = G * GT;
-    complex<double> *tmp;
+    complex<double> *tmp = NULL;
     if(myrank == 0)
     {
         tmp = new complex<double> [nx*ny*nz];
@@ -155,6 +155,6 @@ int main(int argc,char **argv)
     MPI_Barrier(MPI_COMM_WORLD);
     delete [] rhog;
     delete [] rhor;
-    if(myrank == 0) delete [] tmp; 
+    if(tmp!=NULL) delete []tmp;
     return 0;
 }