CUDA_full / P05_MinimalResudualMethod / MinResidual

Постановка задачи

Метод минимальных невязок — это итеративный метод решения системы линейных уравнений Ax = b

Достаточным условием сходимости при любом начальном приближении x0 является положительная определенность матрицы A

Каждая итерация метода уточняет приближение x следующим образом:

rAx b

Ar, r

Ar, Ar

xx r

Задание

Реализовать с использованием CUBLAS метод минимальных невязок для решения системы линейных уравнений с положительно определенной матрицей

Указание: использовать функции cublasScopy, cublasSgemv, cublasSdot, cublasSaxpy.

Дополнительное задание: реализовать выход из метода по условию ||Ax – b|| < ε