//
// Overloaded function for dispatching to
// * netlib-compatible LAPACK backend (the default), and
// * complex<float> value-type.
//
inline std::ptrdiff_t geqrf( const fortran_int_t m, const fortran_int_t n,
std::complex<float>* a, const fortran_int_t lda,
std::complex<float>* tau, std::complex<float>* work,
const fortran_int_t lwork ) {
fortran_int_t info(0);
LAPACK_CGEQRF( &m, &n, a, &lda, tau, work, &lwork, &info );
return info;
}
inline
void geqrf (int const m, int const n,
traits::complex_f* a, int const lda,
traits::complex_f* tau, traits::complex_f* work,
int const lwork, int& info)
{
LAPACK_CGEQRF (&m, &n,
traits::complex_ptr (a), &lda,
traits::complex_ptr (tau),
traits::complex_ptr (work), &lwork, &info );
}
