[Teem-commit] teem/src/ten tendPoint.c,1.8,1.9

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

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");