Mercurial > repos > public > sbplib_julia

view plotDerivative.jl @ 1320:6ae2ec4cef5a tooling/benchmarks

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Rename functions for comparing benchmarks. Refactor their code. Allow comparing the current working directory
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 28 Apr 2023 08:53:39 +0200
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)