Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/volumeops/laplace/laplace.jl @ 1712:b42a38900c9d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Improve type handling in positivity_limits()
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 05 Sep 2024 08:41:02 +0200
parents 84aed3abab94
children b5690ab5f0b8
comparison
equal deleted inserted replaced
1711:84aed3abab94 1712:b42a38900c9d
116 function positivity_limits(Δ::Laplace, g::EquidistantGrid, bc::DirichletCondition) 116 function positivity_limits(Δ::Laplace, g::EquidistantGrid, bc::DirichletCondition)
117 h = spacing(g) 117 h = spacing(g)
118 θ_H = parse_scalar(Δ.stencil_set["H"]["closure"][1]) 118 θ_H = parse_scalar(Δ.stencil_set["H"]["closure"][1])
119 θ_R = parse_scalar(Δ.stencil_set["D2"]["positivity"]["theta_R"]) 119 θ_R = parse_scalar(Δ.stencil_set["D2"]["positivity"]["theta_R"])
120 120
121 τ_H = 1/(h*θ_H) 121 τ_H = one(eltype(Δ))/(h*θ_H)
122 τ_R = 1/(h*θ_R) 122 τ_R = one(eltype(Δ))/(h*θ_R)
123 return τ_H, τ_R 123 return τ_H, τ_R
124 end 124 end
125 125
126 function positivity_limits(Δ::Laplace, g::TensorGrid, bc::DirichletCondition) 126 function positivity_limits(Δ::Laplace, g::TensorGrid, bc::DirichletCondition)
127 τ_H, τ_R = positivity_limits(Δ, g.grids[grid_id(boundary(bc))], bc) 127 τ_H, τ_R = positivity_limits(Δ, g.grids[grid_id(boundary(bc))], bc)