From 372072c004f8c4ea1859c71c7eb17ecf3dd9ccb8 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Wed, 24 Nov 2010 12:08:15 +0100
Subject: Display the matrix before solving

---
 src/post-refinement.c | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

diff --git a/src/post-refinement.c b/src/post-refinement.c
index bc1caa4c..f1c59bff 100644
--- a/src/post-refinement.c
+++ b/src/post-refinement.c
@@ -268,6 +268,20 @@ double mean_partial_dev(struct image *image, struct cpeak *spots, int n,
 }
 
 
+static void show_matrix_eqn(gsl_matrix *M, gsl_vector *v, int r)
+{
+	int i, j;
+
+	for ( i=0; i<r; i++ ) {
+		STATUS("[ ");
+		for ( j=0; j<r; j++ ) {
+			STATUS("%+9.3e ", gsl_matrix_get(M, i, j));
+		}
+		STATUS("][ a%2i ] = [ %+9.3e ]\n", i, gsl_vector_get(v, i));
+	}
+}
+
+
 /* Perform one cycle of post refinement on 'image' against 'i_full' */
 double pr_iterate(struct image *image, double *i_full, const char *sym,
                   struct cpeak **pspots, int *n)
@@ -341,6 +355,7 @@ double pr_iterate(struct image *image, double *i_full, const char *sym,
 		}
 
 	}
+	show_matrix_eqn(M, v, NUM_PARAMS);
 
 	shifts = gsl_vector_alloc(NUM_PARAMS);
 	gsl_linalg_HH_solve(M, v, shifts);
-- 
cgit v1.2.3