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

Introduced new parameter NaN, not a number. It should be assigned for invalid results, such as currently for misorientation of unknown crystal symmetry.

This commit is contained in:
Claudio Zambaldi 2010-04-19 13:40:22 +00:00
parent cab2261693
commit e8ab495a6c
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
code/math.f90
View File

@ -10,6 +10,7 @@
real(pReal), parameter :: pi = 3.14159265358979323846264338327950288419716939937510_pReal real(pReal), parameter :: pi = 3.14159265358979323846264338327950288419716939937510_pReal
real(pReal), parameter :: inDeg = 180.0_pReal/pi real(pReal), parameter :: inDeg = 180.0_pReal/pi
real(pReal), parameter :: inRad = pi/180.0_pReal real(pReal), parameter :: inRad = pi/180.0_pReal
real(pReal), parameter :: NaN = 0.0_pReal/0.0_pReal ! Not a number
! *** 3x3 Identity *** ! *** 3x3 Identity ***
real(pReal), dimension(3,3), parameter :: math_I3 = & real(pReal), dimension(3,3), parameter :: math_I3 = &
reshape( (/ & reshape( (/ &
@ -171,7 +172,8 @@ subroutine math_misorientation(dQ, Q1, Q2, symmetryType)
dQ = 0.0_pReal dQ = 0.0_pReal
if (symmetryType < 1_pInt .or. symmetryType > 2_pInt) then if (symmetryType < 1_pInt .or. symmetryType > 2_pInt) then
call IO_warning(700) dQ=NaN
!call IO_warning(700)
return return
endif endif