Mercurial > repos > public > sbplib_julia

view plotDerivative2d.jl @ 1445:76a121ee7b80 bugfix/sbp_operators/stencil_return_type

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add tests for apply with NestedStencils and differing types in the stencils and arrays
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 24 Nov 2023 21:23:41 +0100
parents d61bfc8cf6a5
children
line wrap: on
line source

include("sbpPlot.jl")

g = sbp.Grid.EquidistantGrid((100,75), (0.0, 0.0), (2pi, 3/2*pi))
op = sbp.readOperator("d2_4th.txt","h_4th.txt")
Laplace = sbp.Laplace(g, 1.0, op)

init(x,y) = sin(x) + cos(y)
v = sbp.Grid.evalOn(g,init)
u = zero(v)

sbp.apply!(Laplace,u,v)

#@show u
#@show u'*u

plotgridfunction(g,u)