- Timestamp:
- Jan 13, 2009, 11:51:22 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/abstract_2d_finite_volumes/test_pmesh2domain.py
r5897 r6145 3 3 4 4 import unittest 5 6 from Numeric import allclose, array7 8 5 9 6 #from anuga.pyvolution.pmesh2domain import * … … 16 13 from anuga.coordinate_transforms.geo_reference import Geo_reference 17 14 from anuga.pmesh.mesh import importMeshFromFile 15 16 import Numeric as num 17 18 18 19 19 class Test_pmesh2domain(unittest.TestCase): … … 66 66 67 67 tags = {} 68 b1 = Dirichlet_boundary(conserved_quantities = array([0.0]))69 b2 = Dirichlet_boundary(conserved_quantities = array([1.0]))70 b3 = Dirichlet_boundary(conserved_quantities = array([2.0]))68 b1 = Dirichlet_boundary(conserved_quantities = num.array([0.0])) 69 b2 = Dirichlet_boundary(conserved_quantities = num.array([1.0])) 70 b3 = Dirichlet_boundary(conserved_quantities = num.array([2.0])) 71 71 tags["1"] = b1 72 72 tags["2"] = b2 … … 80 80 answer = [[0., 8., 0.], 81 81 [0., 10., 8.]] 82 assert allclose(domain.quantities['elevation'].vertex_values,83 answer)82 assert num.allclose(domain.quantities['elevation'].vertex_values, 83 answer) 84 84 85 85 #print domain.quantities['stage'].vertex_values 86 86 answer = [[0., 12., 10.], 87 87 [0., 10., 12.]] 88 assert allclose(domain.quantities['stage'].vertex_values,89 answer)88 assert num.allclose(domain.quantities['stage'].vertex_values, 89 answer) 90 90 91 91 #print domain.quantities['friction'].vertex_values 92 92 answer = [[0.01, 0.04, 0.03], 93 93 [0.01, 0.02, 0.04]] 94 assert allclose(domain.quantities['friction'].vertex_values,95 answer)96 97 #print domain.quantities['friction'].vertex_values 98 assert allclose(domain.tagged_elements['dsg'][0],0)99 assert allclose(domain.tagged_elements['ole nielsen'][0],1)94 assert num.allclose(domain.quantities['friction'].vertex_values, 95 answer) 96 97 #print domain.quantities['friction'].vertex_values 98 assert num.allclose(domain.tagged_elements['dsg'][0],0) 99 assert num.allclose(domain.tagged_elements['ole nielsen'][0],1) 100 100 101 101 self.failUnless( domain.boundary[(1, 0)] == '1', … … 159 159 160 160 tags = {} 161 b1 = Dirichlet_boundary(conserved_quantities = array([0.0]))162 b2 = Dirichlet_boundary(conserved_quantities = array([1.0]))163 b3 = Dirichlet_boundary(conserved_quantities = array([2.0]))161 b1 = Dirichlet_boundary(conserved_quantities = num.array([0.0])) 162 b2 = Dirichlet_boundary(conserved_quantities = num.array([1.0])) 163 b3 = Dirichlet_boundary(conserved_quantities = num.array([2.0])) 164 164 tags["1"] = b1 165 165 tags["2"] = b2 … … 174 174 answer = [[0., 8., 0.], 175 175 [0., 10., 8.]] 176 assert allclose(domain.quantities['elevation'].vertex_values,177 answer)176 assert num.allclose(domain.quantities['elevation'].vertex_values, 177 answer) 178 178 179 179 #print domain.quantities['stage'].vertex_values 180 180 answer = [[0., 12., 10.], 181 181 [0., 10., 12.]] 182 assert allclose(domain.quantities['stage'].vertex_values,183 answer)182 assert num.allclose(domain.quantities['stage'].vertex_values, 183 answer) 184 184 185 185 #print domain.quantities['friction'].vertex_values 186 186 answer = [[0.01, 0.04, 0.03], 187 187 [0.01, 0.02, 0.04]] 188 assert allclose(domain.quantities['friction'].vertex_values,189 answer)190 191 #print domain.quantities['friction'].vertex_values 192 assert allclose(domain.tagged_elements['dsg'][0],0)193 assert allclose(domain.tagged_elements['ole nielsen'][0],1)188 assert num.allclose(domain.quantities['friction'].vertex_values, 189 answer) 190 191 #print domain.quantities['friction'].vertex_values 192 assert num.allclose(domain.tagged_elements['dsg'][0],0) 193 assert num.allclose(domain.tagged_elements['ole nielsen'][0],1) 194 194 195 195 self.failUnless( domain.boundary[(1, 0)] == '1', … … 228 228 #domain.set_tag_dict(tag_dict) 229 229 #Boundary tests 230 b1 = Dirichlet_boundary(conserved_quantities = array([0.0]))231 b2 = Dirichlet_boundary(conserved_quantities = array([1.0]))232 b3 = Dirichlet_boundary(conserved_quantities = array([1.0]))230 b1 = Dirichlet_boundary(conserved_quantities = num.array([0.0])) 231 b2 = Dirichlet_boundary(conserved_quantities = num.array([1.0])) 232 b3 = Dirichlet_boundary(conserved_quantities = num.array([1.0])) 233 233 #test adding a boundary 234 234 tags = {}
Note: See TracChangeset
for help on using the changeset viewer.