Update fenics_sdc.py

Parallel-in-Time · Dec 24, 2023 · 511cca0 · 511cca0
1 parent 850b4ec
commit 511cca0
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/pySDC/playgrounds/FEniCS/fenics_sdc.py b/pySDC/playgrounds/FEniCS/fenics_sdc.py
@@ -43,6 +43,7 @@ def uexact(t, order, V):
 
 # set boundary values
 bc = df.DirichletBC(V, df.Constant(1.0), Boundary)
+bh = df.DirichletBC(V, df.Constant(0.0), Boundary)
 
 # Stiffness term (Laplace)
 u_trial = df.TrialFunction(V)
@@ -117,6 +118,7 @@ def uexact(t, order, V):
         for j in range(nnodes + 1):
             res -= dt * Q[m + 1, j] * (fimpl[j] + fexpl[j])
         res = df.project(res, V)
+        bh.apply(res.vector())
         res_norm.append(df.norm(res, 'L2'))
 
     uex = uexact(t0 + dt, order, V)