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

quat->Eulers now returns always positive angles

This commit is contained in:
Philip Eisenlohr 2018-12-08 11:52:55 -05:00
parent eaf2c77903
commit b777253dd6
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
lib/damask/orientation.py
View File

@ -263,8 +263,8 @@ class Quaternion:
self.normalize() self.normalize()
s = math.sqrt(1. - self.q**2) s = math.sqrt(1. - self.q**2)
x = 2*self.q**2 - 1. x = 2.*self.q**2 - 1.
y = 2*self.q * s y = 2.*self.q * s
angle = math.atan2(y,x) angle = math.atan2(y,x)
if angle < 0.0: if angle < 0.0:
@ -300,6 +300,7 @@ class Quaternion:
math.atan2((P*self.q*self.p[1]+self.p[0]*self.p[2])/chi,( self.p[1]*self.p[2]-P*self.q*self.p[0])/chi), math.atan2((P*self.q*self.p[1]+self.p[0]*self.p[2])/chi,( self.p[1]*self.p[2]-P*self.q*self.p[0])/chi),
]) ])
eulers %= 2.0*math.pi
return np.degrees(eulers) if degrees else eulers return np.degrees(eulers) if degrees else eulers