From a527e38ccc907a64420a2bd73245a0512c8baa87 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Wed, 13 Mar 2019 16:56:54 +0100 Subject: cell_transform_rational: Complain if matrix is singular --- libcrystfel/src/cell.c | 4 ++++ 1 file changed, 4 insertions(+) (limited to 'libcrystfel/src') diff --git a/libcrystfel/src/cell.c b/libcrystfel/src/cell.c index 1e463bf7..ce9d37bb 100644 --- a/libcrystfel/src/cell.c +++ b/libcrystfel/src/cell.c @@ -876,9 +876,13 @@ UnitCell *cell_transform_rational(UnitCell *cell, RationalMatrix *m) gsl_matrix *tm; char ncen; int i, j; + Rational det; if ( m == NULL ) return NULL; + det = rtnl_mtx_det(m); + if ( rtnl_cmp(det, rtnl_zero()) == 0 ) return NULL; + tm = gsl_matrix_alloc(3,3); if ( tm == NULL ) { return NULL; -- cgit v1.2.3