Вопросы с тегом «linear-solver»

У меня есть несколько сложных невыпуклых задач глобальной оптимизации. В настоящее время я использую MATLAB Optimization Toolbox (в частности, fmincon()с алгоритмом = 'sqp'), что довольно эффективно . Тем не менее, большая часть моего кода написана на Python, и я бы тоже хотел провести оптимизацию...

Разреженные линейные системы появляются с возрастающей частотой в приложениях. Для решения этих систем нужно выбирать из множества процедур. На самом высоком уровне существует водораздел между прямым (например, разреженным методом исключения Гаусса или разложением Холецкого, со специальными...

Насколько мне известно, есть 4 способа решения системы линейных уравнений (поправьте меня, если их больше): Если системная матрица является квадратной матрицей полного ранга, вы можете использовать правило Крамера; Вычислить обратную или псевдообратную матрицу системы; Используйте методы матричного...

У меня есть матрицы A и G . A является разреженным и имеет размер n × n с очень большим n (может быть порядка нескольких миллионов). G является матрицей высотой n × m с довольно небольшим m ( 1 < m < 1000 ), и в каждом столбце может быть только один 1 запись с остальным 0 «с, таким образом,...

Существует целый ряд различных библиотек, которые решают разреженную линейную систему уравнений, однако мне трудно понять, в чем различия. Насколько я могу судить, есть три основных пакета: Trilinos , PETSc и Intel MKL . Все они могут выполнять разреженные матричные решения, все они быстры...

Меня интересуют случаи, когда градиент сопряжения работает намного лучше, чем метод GMRES. Вообще, CG является предпочтительным выбором во многих случаях SPD (симметрично-положительно определенного), потому что он требует меньшего объема памяти и теоретическая оценка скорости сходимости для CG...

Я очень впечатлен последовательной производительностью многоуровневых прекондиционеров ILU с обратной связью , особенно для гетерогенного Гельмгольца , но я удивлен, что не смог найти никаких реализаций с открытым исходным кодом. В частности, ILUPACK делает двоичные файлы свободно доступными для...

Предположим, у нас есть линейная система, и мы ничего не знаем о ее обусловленности и не имеем предварительной информации о решении. Мы слепо применяем исключение Гаусса и получаем некоторое решение xxx . Можно ли определить, заслуживает ли доверия данное решение (т. Е. Хорошо ли обусловлена...

В таких методах, как gmres или bicgstab, может быть привлекательным использовать другой метод Крылова в качестве предварительного кондиционера. В конце концов, их легко реализовать без матрицы и в параллельной среде. Например, один кул может использовать несколько (скажем, ~ 5) итераций...

Числа с плавающей запятой одинарной точности занимают половину памяти, и на современных машинах (кажется, даже на графических процессорах) операции могут выполняться с ними почти вдвое быстрее, чем с двойной точностью. Многие коды FDTD, которые я обнаружил, используют исключительно арифметику...

Мне интересно, является ли алгоритм Томаса самым быстрым (доказуемо?) Решением симметричной диагонально доминирующей разреженной трехдиагональной системы с точки зрения алгоритмической сложности (не ища пакетов реализации, таких как LAPACK и т. Д.). Я знаю, что и алгоритм Томаса, и многосетка имеют...

Мультисетка (MG) может использоваться для решения линейной системы путем построения начального предположения x 0 и повторения следующего для i = 0 , 1 .. до сходимости:A x = bAИксзнак равнобAx=bИкс0Икс0x_0я = 0 , 1 ..язнак равно0,1 ..i=0,1.. Вычислить остаток ря= Б - хярязнак равноб-AИксяr_i =...

Я хочу решить нелинейную задачу с нелинейными ограничениями равенства, и я использую расширенный лагранжиан с членом регуляризации штрафа, который, как известно, портит число условий моих линеаризованных систем (на каждой итерации Ньютона, которую я имею в виду) , Чем больше срок штрафа, тем хуже...

Мне нужно решить ту же самую разреженную линейную систему (от 300x300 до 1000x1000) со многими правыми сторонами (от 300 до 1000). В дополнение к этой первой проблеме, я также хотел бы решить различные системы, но с одинаковыми ненулевыми элементами (только с разными значениями), то есть с...

Я использую MATLAB для решения проблемы, которая включает в себя решение на каждом временном шаге, где b изменяется со временем. Прямо сейчас я делаю это, используя MATLAB :A x = bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}бb\mathbf{b}mldivide x = A\b У меня есть возможность делать столько предварительных...

Я портирую существующий код из MATLAB на C ++ и имею линейную систему для решения (вместо более типичной формы )A x = bх = ЬxA=bxA=bA x = bAx=bAx=b Матрица плотная и общего вида, но не больше 1000x1000. Таким образом, в MATLAB решение находится по функции или через косую чертуAAAmrdivide(b,A)x =...

Учитывая систему где A ∈ R n × n , я прочитал, что, если итерация Якоби используется в качестве решателя, метод не будет сходиться, если b имеет ненулевую компоненту в нуль-пространстве A , Итак, как можно формально утверждать, что при условии, что b имеет ненулевой компонент, охватывающий нулевое...

Есть ли надежда на эффективное решение следующей линейной системы итерационным методом? A ∈ Rn × n, x ∈ RN, b ∈ RNс  n > 106A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n > 10^6 A x = bAx=bAx=b с участием A = ( Δ - K)A=(Δ−K)...

Пожалуйста, перечислите пакет Python (petsc4py и т. Д.) И разреженные прямые решатели, которые он поддерживает. Один (Community-Wiki) ответ на пакет,

Методы Ньютона-Крылова без Якобиана (JFNK) и методы Крылова в целом могут быть очень полезными, поскольку они не требуют явного хранения или построения матрицы, а только результатов произведений матрицы-вектора. Если вы действительно формируете разреженную систему, для вас есть много...