Mercurial > repos > public > sbplib_julia

view plotDerivative.jl @ 807:1fc6c4daec44 refactor/sbp_operators_method_signatures

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add type parameter to Vararg to make better error messages when Stencil is called with incompatible weights such as (1,1.,1)
author Jonatan Werpers <jonatan@werpers.com>
date Sat, 18 Sep 2021 23:30:49 +0200
parents 8d505e9bc715
children
line wrap: on
line source

g = sbp.Grid.EquidistantGrid((200,), (0.0,), (2pi,))
op =sbp.readOperator("d2_4th.txt","h_4th.txt")
Laplace = sbp.Laplace(g,1.0,op)

init(x) = cos(x)
v = sbp.Grid.evalOn(g,init)
u = zeros(length(v))

sbp.apply!(Laplace,u,v)

@show u
sbp.Grid.plotgridfunction(g,u)