Changeset 6042 for anuga_work/development/anuga_1d/comp_flux_vel_ext.c

Timestamp:

Dec 2, 2008, 9:46:45 AM (15 years ago)

Author:

steve

Message:

Found bug associated with changing max_speed_array to an integer

File:

• anuga_work/development/anuga_1d/comp_flux_vel_ext.c (modified) (5 diffs)

anuga_work/development/anuga_1d/comp_flux_vel_ext.c

@@ -145 +145 @@
     double max_speed, normal;
     int k, i, ki, n, m, nm=0;
-    
+
+    //printf("Inside _comp\n");    
     
     for (k=0; k<number_of_elements; k++) {
         flux[0] = 0.0;
         flux[1] = 0.0;
-        
+        //printf("k = %d\n",k);
+
         for (i=0; i<2; i++) {
             ki = k*2+i;
@@ -526 +528 @@
                         &velocity)) {
       PyErr_SetString(PyExc_RuntimeError, "comp_flux_vel_ext.c: compute_fluxes_vel_ext could not parse input");
+      printf("comp_flux_vel_ext.c: compute_fluxes_vel_ext could not parse input");
       return NULL;
   }
@@ -560 +563 @@
     number_of_elements = stage_vertex_values -> dimensions[0];
  
+
+    //printf("comp_flux_vel 1 N = %d %15.5e %15.5e %15.5e %15.5e %15.5e \n",number_of_elements,cfl,timestep,epsilon,g,h0);
+
+
+
     // Call underlying flux computation routine and update 
     // the explicit update arrays 
@@ -585 +593 @@
                                        number_of_elements,
                                        (double*) max_speed_array -> data);
-    
+
+    //printf("comp_flux_vel 2 \n");    
     
     Py_DECREF(neighbours);
@@ -605 +614 @@
     Py_DECREF(max_speed_array);
     
-    
+    //printf("comp_flux_vel 3 \n");
+
     // Return updated flux timestep
     return Py_BuildValue("d", timestep);