DLLEXPORT int s_cholesky_factor(int n, float a[]){
int info = clapack_spotrf(CblasColMajor, CblasLower, n, a, n);
for (int i = 0; i < n; ++i)
{
int index = i * n;
for (int j = 0; j < n && i > j; ++j)
{
a[index + j] = 0;
}
}
return info;
}
inline int potrf(
CBLAS_ORDER const Order, CBLAS_UPLO const Uplo,
int const N, float *A, int const lda
) {
return clapack_spotrf(Order, Uplo, N, A, lda);
}
inline std::ptrdiff_t potrf( Order, const UpLo, const int n, float* a,
const int lda ) {
return clapack_spotrf( clapack_option< Order >::value, clapack_option<
UpLo >::value, n, a, lda );
}
