Решение линейной системы с собственным с помощью (почти) диагональной матрицы ⇐ C++

[Anonymous]

мне нужно решить систему типа < /p>
Ax = b< /code> < /p>
, который обычно делается в собственном собственном виде x = a.housholderqr (). Решение (b) < /code>.
Теперь в моем случае матрица A в основном является блочной диагональной матрицей с несколькими вне диагональных элементов.
зависит от системы, которые мне нужны для решения, его Demages Ranges Ranges Ranges. Чтобы дать идею, здесь есть меньшая версия (которая сохраняет шаблон)
< /p>
Числа от основного диагонала. Размеры Там будет просто больше нулей между последовательности чисел не ну-нулей. Модуль?) < /p>
Спасибо за любую помощь по этому вопросу. Они содержит фактические значения задачи, которую я решаю. < /P>
Я добавляю здесь какой -то контекст из комментариев, которые я получил @peter, @chtz, @abcdefg и @reinderien. a заполняется той же положительной определенной матрицей (перспективной), которая повторяется в диагональном и никогда не изменяется во время симуляции < /li>
Я сейчас исследую формулу инверсии Шермана-Моррисон-Вудбери < /li>
< /ul>
Я также заметил следующие: если мы снимаем (концентрированные). вне диагональных компонентов, которые образуют симметричную матрицу, поможет ли это повысить эффективность алгоритма?>

Подробнее здесь: https://stackoverflow.com/questions/796 ... nal-matrix