Решите систему линейных уравнений — метод квази-минимальных невязок
попытки решить систему линейных уравнений x
= qmr(A
,b
)A*x = b
для x
использование Метода квази-минимальных невязок. Когда попытка успешна, qmr
отображает сообщение, чтобы подтвердить сходимость. Если qmr
сбои, чтобы сходиться после максимального количества итераций или остановов по любой причине, это отображает диагностическое сообщение, которое включает относительный остаточный norm(b-A*x)/norm(b)
и номер итерации, в который остановленный метод.
Сходимость большинства итерационных методов зависит от числа обусловленности матрицы коэффициентов, cond(A)
. Можно использовать equilibrate
улучшить число обусловленности A
, и самостоятельно это облегчает для большинства итеративных решателей сходиться. Однако использование equilibrate
также приводит к лучшим качественным матрицам перед формирователем, когда вы впоследствии учитываете уравновешенный матричный B = R*P*A*C
.
Можно использовать матричные функции переупорядочения, такие как dissect
и symrcm
переставить строки и столбцы матрицы коэффициентов и минимизировать количество ненулей, когда матрица коэффициентов учтена, чтобы сгенерировать предварительный формирователь. Это может уменьшать память и время, требуемое впоследствии решить предобусловленную линейную систему.
[1] Барретт, R., М. Берри, Т. Ф. Чан, и др., Шаблоны для Решения Линейных систем: Базовые блоки для Итерационных методов, SIAM, Филадельфия, 1994.
[2] Freund, Роланд В. и Неель М. Нэчтигэл, “QMR: метод квази-минимальных невязок для неэрмитовых линейных систем”, SIAM Journal: Numer. Математика. 60, 1991, стр 315–339.