[Teem-commit] teem/src/ten aniso.c,1.29,1.30

[Gordon K. <kin...@us...>]

Update of /cvsroot/teem/teem/src/ten
In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv8015

Modified Files:
	aniso.c 
Log Message:
fixed bug to avoid NaN generation

Index: aniso.c
===================================================================
RCS file: /cvsroot/teem/teem/src/ten/aniso.c,v
retrieving revision 1.29
retrieving revision 1.30
diff -C2 -d -r1.29 -r1.30
*** aniso.c	22 Apr 2004 04:49:26 -0000	1.29
--- aniso.c	22 Apr 2004 07:13:01 -0000	1.30
***************
*** 102,106 ****
    c[tenAniso_S] = e0*e0 + e1*e1 + e2*e2;
    c[tenAniso_Skew] = R/(FLT_EPSILON + sqrt(2*Q*Q*Q));
!   c[tenAniso_Th] = acos(sqrt(2)*c[tenAniso_Skew])/3;
    c[tenAniso_Cz] = ((e0 + e1)/(FLT_EPSILON + e2) 
  		    + (e1 + e2)/(FLT_EPSILON + e0) 
--- 102,106 ----
    c[tenAniso_S] = e0*e0 + e1*e1 + e2*e2;
    c[tenAniso_Skew] = R/(FLT_EPSILON + sqrt(2*Q*Q*Q));
!   c[tenAniso_Th] = acos(AIR_CLAMP(-1, sqrt(2)*c[tenAniso_Skew], 1))/3;
    c[tenAniso_Cz] = ((e0 + e1)/(FLT_EPSILON + e2) 
  		    + (e1 + e2)/(FLT_EPSILON + e0)