Mercurial > repos > public > sbplib_julia

comparison plotDerivative2d.jl @ 138:d61bfc8cf6a5 cell_based_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change from using PyPlot to Makie. Add ugly workaround in order to enable use of sbpPlot
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 21 Feb 2019 16:51:27 +0100
parents 8d505e9bc715
children
comparison
equal deleted inserted replaced
137:99308f68e548 138:d61bfc8cf6a5
1 include("sbpPlot.jl")
2
1 g = sbp.Grid.EquidistantGrid((100,75), (0.0, 0.0), (2pi, 3/2*pi)) 3 g = sbp.Grid.EquidistantGrid((100,75), (0.0, 0.0), (2pi, 3/2*pi))
2 op = sbp.readOperator("d2_4th.txt","h_4th.txt") 4 op = sbp.readOperator("d2_4th.txt","h_4th.txt")
3 Laplace = sbp.Laplace(g, 1.0, op) 5 Laplace = sbp.Laplace(g, 1.0, op)
4 6
5 init(x,y) = sin(x) + cos(y) 7 init(x,y) = sin(x) + cos(y)
6 v = sbp.Grid.evalOn(g,init) 8 v = sbp.Grid.evalOn(g,init)
7 9 u = zero(v)
8 u = zeros(length(v))
9 10
10 sbp.apply!(Laplace,u,v) 11 sbp.apply!(Laplace,u,v)
11 12
12 #@show u 13 #@show u
13 #@show u'*u 14 #@show u'*u
14 15
15 sbp.Grid.plotgridfunction(g,u) 16 plotgridfunction(g,u)
16