inline
int gesv (MatrA& a, MatrB& b) {
// with 'internal' pivot vector
// gesv() errors:
// if (info == 0), successful
// if (info < 0), the -info argument had an illegal value
// -- we will use -101 if allocation fails
// if (info > 0), U(i-1,i-1) is exactly zero
int info = -101;
traits::detail::array<int> ipiv (traits::matrix_size1 (a));
if (ipiv.valid())
info = gesv (a, ipiv, b);
return info;
}
inline
int lu_solve (MatrA& a, MatrB& b) {
return gesv (a, b);
}