Mercurial > repos > public > sbplib_julia

view plotDerivative.jl @ 1134:667e9c588f23 feature/boundary_conditions

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add type for different kind of boundary data and refactor the BoundaryConditions such that they are general w.r.t the boundary data type
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 11 Oct 2022 18:12:40 +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)