Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/derivatives/dissipation.jl @ 1096:9f0121e465a5 feature/dissipation_operators

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Docs for undivided_dissipation
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 10 May 2022 21:10:31 +0200
parents 423a6442efc3
children 254934aac3f8
comparison
equal deleted inserted replaced
1095:2be39d485221 1096:9f0121e465a5
1 # REVIEW: Add documentation. Also would it be more correct to 1 # REVIEW: Would it be more correct to
2 # call these undivided_differences instead of dissipation? 2 # call these undivided_differences instead of dissipation?
3 # If I understand it correctly, this method simply provides 3 # If I understand it correctly, this method simply provides
4 # the operators required in order to compose a dissipation operator 4 # the operators required in order to compose a dissipation operator
5 # and the dissipation operator are formed by a linear combination 5 # and the dissipation operator are formed by a linear combination
6 # of the products of Dᵀ and D for different orders. 6 # of the products of Dᵀ and D for different orders.
7 """
8 undivided_dissipation(g::EquidistantGrid, p, direction)
9
10 Create undivided difference operators approximating the `p`th derivative for
11 building artificial dissipation. The operators and how they are used to create accurate artifical dissipation is described in
12 "K. Mattsson, M. Svärd, and J. Nordström, “Stable and Accurate Artificial Dissipation,” Journal of Scientific Computing, vol. 21, no. 1, pp. 57–79, Aug. 2004"
13 """
7 function undivided_dissipation(g::EquidistantGrid, p, direction) 14 function undivided_dissipation(g::EquidistantGrid, p, direction)
8 T = eltype(g) 15 T = eltype(g)
9 interior_weights = T.(dissipation_interior_weights(p)) 16 interior_weights = T.(dissipation_interior_weights(p))
10 17
11 D = stencil_operator_distinct_closures( 18 D = stencil_operator_distinct_closures(