223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

just figured out that RK4 integrator is totally broken 

readable code helps ;)
This commit is contained in:
Martin Diehl 2019-01-30 00:17:04 +01:00
parent a09036ff48
commit 77f1f45c23
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/crystallite.f90
View File

@ -1888,6 +1888,7 @@ end subroutine integrateStateAdaptiveEuler
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief integrate stress, state with 4th order explicit Runge Kutta method !> @brief integrate stress, state with 4th order explicit Runge Kutta method
! ToDo: This is totally BROKEN: RK4dotState is never used!!!
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine integrateStateRK4() subroutine integrateStateRK4()
use, intrinsic :: & use, intrinsic :: &
@ -1941,7 +1942,6 @@ subroutine integrateStateRK4()
call update_deltaState call update_deltaState
call update_dependentState call update_dependentState
call update_stress(TIMESTEPFRACTION(n)) call update_stress(TIMESTEPFRACTION(n))
! --- dot state and RK dot state--- ! --- dot state and RK dot state---
first3steps: if (n < 4) then first3steps: if (n < 4) then