Changeset 4200 for anuga_core/source/anuga/shallow_water/shallow_water_ext.c

Timestamp:

Jan 31, 2007, 12:07:32 PM (17 years ago)

Author:

ole

Message:

Added diagnosticts about which triangle is responsible for the largest computed speed.

File:

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

anuga_core/source/anuga/shallow_water/shallow_water_ext.c

@@ -1217 +1217 @@
     *xmom_explicit_update,
     *ymom_explicit_update,
-    *already_computed_flux;//tracks whether the flux across an edge has already been computed
-
+    *already_computed_flux,//tracks whether the flux across an edge has already been computed
+    *max_speed_array; //Keeps track of max speeds for each triangle
 
   //Local variables
@@ -1231 +1231 @@
 
   // Convert Python arguments to C
-  if (!PyArg_ParseTuple(args, "dddOOOOOOOOOOOOOOOOOO",
+  if (!PyArg_ParseTuple(args, "dddOOOOOOOOOOOOOOOOOOO",
                         &timestep,
                         &epsilon,
@@ -1250 +1250 @@
                         &xmom_explicit_update,
                         &ymom_explicit_update,
-                        &already_computed_flux)) {
+                        &already_computed_flux,
+                        &max_speed_array)) {
     PyErr_SetString(PyExc_RuntimeError, "Input arguments failed");
     return NULL;
@@ -1261 +1262 @@
     ((double *) stage_explicit_update -> data)[k]=0.0;
     ((double *) xmom_explicit_update -> data)[k]=0.0;
-    ((double *) ymom_explicit_update -> data)[k]=0.0;
-  }
+    ((double *) ymom_explicit_update -> data)[k]=0.0;  
+  }
+  
   //Loop through neighbours and compute edge flux for each
   for (k=0; k<number_of_elements; k++) {
@@ -1327 +1329 @@
     }
     } // end for i
+    
     //Normalise by area and store for when all conserved
     //quantities get updated
@@ -1332 +1335 @@
     ((double *) xmom_explicit_update -> data)[k] /= ((double *) areas -> data)[k];
     ((double *) ymom_explicit_update -> data)[k] /= ((double *) areas -> data)[k];
+    
+   
+    //Keep track of maximal speeds
+    ((double *) max_speed_array -> data)[k] = max_speed;    
+    
   } //end for k
   return Py_BuildValue("d", timestep);