Add Elasticity/Stokes using biharmonic/Laplace kernels

Author: alexfikl

doc/conf.py

@@ -80,6 +80,7 @@
     # sumpy
     "ExpansionBase": "class:sumpy.expansion.ExpansionBase",
     "ExpansionFactoryBase": "class:sumpy.expansion.ExpansionFactoryBase",
+    "ExpressionKernel": "class:sumpy.kernel.ExpressionKernel",
     "Kernel": "class:sumpy.kernel.Kernel",
     "P2PBase": "class:sumpy.p2p.P2PBase",
     # pytential
@@ -91,6 +92,7 @@
     "KernelArgumentMapping": "obj:pytential.symbolic.primitives.KernelArgumentMapping",
     "Operand": "obj:pytential.symbolic.primitives.Operand",
     "QBXForcedLimit": "obj:pytential.symbolic.primitives.QBXForcedLimit",
+    "VectorExpression": "obj:pytential.symbolic.elasticity.VectorExpression",
     "TargetOrDiscretization": "obj:pytential.target.TargetOrDiscretization",
     "pytential.symbolic.dof_desc.DOFDescriptorLike":
         "data:pytential.symbolic.dof_desc.DOFDescriptorLike",

doc/symbolic.rst

@@ -38,6 +38,11 @@ Maxwell's equations
 
 .. automodule:: pytential.symbolic.pde.maxwell
 
+Elasticity equations
+^^^^^^^^^^^^^^^^^^^^
+
+.. automodule:: pytential.symbolic.elasticity
+
 Stokes' equations
 ^^^^^^^^^^^^^^^^^
 
@@ -48,6 +53,11 @@ Scalar Beltrami equations
 
 .. automodule:: pytential.symbolic.pde.beltrami
 
+Rewriting expressions with layer potentials
+^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+
+.. automodule:: pytential.symbolic.pde.system_utils
+
 Internal affairs
 ----------------
 
@@ -62,3 +72,10 @@ How a symbolic operator gets executed
 .. automodule:: pytential.symbolic.execution
 
 .. automodule:: pytential.symbolic.compiler
+
+Rewriting expressions with ``IntG``\ s internals
+^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+
+.. automethod:: pytential.symbolic.pde.system_utils.convert_target_transformation_to_source
+.. automethod:: pytential.symbolic.pde.system_utils.rewrite_int_g_using_base_kernel
+

pyproject.toml

@@ -145,6 +145,8 @@ extend-ignore-re = [
 
 [tool.typos.default.extend-words]
 "nd" = "nd"
+# short for multi-indices
+mis = "mis"
 
 [tool.basedpyright]
 reportImplicitStringConcatenation = "none"

pytential/qbx/__init__.py

@@ -1047,9 +1047,12 @@ def get_flat_strengths_from_densities(
             places,
             sym.weights_and_area_elements(places.ambient_dim, dofdesc=dofdesc),
             )(actx)
+
+    from numbers import Number
+
     density_dofarrays = [evaluate(density) for density in densities]
     for i, ary in enumerate(density_dofarrays):
-        if not isinstance(ary, DOFArray):
+        if not isinstance(ary, DOFArray | Number):
             raise ValueError(
                 f"DOFArray expected for density '{densities[i]}', "
                 f"{type(ary)} received instead")