Mercurial > repos > public > sbplib_julia

view plotDerivative2d.jl @ 91:c0f33eccd527 cell_based_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Create types StencilIndex, LowerClosureIndex, UpperClosureIndex and InteriorIndex. First attempt at seperating out interior and closure indices from. Not fully implemented.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 29 Jan 2019 14:32:28 +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)