Home Research Papers Posters Manuals Talks Download People Help Internal

libs/ck-libs/pmaf/matrix.h

Go to the documentation of this file.

class Matrix {
 protected:
  double **mtx;
  int order;
  
 public:

  Matrix(int dim) {
    mtx = (double **)malloc(dim*sizeof(double *));
    for (int i=0; i<dim; i++)
      mtx[i] = (double *)malloc(dim*sizeof(double));
    order = dim;
  }
  
  void setElement(int i, int j, double value) { 
    CmiAssert((i < order) && (j < order));
    mtx[i][j] = value; 
  }

  const double &elem(int i, int j) const { 
    CmiAssert((i < order) && (j < order));
    return mtx[i][j]; 
  }
  
  double determinant(void) const {
    double d = 0.0;
    if (order <= 1) return elem(0, 0);
    for (int i=0; i<order; i++) {
      double sign = (i%2) ? -1 : 1;
      d += sign * elem(i, 0) * cofactor(i, 0).determinant();
    }
    return d;
  }


  Matrix cofactor(int aI, int aJ) const {
    CmiAssert((aI < order) && (aJ < order));
    Matrix a(order-1);
    for (int i=0, k=0; i<order; i++)
      if (i != aI) {
        for (int j=0, l=0; j<order; j++) {
          if (j != aJ) {
            a.setElement(k, l, elem(i, j));
            l++;
          }
        }
        k++;
      }
    return a;
  }
};

---