Оптимизация на одном потоке ЦП

#include <iostream>
#include <chrono>
 
void gaussElimination(double** A, double* b, int n) {
    for (int i = 0; i < n; ++i) {
        for (int j = i + 1; j < n; ++j) {
            double factor = A[j][i] / A[i][i];
            b[j] -= factor * b[i];
            for (int k = i; k < n; ++k) {
                A[j][k] -= factor * A[i][k];
            }
        }
    }
 
    double* x = new double[n];
    for (int i = n - 1; i >= 0; --i) {
        double sum = 0.0;
        for (int j = i + 1; j < n; ++j) {
            sum += A[i][j] * x[j];
        }
        x[i] = (b[i] - sum) / A[i][i];
    }
 
    std::cout << "Решение x:\n";
    for (int i = 0; i < n; ++i) {
        std::cout << x[i] << " ";
    }
    std::cout << std::endl;
 
    delete[] x;
}
 
int main() {
    setlocale(LC_ALL, "rus");
    int n = 5000;
 
    double** A = new double* [n];
    double* b = new double[n];
 
    for (int i = 0; i < n; ++i) {
        A[i] = new double[n];
        for (int j = 0; j < n; ++j) {
            A[i][j] = 1.0 + static_cast<double>(rand()) / (RAND_MAX / 9.0);
        }
        b[i] = 1.0 + static_cast<double>(rand()) / (RAND_MAX / 9.0);
    }
 
    auto start = std::chrono::high_resolution_clock::now();
 
    gaussElimination(A, b, n);
 
    auto end = std::chrono::high_resolution_clock::now();
    std::chrono::duration<double> duration = end - start;
    std::cout << "Время выполнения: " << duration.count() << " секунд." << std::endl;
 
    for (int i = 0; i < n; ++i) {
        delete[] A[i];
    }
    delete[] A;
    delete[] b;
 
    return 0;
}