Mercurial > repos > public > sbplib_julia

view plotDerivative2d.jl @ 1284:f19c51ddbd72 tooling/benchmarks

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Update readme with info about running benchmarks from make. Fix typos
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 02 Mar 2023 09:59:42 +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)