Mercurial > repos > public > sbplib_julia

changeset 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 02d816fc70ef fe298044e044
files src/SbpOperators/volumeops/laplace/laplace.jl
diffstat 1 files changed, 2 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/laplace/laplace.jl	Thu Sep 05 08:36:11 2024 +0200
+++ b/src/SbpOperators/volumeops/laplace/laplace.jl	Thu Sep 05 08:41:02 2024 +0200
@@ -118,8 +118,8 @@
     θ_H = parse_scalar(Δ.stencil_set["H"]["closure"][1])
     θ_R = parse_scalar(Δ.stencil_set["D2"]["positivity"]["theta_R"])
 
-    τ_H = 1/(h*θ_H)
-    τ_R = 1/(h*θ_R)
+    τ_H = one(eltype(Δ))/(h*θ_H)
+    τ_R = one(eltype(Δ))/(h*θ_R)
     return τ_H, τ_R
 end