From 850b4ecaff2fdeb8d146ab6949efcd1555414677 Mon Sep 17 00:00:00 2001 From: Robert Speck Date: Sun, 24 Dec 2023 11:48:31 +0100 Subject: [PATCH] Update fenics_sdc.py --- pySDC/playgrounds/FEniCS/fenics_sdc.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/pySDC/playgrounds/FEniCS/fenics_sdc.py b/pySDC/playgrounds/FEniCS/fenics_sdc.py index c920032756..e0419915c5 100644 --- a/pySDC/playgrounds/FEniCS/fenics_sdc.py +++ b/pySDC/playgrounds/FEniCS/fenics_sdc.py @@ -55,6 +55,9 @@ def uexact(t, order, V): M = df.assemble(a_M) K = df.assemble(a_K) +bc.apply(M) +bc.apply(K) + # set forcing term as expression g = df.Expression( '-sin(a*x[0]) * (sin(t) - b*a*a*cos(t))', @@ -64,7 +67,6 @@ def uexact(t, order, V): degree=order, ) - params = {'num_nodes': nnodes, 'quad_type': quad_type, 'node_type': node_type} sweeper = Sweeper(params) Q = sweeper.coll.Qmat