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

IEEE infinite causes problems with older GNU

This commit is contained in:
Martin Diehl 2019-09-22 14:46:27 -07:00
parent 171a93ca30
commit 4ff292ba70
1 changed files with 9 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/rotations.f90
View File

@ -952,10 +952,10 @@ pure function ro2ax(ro) result(ax)
ta = ro(4) ta = ro(4)
if (dEq0(ta)) then if (.not. IEEE_is_finite(ta)) then
ax = [ 0.0, 0.0, 1.0, 0.0 ]
elseif (.not. IEEE_is_finite(ta)) then
ax = [ ro(1), ro(2), ro(3), PI ] ax = [ ro(1), ro(2), ro(3), PI ]
elseif (dEq0(ta)) then
ax = [ 0.0, 0.0, 1.0, 0.0 ]
else else
angle = 2.0*atan(ta) angle = 2.0*atan(ta)
ta = 1.0/norm2(ro(1:3)) ta = 1.0/norm2(ro(1:3))
@ -1210,6 +1210,10 @@ subroutine unitTest
msg = '' msg = ''
#if defined(__GFORTRAN__) && __GNUC__<9
if(i>2 .and. i<7) cycle
#endif
if(i==1) then if(i==1) then
qu = om2qu(math_I3) qu = om2qu(math_I3)
elseif(i==2) then elseif(i==2) then
@ -1219,7 +1223,7 @@ subroutine unitTest
elseif(i==4) then elseif(i==4) then
qu = [0.0_pReal,0.0_pReal,1.0_pReal,0.0_pReal] qu = [0.0_pReal,0.0_pReal,1.0_pReal,0.0_pReal]
elseif(i==5) then elseif(i==5) then
qu = ro2qu([1.0_pReal,0.0_pReal,0.0_pReal,ieee_value(1.0_pReal, IEEE_positive_inf)]) qu = ro2qu([1.0_pReal,0.0_pReal,0.0_pReal,IEEE_value(1.0_pReal, IEEE_positive_inf)])
elseif(i==6) then elseif(i==6) then
qu = ax2qu([1.0_pReal,0.0_pReal,0.0_pReal,0.0_pReal]) qu = ax2qu([1.0_pReal,0.0_pReal,0.0_pReal,0.0_pReal])
else else
@ -1232,7 +1236,7 @@ subroutine unitTest
sin(2.0_pReal*PI*r(1))*A] sin(2.0_pReal*PI*r(1))*A]
if(qu(1)<0.0_pReal) qu = qu * (-1.0_pReal) if(qu(1)<0.0_pReal) qu = qu * (-1.0_pReal)
endif endif
if(dNeq0(norm2(om2qu(qu2om(qu))-qu),1.0e-12_pReal)) msg = trim(msg)//'om2qu/qu2om,' if(dNeq0(norm2(om2qu(qu2om(qu))-qu),1.0e-12_pReal)) msg = trim(msg)//'om2qu/qu2om,'
if(dNeq0(norm2(eu2qu(qu2eu(qu))-qu),1.0e-12_pReal)) msg = trim(msg)//'eu2qu/qu2eu,' if(dNeq0(norm2(eu2qu(qu2eu(qu))-qu),1.0e-12_pReal)) msg = trim(msg)//'eu2qu/qu2eu,'
if(dNeq0(norm2(ax2qu(qu2ax(qu))-qu),1.0e-12_pReal)) msg = trim(msg)//'ax2qu/qu2ax,' if(dNeq0(norm2(ax2qu(qu2ax(qu))-qu),1.0e-12_pReal)) msg = trim(msg)//'ax2qu/qu2ax,'