From fbdab24f4b9285cf98e114138e507e188b617081 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Mon, 20 Nov 2023 13:27:34 +0100 Subject: [PATCH 1/2] bugfix: used wrong definition looks like a copy and paste error --- src/crystal.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/crystal.f90 b/src/crystal.f90 index 32040bce4..66930d488 100644 --- a/src/crystal.f90 +++ b/src/crystal.f90 @@ -512,13 +512,13 @@ function crystal_C66_twin(Ntwin,C66,lattice,CoverA) select case(lattice) case('cF') - coordinateSystem = buildCoordinateSystem(Ntwin,CF_NSLIPSYSTEM,CF_SYSTEMTWIN,& + coordinateSystem = buildCoordinateSystem(Ntwin,CF_NTWINSYSTEM,CF_SYSTEMTWIN,& lattice,0.0_pREAL) case('cI') - coordinateSystem = buildCoordinateSystem(Ntwin,CI_NSLIPSYSTEM,CI_SYSTEMTWIN,& + coordinateSystem = buildCoordinateSystem(Ntwin,CI_NTWINSYSTEM,CI_SYSTEMTWIN,& lattice,0.0_pREAL) case('hP') - coordinateSystem = buildCoordinateSystem(Ntwin,HP_NSLIPSYSTEM,HP_SYSTEMTWIN,& + coordinateSystem = buildCoordinateSystem(Ntwin,HP_NTWINSYSTEM,HP_SYSTEMTWIN,& lattice,cOverA) case default call IO_error(137,ext_msg='crystal_C66_twin: '//trim(lattice)) From 5db4c3ecf6bd68d574c3bee32d085613fb6e25ff Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Mon, 20 Nov 2023 16:19:20 +0100 Subject: [PATCH 2/2] document --- src/crystal.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/crystal.f90 b/src/crystal.f90 index 66930d488..62679e309 100644 --- a/src/crystal.f90 +++ b/src/crystal.f90 @@ -525,7 +525,7 @@ function crystal_C66_twin(Ntwin,C66,lattice,CoverA) end select do i = 1, sum(Ntwin) - call R%fromAxisAngle([coordinateSystem(1:3,2,i),PI],P=1) ! ToDo: Why always 180 deg? + call R%fromAxisAngle([coordinateSystem(1:3,2,i),PI]) ! mirror on habit (twin shear) plane crystal_C66_twin(1:6,1:6,i) = R%rotStiffness(C66) end do