main.cpp

Go to the documentation of this file.

 
#include <iostream>
#include <string>
#include <vector>
#include <cmath> // Included for fabs function
#include <filesystem>
#include <limits> // For std::numeric_limits
#include "utils.h"
#include "methods.h"
#include "interaction.h"
 
using namespace std;
 
int main()
{
    char choice;
    do
    {
        string selected_file = SelectInputFile();
        if (selected_file.empty())
        {
            return 1; // File selection failed
        }
 
        // Start timer after selecting the file
        auto start_time = StartTimer();
 
        vector<vector<double>> matrix;
        int rows, cols;
        if (!InitMatrix(matrix, selected_file, rows, cols))
        {
            return 1; // Matrix initialization failed
        }
 
        ShowEquations(matrix, rows, cols);
        cout << "Starting Gaussian elimination process..." << endl;
        int exchange_count = GaussianElimination(matrix, rows, cols);
        cout << "Gaussian elimination completed." << endl
             << endl;
 
        int rank = DetermineRank(matrix, rows, cols);
        bool consistent = CheckConsistency(matrix, rows, cols);
 
        if (!consistent)
        {
            cout << "The system of equations is inconsistent and has no solution." << endl;
        }
        else if (rank < (cols - 1))
        {
            ShowGeneralSolution(matrix, rows, cols, rank);
        }
        else
        {
            vector<double> solution;
            bool solvable = BackSubstitution(matrix, rows, cols, solution);
            if (solvable)
            {
                DisplaySolution(solution);
            }
            else
            {
                cout << "The system of equations is inconsistent and has no solution." << endl;
            }
        }
 
        // Stop timer after the solution is displayed
        StopTimer(start_time);
        choice = AskRunAgain();
 
    } while (choice == 'y' || choice == 'Y');
 
    WaitForExit();
    return 0;
}