|
From: Gordon K. <kin...@us...> - 2004-03-30 03:44:04
|
Update of /cvsroot/teem/teem/src/ten In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv27237 Modified Files: tendPoint.c Log Message: added more information to output Index: tendPoint.c =================================================================== RCS file: /cvsroot/teem/teem/src/ten/tendPoint.c,v retrieving revision 1.8 retrieving revision 1.9 diff -C2 -d -r1.8 -r1.9 *** tendPoint.c 25 Mar 2004 03:06:07 -0000 1.8 --- tendPoint.c 30 Mar 2004 03:32:30 -0000 1.9 *************** *** 35,39 **** int loc[3], sx, sy, sz, i; Nrrd *nin; ! float *tdata, eval[3], evec[9], c[TEN_ANISO_MAX+1]; hestOptAdd(&hopt, "p", "x y z", airTypeInt, 3, 3, loc, NULL, --- 35,40 ---- int loc[3], sx, sy, sz, i; Nrrd *nin; ! float *tdata, eval[3], evec[9], c[TEN_ANISO_MAX+1], ! angle, axis[3], mat[9]; hestOptAdd(&hopt, "p", "x y z", airTypeInt, 3, 3, loc, NULL, *************** *** 83,86 **** --- 84,97 ---- fprintf(stderr, "% 15.7f : % 15.7f % 15.7f % 15.7f\n", eval[2], evec[6], evec[7], evec[8]); + angle = ell_3m_to_aa_f(axis, evec); + fprintf(stderr, "eigenvector rotation: %g around {%g,%g,%g}\n", + angle, axis[0], axis[1], axis[2]); + ell_aa_to_3m_f(mat, angle, axis); + fprintf(stderr, "% 15.7f % 15.7f % 15.7f\n", + mat[0], mat[1], mat[2]); + fprintf(stderr, "% 15.7f % 15.7f % 15.7f\n", + mat[3], mat[4], mat[5]); + fprintf(stderr, "% 15.7f % 15.7f % 15.7f\n", + mat[6], mat[7], mat[8]); tenAnisoCalc(c, eval); fprintf(stderr, "anisotropies = \n"); |