dvector_functions.h

Go to the documentation of this file.

#ifndef DVECTOR_FUNCTIONS
#define DVECTOR_FUNCTIONS

#include "lapack_functions.h"


//y=x
void DoubleVectorCopy(const tLVectorIndex& nx,
                      const tLVectorIncrement& incx,
                      const double *x,
                      const tLVectorIncrement& incy,
                      double* y);

//y+=alpha.x
void DoubleVectorAdd(const tLVectorIndex& nX,
                     const tLVectorIncrement& incX,const double *X,
                     const double& alpha,
                     const tLVectorIncrement& incY,double* y);

//x*=alpha
void DoubleVectorScale(const tLVectorIndex& nX,
                       const tLVectorIncrement& incX,double *X,const double& alpha);


//<y,x>
double DoubleVectorDot(const tLVectorIndex& nX,
                       const tLVectorIncrement& incX,const double *x,
                       const tLVectorIncrement& incY,const double *y);

void DoubleVectorRotation(const tLVectorIndex& nX,
                          const tLVectorIncrement& incX,double *x,
                          const double& c,const double& s,
                          const tLVectorIncrement& incY,double* y);

void DoubleVectorRotationPlane(double& dx,double& dy,double& c,double& s );
//Y<->X
void DoubleVectorSwap(const tLVectorIndex& nX,
                      const tLVectorIncrement& incX,double* X,
                      const tLVectorIncrement& incY,double *Y);

tLVectorIndex DoubleVectorIndexMax(const tLVectorIndex& nX,
                                  const tLVectorIncrement& incX,const double *x);

double DoubleVectorNorm1(const tLVectorIndex& nX,
                         const tLVectorIncrement& incX,const double *x);


double DoubleVectorNorm2(const tLVectorIndex& nX,
                         const tLVectorIncrement& incX,const double *x);



#endif