gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Revert "INTERPOLATIONS: Make sure MC, VL, and SB limiters do not overshoot."

Browse Source 
This reverts commit 99cc1ca2abcf3e52ecfcf6bc66ee051cf034bbce.
This commit is contained in:
Grzegorz Kowal 2015-12-10 11:22:23 -02:00
parent 9cf87ed959
commit 164f809ee5
1 changed files with 4 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
src/interpolations.F90
View File

@ -2954,15 +2954,10 @@ module interpolations
!
real(kind=8), intent(in) :: x, a, b
real(kind=8) :: c
! local variables
!
real(kind=8) :: y
!
!-------------------------------------------------------------------------------
!
y = x - eps
c = (sign(y, a) + sign(y, b)) * min(abs(a), abs(b), 2.5d-01 * abs(a + b))
c = (sign(x, a) + sign(x, b)) * min(abs(a), abs(b), 2.5d-01 * abs(a + b))
!-------------------------------------------------------------------------------
!
@ -2993,16 +2988,11 @@ module interpolations
!
real(kind=8), intent(in) :: x, a, b
real(kind=8) :: c
! local variables
!
real(kind=8) :: y
!
!-------------------------------------------------------------------------------
!
y = x - eps
c = (sign(y, a) + sign(y, b)) &
* max(min(abs(a), 0.5d+00 * abs(b)), min(0.5d+00 * abs(a), abs(b)))
c = 0.5d+00 * (sign(x, a) + sign(x, b)) &
* max(min(2.0d+00 * abs(a), abs(b)), min(abs(a), 2.0d+00 * abs(b)))
!-------------------------------------------------------------------------------
!
@ -3038,7 +3028,7 @@ module interpolations
!
c = a * b
if (c > 0.0d+00) then
c = 2.0d+00 * (x - eps) * c / (a + b)
c = 2.0d+00 * x * c / (a + b)
else
c = 0.0d+00
end if