Changeset f5586e for pcp/src/init.c

Timestamp:

Apr 22, 2008, 8:43:24 AM (17 years ago)

Author:

Frederik Heber <heber@…>

Children:

90c027

Parents:

cc46b0

Message:

RealBasisQ[] removed

RealBasisQ[] is again replaced by sqrt(RealBasisSQ[]) as we are going to switch to a pure matrix transformation for calculating whether points are still within the periodic cell and so forth instead of "hoping" the simulation box is rectangular.

File:

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

pcp/src/init.c ¶

@@ -148 +148 @@
   for (i=0; i < NDIM_NDIM; i++) Lat->ReciBasis[i] *= factor;            // SM(Lat->ReciBasis, factor, NDIM_NDIM);
   for (i=0; i < NDIM; i++) {
+    Lat->RealBasisSQ[i] = RNORMSQ3(&(Lat->RealBasis[i*NDIM]));
     Lat->ReciBasisSQ[i] = RNORMSQ3(&(Lat->ReciBasis[i*NDIM]));
-    Lat->RealBasisSQ[i] = RNORMSQ3(&(Lat->RealBasis[i*NDIM]));
-    Lat->RealBasisQ[i] = sqrt(Lat->RealBasisSQ[i]);
   }
         // Prepares LatticeLevels
@@ -279 +278 @@
       if (i%NDIM == NDIM-1) fprintf(stderr,"\n"); 
     }
-    fprintf(stderr,"(%i) Lat->RealbasisQ =  ",P->Par.me);
-    for (i=0; i < NDIM; i++)
-      fprintf(stderr,"%lg ",Lat->RealBasisQ[i]);
-    fprintf(stderr,"\n");
   }