MatrixMeschach.h

Go to the documentation of this file.

//   COCOA class header file
//Id:  MatrixMeschach.h
//CAT: Model
//
//   Class for matrices
//
//   History: v1.0
//   Pedro Arce
 
#ifndef _ALIMATRIX_HH
#define _ALIMATRIX_HH
 
#include "Alignment/CocoaUtilities/interface/CocoaGlobals.h"
#include <vector>
#include <iostream>
 
extern "C" {
#include <matrix.h>
#include <matrix2.h>
}
 
// Meschach external (matrix.h) defines C++-incompatible macros
// which break other code (e.g. standard <limits>).
// Since these are not used here, undef them.
#undef max
#undef min
#undef catch
#undef Real
 
class MatrixMeschach {
public:
  MatrixMeschach();
  MatrixMeschach(ALIint NoCol, ALIint NoLin);
  MatrixMeschach(const MatrixMeschach& mat);
  ~MatrixMeschach();
 
  void AddData(ALIuint col, ALIuint lin, ALIdouble data);
  void transpose();
  void inverse();
  void Dump(const ALIstring& mtext);
  void ostrDump(std::ostream& fout, const ALIstring& mtext);
  void EliminateLines(ALIint lin_first, ALIint lin_last);
  void EliminateColumns(ALIint lin_first, ALIint lin_last);
  void SetCorrelation(ALIint i1, ALIint i2, ALIdouble corr);
 
  MatrixMeschach& operator=(const MatrixMeschach& mat);
  void operator*=(const MatrixMeschach& mat);
  void operator+=(const MatrixMeschach& mat);
  void operator*=(const ALIdouble num);
  ALIdouble operator()(int i, int j) const;
 
  //ACCESS PRIVATE DATA MEMBERS
  ALIint NoLines() const { return _NoLines; }
  ALIint NoColumns() const { return _NoColumns; }
  void setNoColumns(ALIint ncol) { _NoColumns = ncol; }
  void setNoLines(ALIint nlin) { _NoLines = nlin; }
  const MAT* Mat() const { return _Mat; }
  void setMat(MAT* mat) { _Mat = mat; }
  MAT* MatNonConst() const { return _Mat; }
 
private:
  // private data members
  ALIint _NoLines;
  ALIint _NoColumns;
  //  vector< ALIdouble> _data;
  MAT* _Mat;
 
  void copy(const MatrixMeschach& mat);
};
 
MatrixMeschach operator*(const MatrixMeschach& mat1, const MatrixMeschach& mat2);
MatrixMeschach operator+(const MatrixMeschach& mat1, const MatrixMeschach& mat2);
MatrixMeschach operator-(const MatrixMeschach& mat1, const MatrixMeschach& mat2);
MatrixMeschach operator*(const ALIdouble doub, const MatrixMeschach& mat);
MatrixMeschach operator*(const MatrixMeschach& mat, const ALIdouble doub);
 
MatrixMeschach* MatrixByMatrix(const MatrixMeschach& mat1, const MatrixMeschach& mat2);
typedef MatrixMeschach ALIMatrix;
 
#endif