I just started to play around with DiffEqOperators. I’d like to extend the example at
https://github.com/JuliaDiffEq/DiffEqOperators.jl/blob/master/examples/poisson.jl to a 2D case.
Here is the adapted code:
using DiffEqOperators, Plots f = 5.0 a = -1.0 b = 2.0 nknots = 10 h = 1.0/(nknots+1) ord_deriv = 2 ord_approx = 2 Δ = CenteredDifference(ord_deriv, ord_approx, h, nknots) bc = DirichletBC(a, b) u = (Δ*bc) \ fill(f, (nknots, nknots)) knots = range(h, step=h, length=nknots) plot(knots, u[nknots÷2, 1:end]) plot!(knots, u[1:end, nknots÷2])
It works somehow - but so to say in 1D
So, my question: How do I define the BC for 2D (or more)? Or maybe there are more things to be changed in the script?