Changeset 6252 for anuga_core/source/anuga/fit_interpolate/fit.py
- Timestamp:
- Jan 31, 2009, 4:02:38 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/fit_interpolate/fit.py
r6244 r6252 267 267 self.mesh_boundary_polygon, 268 268 closed=True, 269 verbose=False) # Too much output if True 270 269 verbose=False) # Suppress output 271 270 272 271 n = len(inside_indices) … … 301 300 AtA[j,k] += sigmas[j]*sigmas[k] 302 301 else: 303 # FIXME(Ole): This is the message referred to in ticket:314304 305 302 flag = is_inside_polygon(x, 306 303 self.mesh_boundary_polygon, 307 304 closed=True, 308 verbose=False) # Too much output if True305 verbose=False) # Suppress output 309 306 msg = 'Point (%f, %f) is not inside mesh boundary' % tuple(x) 310 307 assert flag is True, msg 311 308 309 # FIXME(Ole): This is the message referred to in ticket:314 312 310 minx = min(self.mesh_boundary_polygon[:,0]) 313 311 maxx = max(self.mesh_boundary_polygon[:,0]) … … 317 315 msg += 'Mesh boundary extent is (%.f, %.f), (%.f, %.f)'\ 318 316 % (minx, maxx, miny, maxy) 317 raise RuntimeError, msg 319 318 320 321 raise Exception(msg)322 319 self.AtA = AtA 323 320
Note: See TracChangeset
for help on using the changeset viewer.