From e186f3edf5a3ca217d43bd3d17ee2b2aca2d9bf8 Mon Sep 17 00:00:00 2001 From: Kanav Gupta Date: Sat, 10 Aug 2019 21:38:57 +0530 Subject: [PATCH] fix Brusselator Problem --- src/ode/brusselator_prob.jl | 12 ++++-------- 1 file changed, 4 insertions(+), 8 deletions(-) diff --git a/src/ode/brusselator_prob.jl b/src/ode/brusselator_prob.jl index 3504f78..5b61a17 100644 --- a/src/ode/brusselator_prob.jl +++ b/src/ode/brusselator_prob.jl @@ -97,14 +97,10 @@ prob_ode_brusselator_2d = ODEProblem(brusselator_2d_loop, length(xyd_brusselator))) const N_brusselator_1d = 40 -const D_brusselator_u = DerivativeOperator{Float64}(2,2,1/(N_brusselator_1d-1), - N_brusselator_1d, - :Dirichlet,:Dirichlet; - BC=(1.,1.)) -const D_brusselator_v = DerivativeOperator{Float64}(2,2,1/(N_brusselator_1d-1), - N_brusselator_1d, - :Dirichlet,:Dirichlet; - BC=(3.,3.)) +const D_brusselator_u = CenteredDifference{Float64}(2,2,1/(N_brusselator_1d-1), + N_brusselator_1d) +const D_brusselator_v = CenteredDifference{Float64}(2,2,1/(N_brusselator_1d-1), + N_brusselator_1d) function brusselator_1d(du, u_, p, t) A, B, α, buffer = p u = @view(u_[:, 1])