From 066c598203a3bdb3f5e84185dd7051712fa17c0c Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Tue, 29 Jan 2019 10:52:00 +0100 Subject: [PATCH] wrong dot product in state damper --- src/crystallite.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/crystallite.f90 b/src/crystallite.f90 index 8efe15040..ef898bd77 100644 --- a/src/crystallite.f90 +++ b/src/crystallite.f90 @@ -1744,8 +1744,8 @@ subroutine integrateStateFPI() real(pReal) :: dot_prod12, dot_prod22 - dot_prod12 = dot_product(current - previous, previous - previous2) - dot_prod22 = dot_product(current - previous2, previous - previous2) + dot_prod12 = dot_product(current - previous, previous - previous2) + dot_prod22 = dot_product(previous - previous2, previous - previous2) if (dot_prod22 > 0.0_pReal .and. (dot_prod12 < 0.0_pReal .or. dot_product(current,previous) < 0.0_pReal)) then damper = 0.75_pReal + 0.25_pReal * tanh(2.0_pReal + 4.0_pReal * dot_prod12 / dot_prod22) else