Hello i am trying to calculate the laplacian of a 2D array for the use in diffeq.jl And expand later to 3D hopefully. I am currently implementing it using DiffEqOperators.jl with centereddifference() and neuman0BC() However i dont see an easy way to also implement the neumann BC’s in the y direct…

Hi, What you want to do is: dx = 0.01 axis = collect(-1:dx:1); umesh = [(i, j) for i in axis, j in axis]; sines((x,y)) = sin(pi*x)*sin(pi*y) u = sines.(umesh) Qx, Qy = Dirichlet0BC(Float64, size(u)) Dxx = CenteredDifference{1}(2,2,dx,N) Dyy = CenteredDifference{2}(2,2,dx,N) A = (Dxx + Dyy)*com…

Second order 2D differentiation DiffEqOperators.jl

Specific Domains Numerics

ChrisRackauckas March 17, 2022, 11:42pm 6

MethodOfLines.jl is just a better way to do it. DiffEqOperators.jl has been effectively deprecated for years.

Topic		Replies	Views
Trouble with boundary conditions in DiffEqOperators.jl Numerics diffeq	9	1796	March 6, 2019
Question on boundary conditions in DiffEqOperators.jl for 2D-case Numerics question	0	490	January 7, 2020
Optimize a simple code for performance General Usage	3	686	January 20, 2019
(seems) a more efficient way to use diffEqOperator.jl Modelling & Simulations	5	537	July 2, 2021
Concrete Operators in DiffEqOperators.jl for time-indepedent problems Numerics question	9	975	March 7, 2022

Second order 2D differentiation DiffEqOperators.jl

Related topics