chiark / gitweb /
before try polar decomposition
authorIan Jackson <ian@davenant.relativity.greenend.org.uk>
Sun, 30 Dec 2007 23:47:55 +0000 (23:47 +0000)
committerIan Jackson <ian@davenant.relativity.greenend.org.uk>
Sun, 30 Dec 2007 23:47:55 +0000 (23:47 +0000)
project.c

index 0c4022e..1fb8576 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();