diff --git a/src/phase.f90 b/src/phase.f90
index b7bace81a..098cdebe8 100644
--- a/src/phase.f90
+++ b/src/phase.f90
@@ -436,6 +436,8 @@ subroutine phase_allocateState(state, &
   allocate(state%dotState      (sizeDotState,NEntries), source=0.0_pReal)
 
   allocate(state%deltaState  (sizeDeltaState,NEntries), source=0.0_pReal)
+  state%deltaState2 => state%state(state%offsetDeltaState+1: &
+                                   state%offsetDeltaState+state%sizeDeltaState,:)
 
 end subroutine phase_allocateState
 
diff --git a/src/phase_mechanical.f90 b/src/phase_mechanical.f90
index 04cc4b946..733219d56 100644
--- a/src/phase_mechanical.f90
+++ b/src/phase_mechanical.f90
@@ -81,13 +81,13 @@ submodule(phase) mechanical
 
     module function plastic_dotState(subdt,co,ip,el,ph,en) result(dotState)
       integer, intent(in) :: &
-        co, &                                                                                           !< component-ID of integration point
-        ip, &                                                                                           !< integration point
-        el, &                                                                                           !< element
+        co, &                                                                                       !< constituent
+        ip, &                                                                                       !< integration point
+        el, &                                                                                       !< element
         ph, &
         en
       real(pReal),  intent(in) :: &
-        subdt                                                                                           !< timestep
+        subdt                                                                                       !< timestep
       real(pReal), dimension(plasticState(ph)%sizeDotState) :: &
         dotState
     end function plastic_dotState
diff --git a/src/phase_mechanical_plastic.f90 b/src/phase_mechanical_plastic.f90
index 912aadc03..bb5a9819a 100644
--- a/src/phase_mechanical_plastic.f90
+++ b/src/phase_mechanical_plastic.f90
@@ -415,10 +415,9 @@ module function plastic_deltaState(ph, en) result(broken)
 
       broken = any(IEEE_is_NaN(plasticState(ph)%deltaState(:,en)))
       if (.not. broken) then
-        myOffset = plasticState(ph)%offsetDeltaState
-        mySize   = plasticState(ph)%sizeDeltaState
-        plasticState(ph)%state(myOffset + 1:myOffset + mySize,en) = &
-        plasticState(ph)%state(myOffset + 1:myOffset + mySize,en) + plasticState(ph)%deltaState(1:mySize,en)
+        mySize = plasticState(ph)%sizeDeltaState
+        plasticState(ph)%deltaState2(1:mySize,en) = plasticState(ph)%deltaState2(1:mySize,en) &
+                                                  + plasticState(ph)%deltaState(1:mySize,en)
       end if
 
   end select
diff --git a/src/prec.f90 b/src/prec.f90
index 8de82fee8..61fa141ba 100644
--- a/src/prec.f90
+++ b/src/prec.f90
@@ -45,6 +45,8 @@ module prec
       state, &                                                                                      !< state
       dotState, &                                                                                   !< rate of state change
       deltaState                                                                                    !< increment of state change
+    real(pReal), pointer,     dimension(:,:)  :: &
+      deltaState2
   end type
 
   type, extends(tState) :: tPlasticState