diff --git a/test/interpolation/constrained.jl b/test/interpolation/constrained.jl
index 92191fb..d200e87 100644
--- a/test/interpolation/constrained.jl
+++ b/test/interpolation/constrained.jl
@@ -62,7 +62,6 @@ end
     tri = triangulate(points; boundary_nodes)
     A = get_area(tri)
     D = 6.25e-4
-    refine!(tri; max_area=1e-2A)
     Tf = (x, y) -> let r = sqrt(x^2 + y^2)
         (R₂^2 - r^2) / (4D) + R₁^2 * log(r / R₂) / (2D)
     end