Changeset 8863

Timestamp:

May 9, 2013, 6:23:53 PM (12 years ago)

Author:

steve

Message:

Getting rid of test of nan

Location:

trunk/anuga_work/development/2010-projects/anuga_1d/base

Files:

• generic_domain.py (modified) (2 diffs)
• limiters_python.py (modified) (1 diff)

trunk/anuga_work/development/2010-projects/anuga_1d/base/generic_domain.py

@@ -1059 +1059 @@
         """
 
-        import numpy as np
-        Stage      = self.quantities['stage']
-        print 'w_ex_update', np.any(np.isnan(Stage.explicit_update))
-
         # Compute fluxes across each element edge
         self.compute_fluxes()
 
-        print 'w_ex_update', np.any(np.isnan(Stage.explicit_update))
 
         # Update timestep to fit yieldstep and finaltime
@@ -1094 +1089 @@
         Q^{n+1} = 0.5 Q^n + 0.5 E(h)^2 Q^n
         """
+
 
         # Save initial conserved quantities values

trunk/anuga_work/development/2010-projects/anuga_1d/base/limiters_python.py

@@ -39 +39 @@
     from numpy import fabs, where
 
-    return where((fabs(a) + fabs(b) >= 1.0e-12), (a*fabs(b)+fabs(a)*b)/(fabs(a)+fabs(b)), 0.0)
+    return where((fabs(a) + fabs(b) >= 1.0e-12), (a*fabs(b)+fabs(a)*b)/(fabs(a)+fabs(b)+ 1.0e-50), 0.0)