Changeset 1074cb

Timestamp:

Feb 10, 2010, 2:47:24 PM (16 years ago)

Author:

Frederik Heber <heber@…>

Children:

850e50

Parents:

2bc06b

Message:

NanoCreator - radius in Tube section is now outside of loop.

Signed-off-by: Frederik Heber <heber@…>

File:

• util/src/NanoCreator.c (modified) (2 diffs)

util/src/NanoCreator.c

@@ -1422 +1422 @@
       halfaxis[i] = factors[0]*Norm(Tubevector[axis[0]])/Norm(Tubevector[i]);
 
-    double arg, radius;
+    double arg;
+    const double radius = 1./(2.*M_PI);  // is length of sheet in units of axis vector, divide by pi to get radius (from circumference)
     for (i=0;i<numbersheet;i++) {
       // scan next atom
@@ -1435 +1436 @@
       // roll up (project (x,y,z) on cylindrical coordinates (radius,arg,z))
       arg = 2.*M_PI*x[axis[0]]/(factors[0]) - M_PI;  // is angle
-      radius = 1./(2.*M_PI);  // is length of sheet in units of axis vector, divide by pi to get radius (from circumference)
       // fprintf(stdout, "arg: %5.2f (c%2.2f,s%2.2f)\t",$arg, cos($arg), sin($arg));
       x[axis[0]] = cos(arg)*halfaxis[axis[0]]*(radius+x[axis[2]]/halfaxis[axis[2]]); // as both vectors are not normalized additional betrag has to be taken into account!