agreed on wrong calc of scaled state residual
This commit is contained in:
parent
09495e852d
commit
4718abfc70
|
@ -431,7 +431,6 @@ stress: do ! inner iteration: stress
|
||||||
Rstate = state_new - (state_old+dstate)
|
Rstate = state_new - (state_old+dstate)
|
||||||
RstateS = 0.0_pReal
|
RstateS = 0.0_pReal
|
||||||
forall (i=1:constitutive_Nstatevars(grain,CPFEM_in,cp_en), state_new(i)/=0.0_pReal) &
|
forall (i=1:constitutive_Nstatevars(grain,CPFEM_in,cp_en), state_new(i)/=0.0_pReal) &
|
||||||
! RstateS(i) = Rstress(i)/state_new(i) ???
|
|
||||||
RstateS(i) = Rstate(i)/state_new(i)
|
RstateS(i) = Rstate(i)/state_new(i)
|
||||||
if (maxval(abs(RstateS)) < tol_State) exit state
|
if (maxval(abs(RstateS)) < tol_State) exit state
|
||||||
state_new = state_old+dstate
|
state_new = state_old+dstate
|
||||||
|
|
Loading…
Reference in New Issue