Changeset 1856 for trunk/source/optimize_phymir.cpp

Timestamp:

03/17/08 08:31:47 (10 months ago)

Author:

rporter

Message:

trunk/source/cddefines.h - comment out some of the items in enum colliders because they are not used and just waste memory at this point.

trunk/source/helike_einsta.cpp - add M1 transitions from 3^{3S to 2}3S for he-like sequence. May be important in early universe.

trunk/source/nemala.cpp - make sure this is only done once per coreload.

trunk/source/optimize_phymir.cpp - remove some superfluous realnum casts

trunk/source/parse_interp.cpp - BUGFIX at end of parsing start one higher in zeroing out loop

trunk/source/radius_next.cpp - change punch dr printout from %g to %11.3e so numbers always line up in output, easier to scan file this way.

Files:

• trunk/source/optimize_phymir.cpp (modified) (5 diffs)

trunk/source/optimize_phymir.cpp

@@ -401 +401 @@
                         c2[i] += POW2(a2[i][j]/c1[j]);
                 }
-                c2[i] = 1.f/(realnum)sqrt(c2[i]);
+                c2[i] = 1.f/sqrt(c2[i]);
                 xc[i] = xp[i][jmin];
                 xp[i][0] = xc[i];
@@ -422 +422 @@
                 {
                         /* when we make progress, dmax may become bigger */
-                        r2 = (realnum)sqrt(1.f/F1);
+                        r2 = sqrt(1.f/F1);
                 }
                 else
@@ -428 +428 @@
                         /* when there is no progress force dmax smaller, otherwise there
                          * may never be an end */
-                        r2 = (realnum)sqrt(F1);
+                        r2 = sqrt(F1);
                 }
         }
@@ -445 +445 @@
                 c2[i] = c1[i];
         }
-        diff = (realnum)sqrt(diff);
+        diff = sqrt(diff);
         dmax = dold;
         lgFinish = diff <= toler;
@@ -468 +468 @@
                 for( k=0; k < n; k++ )
                         ip += POW2(a[i][k]);
-                ip = (realnum)sqrt(ip);
+                ip = sqrt(ip);
                 for( k=0; k < n; k++ )
                         a[i][k] /= ip;