diff options
Diffstat (limited to 'bench/btl/libs/C_BLAS/C_BLAS_interface.hh')
-rw-r--r-- | bench/btl/libs/C_BLAS/C_BLAS_interface.hh | 10 |
1 files changed, 10 insertions, 0 deletions
diff --git a/bench/btl/libs/C_BLAS/C_BLAS_interface.hh b/bench/btl/libs/C_BLAS/C_BLAS_interface.hh index 21d9bde1a..d0148a29d 100644 --- a/bench/btl/libs/C_BLAS/C_BLAS_interface.hh +++ b/bench/btl/libs/C_BLAS/C_BLAS_interface.hh @@ -253,6 +253,16 @@ public : #endif } + static inline void trisolve_lower_matrix(const gene_matrix & L, const gene_matrix& B, gene_matrix & X, int N){ + #ifdef PUREBLAS + scopy_(&N, B, &intone, X, &intone); + strsv_(&lower, ¬rans, &nonunit, &N, L, &N, X, &intone); + #else + cblas_scopy(N, B, 1, X, 1); + cblas_strsm(CblasColMajor, CblasRight, CblasLower, CblasNoTrans, CblasNonUnit, N, N, 1, L, N, X, N); + #endif + } + }; |