Решение Ax =B с помощью UMFPACK

Question

Решение Ax =B с помощью UMFPACK

Я использую UMFPACK для решения Ax = B, где A и B - матрицы. В настоящее время в моей реализации я решаю x(:,i) = A\b(:,i) один за другим, чтобы получить x. Есть ли способ напрямую получить x без необходимости циклически перебирать каждый столбец B, как это делает LAPACK?

0

c++ linear-algebra numerical-methods umfpack

Источник

user4451072 11 фев '16 в 17:52

1 ответ

Решение

Другие вопросы по тегам c++ linear-algebra numerical-methods umfpack

user1996490 21 фев '16 в 18:27 2016-02-21 18:27 · Accepted Answer · 2016-02-21 18:27

Боюсь, что нет. Руководство пользователя UMFPACK p59, о umfpack_*_ решить рутину (ы):

Эта процедура не решает несколько правых частей и не позволяет хранить b в виде разреженных столбцов.

Также в разделе "Будущая работа":

Процедуры решения могут обрабатывать несколько правых и редких правых сторон.

Что скорее подразумевает, что это не поддерживается в настоящее время.