From 4f7c412d650e8493e0a356f4159be650cffb4d0b Mon Sep 17 00:00:00 2001 From: Thomas White Date: Fri, 7 Jul 2023 15:02:45 +0200 Subject: Factorise matrix operations This makes the code much clearer. Note that two opposing sign errors have been fixed in the gradient calculation. --- libcrystfel/src/utils.h | 3 +++ 1 file changed, 3 insertions(+) (limited to 'libcrystfel/src/utils.h') diff --git a/libcrystfel/src/utils.h b/libcrystfel/src/utils.h index 82a9aa1c..6d2ff253 100644 --- a/libcrystfel/src/utils.h +++ b/libcrystfel/src/utils.h @@ -77,6 +77,9 @@ extern void show_matrix(gsl_matrix *M); extern void show_vector(gsl_vector *M); extern gsl_vector *solve_svd(gsl_vector *v, gsl_matrix *M, int *n_filt, int verbose); +extern gsl_matrix *matrix_mult2(gsl_matrix *A, gsl_matrix *B); +extern gsl_matrix *matrix_mult3(gsl_matrix *A, gsl_matrix *B, gsl_matrix *C); +extern gsl_matrix *matrix_invert(gsl_matrix *m); extern size_t notrail(char *s); extern int convert_int(const char *str, int *pval); -- cgit v1.2.3