From c36968fc68c58a48b4ed9f972e859f4070655829 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Wed, 30 Mar 2011 11:57:41 +0000 Subject: [PATCH] distance for periodic images in internal stress calculation was incorrect --- code/constitutive_nonlocal.f90 | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90 index ce6663e19..7c738a589 100644 --- a/code/constitutive_nonlocal.f90 +++ b/code/constitutive_nonlocal.f90 @@ -984,9 +984,8 @@ if (.not. phase_localConstitution(phase)) then .and. deltaX == 0 .and. deltaY == 0 .and. deltaZ == 0) then cycle ! this is myself endif - neighboring_ipCoords = mesh_ipCenterOfGravity(1:3,neighboring_ip,neighboring_el) + dble(deltaX) * meshSize(1) & - + dble(deltaY) * meshSize(2) & - + dble(deltaZ) * meshSize(3) + neighboring_ipCoords = mesh_ipCenterOfGravity(1:3,neighboring_ip,neighboring_el) & + + (/dble(deltaX), dble(deltaY), dble(deltaZ)/) * meshSize connection = neighboring_ipCoords - ipCoords distance = sqrt(sum(connection ** 2.0_pReal)) if (.not. phase_localConstitution(neighboring_phase) &