ocmatrix.h

Go to the documentation of this file.

#ifndef ocmatrix_h
#define ocmatrix_h

#ifndef MATRIXH
#define MAT void
#define SPMAT void
#define PERM void
#endif

struct Object;
class IvocVect;
class OcFullMatrix;
#define Vect IvocVect
#define Matrix OcMatrix

class OcMatrix {
public:
   enum {MFULL=1, MSPARSE, MBAND};
   static OcMatrix* instance(int nrow, int ncol, int type = MFULL);
   virtual ~OcMatrix();

   virtual double* mep(int i, int j) {unimp(); return nil;} //matrix element pointer
   virtual double getval(int i, int j) { unimp(); return 0.; }
   virtual int nrow() {unimp(); return 0;}
   virtual int ncol() {unimp(); return 0;}
   virtual void resize(int, int) {unimp();}

   OcFullMatrix* full();
   
   virtual void mulv(Vect* in, Vect* out){unimp();}
   virtual void mulm(Matrix* in, Matrix* out){unimp();}
   virtual void muls(double, Matrix* out){unimp();}
   virtual void add(Matrix*, Matrix* out){unimp();}
   virtual void getrow(int, Vect* out){unimp();}
   virtual void getcol(int, Vect* out){unimp();}
   virtual void getdiag(int, Vect* out){unimp();}
   virtual void setrow(int, Vect* in){unimp();}
   virtual void setcol(int, Vect* in){unimp();}
   virtual void setdiag(int, Vect* in){unimp();}
   virtual void setrow(int, double in){unimp();}
   virtual void setcol(int, double in){unimp();}
   virtual void setdiag(int, double in){unimp();}
   virtual void zero(){unimp();}
   virtual void ident(){unimp();}
   virtual void exp(Matrix* out){unimp();}
   virtual void pow(int, Matrix* out){unimp();}
   virtual void inverse(Matrix* out){unimp();}
   virtual void solv(Vect* vin, Vect* vout, boolean use_lu) {unimp();}
   virtual void copy(Matrix* out){unimp();}
   virtual void bcopy(Matrix* mout, int i0, int j0, int n0, int m0, int i1, int j1){unimp();}
   virtual void transpose(Matrix* out){unimp();}
   virtual void symmeigen(Matrix* mout, Vect* vout){unimp();}
   virtual void svd1(Matrix* u, Matrix* v, Vect* d){unimp();}
   virtual double det(int* e){unimp(); return 0.0;}
   virtual int sprowlen(int){unimp(); return 0;}
   virtual double spgetrowval(int i, int jindx, int* j){unimp(); return 0.;}

   void unimp();
   Object** temp_objvar();
protected:
   OcMatrix(int type);
public:
   Object* obj_;
private:
   int type_;
};

extern "C" {
   extern Matrix* matrix_arg(int);
}

class OcFullMatrix : public OcMatrix { // type 1
public:
   OcFullMatrix(int, int);
   virtual ~OcFullMatrix();

   virtual double* mep(int, int);
   virtual double getval(int i, int j);
   virtual int nrow();
   virtual int ncol();
   virtual void resize(int, int);

   virtual void mulv(Vect* in, Vect* out);
   virtual void mulm(Matrix* in, Matrix* out);
   virtual void muls(double, Matrix* out);
   virtual void add(Matrix*, Matrix* out);
   virtual void getrow(int, Vect* out);
   virtual void getcol(int, Vect* out);
   virtual void getdiag(int, Vect* out);
   virtual void setrow(int, Vect* in);
   virtual void setcol(int, Vect* in);
   virtual void setdiag(int, Vect* in);
   virtual void setrow(int, double in);
   virtual void setcol(int, double in);
   virtual void setdiag(int, double in);
   virtual void zero();
   virtual void ident();
   virtual void exp(Matrix* out);
   virtual void pow(int, Matrix* out);
   virtual void inverse(Matrix* out);
   virtual void solv(Vect* vin, Vect* vout, boolean use_lu);
   virtual void copy(Matrix* out);
   virtual void bcopy(Matrix* mout, int i0, int j0, int n0, int m0, int i1, int j1);
   virtual void transpose(Matrix* out);
   virtual void symmeigen(Matrix* mout, Vect* vout);
   virtual void svd1(Matrix* u, Matrix* v, Vect* d);
   virtual double det(int* exponent);
private:
   MAT* m_;
   MAT* lu_factor_;
   PERM* lu_pivot_;
};

class OcSparseMatrix : public OcMatrix {  // type 2
public:
   OcSparseMatrix(int, int);
   virtual ~OcSparseMatrix();

   virtual double* mep(int, int);
   virtual double* pelm(int, int); // nil if element does not exist
   virtual int nrow();
   virtual int ncol();
   virtual double getval(int, int);
   virtual void mulv(Vect* in, Vect* out);
   virtual void solv(Vect* vin, Vect* vout, boolean use_lu);

   virtual void setrow(int, Vect* in);
   virtual void setcol(int, Vect* in);
   virtual void setdiag(int, Vect* in);
   virtual void setrow(int, double in);
   virtual void setcol(int, double in);
   virtual void setdiag(int, double in);

   virtual int sprowlen(int); // how many elements in row
   virtual double spgetrowval(int i, int jindx, int* j); 
private:
   SPMAT* m_;
   SPMAT* lu_factor_;
   PERM* lu_pivot_;
};

#ifndef MATRIXH
#undef MAT
#undef SPMAT
#endif

#endif