From: Ian Jackson <ian@davenant.relativity.greenend.org.uk>
Date: Sun, 30 Dec 2007 23:47:55 +0000 (+0000)
Subject: before try polar decomposition
X-Git-Url: http://www.chiark.greenend.org.uk/ucgi/~ian/git?p=moebius2.git;a=commitdiff_plain;h=d7854bc606c02084a6287323cb6c994da977d121;ds=sidebyside

before try polar decomposition
---

diff --git a/project.c b/project.c
index 0c4022e..1fb8576 100644
--- 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();