- Timestamp:
- Oct 27, 2010, 5:22:43 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_work/development/2010-projects/kv/test_kinematic_viscosity.py
r8051 r8052 1 from anuga.abstract_2d_finite_volumes. domain importDomain2 from anuga .abstract_2d_finite_volumes.generic_boundary_conditionsimport Dirichlet_boundary1 from anuga.abstract_2d_finite_volumes.generic_domain import Generic_Domain 2 from anuga import Dirichlet_boundary 3 3 from kinematic_viscosity import Kinematic_Viscosity_Operator 4 4 import numpy as num … … 21 21 boundary_map[(0,1)] = Dirichlet_boundary([1,1,2]) 22 22 boundary_map[(0,2)] = Dirichlet_boundary([1,1,0]) 23 domain = Domain(source=points,triangles=elements,boundary=boundary_map)23 domain = Generic_Domain(source=points,triangles=elements,boundary=boundary_map) 24 24 return Kinematic_Viscosity_Operator(domain) 25 25 … … 33 33 boundary_map[(1,0)] = Dirichlet_boundary([1,1,0]) 34 34 boundary_map[(1,2)] = Dirichlet_boundary([1,2,1]) 35 domain = Domain(source=points,triangles=elements,boundary=boundary_map)35 domain = Generic_Domain(source=points,triangles=elements,boundary=boundary_map) 36 36 return Kinematic_Viscosity_Operator(domain) 37 37 … … 146 146 U_mod[0,:] = U_new 147 147 assert num.allclose(U_new - operator1.dt * 2 * num.mat(A)*num.mat(U_mod), U[0,:]) 148 148 149 ################################################################################ 149 150 150 151 if __name__ == "__main__": 151 suite = unittest.makeSuite(Test_Kinematic_Viscosity, 'test' )152 suite = unittest.makeSuite(Test_Kinematic_Viscosity, 'test', verbose=True) 152 153 runner = unittest.TextTestRunner() 153 154 runner.run(suite)
Note: See TracChangeset
for help on using the changeset viewer.