From c195f316e1a448222fe7a681869438c00dec37e4 Mon Sep 17 00:00:00 2001 From: Franz Roters Date: Mon, 21 Jan 2013 08:17:43 +0000 Subject: [PATCH] corrected calculation of twin stiffness matrix --- code/constitutive_dislotwin.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/code/constitutive_dislotwin.f90 b/code/constitutive_dislotwin.f90 index a9be85065..62bb872de 100644 --- a/code/constitutive_dislotwin.f90 +++ b/code/constitutive_dislotwin.f90 @@ -648,7 +648,7 @@ do i = 1_pInt,maxNinstance lattice_Qtwin(o,s,index_otherFamily+j,myStructure) enddo ; enddo ; enddo ; enddo enddo ; enddo ; enddo ; enddo - constitutive_dislotwin_Ctwin_66(:,:,k,i) = math_Mandel3333to66(constitutive_dislotwin_Ctwin_3333(:,:,:,:,k,i)) + constitutive_dislotwin_Ctwin_66(:,:,index_myFamily+j,i) = math_Mandel3333to66(constitutive_dislotwin_Ctwin_3333(:,:,:,:,index_myFamily+j,i)) !* Interaction matrices