From d7854bc606c02084a6287323cb6c994da977d121 Mon Sep 17 00:00:00 2001 From: Ian Jackson Date: Sun, 30 Dec 2007 23:47:55 +0000 Subject: [PATCH] before try polar decomposition --- project.c | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) 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(); -- 2.30.2