Changeset f91abc for pcp

Timestamp:

Apr 22, 2008, 11:44:28 AM (17 years ago)

Author:

Frederik Heber <heber@…>

Children:

519b83

Parents:

1d77026

Message:

CalculatePerturbationOperator_RxP(): arrays for N, NUp and Wcentre transformed to pointers referring to orginal values

File:

• pcp/src/perturbed.c (modified) (1 diff)

pcp/src/perturbed.c

@@ -1206 +1206 @@
   fftw_complex *PsiC = Dens0->DensityCArray[ActualPsiDensity];
   fftw_real *PsiCR = (fftw_real *) PsiC;
-  double x[NDIM], X[NDIM], fac[NDIM], Wcentre[NDIM];
+  double x[NDIM], X[NDIM], fac[NDIM], *Wcentre;
   int n[NDIM], n0, g, Index, iS, i0; //pos, 
-  int N[NDIM], NUp[NDIM];
+  const int *N, *NUp;
   const int N0 = LevS->Plan0.plan->local_nx;
-  N[0] = LevS->Plan0.plan->N[0];
-  N[1] = LevS->Plan0.plan->N[1];
-  N[2] = LevS->Plan0.plan->N[2];
-  NUp[0] = LevS->NUp[0]; 
-  NUp[1] = LevS->NUp[1]; 
-  NUp[2] = LevS->NUp[2]; 
-  Wcentre[0] = Lat->Psi.AddData[phi0nr].WannierCentre[0];
-  Wcentre[1] = Lat->Psi.AddData[phi0nr].WannierCentre[1];
-  Wcentre[2] = Lat->Psi.AddData[phi0nr].WannierCentre[2];
+  N = LevS->Plan0.plan->N;
+  NUp = LevS->NUp; 
+  Wcentre = Lat->Psi.AddData[phi0nr].WannierCentre;
   // init pointers and values
   const int myPE = P->Par.me_comm_ST_Psi;