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

EQUATIONS: Prepare the initial guess in the Newton-Raphson subroutines.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2015-02-20 13:29:48 -02:00
parent 22a8711259
commit 8c7d891a4d
1 changed files with 8 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
src/equations.F90
View File

@ -4235,17 +4235,6 @@ module equations
en = u(ien,i) + u(idn,i)
dn = u(idn,i)
! prepare the initial guess using velocity and pressure from the previous step
!
vv = sum(q(ivx:ivz,i) * q(ivx:ivz,i))
vb = sum(q(ivx:ivz,i) * u(ibx:ibz,i))
vm = 1.0d+00 + vv
wm = en + 0.5d+00 * (vb * vb - vm * bb) + q(ipr,i)
! set the initial W to the minimum value
!
w = wm
! find the exact W using an Newton-Ralphson interative method
!
call nr_iterate(mm, bb, mb, en, dn, wm, w, vv)
@ -4700,6 +4689,10 @@ module equations
it = nmax
cn = next
! get the initial guess
!
w = en + pmin
! iterate using the Newton-Raphson method in order to find a root w of the
! function
!
@ -4848,6 +4841,10 @@ module equations
it = nmax
cn = next
! get the initial guess
!
w = en + pmin
! iterate using the Newton-Raphson method in order to find a root w of the
! function
!