|
From: Gordon K. <kin...@us...> - 2004-04-22 07:59:55
|
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) |