I have tried the demo in
And I find there are some problem in “implementation” chapter
1.In Defining the variational problem, the code should be
a = ufl.inner(ufl.grad(u), ufl.grad(v)) * ufl.dx
L = ufl.inner(p, v) * ufl.dx
2.In Plotting the solution over a line,when the debuger run into
plotter.show()
Then going to “connected”.
I’m not sure whether there is any wrong. I just comment out these lines and go on.