Test Stokes Identities

pytential/__init__.py

@@ -100,11 +100,8 @@ def norm(discr, x, p=2):
     if isinstance(x, MultiVector):
         x = x.as_vector(object)
 
-    from meshmode.dof_array import DOFArray
     num_components = None
-    if (isinstance(x, np.ndarray)
-            and x.dtype.char == "O"
-            and not isinstance(x, DOFArray)):
+    if isinstance(x, np.ndarray) and x.dtype.char == "O":
         num_components, = x.shape
 
     if p == 2:

pytential/symbolic/primitives.py

@@ -1056,6 +1056,17 @@ def h_max(ambient_dim, dim=None, dofdesc=None):
             cse_scope.DISCRETIZATION)
 
 
+def h_min(ambient_dim, dim=None, dofdesc=None):
+    """Yields an approximate minimum element size in the discretization."""
+
+    dofdesc = as_dofdesc(dofdesc).copy(granularity=GRANULARITY_ELEMENT)
+    r = _quad_resolution(ambient_dim, dim=dim, dofdesc=dofdesc)
+
+    return cse(NodeMin(r),
+            "h_min",
+            cse_scope.DISCRETIZATION)
+
+
 def weights_and_area_elements(ambient_dim, dim=None, dofdesc=None):
     """Combines :func:`area_element` and :class:`QWeight`."""
 

test/extra_int_eq_data.py

@@ -291,6 +291,16 @@ class CircleTestCase(EllipseTestCase):
     aspect_ratio = 1.0
     radius = 1.0
 
+
+class StarfishTestCase(CurveTestCase):
+    name = "starfish"
+    narms = 5
+    amplitude = 0.25
+
+    def _curve_fn(self, t):
+        from meshmode.mesh.generation import NArmedStarfish
+        return NArmedStarfish(self.narms, self.amplitude)(t)
+
 # }}}
 
 
@@ -356,17 +366,32 @@ class SphereTestCase(IntegralEquationTestCase):
 
     # test case
     resolutions = [1, 2]
-    inner_radius = 0.4
-    outer_radius = 5.0
     check_gradient = False
     check_tangential_deriv = False
 
+    radius = 1.0
+    inner_radius = 0.4
+    outer_radius = 5.0
+
     def get_mesh(self, resolution, mesh_order):
         from meshmode.mesh.generation import generate_sphere
-        return generate_sphere(1.0, mesh_order,
+        return generate_sphere(self.radius, mesh_order,
                 uniform_refinement_rounds=resolution)
 
 
+class SpheroidTestCase(SphereTestCase):
+    name = "spheroid"
+    aspect_ratio = 2.0
+
+    def get_mesh(self, resolution, mesh_order):
+        mesh = super().get_mesh(resolution, mesh_order)
+
+        from meshmode.mesh.processing import affine_map
+        return affine_map(mesh, A=np.diag([
+            self.radius, self.radius, self.radius / self.aspect_ratio,
+            ]))
+
+
 class TorusTestCase(IntegralEquationTestCase):
     ambient_dim = 3
     name = "torus"