From eb8265b9149bd6e5c164a23140d0b4532973e450 Mon Sep 17 00:00:00 2001 From: Franz Roters Date: Tue, 8 May 2012 07:16:00 +0000 Subject: [PATCH] fixed potential division by zero error discovered by our American friends --- code/constitutive_nonlocal.f90 | 1 + 1 file changed, 1 insertion(+) diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90 index 87dd3d228..22db3b773 100644 --- a/code/constitutive_nonlocal.f90 +++ b/code/constitutive_nonlocal.f90 @@ -1690,6 +1690,7 @@ endif do s = 1_pInt,ns ! loop over slip systems sLattice = constitutive_nonlocal_slipSystemLattice(s,myInstance) tau(s) = math_mul6x6(Tstar_v, lattice_Sslip_v(1:6,sLattice,myStructure)) + tauBack(s) + if (abs(tau(s)) < 1.0e-15_pReal) tau(s) = 1.0e-15_pReal enddo dLower = constitutive_nonlocal_minimumDipoleHeight(1:ns,1:2,myInstance)