Матричное произведение

[Andrey_Kireew]

Возникла необходимость вычислить матричное произведение вида

$$p=r^T\cdot M^{-1}\cdot r$$

Для снижения вычислительной сложности, возникла идея использовать разложение Холецкого

$$M=L\cdot L^T$$

С его помощью это произведение сводится к вычислению нормы вектора

$$p=v^T\cdot v$$

где $v$ решение СЛАУ с нижнетреугольной матрицей

$$L\cdot v=r$$

Нет ли здесь ошибок?

[nnosipov]

Andrey_Kireew
А Вы сами не в состоянии это проверить? Взять одну формулу, подставить в другую и убедиться, что все сошлось. У меня это заняло секунд 20 от силы. Нельзя же так себе не доверять.

[Andrey_Kireew]

nnosipov в сообщении #1429028 писал(а):

Andrey_Kireew
А Вы сами не в состоянии это проверить?

Я то в состоянии, но никто не застрахован от ошибок. Обнаружить их потом бывает очень нелегко. Поэтому, при малейшем сомнении, лучше всё перепроверить. Если не хотите отвечать по существу - не отвечайте, я не настаиваю. А что касается недоверия - так это моё сугубо личное дело. Никто не обязан никому доверять, даже самому себе.

[arseniiv]

Я тоже проверил, нормально всё. И я не думаю, что nnosipov не дал добро. Если бы он нашёл ошибку, написал бы совсем другое, не «убедиться, что всё сошлось».

-- Пт дек 06, 2019 11:20:03 --

Не надо нервничать, это не полезно.