Разработать программу для решения разреженных СЛАУ на С++

Мне нужно реализовать программу: Ускоренное решение разреженной СЛАУ методом Гаусса с использованием структуры данных биноминальной кучи.
(Ax=B) Сгенерировать матрицы A и B. Матрица А должна быть разреженной матрице размера n на n (в итоге программа должна производить вычисления для матриц n=150000) и записать эти матрицы в файл . Потом из файла уже считать эти матрицы и начать хранить хранить разреженную матрицу в Формате CRS и с помощью структуры данных биноминальная куча.
Дальше произвести вычисление СЛАУ методом Гаусса. Необходимо при выполнении всех этих процессах вычислений использовать технологию для распараллеливания программ Open MP. Во вложении пример работы, на основе которой нужно сделать программу.
И вот примерно какую таблицу хочется получить. Чтобы можно было показать результаты работы программы для матриц разного размера.