Extend interpolate_geometry tests to GLL (#4223)

* Also lift to gll geometry

* gll_isaac - what does it do on quads?

* ruff

Author: jhale

python/test/unit/fem/test_interpolate_geometry.py

@@ -80,7 +80,7 @@ def test_interpolate_geometry_p1_roundtrip(dtype):
         np.testing.assert_allclose(x_new[dm_new[c]], x_old[dm_old[c]], atol=atol, rtol=0.0)
 
 
-def _curve_mesh_errors(N, degree, dtype, R, cell_type):
+def _curve_mesh_errors(N, degree, dtype, R, cell_type, lagrange_variant):
     """Return (area_error, circumference_error) for a degree-p curved disk mesh with N cells."""
     mesh = create_rectangle(
         MPI.COMM_WORLD,
@@ -98,7 +98,7 @@ def transform(x):
         x_c[:, 1] = R * x[:, 1] * np.sqrt(1.0 - (x[:, 0] ** 2 / (2.0)))
         return x_c
 
-    cmap = coordinate_element(cell_type, degree, variant=LagrangeVariant.equispaced, dtype=dtype)
+    cmap = coordinate_element(cell_type, degree, variant=lagrange_variant, dtype=dtype)
     curved_mesh = interpolate_geometry(mesh, cmap)
     curved_mesh.geometry.x[:] = transform(curved_mesh.geometry.x)
 
@@ -131,11 +131,14 @@ def transform(x):
 
 @pytest.mark.parametrize("dtype", [np.float32, np.float64])
 @pytest.mark.parametrize("degree", [1, 2, 3])
-@pytest.mark.parametrize("R", [0.1, 1])
+@pytest.mark.parametrize("R", [0.1])
 @pytest.mark.parametrize("cell_type", [CellType.triangle, CellType.quadrilateral])
-def test_curve_mesh(degree, dtype, R, cell_type):
+@pytest.mark.parametrize(
+    "lagrange_variant", [LagrangeVariant.equispaced, LagrangeVariant.gll_isaac]
+)
+def test_curve_mesh(degree, dtype, R, cell_type, lagrange_variant):
     Ns = [4, 8, 16, 32]
-    errors = [_curve_mesh_errors(N, degree, dtype, R, cell_type) for N in Ns]
+    errors = [_curve_mesh_errors(N, degree, dtype, R, cell_type, lagrange_variant) for N in Ns]
 
     area_errors = np.array([e[0] for e in errors])
     circ_errors = np.array([e[1] for e in errors])