Skip to content

Commit

Permalink
Update fenics_sdc.py
Browse files Browse the repository at this point in the history
  • Loading branch information
pancetta committed Dec 24, 2023
1 parent 047291c commit 850b4ec
Showing 1 changed file with 3 additions and 1 deletion.
4 changes: 3 additions & 1 deletion pySDC/playgrounds/FEniCS/fenics_sdc.py
Original file line number Diff line number Diff line change
Expand Up @@ -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))',
Expand All @@ -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
Expand Down

0 comments on commit 850b4ec

Please sign in to comment.