Changeset e24bbb for pcp/src/output.c

Timestamp:

Apr 21, 2008, 3:33:52 PM (17 years ago)

Author:

Frederik Heber <heber@…>

Children:

48cbc9

Parents:

d8bb59

Message:

PlotRealDensity(): nodes now calculated via ReciBasis-matrix multiplication

We introcude cartesian coordinate in the config file which is then transformed to integer cell (node) coordinate.

File:

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

pcp/src/output.c

@@ -1735 +1735 @@
   int PE, zahl;
   double *buffer;
-  int node = ceil(n_orth/Lat->RealBasisQ[index]*N[index]);
   MPI_Status status;
   int sizes[P->Par.Max_me_comm_ST_Psi], c0, c1;
+  double nodes[NDIM], node[NDIM];
+  
+  for(i=0;i<NDIM;i++) {
+    nodes[i] = (i == index) ? n_orth : 0.;
+    node[i] = 0.;
+  }
+  RMat33Vec3(node, Lat->ReciBasis, nodes);  // transform cartesian coordinates into cell coordinates [0,1]^3
+  for(i=0;i<NDIM;i++) // now N^3 within node range of discrete grid
+    node[i] = (int)(node[i]*N[i]/(2.*PI));
+  fprintf(stderr,"(%i) n_orth %lg, index %i converted to plane offset vector (%lg, %lg, %lg).\n", P->Par.me, n_orth, index, node[0], node[1], node[2]);
   
   switch (index) {
@@ -1750 +1759 @@
       break;
   } 
+  fprintf(stderr,"(%i) buffer size %i\n", P->Par.me, zahl);
   buffer = Malloc(sizeof(double)*zahl,"PlotRealDensity: buffer");
   
@@ -1755 +1765 @@
   c1 = cross(index,1);
   // then for every point on the grid in real space ...
+  i=0;
   for (n0=0;n0<N0;n0++)  // only local points on x axis
     for (n[1]=0;n[1]<N[1];n[1]++)
       for (n[2]=0;n[2]<N[2];n[2]++) {
         n[0]=n0 + N0*myPE; // global relative coordinate: due to partitoning of x-axis in PEPGamma>1 case
-        if (n[index] == node) { // only on the correct plane orthogonal to desired axis and at desired node ...
+        if (n[index] == (int)node[index]) { // only on the correct plane orthogonal to desired axis and at desired node ...
           fac[0] = (double)n[0]/(double)N[0];
           fac[1] = (double)n[1]/(double)N[1];