Доброго времени суток. Такой вопрос: у меня есть разреженная матрица c заполненными 7 диагоналями (не ленточная матрица). Пытаясь воспроизвести что-то похожее на метод Шварца выделяю 3 (например) матрицы с пересечением и решаю параллельно методом GMRES. На каждой итерации значения, не входящие в матрицы и расположенные в тех же строках, умножаются на текущий вектор результата и вычитаются из вектора правых частей. Так вот решается медленнее, чем в одном потоке - может у вас есть комментарии/рекомендации?
p.s. Хотел загрузить рисунок, чтобы было понятнее, но на этом форуме его вроде не прикрепить?
|