Changeset 274 for inundation/ga/storm_surge/pyvolution/test_quantity.py

Timestamp:

Sep 6, 2004, 1:57:41 PM (21 years ago)

Author:

ole

Message:

Cleaned up distribute

File:

• inundation/ga/storm_surge/pyvolution/test_quantity.py (modified) (1 diff)

inundation/ga/storm_surge/pyvolution/test_quantity.py

@@ -254 +254 @@
 
 
-#     def test_limiter(self):
-
-#         initialise_consecutive_datastructure(points=6+4, elements=4)          
-        
-#         a = Point (0.0, 0.0)
-#         b = Point (0.0, 2.0)
-#         c = Point (2.0, 0.0)
-#         d = Point (0.0, 4.0)
-#         e = Point (2.0, 2.0)
-#         f = Point (4.0, 0.0)
-
-#         #Set up for a gradient of (3,1), f(x) = 3x+y
-#         v1 = Volume(b,a,c,array([0.0,0,0]))        
-#         v2 = Volume(b,c,e,array([1.0,0,0]))
-#         v3 = Volume(e,c,f,array([10.0,0,0]))
-#         v4 = Volume(d,b,e,array([0.0,0,0]))
-
-#         #Setup neighbour structure
-#         domain = Domain([v1,v2,v3,v4])
-#         domain.precompute()
-        
-#         #Lets's check first order first, hey
-#       domain.order = 1
-#       domain.limiter = None
-#         distribute_to_vertices_and_edges(domain)
-#         assert allclose(v2.conserved_quantities_vertex0,
-#                         v2.conserved_quantities_centroid)
-#         assert allclose(v2.conserved_quantities_vertex1,
-#                         v2.conserved_quantities_centroid)
-#         assert allclose(v2.conserved_quantities_vertex2,
-#                         v2.conserved_quantities_centroid)
-
-
-#         #Gradient of fitted pwl surface
-#         a, b = compute_gradient(v2.id)
-        
-
-#         assert abs(a[0] - 5.0) < epsilon
-#         assert abs(b[0]) < epsilon        
-#         #assert qminr[0] == 0.0
-#         #assert qmaxr[0] == 10.0
-        
-#         #And now for the second order stuff        
-#         # - the full second order extrapolation
-#       domain.order = 2        
-#         distribute_to_vertices_and_edges(domain)
-
-
-#         qmin = qmax = v2.conserved_quantities_centroid
-
-#         qmin = minimum(qmin, v1.conserved_quantities_centroid)
-#         qmax = maximum(qmax, v1.conserved_quantities_centroid)
-
-#         qmin = minimum(qmin, v3.conserved_quantities_centroid)
-#         qmax = maximum(qmax, v3.conserved_quantities_centroid)
-
-#         qmin = minimum(qmin, v4.conserved_quantities_centroid)
-#         qmax = maximum(qmax, v4.conserved_quantities_centroid)                    
-#         #assert qminr == qmin
-#         #assert qmaxr == qmax
-
-#         assert v2.conserved_quantities_vertex0 <= qmax
-#         assert v2.conserved_quantities_vertex0 >= qmin
-#         assert v2.conserved_quantities_vertex1 <= qmax
-#         assert v2.conserved_quantities_vertex1 >= qmin
-#         assert v2.conserved_quantities_vertex2 <= qmax
-#         assert v2.conserved_quantities_vertex2 >= qmin                    
-
-
-#         #Check that volume has been preserved    
-
-#         q = v2.conserved_quantities_centroid[0]
-#         w = (v2.conserved_quantities_vertex0[0] +
-#              v2.conserved_quantities_vertex1[0] +
-#              v2.conserved_quantities_vertex2[0])/3
-
-#         assert allclose(q, w)
-
-
-        
-
-
     def test_first_order_extrapolator(self):
         quantity = Conserved_quantity(self.mesh4)