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