Mercurial > repos > public > sbplib_julia

annotate plotDerivative2d.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 d61bfc8cf6a5
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
138
d61bfc8cf6a5 Change from using PyPlot to Makie. Add ugly workaround in order to enable use of sbpPlot
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 85
diff changeset
1 include("sbpPlot.jl")
d61bfc8cf6a5 Change from using PyPlot to Makie. Add ugly workaround in order to enable use of sbpPlot
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 85
diff changeset
2
76
81d9510cb2d0 Make Laplace take dimension as a parameter
Jonatan Werpers <jonatan@werpers.com>
parents: 54
diff changeset
3 g = sbp.Grid.EquidistantGrid((100,75), (0.0, 0.0), (2pi, 3/2*pi))
50
2a7d0ed7ac10 Add script for differentiation and plotting of 2d function
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 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
5 Laplace = sbp.Laplace(g, 1.0, op)
50
2a7d0ed7ac10 Add script for differentiation and plotting of 2d function
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6
2a7d0ed7ac10 Add script for differentiation and plotting of 2d function
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 init(x,y) = sin(x) + cos(y)
54
4300a3fbd818 switch grid to Grid in diffOp
Ylva Rydin <ylva.rydin@telia.com>
parents: 50
diff changeset
8 v = sbp.Grid.evalOn(g,init)
138
d61bfc8cf6a5 Change from using PyPlot to Makie. Add ugly workaround in order to enable use of sbpPlot
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 85
diff changeset
9 u = zero(v)
50
2a7d0ed7ac10 Add script for differentiation and plotting of 2d function
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10
2a7d0ed7ac10 Add script for differentiation and plotting of 2d function
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11 sbp.apply!(Laplace,u,v)
2a7d0ed7ac10 Add script for differentiation and plotting of 2d function
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12
56
27a8d3021a1c Convert apply functions to cell-based
Ylva Rydin <ylva.rydin@telia.com>
parents: 54
diff changeset
13 #@show u
27a8d3021a1c Convert apply functions to cell-based
Ylva Rydin <ylva.rydin@telia.com>
parents: 54
diff changeset
14 #@show u'*u
50
2a7d0ed7ac10 Add script for differentiation and plotting of 2d function
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
15
138
d61bfc8cf6a5 Change from using PyPlot to Makie. Add ugly workaround in order to enable use of sbpPlot
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents: 85
diff changeset
16 plotgridfunction(g,u)