check for NaN in math_QuaternionInSST

2012-05-11 07:37:32 +00:00 · 2012-05-11 07:37:32 +00:00 · f8f9c089d6
parent babf2647fb
commit f8f9c089d6
1 changed files with 16 additions and 12 deletions
--- a/code/math.f90
+++ b/code/math.f90
@ -1804,6 +1804,9 @@ pure function math_QuaternionInSST(Q, symmetryType)
  real(pReal), dimension(3) ::                  Rodrig                      ! Rodrigues vector of Q
 
  Rodrig = math_QuaternionToRodrig(Q)
+  if (any(Rodrig/=Rodrig)) then
+    math_QuaternionInSST = .false.
+  else
    select case (symmetryType)
      case (1_pInt)
        math_QuaternionInSST = Rodrig(1) > Rodrig(2) .and. &
@ -1816,6 +1819,7 @@ pure function math_QuaternionInSST(Q, symmetryType)
      case default
        math_QuaternionInSST = .true.
    end select
+  endif
  
 end function math_QuaternionInSST