Changeset 3757
- Timestamp:
- Oct 11, 2006, 5:05:39 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_validation/automated_validation_tests/okushiri_tank_validation/validate_okushiri.py
r3722 r3757 22 22 import project 23 23 24 try: 25 import pylab 26 except: 27 plotting = False 28 else: 29 plotting = True 30 31 24 32 finaltime = 22.5 25 33 timestep = 0.05 … … 28 36 # Create Domain 29 37 #------------------------- 30 #N = 5031 N = 1038 N = 50 39 #N = 10 32 40 points, vertices, boundary = rectangular_cross(N, N/5*3, 33 41 len1=5.448, len2=3.402) … … 110 118 111 119 112 predicted_gauge_values = [[], [], [],[]]120 predicted_gauge_values = [[], [], [], []] 113 121 114 122 … … 126 134 stage = domain.get_quantity('stage') 127 135 for i, (x, y) in enumerate(gauge_locations): 128 w = stage.get_values(interpolation_points=[[x,y]]) 136 w = stage.get_values(interpolation_points=[[x,y]])[0] 129 137 predicted_gauge_values[i].append(w) 130 138 131 print 'difference', x, y, w - reference_gauge_values[i][j] 139 #print 'difference', x, y,\ 140 # predicted_gauge_values[i][j] - reference_gauge_values[i][j] 132 141 133 142 … … 138 147 predicted_gauge_values[i] = ensure_numeric(predicted_gauge_values[i]) 139 148 149 print predicted_gauge_values[i].shape 150 print reference_gauge_values[i].shape 151 #print predicted_gauge_values.shape 152 #print reference_gauge_values.shape 140 153 141 sqsum = sum((predicted_gauge_values[i] - reference_gauge_values)**2) 142 denom = sum(reference_gauge_values**2) 154 sqsum = sum((predicted_gauge_values[i].flat - \ 155 reference_gauge_values[i].flat)**2) 156 denom = sum(reference_gauge_values[i].flat**2) 143 157 144 158 print i, sqsum, sqsum/denom 145 159 146 msg 'Rms error is %f' %(sqsum/denom) 160 msg = 'Rms error is %f' %(sqsum/denom) 161 print msg 162 147 163 assert sqsum/denom < 0.01, msg 164 148 165 149 try: 166 if plotting is True: 167 150 168 from pylab import * 151 except: 152 pass 153 else: 169 154 170 ion() 155 171 hold(False) … … 165 181 166 182 167 try: 168 from pylab import * 169 except: 170 pass 171 else: 183 if plotting is True: 172 184 show() 173 185
Note: See TracChangeset
for help on using the changeset viewer.