Mercurial > repos > public > sbplib_julia

comparison plotDerivative.jl @ 46:50c6c252d954

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix issue in plotgridfunction in 2d. Remove 3d case. Will be added later
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 10 Jan 2019 17:39:42 +0100
parents 56722ba36604
children 614b56a017b9
comparison
equal deleted inserted replaced
45:56722ba36604 46:50c6c252d954
1 g = sbp.grid.EquidistantGrid(200,(0, 1)) 1 g = sbp.grid.EquidistantGrid(200,(0, 2pi))
2 op =sbp.readOperator("d2_4th.txt","h_4th.txt") 2 op =sbp.readOperator("d2_4th.txt","h_4th.txt")
3 Laplace = sbp.Laplace1D(g,1,op) 3 Laplace = sbp.Laplace1D(g,1,op)
4 4
5 init(x) = sin(x) 5 init(x) = sin(x)
6 v = sbp.grid.evalOn(g,init) 6 v = sbp.grid.evalOn(g,init)
8 8
9 sbp.apply!(Laplace,u,v) 9 sbp.apply!(Laplace,u,v)
10 10
11 @show u 11 @show u
12 sbp.grid.plotgridfunction(g,u) 12 sbp.grid.plotgridfunction(g,u)
13