From 4515920b695e01aa6b0a343eb3ffa41887e7fc4d Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Thu, 11 Feb 2021 23:54:43 +0100 Subject: [PATCH] not needed --- src/homogenization.f90 | 20 -------------------- 1 file changed, 20 deletions(-) diff --git a/src/homogenization.f90 b/src/homogenization.f90 index 62ac52f06..14dd26726 100644 --- a/src/homogenization.f90 +++ b/src/homogenization.f90 @@ -311,26 +311,6 @@ 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) elementLooping3: do el = FEsolving_execElem(1),FEsolving_execElem(2) ho = material_homogenizationAt(el)