Mercurial > repos > public > sbplib_julia

view plotDerivative2d.jl @ 92:b8c9e2db126f stencil_index

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
New branch for implementing specialized indices used by the diffOps apply functions.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 04 Feb 2019 09:11:53 +0100
parents 8d505e9bc715
children d61bfc8cf6a5
line wrap: on
line source

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 = zeros(length(v))

sbp.apply!(Laplace,u,v)

#@show u
#@show u'*u

sbp.Grid.plotgridfunction(g,u)