From ea21ba73273891261ba2e4d0d85729f308c54d72 Mon Sep 17 00:00:00 2001 From: Pjotr Prins Date: Sat, 26 Aug 2017 07:41:48 +0000 Subject: Tests and enforces added related to https://github.com/genetics-statistics/GEMMA/issues/78 --- src/lapack.h | 29 ++--------------------------- 1 file changed, 2 insertions(+), 27 deletions(-) (limited to 'src/lapack.h') diff --git a/src/lapack.h b/src/lapack.h index 9596667..9115201 100644 --- a/src/lapack.h +++ b/src/lapack.h @@ -23,47 +23,22 @@ using namespace std; -void lapack_float_cholesky_decomp(gsl_matrix_float *A); void lapack_cholesky_decomp(gsl_matrix *A); -void lapack_float_cholesky_solve(gsl_matrix_float *A, const gsl_vector_float *b, - gsl_vector_float *x); void lapack_cholesky_solve(gsl_matrix *A, const gsl_vector *b, gsl_vector *x); -void lapack_sgemm(char *TransA, char *TransB, float alpha, - const gsl_matrix_float *A, const gsl_matrix_float *B, - float beta, gsl_matrix_float *C); void lapack_dgemm(char *TransA, char *TransB, double alpha, const gsl_matrix *A, const gsl_matrix *B, double beta, gsl_matrix *C); -void lapack_float_eigen_symmv(gsl_matrix_float *A, gsl_vector_float *eval, - gsl_matrix_float *evec, - const size_t flag_largematrix); void lapack_eigen_symmv(gsl_matrix *A, gsl_vector *eval, gsl_matrix *evec, const size_t flag_largematrix); - 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); - double CholeskySolve(gsl_matrix *Omega, gsl_vector *Xty, gsl_vector *OiXty); -double CholeskySolve(gsl_matrix_float *Omega, gsl_vector_float *Xty, - gsl_vector_float *OiXty); - void LUDecomp(gsl_matrix *LU, gsl_permutation *p, int *signum); -void LUDecomp(gsl_matrix_float *LU, gsl_permutation *p, int *signum); -void LUInvert(const gsl_matrix *LU, const gsl_permutation *p, - gsl_matrix *inverse); -void LUInvert(const gsl_matrix_float *LU, const gsl_permutation *p, - gsl_matrix_float *inverse); -double LULndet(gsl_matrix *LU); -double LULndet(gsl_matrix_float *LU); +void LUInvert(const gsl_matrix *LU, const gsl_permutation *p, gsl_matrix *ret_inverse); +double LULndet(const gsl_matrix *LU); void LUSolve(const gsl_matrix *LU, const gsl_permutation *p, const gsl_vector *b, gsl_vector *x); -void LUSolve(const gsl_matrix_float *LU, const gsl_permutation *p, - const gsl_vector_float *b, gsl_vector_float *x); - bool lapack_ddot(vector &x, vector &y, double &v); -bool lapack_sdot(vector &x, vector &y, double &v); #endif -- cgit 1.4.1