diff --git a/src/homogenization.f90 b/src/homogenization.f90
index 48919b46d..9b8c33c2e 100644
--- a/src/homogenization.f90
+++ b/src/homogenization.f90
@@ -313,25 +313,25 @@ subroutine materialpoint_stressAndItsTangent(dt,FEsolving_execIP,FEsolving_execE
     enddo
     !$OMP END DO
 
-    !$OMP DO PRIVATE(ho,ph,ce)
-    do el = FEsolving_execElem(1),FEsolving_execElem(2)
-      if (terminallyIll) continue
-      ho = material_homogenizationAt(el)
-      do ip = FEsolving_execIP(1),FEsolving_execIP(2)
-        ce = (el-1)*discretization_nIPs + ip
-      ! call damage_partition(ce)
-      !  do co = 1, homogenization_Nconstituents(ho)
-      !    ph = material_phaseAt(co,el)
-      !    if (.not. thermal_stress(dt,ph,material_phaseMemberAt(co,ip,el))) then
-      !      if (.not. terminallyIll) &                                                           ! so first signals terminally ill...
-      !        print*, ' Integration point ', ip,' at element ', el, ' terminally ill'
-      !      terminallyIll = .true.                                                                  ! ...and kills all others
-      !   endif
-      !   call thermal_homogenize(ip,el)
-      !  enddo
-      enddo
-    enddo
-    !$OMP END DO
+!    !$OMP DO PRIVATE(ho,ph,ce)
+!    do el = FEsolving_execElem(1),FEsolving_execElem(2)
+!      if (terminallyIll) continue
+!      ho = material_homogenizationAt(el)
+!      do ip = FEsolving_execIP(1),FEsolving_execIP(2)
+!        ce = (el-1)*discretization_nIPs + ip
+!        call damage_partition(ce)
+!        do co = 1, homogenization_Nconstituents(ho)
+!          ph = material_phaseAt(co,el)
+!          if (.not. thermal_stress(dt,ph,material_phaseMemberAt(co,ip,el))) then
+!            if (.not. terminallyIll) &                                                           ! so first signals terminally ill...
+!              print*, ' Integration point ', ip,' at element ', el, ' terminally ill'
+!            terminallyIll = .true.                                                                  ! ...and kills all others
+!         endif
+!         call thermal_homogenize(ip,el)
+!        enddo
+!      enddo
+!    enddo
+!    !$OMP END DO
 
     !$OMP DO PRIVATE(ho)
     elementLooping3: do el = FEsolving_execElem(1),FEsolving_execElem(2)