Mercurial > repos > public > sbplib_julia

annotate plotDerivative.jl @ 637:4a81812150f4 feature/volume_and_boundary_operators

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change qudrature closure from tuple of reals to tuple of Stencils. Also remove parametrization of stencil width in D2 since this was illformed for the 2nd order case.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Sun, 03 Jan 2021 18:15:14 +0100
parents 8d505e9bc715
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
76
81d9510cb2d0 Make Laplace take dimension as a parameter
Jonatan Werpers <jonatan@werpers.com>
parents: 54
diff changeset
1 g = sbp.Grid.EquidistantGrid((200,), (0.0,), (2pi,))
45
56722ba36604 ass script that plots the derivative
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
2 op =sbp.readOperator("d2_4th.txt","h_4th.txt")
76
81d9510cb2d0 Make Laplace take dimension as a parameter
Jonatan Werpers <jonatan@werpers.com>
parents: 54
diff changeset
3 Laplace = sbp.Laplace(g,1.0,op)
45
56722ba36604 ass script that plots the derivative
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
4
56
27a8d3021a1c Convert apply functions to cell-based
Ylva Rydin <ylva.rydin@telia.com>
parents: 54
diff changeset
5 init(x) = cos(x)
51
614b56a017b9 Split grid.jl into AbstractGrid.jl and EquidistantGrid.jl
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 46
diff changeset
6 v = sbp.Grid.evalOn(g,init)
45
56722ba36604 ass script that plots the derivative
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
7 u = zeros(length(v))
56722ba36604 ass script that plots the derivative
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
8
56722ba36604 ass script that plots the derivative
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
9 sbp.apply!(Laplace,u,v)
56722ba36604 ass script that plots the derivative
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
10
56722ba36604 ass script that plots the derivative
Ylva Rydin <ylva.rydin@telia.com>
parents:
diff changeset
11 @show u
54
4300a3fbd818 switch grid to Grid in diffOp
Ylva Rydin <ylva.rydin@telia.com>
parents: 51
diff changeset
12 sbp.Grid.plotgridfunction(g,u)