Try forcing manual integration

Author: jorgensd

tests/test_average_coefficient.py

@@ -158,23 +158,20 @@ def f(x):
     J_gamma = assemble_scalar(J)
 
     # Symbolic integration of cylinder
-    if domain.comm.rank == 0:
-        x, y, z, R_s, Z_lower, Z_upper = sp.symbols("x y z R Z_lower Z_upper")
-        integral_result = sp.integrate(
-            1 / (sp.pi * R_s**2) * f([x, y, z]),
-            (z, Z_lower, Z_upper),
-            (
-                y,
-                y_c - sp.sqrt(R_s**2 - (x - x_c) ** 2),
-                y_c + sp.sqrt(R_s**2 - (x - x_c) ** 2),
-            ),
-            (x, x_c - R_s, x_c + R_s),
-        )
-        float_result = complex(
-            integral_result.subs({R_s: R, Z_lower: c_min, Z_upper: c_max})
-        ).real
-    reference = domain.comm.bcast(
-        float_result if domain.comm.rank == 0 else None, root=0
+    x, y, z, R_s, Z_lower, Z_upper = sp.symbols("x y z R Z_lower Z_upper")
+    integral_result = sp.integrate(
+        1 / (sp.pi * R_s**2) * f([x, y, z]),
+        (z, Z_lower, Z_upper),
+        (
+            y,
+            y_c - sp.sqrt(R_s**2 - (x - x_c) ** 2),
+            y_c + sp.sqrt(R_s**2 - (x - x_c) ** 2),
+        ),
+        (x, x_c - R_s, x_c + R_s),
+        manual=True,
     )
+    float_result = complex(
+        integral_result.subs({R_s: R, Z_lower: c_min, Z_upper: c_max})
+    ).real
 
-    np.testing.assert_allclose(reference, J_gamma)
+    np.testing.assert_allclose(float_result, J_gamma)