#include "Solver.h"

#include "test_assert.h"
#include <cstdlib>
#include <future>
#include <iostream>
#include <vector>

static constexpr int EXECUTION_COUNT = 100;
static constexpr int KERNEL_CHECKS = 100;
static constexpr int MATRIX_MAX_SIZE = 100;

static const Solver solver;

static unsigned int GetRandomInt() {
	return rand() % MATRIX_MAX_SIZE + 1;
}

static Matrix GetRandomMatrix(std::size_t a_Raw, std::size_t a_Column) {
	Matrix matrix {a_Raw, a_Column};

	for (std::size_t i = 0; i < matrix.GetRawCount(); i++) {
		for (std::size_t j = 0; j < matrix.GetColumnCount(); j++) {
			matrix.at(i, j) = GetRandomInt();
		}
	}

	return matrix;
}

static void Test() {
	Matrix matrix = GetRandomMatrix(GetRandomInt(), GetRandomInt());

	for (std::size_t i = 0; i < matrix.GetRawCount(); i++) {
		for (std::size_t j = 0; j < matrix.GetColumnCount(); j++) {
			matrix.at(i, j) = GetRandomInt();
		}
	}

	Vect kernel = solver.Kernel(matrix);

	Matrix nullVector {matrix.GetRawCount(), 1};
	nullVector.Fill(0.0);

	for (std::size_t i = 0; i < kernel.GetCardinal(); i++) {
		test_assert(matrix * kernel.GetVector(i) == nullVector);
	}

	for (std::size_t i = 0; i < KERNEL_CHECKS; i++) {
		Matrix vector = GetRandomMatrix(kernel.GetDimension(), 1);

		test_assert(kernel.IsElementOf(vector) == (matrix * vector == nullVector));
	}

	Vect kernel2 = solver.Kernel(kernel.GetLinearSystem());

	test_assert(kernel == kernel2);
}

int main() {
	srand(time(0));

	std::vector<std::future<void>> results;

	// appelle la fonction Test() en parallèle
	for (int i = 0; i < EXECUTION_COUNT; i++) {
		auto handle = std::async(std::launch::async, &Test);
		results.push_back(std::move(handle));
	}

	return EXIT_SUCCESS;
}