diff --git a/src/crystallite.f90 b/src/crystallite.f90
index 0595aa197..53ddca1c7 100644
--- a/src/crystallite.f90
+++ b/src/crystallite.f90
@@ -1654,29 +1654,7 @@ subroutine integrateStateFPI()
    !$OMP END PARALLEL DO
 
    call update_dependentState
-   !$OMP PARALLEL
-
-   ! --- STRESS INTEGRATION ---
-
-#ifdef DEBUG
-   if (iand(debug_level(debug_crystallite), debug_levelExtensive) /= 0_pInt) &
-     write(6,'(a,i8,a)') '<< CRYST >> ', count(crystallite_todo(:,:,:)),' grains todo before stress integration'
-#endif
-   !$OMP DO
-     do e = eIter(1),eIter(2); do i = iIter(1,e),iIter(2,e); do g = gIter(1,e),gIter(2,e)                  ! iterate over elements, ips and grains
-       !$OMP FLUSH(crystallite_todo)
-       if (crystallite_todo(g,i,e) .and. .not. crystallite_converged(g,i,e)) then
-         crystallite_todo(g,i,e) = integrateStress(g,i,e)
-         !$OMP FLUSH(crystallite_todo)
-         if (.not. crystallite_todo(g,i,e) .and. .not. crystallite_localPlasticity(g,i,e)) then            ! broken non-local...
-           !$OMP CRITICAL (checkTodo)
-             crystallite_todo = crystallite_todo .and. crystallite_localPlasticity                         ! ... then all non-locals skipped
-           !$OMP END CRITICAL (checkTodo)
-         endif
-       endif
-     enddo; enddo; enddo
-   !$OMP ENDDO
-!$OMP END PARALLEL
+   call update_stress(1.0_pReal)
    call update_dotState(1.0_pReal)
 !$OMP PARALLEL
    ! --- UPDATE STATE  ---
@@ -1923,36 +1901,14 @@ eIter = FEsolving_execElem(1:2)
  call update_State(1.0_pReal)
  call update_deltaState
  call update_dependentState
+ call update_stress(1.0_pReal)
 
 
- !$OMP PARALLEL
- ! --- STRESS INTEGRATION ---
-
- !$OMP DO
-   do e = eIter(1),eIter(2); do i = iIter(1,e),iIter(2,e); do g = gIter(1,e),gIter(2,e)                    ! iterate over elements, ips and grains
-     !$OMP FLUSH(crystallite_todo)
-     if (crystallite_todo(g,i,e) .and. .not. crystallite_converged(g,i,e)) then
-       crystallite_todo(g,i,e) = integrateStress(g,i,e)
-       !$OMP FLUSH(crystallite_todo)
-       if (.not. crystallite_todo(g,i,e) .and. .not. crystallite_localPlasticity(g,i,e)) then                  ! if broken non-local...
-         !$OMP CRITICAL (checkTodo)
-           crystallite_todo = crystallite_todo .and. crystallite_localPlasticity                           ! ...all non-locals skipped
-         !$OMP END CRITICAL (checkTodo)
-       endif
-     endif
-   enddo; enddo; enddo
- !$OMP ENDDO
-
-
- ! --- SET CONVERGENCE FLAG ---
-
- !$OMP DO
+ !$OMP PARALLEL DO 
    do e = eIter(1),eIter(2); do i = iIter(1,e),iIter(2,e); do g = gIter(1,e),gIter(2,e)                    ! iterate over elements, ips and grains
      crystallite_converged(g,i,e) = crystallite_todo(g,i,e) .or. crystallite_converged(g,i,e)              ! if still "to do" then converged per definitionem
    enddo; enddo; enddo
- !$OMP ENDDO
-
- !$OMP END PARALLEL
+ !$OMP END PARALLEL DO 
 
 
  ! --- CHECK NON-LOCAL CONVERGENCE ---
@@ -2088,27 +2044,8 @@ subroutine integrateStateAdaptiveEuler()
  !$OMP END PARALLEL
   call update_deltaState
   call update_dependentState
