MATH_Vector.h

Go to the documentation of this file.

#ifndef MATH_Vector_H
#define MATH_Vector_H

#include "CORE_Object.h"



DEFINE_SPTR(MATH_Vector);

class MATH_Vector : public CORE_Object {
    SP_OBJECT(MATH_Vector);
    
    
protected:
    MATH_Vector(void) {
        
    }
    
    
    
    // DESTRUCTORS 
    
  
    virtual ~MATH_Vector(void) {
    }

public:
    virtual SP::MATH_Vector NewInstance() const=0;

    //operators
    //========

   
     
    virtual tReal& operator[](const tUIndex& i)=0;
    
    virtual const tReal& operator[](const tUIndex& i) const=0;
    
    virtual  const tReal& operator()(const tUIndex& i) const {
        return (*this)[i];
       
    }
    virtual  tReal& operator()(const tUIndex& i)  {
        return (*this)[i];
       
    }
    //set methods
    //===========
    /* \brief copy method
     * @param c: vector to copy
     */
    virtual void copy(const MATH_Vector& v) {
    }


    virtual void setCapacity(const tUIndex& n)=0;
    virtual void setSize(const tUIndex& n)=0;
    
    virtual void init(const tReal& alpha);
    
    virtual void init(const tReal& alpha,const MATH_Vector& y);

    inline void setValue(const tUIndex& i,const tReal& v) {
        getValues()[i]=v;
    }
    
    virtual void setValuesByReference(const tUIndex& n,tReal *Vs )=0;
    
    //get methods
    //============
    virtual tULLInt getMemorySize() const {
        return 0;
    }
    virtual const tReal* getValues() const=0;

    virtual tReal* getValues()=0;

    virtual tUIndex getSize() const=0;

    virtual tUIndex getCapacity() const=0;

    
    
    
    //arithmetic operations
    //========================
   

    virtual void add(const MATH_Vector& y);
    
    virtual void divide(const MATH_Vector& x,MATH_Vector& y) const; 
    inline void divide(const MATH_Vector& x) {
        divide(x,*this);
    }
    
    virtual void add(const tReal& alpha,const MATH_Vector& y);
    
    
    virtual void add(const tReal& beta,const MATH_Vector& x,const tReal& alpha);
    
    virtual void add(const tReal& beta,const MATH_Vector& x,const tReal& alpha,MATH_Vector& R);
    
    virtual void add(const tReal& beta,const MATH_Vector& x,MATH_Vector& R) {
        add(beta,x,1,R);
    }
    
    virtual tReal  dot(const MATH_Vector& x) const;

    virtual void dot(const tReal& alpha);

    inline tReal norm2() const {
        return dot(*this);
    }
    inline tReal norm() const {
        return sqrt(dot(*this));
    }
    
    
  

    
};
#endif