diff --git a/VERSION b/VERSION index 5190c13c9..62377fd35 100644 --- a/VERSION +++ b/VERSION @@ -1 +1 @@ -v3.0.0-alpha2-329-gda93031cb +v3.0.0-alpha2-355-gbd9bb5c7c 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)