Changeset 88e890 for pcp/src/init.c

Timestamp:

Apr 22, 2008, 8:52:36 AM (18 years ago)

Author:

Frederik Heber <heber@…>

Children:

53b5b6

Parents:

90c027

Message:

RealBasisCenter[] is present and initialised.

RealBasisCenter[] is part of Lattice structure (fixed NDIM_NDIM) and initialised during InitLattice() by vector (.5, .5, .5) and MatrixTrafo with RealBasis[].

File:

• pcp/src/init.c (modified) (3 diffs)

pcp/src/init.c

@@ -133 +133 @@
   int i,j,UpDummy;
   int MaxG[NDIM], N[NDIM], NFactor, Lev1N[NDIM];
+  double x[NDIM];
   struct LatticeLevel *Lev0;
   int RL=-1;
@@ -142 +143 @@
     return(1);
   }
+  for (i=0; i < NDIM; i++) // create center point 
+    x[i] = .5;
+  RMat33Vec3(Lat->RealBasisCenter, Lat->RealBasis, x);
+  if (P->Call.out[ValueOut]) fprintf(stderr,"Unit cell center = (%lg, %lg, %lg)\n", Lat->RealBasisCenter[0], Lat->RealBasisCenter[1], Lat->RealBasisCenter[2]);
   // square-rooted and normal reciprocal (transposed) Base
   for (i=0; i < NDIM_NDIM; i++) h[i] = Lat->RealBasis[i];
@@ -148 +153 @@
   for (i=0; i < NDIM_NDIM; i++) Lat->ReciBasis[i] *= factor;            // SM(Lat->ReciBasis, factor, NDIM_NDIM);
   for (i=0; i < NDIM; i++) {
+    Lat->ReciBasisSQ[i] = RNORMSQ3(&(Lat->ReciBasis[i*NDIM]));
     Lat->RealBasisSQ[i] = RNORMSQ3(&(Lat->RealBasis[i*NDIM]));
-    Lat->ReciBasisSQ[i] = RNORMSQ3(&(Lat->ReciBasis[i*NDIM]));
   }
         // Prepares LatticeLevels