Mercurial > repos > public > sbplib_julia

view plotDerivative2d.jl @ 1060:ff0e819f2075 feature/nested_stencils

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Refactor code for regular stencils to use fewer type parameters and allow promotion (grafted from 737cd68318c7118f9f20a5fbc31431c48962bc71)
author Jonatan Werpers <jonatan@werpers.com>
date Sat, 12 Feb 2022 22:02:06 +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)