223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

left over "if" caused wrong coordinate systems for hex

This commit is contained in:
Martin Diehl 2018-12-14 06:39:49 +01:00
parent 4edaab6da6
commit 25bd6faf7c
1 changed files with 0 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/lattice.f90
View File

@ -2657,12 +2657,9 @@ function buildCoordinateSystem(active,complete,system,structure,cOverA)
normal = system(4:6,c) normal = system(4:6,c)
case ('hex') case ('hex')
if (cOverA < 1.0_pReal .or. cOverA > 2.0_pReal) &
direction = [ system(1,c)*1.5_pReal, & direction = [ system(1,c)*1.5_pReal, &
(system(1,c)+2.0_pReal*system(2,c))*sqrt(0.75_pReal), & (system(1,c)+2.0_pReal*system(2,c))*sqrt(0.75_pReal), &
system(4,c)*cOverA ] ! direction [uvtw]->[3u/2 (u+2v)*sqrt(3)/2 w*(c/a)]) system(4,c)*cOverA ] ! direction [uvtw]->[3u/2 (u+2v)*sqrt(3)/2 w*(c/a)])
normal = [ system(5,c), & normal = [ system(5,c), &
(system(5,c)+2.0_pReal*system(6,c))/sqrt(3.0_pReal), & (system(5,c)+2.0_pReal*system(6,c))/sqrt(3.0_pReal), &
system(8,c)/cOverA ] ! plane (hkil)->(h (h+2k)/sqrt(3) l/(c/a)) system(8,c)/cOverA ] ! plane (hkil)->(h (h+2k)/sqrt(3) l/(c/a))