#include <iostream>
#include <fstream>

using namespace std;

void bezwyboru(unsigned rozmiarmacierzy, double**A, double*b, double*x)
{    double* p;

                p = new double[rozmiarmacierzy];

            for(int zmienna=0; zmienna < rozmiarmacierzy-1; zmienna++)
            {
                for(int i=zmienna; i<rozmiarmacierzy; i++)
                    p[i]=A[i][zmienna]/A[zmienna][zmienna];
                    
                for(int j=zmienna+1; j<rozmiarmacierzy; j++)
                {
                    for(int k=zmienna; k < rozmiarmacierzy; k++)
                    {
                        A[j][k]=A[j][k]-p[j]*A[zmienna][k];
                    }

                    b[j]=b[j]-p[j]*b[zmienna];
                    cout << b[j]<<endl;
                }
            }

}

int main()
{
    unsigned rozmiarmacierzy;
    double** A;
    double* b;
    double* x;

    while(true)
    {
        int wybor;
        cout << "\nMENU:" <<endl;
        cout <<"1 - wczytanie danych" << endl;
        cout <<"2 - metoda eliminacji Gaussa bez wyboru elementu"<< endl;
        cout <<"3 - metoda eliminacji Gaussa z wyborem elementu w wierszach"<<endl;
        cout <<"4 - metoda eliminacji Gaussa z wyborem elementu w kolumnach"<<endl;
        cout <<"5 - metoda eliminacji Gaussa z pelnym wyborem elementu"<< endl;
        cout <<"6 - wyjscie z programu"<<endl;

        cin>>wybor;

        switch(wybor)
        {
        case 1:
{
            ifstream plik_zrodlowy("test.csv");
            if (!plik_zrodlowy.is_open())
            {
                cout <<"Plik nie zosstal otwarty!"<<endl;
                return -1;
            }

            plik_zrodlowy >> rozmiarmacierzy;

            A = new double*[rozmiarmacierzy];

            A[0] = new double[rozmiarmacierzy*rozmiarmacierzy];

            for(unsigned i = 1; i< rozmiarmacierzy; i++)
                A[i] = A[i-1] + rozmiarmacierzy;

            b = new double[rozmiarmacierzy];

            char pol_kolumny;

            for (unsigned i = 0; i < rozmiarmacierzy+1; i++)
                plik_zrodlowy >> pol_kolumny;

            for (unsigned i = 0; i < rozmiarmacierzy; i++)
            {
                for (unsigned j = 0; j < rozmiarmacierzy; j++)
                {
                    plik_zrodlowy >> A[i][j];
                    plik_zrodlowy >> pol_kolumny;
                }
                plik_zrodlowy >> pol_kolumny;
                plik_zrodlowy >> b[i];
            }
            plik_zrodlowy.close();

            cout << "\nPobrano dane"<<endl;
            for(int i=0; i<rozmiarmacierzy; i++)
            {
                for(int j=0; j<rozmiarmacierzy; j++)
                {
                    cout <<"\t"<< A[i][j];
                }
                cout<<endl;
            }
            cout << rozmiarmacierzy<<endl;

            break;}

        case 2:

       bezwyboru(rozmiarmacierzy,A,b,x);
                break;
        case 3:
            {
            double suma=0.0;
            for(int k=rozmiarmacierzy-1; k>=0;k--)
            {
                suma=0.0;
                for(int j=0;j<rozmiarmacierzy;j++)
                {
                    if(k==j)
                        suma=suma+0;
                    else
                        suma=suma+(A[k][j]*b[j]);

                b[k]=(A[k][rozmiarmacierzy]-suma)/A[k][k];
                cout <<b[j]<<endl;
                }
            }
            }

            break;

        case 4:
            break;
        case 5:
            break;
        case 6:
            return 0;

    }
}
}