inline std::ptrdiff_t hptrf( const UpLo, const fortran_int_t n, double* ap, fortran_int_t* ipiv ) { fortran_int_t info(0); LAPACK_DSPTRF( &lapack_option< UpLo >::value, &n, ap, ipiv, &info ); return info; }
inline void sptrf (char const uplo, int const n, double* ap, int* ipiv, int* info) { LAPACK_DSPTRF (&uplo, &n, ap, ipiv, info); }