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

no need to store relative residual pointwise

This commit is contained in:
Martin Diehl 2019-01-30 11:04:49 +01:00
parent fd069a96cd
commit 3dd21177a0
1 changed files with 38 additions and 38 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

76
src/crystallite.f90
View File

@ -1890,8 +1890,6 @@ end subroutine integrateStateAdaptiveEuler
! ToDo: This is totally BROKEN: RK4dotState is never used!!! ! ToDo: This is totally BROKEN: RK4dotState is never used!!!
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine integrateStateRK4() subroutine integrateStateRK4()
use, intrinsic :: &
IEEE_arithmetic
use mesh, only: & use mesh, only: &
mesh_element mesh_element
use material, only: & use material, only: &
@ -1960,8 +1958,8 @@ end subroutine integrateStateRK4
!> adaptive step size (use 5th order solution to advance = "local extrapolation") !> adaptive step size (use 5th order solution to advance = "local extrapolation")
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine integrateStateRKCK45() subroutine integrateStateRKCK45()
use, intrinsic :: & use prec, only: &
IEEE_arithmetic dNeq0
use numerics, only: & use numerics, only: &
rTol_crystalliteState rTol_crystalliteState
use mesh, only: & use mesh, only: &
@ -2005,26 +2003,25 @@ subroutine integrateStateRKCK45()
i, & ! integration point index in ip loop i, & ! integration point index in ip loop
g, & ! grain index in grain loop g, & ! grain index in grain loop
stage, & ! stage index in integration stage loop stage, & ! stage index in integration stage loop
u, & ! state index
n, & n, &
p, & p, &
cc, & cc, &
s, & s, &
sizeDotState sizeDotState
! ToDo: MD: once all constitutives use allocate state, attach residuum arrays to the state in case of adaptive Euler ! ToDo: MD: once all constitutives use allocate state, attach residuum arrays to the state in case of RKCK45
! ToDo: MD: rel residuu don't have to be pointwise
real(pReal), dimension(constitutive_plasticity_maxSizeDotState, & real(pReal), dimension(constitutive_plasticity_maxSizeDotState, &
homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems) :: & homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems) :: &
residuum_plastic, & ! residuum from evolution in microstructure residuum_plastic ! relative residuum from evolution in microstructure
residuum_plastic_rel ! relative residuum from evolution in microstructure
real(pReal), dimension(constitutive_source_maxSizeDotState, & real(pReal), dimension(constitutive_source_maxSizeDotState, &
maxval(phase_Nsources), & maxval(phase_Nsources), &
homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems) :: & homogenization_maxNgrains,mesh_maxNips,mesh_NcpElems) :: &
residuum_source, & ! residuum from evolution in microstructure residuum_source ! relative residuum from evolution in microstructure
residuum_source_rel ! relative residuum from evolution in microstructure real(pReal), dimension(constitutive_plasticity_maxSizeDotState) :: &
residuum_plastic_rel
real(pReal), dimension(constitutive_source_maxSizeDotState) :: &
residuum_source_rel
call update_dotState(1.0_pReal) call update_dotState(1.0_pReal)
@ -2076,8 +2073,6 @@ subroutine integrateStateRKCK45()
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! --- STATE UPDATE WITH ERROR ESTIMATE FOR STATE --- ! --- STATE UPDATE WITH ERROR ESTIMATE FOR STATE ---
residuum_plastic_rel = 0.0_pReal
residuum_source_rel = 0.0_pReal
!$OMP PARALLEL DO PRIVATE(sizeDotState,p,cc) !$OMP PARALLEL DO PRIVATE(sizeDotState,p,cc)
do e = FEsolving_execElem(1),FEsolving_execElem(2) do e = FEsolving_execElem(1),FEsolving_execElem(2)
do i = FEsolving_execIP(1,e),FEsolving_execIP(2,e) do i = FEsolving_execIP(1,e),FEsolving_execIP(2,e)
@ -2116,43 +2111,48 @@ subroutine integrateStateRKCK45()
call update_state(1.0_pReal) call update_state(1.0_pReal)
!$OMP PARALLEL
! --- relative residui and state convergence --- ! --- relative residui and state convergence ---
!$OMP DO PRIVATE(sizeDotState,p,cc,u) !$OMP PARALLEL DO PRIVATE(sizeDotState,p,cc)
do e = FEsolving_execElem(1),FEsolving_execElem(2) do e = FEsolving_execElem(1),FEsolving_execElem(2)
do i = FEsolving_execIP(1,e),FEsolving_execIP(2,e) do i = FEsolving_execIP(1,e),FEsolving_execIP(2,e)
do g = 1,homogenization_Ngrains(mesh_element(3,e)) do g = 1,homogenization_Ngrains(mesh_element(3,e))
if (crystallite_todo(g,i,e)) then if (crystallite_todo(g,i,e)) then
p = phaseAt(g,i,e); cc = phasememberAt(g,i,e) p = phaseAt(g,i,e); cc = phasememberAt(g,i,e)
sizeDotState = plasticState(p)%sizeDotState sizeDotState = plasticState(p)%sizeDotState
forall (u = 1_pInt:sizeDotState, abs(plasticState(p)%state(u,cc)) > 0.0_pReal) & where(dNeq0(plasticState(p)%dotState(1:sizeDotState,cc)))
residuum_plastic_rel(u,g,i,e) = & residuum_plastic_rel(1:sizeDotState) = residuum_plastic(1:sizeDotState,g,i,e) &
residuum_plastic(u,g,i,e) / plasticState(p)%state(u,cc) / plasticState(p)%state(1:sizeDotState,cc)
else where
residuum_plastic_rel(1:sizeDotState) = 0.0_pReal
end where
crystallite_todo(g,i,e) = all(abs(residuum_plastic_rel(1:sizeDotState,g,i,e)) < &
rTol_crystalliteState .or. &
abs(residuum_plastic(1:sizeDotState,g,i,e)) < &
plasticState(p)%aTolState(1:sizeDotState))
do s = 1_pInt, phase_Nsources(p) crystallite_todo(g,i,e) = all(abs(residuum_plastic_rel(1:sizeDotState)) < &
sizeDotState = sourceState(p)%p(s)%sizeDotState rTol_crystalliteState .or. &
forall (u = 1_pInt:sizeDotState,abs(sourceState(p)%p(s)%state(u,cc)) > 0.0_pReal) & abs(residuum_plastic(1:sizeDotState,g,i,e)) < &
residuum_source_rel(u,s,g,i,e) = & plasticState(p)%aTolState(1:sizeDotState))
residuum_source(u,s,g,i,e) / sourceState(p)%p(s)%state(u,cc)
do s = 1_pInt, phase_Nsources(p)
sizeDotState = sourceState(p)%p(s)%sizeDotState
where(dNeq0(sourceState(p)%p(s)%dotState(1:sizeDotState,cc)))
residuum_source_rel(1:sizeDotState) = residuum_source(1:sizeDotState,s,g,i,e) &
/ sourceState(p)%p(s)%state(1:sizeDotState,cc)
else where
residuum_source_rel(1:SizeDotState) = 0.0_pReal
end where
sizeDotState = sourceState(p)%p(s)%sizeDotState
crystallite_todo(g,i,e) = crystallite_todo(g,i,e) .and. & crystallite_todo(g,i,e) = crystallite_todo(g,i,e) .and. &
all(abs(residuum_source_rel(1:sizeDotState,s,g,i,e)) < & all(abs(residuum_source_rel(1:sizeDotState)) < &
rTol_crystalliteState .or. & rTol_crystalliteState .or. &
abs(residuum_source(1:sizeDotState,s,g,i,e)) < & abs(residuum_source(1:sizeDotState,s,g,i,e)) < &
sourceState(p)%p(s)%aTolState(1:sizeDotState)) sourceState(p)%p(s)%aTolState(1:sizeDotState))
enddo enddo
endif endif
enddo; enddo; enddo enddo; enddo; enddo
!$OMP ENDDO !$OMP END PARALLEL DO
!$OMP END PARALLEL
call update_deltaState call update_deltaState
call update_dependentState call update_dependentState