void complexVector::axpy(const Real alpha, const complexVector& x)
{
assert( x.length() == length() );
int N=2*length();
BLASRPFX(axpy,BLASINT(N),BLASREAL(alpha),BLASPREAL(x.data()),BLASINT(mnlIntOne),BLASPREAL(v_data),BLASINT(mnlIntOne));
}
void complexVector::operator -=(const complexVector& v2)
{
assert( v2.length()==length() );
BLASCPFX(axpy,BLASINT(m_N),BLASCOMPLEX(mnlComplexMinusOne),BLASPCOMPLEX(v2.data()),BLASINT(mnlIntOne),BLASPCOMPLEX(v_data),BLASINT(mnlIntOne));
}
void complexVector::axpy(const std::complex<Real> alpha, const complexVector& x)
{
assert( x.length() == length() );
BLASCPFX(axpy,BLASINT(m_N),BLASCOMPLEX(alpha),BLASPCOMPLEX(x.data()),BLASINT(mnlIntOne),BLASPCOMPLEX(v_data),BLASINT(mnlIntOne));
}
complexVector::complexVector(const string& n_name, const complexVector& vector) :
basicVector< std::complex<Real> >(n_name,vector.length(),false,const_cast<std::complex<Real>*>(vector.data()))
{
}