From 04df1885ff18be3d87fdb72b0b2145aa627b0588 Mon Sep 17 00:00:00 2001 From: Johnny Vogels <35307256+jmv2009@users.noreply.github.com> Date: Sun, 15 Dec 2024 12:37:27 +0100 Subject: [PATCH] Update demo_tnt-elements.py --- python/demo/demo_tnt-elements.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/python/demo/demo_tnt-elements.py b/python/demo/demo_tnt-elements.py index 071c29fe01..ee67a0e25a 100644 --- a/python/demo/demo_tnt-elements.py +++ b/python/demo/demo_tnt-elements.py @@ -168,7 +168,7 @@ def create_tnt_quad(degree): - assert degree > 1 + assert degree > 0 # Polyset ndofs = (degree + 1) ** 2 + 4 npoly = (degree + 2) ** 2 @@ -309,7 +309,7 @@ def poisson_error(V: fem.FunctionSpace): tnt_ndofs.append(V.dofmap.index_map.size_global) tnt_errors.append(poisson_error(V)) print(f"TNT degree 1 error: {tnt_errors[-1]}") -for degree in range(2, 9): +for degree in range(1, 9): V = fem.functionspace(msh, create_tnt_quad(degree)) tnt_degrees.append(degree + 1) tnt_ndofs.append(V.dofmap.index_map.size_global)