From 15963e99d08aeb8a5f2ed5f3fbfed3aaf98c858c Mon Sep 17 00:00:00 2001 From: "Jack S. Hale" Date: Mon, 13 Jan 2025 09:43:03 +0100 Subject: [PATCH] interpolation_points exposed as array --- demo/demo_nonlinear-naghdi-clamped-semicylinder.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/demo/demo_nonlinear-naghdi-clamped-semicylinder.py b/demo/demo_nonlinear-naghdi-clamped-semicylinder.py index 9ccc1e1..2adda06 100644 --- 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)