Mercurial > repos > public > sbplib_julia

view plotDerivative.jl @ 46:50c6c252d954

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix issue in plotgridfunction in 2d. Remove 3d case. Will be added later
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 10 Jan 2019 17:39:42 +0100
parents 56722ba36604
children 614b56a017b9
line wrap: on
line source

g = sbp.grid.EquidistantGrid(200,(0, 2pi))
op =sbp.readOperator("d2_4th.txt","h_4th.txt")
Laplace = sbp.Laplace1D(g,1,op)

init(x) = sin(x)
v = sbp.grid.evalOn(g,init)
u = zeros(length(v))

sbp.apply!(Laplace,u,v)

@show u
sbp.grid.plotgridfunction(g,u)