//
// Overloaded function for dispatching to
// * netlib-compatible BLAS backend (the default), and
// * complex<double> value-type.
//
inline void axpy( const fortran_int_t n, const std::complex<double> a,
const std::complex<double>* x, const fortran_int_t incx,
std::complex<double>* y, const fortran_int_t incy ) {
BLAS_ZAXPY( &n, &a, x, &incx, y, &incy );
}
inline void axpy( integer_t const n, traits::complex_d const a,
traits::complex_d* x, integer_t const incx, traits::complex_d* y,
integer_t const incy ) {
BLAS_ZAXPY( &n, traits::complex_ptr(&a), traits::complex_ptr(x),
&incx, traits::complex_ptr(y), &incy );
}