#include "TMatrixD.h"
#include "fstream"
#include "iostream"
#include "math.h"
#include "TVectorD.h"
#include "TArrayD.h"
#include "iomanip"
#include "TDecompLU.h"


using namespace std;


int main()
{
    TMatrixD matrixD(3,3);
    

    TArrayD data1(9);
    data1[0]=1;
    data1[1]=1;
    data1[2]=1;
    data1[3]=1;
    data1[4]=2;
    data1[5]=2;
    data1[6]=2;
    data1[7]=1;
    data1[8]=2;

    matrixD.SetMatrixArray(data1.GetArray());
    double DetD=matrixD.Determinant();
    
   
    TVectorD b1(0,2,10.,12.,19.,"END");
    TVectorD solution1(0,2);
    TDecompLU lu1(matrixD);
    bool ok;
    solution1=lu1.Solve(b1,ok);
    cout<<"bbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbb"<<endl;
    solution1.Print();
    cout<<"bbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbb"<<endl;


// ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo
// Fazekas equations
// ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo

   ifstream ifile1("Aij86_matrix_Fazekas.txt");
   ifstream ifile2("Bi86_vector_Fazekas.txt");
   ofstream ofileFazekas("Solution_Fazekas.txt");   
   TVectorD b2(0,8);
   TVectorD solution2(0,8);
   TMatrixD matrixFazekas(9,9);
   
   for(int i=0;i<9;i++)
   {
       for(int j=0;j<9;j++)
       {
           ifile1>>TMatrixDRow(matrixFazekas,i)(j);
       }
   }
   double temp1;
   for(int i=0;i<9;i++)
   {
       ifile2>>b2(i);
       ifile2>>temp1;
   }
   ofstream ofile1("LargeMatrix.txt");
   for(int i=0;i<9;i++)
   {
       for(int j=0;j<9;j++)
       {
           ofile1<<setw(15)<<scientific<<TMatrixDRow(matrixFazekas,i)(j);
       }
       ofile1<<setw(15)<<b2(i)<<endl;
   }
   TDecompLU lu2(matrixFazekas);
   solution2=lu2.Solve(b2,ok);

   cout<<"ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo";
   solution2.Print();
   cout<<"ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo";
   
   for(int i=0;i<9;i++)
   {
       ofileFazekas<<scientific<<solution2(i)<<endl;
   }
    


    return 100000000;

}