#include "TMatrixDSym.h"

#define POINTER
#define SETMATARR

int MatrixFillTest(){
  const Int_t N=3;
  #ifdef POINTER
  printf("Pointer version\n");
  TMatrixDSym* M = new TMatrixDSym(N);
  #else
  TMatrixDSym M(N);
  #endif

  Double_t *data = new Double_t[N * N];
  for(int i=0; i<N*N; i++) data[i]=0;
  
  //Fill the matrix
  for(int i=0; i<N; i++)
    for(int j=i; j<N; j++)
      {
	#ifdef SETMATARR
	data[i + N*j] = (i+1)*(j+3); data[j + N*i]=data[i + N*j]; 
        #else
	M[i][j] = (i+1)*(j+3); M[j][i]=M[i][j]; 
	#endif
      }
  
  #ifdef SETMATARR
  #ifdef POINTER
  M->SetMatrixArray(data);
  #else
  M.SetMatrixArray(data); 
  #endif
  #endif

  #ifdef POINTER
  M->Print();
  #else
  M.Print();
  #endif
  
  return 0;
}