MatematicaMatrix.h

#ifndef MATEMATICAMATRIX_H
#define MATEMATICAMATRIX_H
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "MatematicaStruct.h"
//---------------------Matrice--------------------
class Matrice{
 private:
  void allocate(){
    delete[] data;
    data = new double [NSize*NSize];
  };
  int NSize;
  int NCol;
  int NRow;
  int NZed;
  int NDim;
 public:
  Matrice(int newNSize);// {Matrice(newNSize,newNSize);};
  //Matrice(int newNSize, int newactualsize);
  Matrice(int NRow, int NCol);
  Matrice(int NRow, int NCol,int NZed);
  Matrice(int ExtNRow, int ExtNCol, double *pointer);
  Matrice(SPLINE Wg);
  Matrice(SPLINE Wg,int Dim);
  Matrice(Quadri q,int dim);
  Matrice(double *M,int Nr,int Nc);
  Matrice(double Roll,double Pitch,double Yaw,int NDim);
  Matrice(double *Axis,double Angle,int NDim);
  ~Matrice();
  bool Set(int row, int column,double Val);
  bool Add(int row,int col,double Val);
  int Size(){return NSize;};
  int pNCol(){return NCol;};
  int pNRow(){return NRow;};
  int pNZed(){return NZed;};
  void Apply(double *Known,double *UnKnown);
  int Solve(double *Known,double *UnKnown);
  int getNRow();
  void Shout(const char *s, ... );
  void comparetoidentity();
  void settoproduct(Matrice& left, Matrice& right);
  void copymatrix(Matrice&  source);
  void setNRow(int newNRow);
  void getvalue(int row, int column, double& returnvalue, bool& success);
  void Invert();
  void Mult(Matrice &A,Matrice &B);
  void Mult(Matrice &A);
  Vettore Mult(Matrice &A,Vettore &v);
  Vettore Mult(Vettore &v);
  void Mult(Vettore &v,Vettore &u);
  void RandomFill(double Max);
  void FillDiffOperator(SPLINE Wg,int NDim);
  void FillCanny();
  void FillGaussian(double Sigma,double CutOff);
  void FillGaussian5();
  void Transpose();
  void Normalize();
  void Clear();
  void Multiply(double Val);
  void CopyOn(Matrice *B);
  void Print();
  double Val(int row);
  double Val(int row,int col);
  double Val(int row,int col,int zed);
  double Det();
  Matrice operator+(Matrice&);
  //Matrice operator*(Matrice&) ;
  Matrice operator*(Matrice &A);
  //  Matrice &operator*(const Matrice&) const;
  Matrice operator*(const double&) const;
  Vettore operator*(const Vettore&) const;
  Matrice operator=(Matrice&);
  //  Matrice operator[][](int row,int col);
  Matrice operator^(Matrice&) const;
  double* data;
  void ConvoluteMatrix(double *Plot,int NGrid,int NDim,int IfMinImConv);
  void ConvoluteMatrix1(double *Plot,int NGrid);
  void ConvoluteMatrix1MinImConv(double *Plot,int NGrid);
  void ConvoluteMatrix2(double *Plot,int NGrid);
  void ConvoluteMatrix2MinImConv(double *Plot,int NGrid);
  void ConvoluteMatrix3(double *Plot,int NGrid);
};
#endif //MATEMATICA_H