diff options
author | Pjotr Prins | 2017-08-02 08:46:58 +0000 |
---|---|---|
committer | Pjotr Prins | 2017-08-02 08:46:58 +0000 |
commit | 3935ba39d30666dd7d4a831155631847c77b70c4 (patch) | |
tree | c45fc682b473618a219e324d5c85b5e1f9361d0c /src/lapack.h | |
parent | 84360c191f418bf8682b35e0c8235fcc3bd19a06 (diff) | |
download | pangemma-3935ba39d30666dd7d4a831155631847c77b70c4.tar.gz |
Massive patch using LLVM coding style. It was generated with:
clang-format -style=LLVM -i *.cpp *.h
Please set your editor to replace tabs with spaces and use indentation of 2 spaces.
Diffstat (limited to 'src/lapack.h')
-rw-r--r-- | src/lapack.h | 66 |
1 files changed, 32 insertions, 34 deletions
diff --git a/src/lapack.h b/src/lapack.h index 5e1db35..ff02b96 100644 --- a/src/lapack.h +++ b/src/lapack.h @@ -23,45 +23,43 @@ 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); +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 (gsl_matrix_float *G, gsl_matrix_float *U, - gsl_vector_float *eval, const size_t flag_largematrix); +double EigenDecomp(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); + 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 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); +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 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<double> &x, vector<double> &y, double &v); bool lapack_sdot(vector<float> &x, vector<float> &y, double &v); |