From 302cf9b6c2614df12915d5e71d031901b02500c7 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sun, 7 Oct 2018 22:16:18 +0200 Subject: [PATCH] bugfix, now inline with code from initializeStructure --- src/lattice.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/lattice.f90 b/src/lattice.f90 index 6e1f90752..401172428 100644 --- a/src/lattice.f90 +++ b/src/lattice.f90 @@ -2338,7 +2338,7 @@ function lattice_nonSchmidMatrix(Nslip,nonSchmidCoefficients,sense) result(nonSc + nonSchmidCoefficients(4) * math_tensorproduct33(normal, normal) if (size(nonSchmidCoefficients)>4) nonSchmidMatrix(1:3,1:3,i) = nonSchmidMatrix(1:3,1:3,i) & + nonSchmidCoefficients(5) * math_tensorproduct33(math_crossproduct(normal, direction), & - math_crossproduct(normal, normal)) + math_crossproduct(normal, direction)) if (size(nonSchmidCoefficients)>5) nonSchmidMatrix(1:3,1:3,i) = nonSchmidMatrix(1:3,1:3,i) & + nonSchmidCoefficients(6) * math_tensorproduct33(direction, direction) enddo