-
-
- ! --- STRESS INTEGRATION (EULER INTEGRATION) ---
-
- !$OMP PARALLEL DO
-   do e = eIter(1),eIter(2); do i = iIter(1,e),iIter(2,e); do g = gIter(1,e),gIter(2,e)                    ! iterate over elements, ips and grains
-     !$OMP FLUSH(crystallite_todo)
-     if (crystallite_todo(g,i,e)) then
-       crystallite_todo(g,i,e) = integrateStress(g,i,e)
-       !$OMP FLUSH(crystallite_todo)
-       if (.not. crystallite_todo(g,i,e) .and. .not. crystallite_localPlasticity(g,i,e)) then              ! if broken non-local...
-         !$OMP CRITICAL (checkTodo)
-           crystallite_todo = crystallite_todo .and. crystallite_localPlasticity                           ! ...all non-locals skipped
-         !$OMP END CRITICAL (checkTodo)
-       endif
-     endif
-   enddo; enddo; enddo
- !$OMP END PARALLEL DO
-
-
-   call update_dotState(1.0_pReal)
+  call update_stress(1.0_pReal)
+  call update_dotState(1.0_pReal)
    
    !$OMP PARALLEL
    !$OMP DO PRIVATE(p,c,NaN)
@@ -2325,26 +2262,7 @@ subroutine integrateStateRK4()
    call update_state(TIMESTEPFRACTION(n))
    call update_deltaState
    call update_dependentState
-   !$OMP PARALLEL
-
-
-   ! --- stress integration ---
-
-   !$OMP DO
-     do e = eIter(1),eIter(2); do i = iIter(1,e),iIter(2,e); do g = gIter(1,e),gIter(2,e)                  ! iterate over elements, ips and grains
-       !$OMP FLUSH(crystallite_todo)
-       if (crystallite_todo(g,i,e)) then
-         crystallite_todo(g,i,e) = integrateStress(g,i,e,timeStepFraction(n))                  ! fraction of original times step
-         !$OMP FLUSH(crystallite_todo)
-         if (.not. crystallite_todo(g,i,e) .and. .not. crystallite_localPlasticity(g,i,e)) then            ! if broken non-local...
-           !$OMP CRITICAL (checkTodo)
-             crystallite_todo = crystallite_todo .and. crystallite_localPlasticity                         ! ...all non-locals skipped
-           !$OMP END CRITICAL (checkTodo)
-         endif
-       endif
-     enddo; enddo; enddo
-   !$OMP ENDDO
-!$OMP END PARALLEL
+   call update_stress(TIMESTEPFRACTION(n))
 
    ! --- dot state and RK dot state---
 
@@ -2526,26 +2444,8 @@ subroutine integrateStateRKCK45()
     call update_state(1.0_pReal) !MD: 1.0 correct?
     call update_deltaState
     call update_dependentState
-   !$OMP PARALLEL
-
-   ! --- stress integration ---
-
-   !$OMP DO
-     do e = eIter(1),eIter(2); do i = iIter(1,e),iIter(2,e); do g = gIter(1,e),gIter(2,e)                  ! iterate over elements, ips and grains
-       !$OMP FLUSH(crystallite_todo)
-       if (crystallite_todo(g,i,e)) then
-         crystallite_todo(g,i,e) =  integrateStress(g,i,e,C(stage))                            ! fraction of original time step
-         !$OMP FLUSH(crystallite_todo)
-         if (.not. crystallite_todo(g,i,e) .and. .not. crystallite_localPlasticity(g,i,e)) then            ! if broken non-local...
-           !$OMP CRITICAL (checkTodo)
-             crystallite_todo = crystallite_todo .and. crystallite_localPlasticity                         ! ...all non-locals skipped
-           !$OMP END CRITICAL (checkTodo)
-         endif
-       endif
-     enddo; enddo; enddo
-   !$OMP ENDDO
-   !$OMP END PARALLEL
-   call update_dotState(C(stage))
+    call update_stress(C(stage))
+    call update_dotState(C(stage))
 
  enddo
 
