Changeset 4bb871a for util/src

Timestamp:

Feb 2, 2010, 12:22:06 PM (16 years ago)

Author:

Tillmann Crueger <crueger@…>

Children:

45cd358

Parents:

adcdf8 (diff), 6d0fcaa (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.

Message:

Merge commit 'heber/Analysis_PairCorrelation' into MenuRefactoring

Location:

util/src

Files:

• NanoCreator.c (modified) (11 diffs)
• NanoCreator.h (modified) (1 diff)

util/src/NanoCreator.c

@@ -520 +520 @@
 {
   int c;
+  if ((a == 0) && (b == 0))
+    return 1;
+  else if (a == 0)
+    return b;
+  else if (b == 0)
+    return a;
   do {
     c = a % b;        /* Rest of integer divison */
@@ -672 +678 @@
   double *x = NULL, *coord = NULL, *tempvector = NULL, *offset = NULL;
   //double *cog = NULL, *cog_projected = NULL;
-  char stage[6];
+  char *stage = NULL;
   int axis[NDIM] = {0,1,2}, chiral[2] = {1,1}, factors[NDIM] = {1,1,1};
   char name[255], line[255], input = 'n';
   char *CellBuffer = NULL, *SheetBuffer = NULL, *TubeBuffer = NULL, *bufptr = NULL;
   double *randomness = NULL, percentage; // array with percentages for presence in sheet and beyond
-  int i,j, ggT;
+  int i,j, ggT, ggTsmall;
   int length;
 
@@ -688 +694 @@
         } else {
                 // store variables
-                filename = argv[1];
-                strncpy(stage, argv[2], 5);
+                filename = Malloc(sizeof(char)*(strlen(argv[1])+1), "Main: filename");
+    stage = Malloc(sizeof(char)*(strlen(argv[2])+1), "Main: stage");
+                strcpy(filename, argv[1]);
+                strcpy(stage, argv[2]);
     fprintf(stdout, "I will begin at stage %s.\n", stage);
 
                 // build filenames
+    length = strlen(filename);
     char *ptr = strrchr(filename, '.');
     if (ptr == NULL) {
       ptr = filename;
     } else {
-      length = strlen(filename);
       if (ptr != NULL) {
         length = ((int)(ptr-filename) >= length-5) ? (int)(ptr-filename) : length;
@@ -841 +849 @@
 
   int numbersheet, biggestdiameter, sheetnr[NDIM], tmp, seed;
+  double dfactors[2];
   double x1,x2,x3, angle;
   char flag = 'n';
@@ -874 +883 @@
 
     do {
-      fprintf(stdout, "\nNow specify the two natural numbers (m n) defining the chiral angle, \nif the result is crap, try flipping to (m,n): ");
+      fprintf(stdout, "\nNow specify the two natural numbers (n m) defining the chiral angle, \nif the result is crap, try flipping to (n,m): ");
       fscanf(stdin, "%d %d", &chiral[0], &chiral[1]);
-      ggT = GCD(2*chiral[1]+chiral[0],2*chiral[0]+chiral[1]);
+      ggTsmall = GCD(chiral[0],chiral[1]);
+      ggT = GCD(2*chiral[0]+chiral[1],2*chiral[1]+chiral[0]);
+      fprintf(stdout, "Greatest Common Denominator of (n, m) is %d\n", ggTsmall);
       fprintf(stdout, "Greatest Common Denominator of (2n+m, 2m+n) is %d\n", ggT);
       fprintf(stdout, "chiral0: %d\tchiral1: %d\n", chiral[0], chiral[1]);
       for (i=0;i<NDIM;i++) {
-        Tubevector[axis[0]][i] = (double)chiral[0] * Vector[axis[0]][i] + (double)chiral[1] * Vector[axis[1]][i];
-        //Tubevector[axis[0]][i] = chiral[0] * Vector[axis[0]][i] + chiral[1] * Vector[axis[1]][i];
+//        Tubevector[axis[0]][i] = (double)chiral[0] * Vector[axis[0]][i] + (double)chiral[1] * Vector[axis[1]][i];
+        Tubevector[axis[0]][i] = chiral[0] * Vector[axis[0]][i] + chiral[1] * Vector[axis[1]][i];
         //Tubevector[axis[0]][i] = (2.*chiral[0]+chiral[1])/(double)ggT * Vector[axis[0]][i] + (-chiral[0]-2.*chiral[1])/(double)ggT * Vector[axis[1]][i];
         //Tubevector[axis[1]][i] = -chiral[1] * Vector[axis[0]][i] + chiral[0] * Vector[axis[1]][i];
-        Tubevector[axis[1]][i] = (double)chiral[0] * OrthoVector[axis[0]][i] - (double)chiral[1] * OrthoVector[axis[1]][i];
+//        Tubevector[axis[1]][i] = (double)chiral[0] * OrthoVector[axis[0]][i] - (double)chiral[1] * OrthoVector[axis[1]][i];
+        Tubevector[axis[1]][i] = (2.*chiral[0]+chiral[1])/(double)ggT * Vector[axis[0]][i] + (-chiral[0]-2.*chiral[1])/(double)ggT * Vector[axis[1]][i];
         //Tubevector[axis[1]][i] = (-chiral[0]-2.*chiral[1])/(double)ggT * Vector[axis[0]][i] + (2.*chiral[0]+chiral[1])/(double)ggT * Vector[axis[1]][i];
 //        fprintf(stderr, "Tubevector[axis[0]][i] = (double)chiral[0] * Vector[axis[0]][i] + (double)chiral[1] * Vector[axis[1]][i]\n = %lg * %lg + %lg * %lg = %lg + %lg = %lg\n\n",
@@ -924 +936 @@
         x1 = gsl_vector_get(u,0)*(double)i;
         x2 = gsl_vector_get(u,1)*(double)i;
-        x3 =
         fprintf(stdout, "%d: %d\t%d vs. %lg\t%lg\n",i, ((int)(x1+x1/fabs(x1)*.5)), ((int)(x2+x2/fabs(x2)*.5)), (x1), (x2));
         if (( fabs( ((int)(x1+x1/fabs(x1)*.5)) - (x1) ) < 1e-6) && ( fabs( ((int)(x2+x2/fabs(x2)*.5)) - (x2) ) < 1e-6 )) {
@@ -931 +942 @@
         }
       }
-      fprintf(stdout, "(c,d) = (%lg,%lg)\n",gsl_vector_get(u,0), gsl_vector_get(u,1));
-
-      // get length
+      x1 = fabs(gsl_vector_get(u,0));
+      x2 = fabs(gsl_vector_get(u,1));
+      fprintf(stdout, "(c,d) = (%d,%d).\n",((int)(x1+x1/fabs(x1)*.5)), ((int)(x2+x2/fabs(x2)*.5)));
+      j = GCD(((int)(x1+x1/fabs(x1)*.5)), ((int)(x2+x2/fabs(x2)*.5)));
+      fprintf(stdout, "GCD(%d,%d) = %i.\n", ((int)(x1+x1/fabs(x1)*.5)), ((int)(x2+x2/fabs(x2)*.5)), j);
+      for (i=0;i<2;i++)
+        gsl_vector_set(u,i, gsl_vector_get(u,i)/(double)j);
+        // get length
       double x[NDIM];
       for (i=0;i<NDIM;i++)
         x[i] = gsl_vector_get(u,0) * Vector[axis[0]][i] + gsl_vector_get(u,1) * Vector[axis[1]][i];
       angle = Norm(x)/Norm(Tubevector[axis[1]]) ;//ScalarProduct(x,Tubevector[axis[1]])/Norm(Tubevector[axis[1]]);
-      fprintf(stdout, "angle is %lg\n", angle);
+      fprintf(stdout, "Angle is %lg.\n", angle);
       for (i=0;i<NDIM;i++) {
         Tubevector[axis[1]][i] = gsl_vector_get(u,0) * Vector[axis[0]][i] + gsl_vector_get(u,1) * Vector[axis[1]][i];
@@ -994 +1010 @@
         );
       fprintf(stdout, "\nGive integer factors for length and radius of tube (multiple of %d suggested) : ", ggT);
-      fscanf(stdin, "%d %d", &factors[1], &factors[0]);
+      fscanf(stdin, "%lg %lg", &dfactors[1], &dfactors[0]);
+      for(i=0;i<2;i++)
+        factors[i] = (int)(dfactors[i]);
+      fprintf(stdout, "\nI got %d and %d.\n", factors[1], factors[0]);
       fprintf(stdout, "\nThe chiral angle then is %5.5f degrees with tube radius %5.5f A and length %5.5f A, i.e. torus radius of %5.5f A.\n",
         acos(Projection(Vector[axis[0]], Tubevector[axis[0]]))/M_PI*180.,
@@ -1323 +1342 @@
               fprintf(SheetFile, "%s\t%5.5lg\t%5.5lg\t%5.5lg\n", atombuffer[nr].name, coord[0], coord[1], coord[2]);
               // rotate to align the sheet in xy plane
-              x1 = coord[0]*cos(-angle) + coord[1] * sin(-angle);
-              x2 = coord[0]*(-sin(-angle)) + coord[1] * cos(-angle);
+              x1 = coord[0]*cos(angle) + coord[1] * sin(angle);
+              x2 = coord[0]*(-sin(angle)) + coord[1] * cos(angle);
               x3 = coord[2];
               fprintf(SheetFileAligned, "%s\t%5.5lg\t%5.5lg\t%5.5lg\n", atombuffer[nr].name, x1, x2, x3);
@@ -1424 +1443 @@
       fprintf(TubeFile, "%s\t%lg\t%lg\t%lg\n", name, atom_transformed[0], atom_transformed[1], atom_transformed[2]);
       // rotate and flip to align tube in z-direction
-      x1 = atom_transformed[0]*cos(-angle) + atom_transformed[1] * sin(-angle);
-      x2 = atom_transformed[0]*(-sin(-angle)) + atom_transformed[1] * cos(-angle);
+      x1 = atom_transformed[0]*cos(angle) + atom_transformed[1] * sin(angle);
+      x2 = atom_transformed[0]*(-sin(angle)) + atom_transformed[1] * cos(angle);
       x3 = atom_transformed[2];
       fprintf(TubeFileAligned, "%s\t%lg\t%lg\t%lg\n", name, x3, x2, x1);  // order so that symmetry is along z axis
@@ -1537 +1556 @@
   if (TubeBuffer != NULL) Free(TubeBuffer, "Main: end of stages - TubeBuffer");
 
-  Free(CellFilename, "Main: end of stafes - CellFilename");
-  Free(SheetFilename, "Main: end of stafes - CellFilename");
-  Free(TubeFilename, "Main: end of stafes - CellFilename");
-  Free(TorusFilename, "Main: end of stafes - CellFilename");
-  Free(SheetFilenameAligned, "Main: end of stafes - CellFilename");
-  Free(TubeFilenameAligned, "Main: end of stafes - CellFilename");
+  Free(CellFilename, "Main: end of stages - CellFilename");
+  Free(SheetFilename, "Main: end of stages - CellFilename");
+  Free(TubeFilename, "Main: end of stages - CellFilename");
+  Free(TorusFilename, "Main: end of stages - CellFilename");
+  Free(SheetFilenameAligned, "Main: end of stages - CellFilename");
+  Free(TubeFilenameAligned, "Main: end of stages - CellFilename");
+  Free(filename, "Main: end of stages - filename");
+  Free(stage, "Main: end of stages - stage");
 
   // exit

util/src/NanoCreator.h

@@ -2 +2 @@
 #define NANOCREATOR_H_
 
-#define MYEPSILON 1e-13
+#define MYEPSILON 2e-1
 
 struct Atoms