diff --git a/grudge/symbolic/primitives.py b/grudge/symbolic/primitives.py
index fd963eba539436c881f70b6e6b491c048217c523..ec75087e595cf58c75fc760f3d6876e3b4a656a4 100644
--- a/grudge/symbolic/primitives.py
+++ b/grudge/symbolic/primitives.py
@@ -673,7 +673,7 @@ def inverse_surface_metric_derivative(rst_axis, xyz_axis,
     else:
         inv_form1 = inverse_first_fundamental_form(ambient_dim, dim=dim, dd=dd)
         imd = sum(
-                inv_form1[rst_axis, d] * forward_metric_derivative(xyz_axis, d, dd=dd)
+                inv_form1[rst_axis, d]*forward_metric_derivative(xyz_axis, d, dd=dd)
                 for d in range(dim))
 
         return cse(imd,
diff --git a/test/mesh_data.py b/test/mesh_data.py
index d1d04dacf9794387e4289a3dab1cca6762c7ac65..f508b520b32f8d051459a6511eb3d8597065bbdd 100644
--- a/test/mesh_data.py
+++ b/test/mesh_data.py
@@ -91,11 +91,11 @@ class SpheroidMeshBuilder(MeshBuilder):
     def get_mesh(self, resolution, mesh_order):
         from meshmode.mesh.io import ScriptSource
         source = ScriptSource("""
-        SetFactory("OpenCASCADE");
-        Sphere(1) = {{0, 0, 0, {r}}};
-        Dilate {{ {{0, 0, 0}}, {{ {r}, {r}, {rr} }} }} {{ Volume{{1}}; }}
-        """.format(r=self.diameter, rr=self.aspect_ratio * self.diameter),
-        "geo"
+            SetFactory("OpenCASCADE");
+            Sphere(1) = {{0, 0, 0, {r}}};
+            Dilate {{ {{0, 0, 0}}, {{ {r}, {r}, {rr} }} }} {{ Volume{{1}}; }}
+            """.format(r=self.diameter, rr=self.aspect_ratio * self.diameter),
+            "geo"
         )
 
         from meshmode.mesh.io import generate_gmsh
diff --git a/test/test_grudge.py b/test/test_grudge.py
index e0d08dfc12253561a5d5bd5ff9e2e944f389891c..366494ec9b6440c666b9f037fa084c8bbc1f0b54 100644
--- a/test/test_grudge.py
+++ b/test/test_grudge.py
@@ -275,11 +275,9 @@ def test_surface_mass_operator_inverse(ctx_factory, name):
     if name == "2-1-ellipse":
         from mesh_data import EllipseMeshBuilder
         builder = EllipseMeshBuilder(radius=3.1, aspect_ratio=2.0)
-        surface_area = _ellipse_surface_area(builder.radius, builder.aspect_ratio)
     elif name == "spheroid":
         from mesh_data import SpheroidMeshBuilder
         builder = SpheroidMeshBuilder()
-        surface_area = _spheroid_surface_area(builder.radius, builder.aspect_ratio)
     else:
         raise ValueError("unknown geometry name: %s" % name)
 
@@ -601,7 +599,7 @@ def test_surface_divergence_theorem(ctx_factory, mesh_name, visualize=False):
     for i, resolution in enumerate(builder.resolutions):
         from meshmode.mesh.processing import affine_map
         mesh = builder.get_mesh(resolution, builder.mesh_order)
-        # mesh = affine_map(mesh, A=mesh_rotation, b=mesh_offset)
+        mesh = affine_map(mesh, A=mesh_rotation, b=mesh_offset)
 
         from meshmode.discretization.poly_element import \
                 QuadratureSimplexGroupFactory