distance for periodic images in internal stress calculation was incorrect

This commit is contained in:
Christoph Kords 2011-03-30 11:57:41 +00:00
parent 08da2b5e4d
commit c36968fc68
1 changed files with 2 additions and 3 deletions

View File

@ -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) &