MMSD_DoubleVector.h

Go to the documentation of this file.

#ifndef MMSD_DoubleVector_H
#define MMSD_DoubleVector_H

#include "LAP_DoubleVector.h"

#include "MMSD_ConstDoubleVector.h"

DEFINE_SPTR(MMSD_DoubleVector);
DEFINE_SVPTR(MMSD_DoubleVector);

class MMSD_DoubleVector : public virtual LAP_DoubleVector {
    SP_OBJECT(MMSD_DoubleVector);
    // ATTRIBUTES
    
private:
 
  
    // ASSOCIATIONS
  
  
public:
    // METHODS
  
    // CONSTRUCTORS 
  
    MMSD_DoubleVector(void);
    MMSD_DoubleVector(const int& n) ;
    MMSD_DoubleVector(const MMSD_DoubleVector& s); 
  
  
    // DESTRUCTORS 
  
  
    virtual ~MMSD_DoubleVector(void);
  
  
  
public:
    //----------------
    // NEW constructors
    // ---------------
    inline static SP::MMSD_DoubleVector New() {
        SP::MMSD_DoubleVector p(new MMSD_DoubleVector(),
                          MMSD_DoubleVector::Delete());
        p->setThis(p);
        return p;
    };
    inline static SP::MMSD_DoubleVector New(const int& nElements) {
        SP::MMSD_DoubleVector p(new MMSD_DoubleVector(nElements),
                          MMSD_DoubleVector::Delete());
        p->setThis(p);
        return p;
    };
    inline MMSD_DoubleVector& operator=(const lapack_real& s)  {
        LAP_DoubleVector::init(s);
        return (*this);
    };
    inline MMSD_DoubleVector& operator=(const MMSD_DoubleVector& s)  {
        copy(s);
        return (*this);
    };
    inline MMSD_DoubleVector& operator=(const MMSD_ConstDoubleVector& s)  {
        copy(s);
        return (*this);
    };


    static inline double triProduct(const tArrayIndex& n,
                                    const double* a,
                                    const double* b,
                                    const double* c) {

      
        
        const double *va=a;
        const double *vb=b;
        const double *vc=c;
        
        double s=0;
        for (tArrayIndex i=0;i<n;i++) {
            s+=(*va)*(*vb)*(*vc);
            va++;
            vb++;
            vc++;
        }
        return s;
        
        
        
    }
    static inline double quadProduct(const tArrayIndex& n,
                                     const double* a,
                                     const double* b,
                                     const double* c,
                                     const double* d) {
       
      
        const double *va=a,*vb=b,*vc=c,*vd=d;
        double s=0;
        for (tArrayIndex i=0;i<n;i++) {
            s+=(*va)*(*vb)*(*vc)*(*vd);
            va++;
            vb++;
            vc++;
            vd++;
        }
        return s;
        
        
        
    }
    
    inline double triProduct(const double * b,
                             const double * c) const {
        return triProduct(getSize(),getValues(),b,c);
        
    }
    inline double quadProduct(const double* a,
                              const double* b,
                              const double* c) const {

        return quadProduct(getSize(),getValues(),a,b,c);
        
    }
                             
   
    
};

#endif