before try polar decomposition

diff --git a/project.c b/project.c
index 0c4022eb141d72f755ed3fc776d8833668967676..1fb857630684fdaf4faf62b73dc42d3d1a7614fc 100644 (file)
--- a/project.c
+++ b/project.c
@@ -262,14 +262,17 @@ static void drag_rotate_delta(double dx, double dy) {
   GA( gsl_blas_dgemm(CblasNoTrans,CblasNoTrans, 1.0,
                     &rotateby_gsl,&transform_gsl, 0.0,&qr_gsl) );
 
-  PMATRIX(rotateby);
+  pmatrix("input", qr);
 
   GA( gsl_linalg_QR_decomp(&qr_gsl, &tau_gsl) );
+
+  pmatrix("mangled", qr);
+
   GA( gsl_linalg_QR_unpack(&qr_gsl, &tau_gsl,
                           &transform_gsl, &rotateby_gsl /*dummy*/) );
 
-  PMATRIX(transform);
-  PMATRIX(rotateby);
+  pmatrix("Q", transform);
+  pmatrix("R", rotateby);
 
   printf("drag_rotate_delta...\n");
   show();