Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/derivatives/dissipation.jl @ 1344:760a4a1ec4b7 refactor/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add 2D tests for dissipation operators and fix bug
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 15 May 2023 22:37:55 +0200
parents e94ddef5e72f
children 08f06bfacd5c
comparison
equal deleted inserted replaced
1343:fa3695f634de 1344:760a4a1ec4b7
12 57–79, Aug. 2004" 12 57–79, Aug. 2004"
13 """ 13 """
14 function undivided_skewed04 end 14 function undivided_skewed04 end
15 15
16 function undivided_skewed04(g::TensorGrid, p, direction) 16 function undivided_skewed04(g::TensorGrid, p, direction)
17 op = undivided_skewed04(g.grids[direction], p) 17 D,Dᵀ = undivided_skewed04(g.grids[direction], p)
18 return LazyTensors.inflate(op, size(g), direction) 18 return (
19 LazyTensors.inflate(D, size(g), direction),
20 LazyTensors.inflate(Dᵀ, size(g), direction),
21 )
19 end 22 end
20 23
21 function undivided_skewed04(g::EquidistantGrid, p) 24 function undivided_skewed04(g::EquidistantGrid, p)
22 T = eltype(g) 25 T = eltype(g)
23 interior_weights = T.(dissipation_interior_weights(p)) 26 interior_weights = T.(dissipation_interior_weights(p))