/**
* C++ version of gsl_linalg_complex_householder_hv().
* @param tau A scalar
* @param v A vector
* @param w Another vector
* @return Error code on failure
*/
inline int complex_householder_hv( complex& tau, vector_complex const& v, vector_complex& w ){
return gsl_linalg_complex_householder_hv( tau.get(), v.get(), w.get() ); }
/**
* C++ version of gsl_linalg_complex_householder_hm().
* @param tau A scalar
* @param v A vector
* @param A A matrix
* @return Error code on failure
*/
inline int complex_householder_hm( complex& tau, vector_complex const& v, matrix_complex& A ){
return gsl_linalg_complex_householder_hm( tau.get(), v.get(), A.get() ); }
/**
* C++ version of gsl_blas_zgerc().
* @param alpha A constant
* @param X A vector
* @param Y A vector
* @param A A matrix
* @return Error code on failure
*/
int zgerc( complex const& alpha, vector_complex const& X, vector_complex const& Y,
matrix_complex& A ){
return gsl_blas_zgerc( alpha.get(), X.get(), Y.get(), A.get() ); }
/**
* C++ version of gsl_blas_zher2().
* @param Uplo Upper or lower triangular
* @param alpha A constant
* @param X A vector
* @param Y A vector
* @param A A matrix
* @return Error code on failure
*/
int zher2( CBLAS_UPLO_t Uplo, complex const alpha, vector_complex const& X,
vector_complex const& Y, matrix_complex& A ){
return gsl_blas_zher2( Uplo, alpha.get(), X.get(), Y.get(), A.get() ); }
/**
* C++ version of gsl_blas_zgemv().
* @param TransA Transpose type
* @param alpha A constant
* @param A A matrix
* @param X A vector
* @param beta Another constant
* @param Y A vector
* @return Error code on failure
*/
int zgemv( CBLAS_TRANSPOSE_t TransA, complex const& alpha, matrix_complex const& A,
vector_complex const& X, complex const& beta, vector_complex& Y ){
return gsl_blas_zgemv( TransA, alpha.get(), A.get(), X.get(), beta.get(), Y.get() ); }
/**
* C++ version of gsl_blas_zhemv().
* @param Uplo Upper or lower triangular
* @param alpha A constant
* @param A A matrix
* @param X A vector
* @param beta Another constant
* @param Y A vector
* @return Error code on failure
*/
int zhemv( CBLAS_UPLO_t Uplo, complex const& alpha, matrix_complex const& A,
vector_complex const& X, complex const& beta, vector_complex& Y ){
return gsl_blas_zhemv( Uplo, alpha.get(), A.get(), X.get(), beta.get(), Y.get() ); }
/**
* C++ version of gsl_blas_zaxpy().
* @param alpha A constant
* @param X A vector
* @param Y A vector
* @return Error code on failure
*/
int zaxpy( complex const& alpha, vector_complex const& X, vector_complex& Y ){
return gsl_blas_zaxpy( alpha.get(), X.get(), Y.get() ); }
/**
* C++ version of gsl_blas_zscal().
* @param alpha A constant
* @param X A vector
*/
void zscal( complex const& alpha, vector_complex& X ){ gsl_blas_zscal( alpha.get(), X.get() ); }
/**
* C++ version of gsl_blas_zher2k().
* @param Uplo Upper or lower triangular
* @param Trans Transpose type
* @param alpha A constant
* @param A A matrix
* @param B Another matrix
* @param beta Another constant
* @param C Another matrix
* @return Error code on failure
*/
int zher2k( CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t Trans, complex const& alpha,
matrix_complex const& A, matrix_complex const& B, double beta,
matrix_complex& C ){
return gsl_blas_zher2k( Uplo, Trans, alpha.get(), A.get(), B.get(), beta, C.get() ); }
/**
* C++ version of gsl_blas_zhemm().
* @param Side Side to apply operation to
* @param Uplo Upper or lower triangular
* @param alpha A constant
* @param A A matrix
* @param B Another matrix
* @param beta Another constant
* @param C Another matrix
* @return Error code on failure
*/
int zhemm( CBLAS_SIDE_t Side, CBLAS_UPLO_t Uplo, complex const& alpha,
matrix_complex const& A, matrix_complex const& B, complex const& beta,
matrix_complex& C ){
return gsl_blas_zhemm( Side, Uplo, alpha.get(), A.get(), B.get(), beta.get(), C.get() ); }
/**
* C++ version of gsl_blas_ztrsm().
* @param Side Side to apply operation to
* @param Uplo Upper or lower triangular
* @param TransA Transpose type
* @param Diag Diagonal type
* @param alpha A constant
* @param A A matrix
* @param B Another matrix
* @return Error code on failure
*/
int ztrsm( CBLAS_SIDE_t Side, CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA,
CBLAS_DIAG_t Diag, complex const& alpha, matrix_complex const& A,
matrix_complex& B ){
return gsl_blas_ztrsm( Side, Uplo, TransA, Diag, alpha.get(), A.get(), B.get() ); }
/**
* C++ version of gsl_blas_zsyrk().
* @param Uplo Upper or lower triangular
* @param Trans Transpose type
* @param alpha A constant
* @param A A matrix
* @param beta Another constant
* @param C Another matrix
* @return Error code on failure
*/
int zsyrk( CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t Trans, complex const& alpha,
matrix_complex const& A, complex const& beta, matrix_complex& C ){
return gsl_blas_zsyrk( Uplo, Trans, alpha.get(), A.get(), beta.get(), C.get() ); }
/**
* C++ version of gsl_blas_zgemm().
* @param TransA Transpose type
* @param TransB Transpose type for B
* @param alpha A constant
* @param A A matrix
* @param B Another matrix
* @param beta Another constant
* @param C Another matrix
* @return Error code on failure
*/
int zgemm( CBLAS_TRANSPOSE_t TransA, CBLAS_TRANSPOSE_t TransB, complex const& alpha,
matrix_complex const& A, matrix_complex const& B, complex const& beta,
matrix_complex& C ){
return gsl_blas_zgemm( TransA, TransB, alpha.get(), A.get(), B.get(), beta.get(), C.get() ); }