inline std::ptrdiff_t hesv( const UpLo, const fortran_int_t n,
const fortran_int_t nrhs, std::complex<double>* a,
const fortran_int_t lda, fortran_int_t* ipiv, std::complex<double>* b,
const fortran_int_t ldb, std::complex<double>* work,
const fortran_int_t lwork ) {
fortran_int_t info(0);
LAPACK_ZHESV( &lapack_option< UpLo >::value, &n, &nrhs, a, &lda, ipiv, b,
&ldb, work, &lwork, &info );
return info;
}
inline
void hesv (char const uplo, int const n, int const nrhs,
traits::complex_d* a, int const lda, int* ipiv,
traits::complex_d* b, int const ldb,
traits::complex_d* w, int const lw, int* info)
{
LAPACK_ZHESV (&uplo, &n, &nrhs,
traits::complex_ptr (a), &lda, ipiv,
traits::complex_ptr (b), &ldb,
traits::complex_ptr (w), &lw, info);
}