interpolation_points exposed as array

FEniCS-Shells · Jan 13, 2025 · 15963e9 · 15963e9
1 parent a2e7273
commit 15963e9
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/demo/demo_nonlinear-naghdi-clamped-semicylinder.py b/demo/demo_nonlinear-naghdi-clamped-semicylinder.py
@@ -455,7 +455,7 @@ def gamma(F, d):
 # Calculate the factor alpha as a function of the mesh size h
 h = ufl.CellDiameter(mesh)
 alpha_FS = functionspace(mesh, element("DG", cell, 0))
-alpha_expr = Expression(t**2 / h**2, alpha_FS.element.interpolation_points())
+alpha_expr = Expression(t**2 / h**2, alpha_FS.element.interpolation_points)
 alpha = Function(alpha_FS)
 alpha.interpolate(alpha_expr)
 
@@ -693,7 +693,7 @@ def F(self, x: PETSc.Vec, b: PETSc.Vec) -> None:
 
 # Interpolate phi in the [P2]³ Space
 phi_FS = functionspace(mesh, blocked_element(P2, shape=(3,)))
-phi_expr = Expression(phi0_ufl + u_P2B3, phi_FS.element.interpolation_points())
+phi_expr = Expression(phi0_ufl + u_P2B3, phi_FS.element.interpolation_points)
 phi_func = Function(phi_FS)
 phi_func.interpolate(phi_expr)