Я хочу знать , какие из классических линейных решателей (например , Гаусс-Зейделя, Jacobi, SOR) гарантированно сходятся для задачи , где положительно полу определена и, конечноA b ∈ i m ( A )A x = bAx=bAx=bAAAb ∈ i m ( A )b∈im(A)b \in im(A) (Примечание является полуопределенным и не...
В настоящее время я работаю над решением очень больших симметричных (но не положительно определенных) систем, порожденных некоторыми определенными алгоритмами. Эти матрицы имеют хороший разреженный блок, который можно использовать для параллельного решения. Но я не могу решить, должен ли я...
в Matlab и Linsolve, и Mldivide используются для решения системы линейных уравнений во всех детерминированных, переопределенных и недоопределенных случаях. Читая их документы, мне было интересно, какие различия между ними? Используют ли они почти одинаковые алгоритмы матричной факторизации и...
Я уже разработал рабочее решение метода конечных элементов для решения задач теплопередачи, используя GPU и OpenCL, используя метод сопряженных градиентов. Основным недостатком этого метода является высокий спрос на память. Более того, в случае видеокарт память часто очень ограничена. Я вижу два...
Мне нужно решить Ax = b, но я понимаю, что, даже если оно редкое, сохранение матричных коэффициентов моей задачи займет слишком много памяти. Поэтому сейчас я рассматриваю возможность использования метода без матрицы, потому что одни и те же коэффициенты появляются в матрице много раз, поэтому я...
Я намереваюсь решить Ax = b, где A - сложная, разреженная, несимметричная и крайне плохо обусловленная (номер условия ~ 1E + 20) квадратная или прямоугольная матрица. Я смог точно решить систему с ZGELSS в LAPACK. Но по мере того, как степени свободы в моей системе растут, требуется много времени...
Есть ли простое эмпирическое правило, чтобы сказать, стоит ли делать SOR вместо Gauss-Seidel? (и возможный способ, как оценить параметр перехвата )ωω\omega Я имею в виду, просто глядя на матрицу , или знание конкретной проблемы, которую представляет матрица? Я читал ответ на этот вопрос: есть ли...
Большая часть моего программирования - это одноразовые исследовательские коды на Си для моего собственного использования. Я никогда не распространял никакого кода, кроме близких соавторов. Я разработал алгоритм, который я публикую в научном журнале. Я хочу предоставить исходный код и, возможно,...
Предположим, что задана следующая линейная система Lx=c,(1)(1)Lx=c,Lx=c,\tag1 где представляет собой взвешенное лапласиан , как известно, положительно определенной с одномерным нуль - пространство , натянутое на , а перевод дисперсия , т. е. не меняет значение функции (производная которой )....
Я ищу быстрое (смею сказать, оптимальное?) Явное решение линейной вещественной задачи 3x3, A x = bAИксзнак равноб\mathbf{A}\mathbf{x} = \mathbf{b}, A ∈р3 × 3, b ∈р3A∈р3×3,б∈р3\mathbf{A} \in \mathbf{R}^{3 \times 3}, \mathbf{b} \in \mathbf{R}^{3}, матрица AA\mathbf{A} является общим, но близким к...
Мне очень интересно оптимизировать решение линейных систем для маленьких матриц (10x10), которые иногда называют крошечными матрицами. Есть ли готовое решение для этого? Матрицу можно считать неособой. Этот решатель должен выполняться более 1 000 000 раз в микросекундах на процессоре Intel. Я...