From 99527865c00b74a3a48daa2e1e5eb7c71bd861b5 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Tue, 22 Aug 2017 07:56:07 +0000
Subject: Fixes matrix checks - looking saner now

- Matrix checks as described in https://github.com/genetics-statistics/GEMMA/issues/72
- introduces -strict switch which will exit on certain conditions
- zero small eigenvalues in EigenDecomp_Zeroed which also checks for negative values
- commented out float versions of functions in lapack.cpp (pre-removal)
- reverted on disabled regression tests (GEMMA shows its previous behaviour now)
---
 src/lapack.h | 2 ++
 1 file changed, 2 insertions(+)

(limited to 'src/lapack.h')

diff --git a/src/lapack.h b/src/lapack.h
index ff02b96..9596667 100644
--- a/src/lapack.h
+++ b/src/lapack.h
@@ -41,6 +41,8 @@ void lapack_eigen_symmv(gsl_matrix *A, gsl_vector *eval, gsl_matrix *evec,
 
 double EigenDecomp(gsl_matrix *G, gsl_matrix *U, gsl_vector *eval,
                    const size_t flag_largematrix);
+double EigenDecomp_Zeroed(gsl_matrix *G, gsl_matrix *U, gsl_vector *eval,
+                   const size_t flag_largematrix);
 double EigenDecomp(gsl_matrix_float *G, gsl_matrix_float *U,
                    gsl_vector_float *eval, const size_t flag_largematrix);
 
-- 
cgit v1.2.3