Mercurial > repos > public > sbplib_julia

view plotDerivative2d.jl @ 1346:d2ffb412a0be feature/variable_derivatives

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Pad closure stencils with zero rows to ensure that SecondDerivativeVarible can infer the size of the closure stencils
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 03 May 2023 08:51:29 +0200
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)