Home Research Papers Posters Manuals Talks Download People Help Internal

libs/ck-libs/ifem/ilsi_vec.h

Go to the documentation of this file.

#ifndef __UIUC_CHARM_ILSI_VEC_H
#define __UIUC_CHARM_ILSI_VEC_H

typedef double real;

//Vector operations:
inline void copy(int n,real *c,const real *a)
{ for (int i=0;i<n;i++) c[i]=a[i]; }

inline void sub(int n,real *c,const real *a,const real *b)
{ for (int i=0;i<n;i++) c[i]=a[i]-b[i]; }

inline void fma(int n,real *c,const real *a,double k,const real *b)
{ 
    for (int i=0;i<n;i++) 
        c[i]=k*b[i]+a[i]; 
}

inline real dot(int n,const real *a,const real *b)
{ 
    real sum=(real)0;
    for (int i=0;i<n;i++) 
        sum+=a[i]*b[i]; 
    return sum;
}


class allocVector {
    real *sto;
public:
    allocVector(int n) {sto=new real[n];}
    ~allocVector() {delete[] sto;}
    
    operator real *() {return sto;}
    operator const real *() const {return sto;}
};


#endif

---