- Timestamp:
- Nov 7, 2006, 1:13:03 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/abstract_2d_finite_volumes/quantity.py
r3918 r3928 29 29 30 30 if vertex_values is None: 31 N = domain.number_of_elements31 N = len(domain) # number_of_elements 32 32 self.vertex_values = zeros((N, 3), Float) 33 33 else: … … 42 42 %N 43 43 msg += 'number of elements in specified domain (%d).'\ 44 % domain.number_of_elements45 46 assert N == domain.number_of_elements, msg44 %len(domain) 45 46 assert N == len(domain), msg 47 47 48 48 self.domain = domain … … 1135 1135 """ 1136 1136 integral = 0 1137 for k in range( self.domain.number_of_elements):1137 for k in range(len(self.domain)): 1138 1138 area = self.domain.areas[k] 1139 1139 qc = self.centroid_values[k] … … 1165 1165 #flux calculations and forcing functions 1166 1166 1167 N = domain.number_of_elements1167 N = len(domain) # number_of_triangles 1168 1168 self.explicit_update = zeros(N, Float ) 1169 1169 self.semi_implicit_update = zeros(N, Float ) … … 1293 1293 1294 1294 #Check each triangle 1295 for k in range( quantity.domain.number_of_elements):1295 for k in range(len(quantity.domain)): 1296 1296 #Centroid coordinates 1297 1297 x, y = quantity.domain.centroid_coordinates[k] … … 1389 1389 from Numeric import zeros, Float 1390 1390 1391 N = quantity.domain.number_of_elements1391 N = len(quantity.domain) 1392 1392 1393 1393 beta_w = quantity.domain.beta_w
Note: See TracChangeset
for help on using the changeset viewer.