Changeset 4829
 Timestamp:
 Nov 19, 2007, 5:05:50 PM (16 years ago)
 Location:
 anuga_core/source/anuga
 Files:

 2 edited
Legend:
 Unmodified
 Added
 Removed

anuga_core/source/anuga/abstract_2d_finite_volumes/domain.py
r4827 r4829 8 8 """ 9 9 10 from Numeric import allclose, argmax 10 from Numeric import allclose, argmax, zeros, Float 11 11 from anuga.config import epsilon 12 12 … … 230 230 # Storage for maximal speeds computed for each triangle by 231 231 # compute_fluxes 232 # This is used for diagnostics only 232 # This is used for diagnostics only (reset at every yieldstep) 233 233 self.max_speed = zeros(N, Float) 234 234 … … 726 726 msg += 'area = %.4f and radius = %.4f ' %(area, radius) 727 727 msg += 'had the largest computed speed: %.6f m/s ' %(max_speed) 728 msg += 'during last time interval. Quantities below ' 729 msg += 'are reported at their present value, and not what ' 730 msg += 'they were at the time the maximal speed was attained.' 731 msg += 'To see this, rerun the model with yieldsteps smaller ' 732 msg += 'than the smallest internal timestep reported.' 728 733 if max_speed > 0.0: 729 734 msg += '(timestep=%.6f)\n' %(radius/max_speed) … … 1055 1060 1056 1061 1057 1062 N = len(self) # Number of triangles 1058 1063 self.yieldtime = 0.0 # Track time between 'yields' 1059 1064 … … 1063 1068 self.number_of_steps = 0 1064 1069 self.number_of_first_order_steps = 0 1070 1065 1071 1066 1072 # Update ghosts … … 1138 1144 self.number_of_steps = 0 1139 1145 self.number_of_first_order_steps = 0 1140 1146 self.max_speed = zeros(N, Float) 1141 1147 1142 1148 def evolve_one_euler_step(self, yieldstep, finaltime): 
anuga_core/source/anuga/shallow_water/shallow_water_ext.c
r4826 r4829 1886 1886 1887 1887 1888 // Keep track of maximal speeds 1889 max_speed_array[k] = max_speed; 1888 // Keep track of maximal speeds so far 1889 if (max_speed > max_speed_array[k]) { 1890 max_speed_array[k] = max_speed; 1891 } 1890 1892 1891 1893 } // End triangle k
Note: See TracChangeset
for help on using the changeset viewer.