Mercurial > repos > public > sbplib_julia

view plotDerivative.jl @ 121:3560f54e3eb3 cell_based_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change implementation of getrange to be more julia(???), mostly as an experiment to see if it breakes stuff
author Jonatan Werpers <jonatan@werpers.com>
date Sun, 10 Feb 2019 19:16:14 +0100
parents 8d505e9bc715
children
line wrap: on
line source

g = sbp.Grid.EquidistantGrid((200,), (0.0,), (2pi,))
op =sbp.readOperator("d2_4th.txt","h_4th.txt")
Laplace = sbp.Laplace(g,1.0,op)

init(x) = cos(x)
v = sbp.Grid.evalOn(g,init)
u = zeros(length(v))

sbp.apply!(Laplace,u,v)

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