@@ -2652,31 +2552,13 @@ subroutine integrateStateRKCK45()
      endif
    enddo; enddo; enddo
  !$OMP ENDDO
-
 !$OMP END PARALLEL
 
  call update_deltaState
  call update_dependentState
-
-   !$OMP PARALLEL
-!--------------------------------------------------------------------------------------------------
-! --- FINAL STRESS INTEGRATION STEP IF RESIDUUM BELOW TOLERANCE ---
- !$OMP DO
-   do e = eIter(1),eIter(2); do i = iIter(1,e),iIter(2,e); do g = gIter(1,e),gIter(2,e)                    ! iterate over elements, ips and grains
-     !$OMP FLUSH(crystallite_todo)
-     if (crystallite_todo(g,i,e)) then
-       crystallite_todo(g,i,e) = integrateStress(g,i,e)
-       !$OMP FLUSH(crystallite_todo)
-       if (.not. crystallite_todo(g,i,e) .and. .not. crystallite_localPlasticity(g,i,e)) then              ! if broken non-local...
-         !$OMP CRITICAL (checkTodo)
-           crystallite_todo = crystallite_todo .and. crystallite_localPlasticity                           ! ...all non-locals skipped
-         !$OMP END CRITICAL (checkTodo)
-       endif
-     endif
-   enddo; enddo; enddo
- !$OMP ENDDO
-
-
+ call update_stress(1.0_pReal)
+ 
+  !$OMP PARALLEL
 !--------------------------------------------------------------------------------------------------
 ! --- SET CONVERGENCE FLAG ---
  !$OMP DO
@@ -2699,6 +2581,43 @@ subroutine integrateStateRKCK45()
 end subroutine integrateStateRKCK45
 
 
+!--------------------------------------------------------------------------------------------------
+!> @brief Standard forwarding of state as state = state0 + dotState * (delta t)
+!--------------------------------------------------------------------------------------------------
+subroutine update_stress(timeFraction)
+ use material, only: &
+   plasticState, &
+   sourceState, &
+   phase_Nsources, &
+   phaseAt, phasememberAt
+
+ implicit none
+ real(pReal), intent(in) :: &
+   timeFraction
+ integer(pInt) :: &
+   e, &                                                                                             !< element index in element loop
+   i, &                                                                                             !< integration point index in ip loop
+   g
+
+ !$OMP PARALLEL DO
+   do e = FEsolving_execElem(1),FEsolving_execElem(2)
+     do i = FEsolving_execIP(1,e),FEsolving_execIP(2,e)
+       do g = 1,homogenization_Ngrains(mesh_element(3,e))
+     !$OMP FLUSH(crystallite_todo)
+     if (crystallite_todo(g,i,e) .and. .not. crystallite_converged(g,i,e)) then
+       crystallite_todo(g,i,e) = integrateStress(g,i,e,timeFraction)
+       !$OMP FLUSH(crystallite_todo)
+       if (.not. crystallite_todo(g,i,e) .and. .not. crystallite_localPlasticity(g,i,e)) then                  ! if broken non-local...
+         !$OMP CRITICAL (checkTodo)
+           crystallite_todo = crystallite_todo .and. crystallite_localPlasticity                           ! ...all non-locals skipped
+         !$OMP END CRITICAL (checkTodo)
+       endif
+     endif
+   enddo; enddo; enddo
+ !$OMP END PARALLEL DO
+
+end subroutine update_stress
+
 !--------------------------------------------------------------------------------------------------
 !> @brief tbd
 !--------------------------------------------------------------------------------------------------
@@ -2725,6 +2644,7 @@ subroutine update_dependentState()
 
 end subroutine update_dependentState
 
+
 !--------------------------------------------------------------------------------------------------
 !> @brief Standard forwarding of state as state = state0 + dotState * (delta t)
 !--------------------------------------------------------------------------------------------------