From 97a014667204ff477a676d9bac223a9e14798475 Mon Sep 17 00:00:00 2001 From: Franz Roters Date: Fri, 11 Mar 2011 10:05:30 +0000 Subject: [PATCH] deactivated one debug statement at the end of crystallite_updateState(g,i,e) as it leads to division by zero if constitutive_state(g,i,e)%p(1:mySize) contains zero values statement should be changed to check for zero values --- code/crystallite.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/code/crystallite.f90 b/code/crystallite.f90 index 4d7488e7d..e29cf445a 100644 --- a/code/crystallite.f90 +++ b/code/crystallite.f90 @@ -2243,9 +2243,9 @@ if (verboseDebugger .and. (e == debug_e .and. i == debug_i .and. g == debug_g) . write(6,*) write(6,'(a,/,12(e12.5,x))') 'updateState: new state',constitutive_state(g,i,e)%p(1:mySize) write(6,*) - write(6,'(a,/,12(f12.1,x))') 'updateState: relative residuum tolerance', abs(residuum / rTol_crystalliteState & - / constitutive_state(g,i,e)%p(1:mySize)) - write(6,*) +! write(6,'(a,/,12(f12.1,x))') 'updateState: relative residuum tolerance', abs(residuum / rTol_crystalliteState & +! / constitutive_state(g,i,e)%p(1:mySize)) +! write(6,*) !$OMP END CRITICAL (write2out) endif