Changeset 7097

Timestamp:

May 27, 2009, 12:57:54 PM (16 years ago)

Author:

ole

Message:

Played with inlining and profiled, but no speedup was achieved, so I rolled it back and commented.

File:

• anuga_core/source/anuga/shallow_water/shallow_water_ext.c (modified) (4 diffs)

anuga_core/source/anuga/shallow_water/shallow_water_ext.c

@@ -27 +27 @@
 
 // Computational function for rotation
-// FIXME: Perhaps inline this and profile
+// Tried to inline, but no speedup was achieved 27th May 2009 (Ole)
+// static inline int _rotate(double *q, double n1, double n2) {
 int _rotate(double *q, double n1, double n2) {
   /*Rotate the momentum component q (q[1], q[2])
@@ -168 +169 @@
 // This is used by flux functions
 // Input parameters uh and h may be modified by this function.
-
-// FIXME: Perhaps inline this and profile
+// Tried to inline, but no speedup was achieved 27th May 2009 (Ole)
+//static inline double _compute_speed(double *uh, 
 double _compute_speed(double *uh, 
-              double *h, 
-              double epsilon, 
-              double h0) {
+                      double *h, 
+                      double epsilon, 
+                      double h0) {
   
   double u;
@@ -280 +281 @@
   double h0 = H0*H0; // This ensures a good balance when h approaches H0.
                      // But evidence suggests that h0 can be as little as
-             // epsilon!
+                     // epsilon!
   
   // Copy conserved quantities to protect from modification
@@ -295 +296 @@
 
   z = 0.5*(z_left + z_right); // Average elevation values. 
-                            // Even though this will nominally allow for discontinuities 
-                            // in the elevation data, there is currently no numerical 
-                            // support for this so results may be strange near jumps in the bed.
+                              // Even though this will nominally allow 
+                              // for discontinuities in the elevation data, 
+                              // there is currently no numerical support for
+                              // this so results may be strange near 
+                              // jumps in the bed.
 
   // Compute speeds in x-direction