Mercurial > repos > public > sbplib_julia

view plotDerivative2d.jl @ 138:d61bfc8cf6a5 cell_based_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change from using PyPlot to Makie. Add ugly workaround in order to enable use of sbpPlot
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 21 Feb 2019 16:51:27 +0100
parents 8d505e9bc715
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)