setting todo and converged flag only for those e,i,g that are to be executed (FE_elemExec etc)

used to spoil convergence checking in non ping pong mode and lead to slow performance
This commit is contained in:
Christoph Kords 2014-01-21 16:09:00 +00:00
parent 1d26cff024
commit 351aa6ba29
1 changed files with 11 additions and 8 deletions

View File

@ -1260,12 +1260,16 @@ subroutine crystallite_stressAndItsTangent(updateJaco,rate_sensitivity)
! --- PERTURB EITHER FORWARD OR BACKWARD --- ! --- PERTURB EITHER FORWARD OR BACKWARD ---
crystallite_subF = F_backup do e = FEsolving_execElem(1),FEsolving_execElem(2)
crystallite_subF(k,l,:,:,:) = crystallite_subF(k,l,:,:,:) + myPert myNgrains = homogenization_Ngrains(mesh_element(3,e))
do i = FEsolving_execIP(1,e),FEsolving_execIP(2,e)
crystallite_converged = convergenceFlag_backup do g = 1,myNgrains
crystallite_todo = crystallite_requested .and. crystallite_converged crystallite_subF(1:3,1:3,g,i,e) = F_backup(1:3,1:3,g,i,e)
where (crystallite_todo) crystallite_converged = .false. ! start out non-converged crystallite_subF(k,l,g,i,e) = crystallite_subF(k,l,g,i,e) + myPert
crystallite_todo(g,i,e) = crystallite_requested(g,i,e) &
.and. convergenceFlag_backup(g,i,e)
if (crystallite_todo(g,i,e)) crystallite_converged(g,i,e) = .false. ! start out non-converged
enddo; enddo; enddo
select case(numerics_integrator(numerics_integrationMode)) select case(numerics_integrator(numerics_integrationMode))
case(1_pInt) case(1_pInt)
@ -3035,8 +3039,7 @@ logical function crystallite_integrateStress(&
dR_dLp2 ! working copy of dRdLp dR_dLp2 ! working copy of dRdLp
real(pReal), dimension(3,3,3,3):: dT_dFe3333, & ! partial derivative of 2nd Piola-Kirchhoff stress real(pReal), dimension(3,3,3,3):: dT_dFe3333, & ! partial derivative of 2nd Piola-Kirchhoff stress
dFe_dLp3333 ! partial derivative of elastic deformation gradient dFe_dLp3333 ! partial derivative of elastic deformation gradient
real(pReal) p_hydro, & ! volumetric part of 2nd Piola-Kirchhoff Stress real(pReal) det, & ! determinant
det, & ! determinant
steplength0, & steplength0, &
steplength, & steplength, &
dt, & ! time increment dt, & ! time increment