Вопросы с тегом «sparse-matrix»

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

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

РЕДАКТИРОВАТЬ: я проверяю, если какие-либо собственные значения имеют величину один или больше. Мне нужно найти наибольшее абсолютное собственное значение большой разреженной несимметричной матрицы. Я использовал eigen()функцию R , которая использует алгоритм QR из EISPACK или LAPACK, чтобы найти...

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

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

Мне было интересно, существует ли быстрый и эффективный способ заранее определить количество ненулевых элементов для операции умножения разреженных матриц при условии, что обе матрицы находятся в формате CSC или CSR. Я знаю, что есть один в пакете smmp, но мне нужно что-то, что уже реализовано в C...

Я хотел бы знать, есть ли быстрый способ вычисления евклидова расстояния двух векторов в октаве. Кажется, что для этого нет специальной функции, так что я должен просто использовать формулу с...

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

Я пытаюсь решить двумерное уравнение Пуассона с помощью конечных разностей. В процессе, я получаю разреженную матрицу только с переменными в каждом уравнении. Например, если переменные были U , то дискретизация даст:555UUU...

У меня есть два графика с почти n ~ 100000 узлов каждый. На обоих графиках каждый узел связан ровно с 3 другими узлами, поэтому матрица смежности является симметричной и очень разреженной. Сложность в том, что мне нужны все собственные значения матрицы смежности, но не собственные векторы. Чтобы...

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

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

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

Умножается ли матричное умножение (как Mat * Mat, так и Mat * Vec) на количество ненулевых элементов или на размер матрицы? Или какая-то комбинация двух. Как насчет формы. Например, у меня есть матрица 100 x 100 с 100 значениями в ней или матрица 1000 x 1000 с 100 значениями в ней. При возведении в...

Влияет ли компиляция PETSc с внешней библиотекой BLAS / LAPACK существенным образом на производительность разреженных матриц, или он использует эти библиотеки только для математики с плотной...

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

Я работаю над некой обусловленной большой разреженной линейной системой уравнений. Я хочу использовать двойную арифметику или квадро-двойную арифметику для их решения. Я знаю, что существует пакет под названием MPACK, разработанный Накатой, Махо, который может выполнять числовые линейные...

Например, библиотеки разреженных матриц C ++, которые я использовал - Eigen и SuiteSparse, похоже, что они не имеют какой-либо функциональности SVD для разреженных матриц. Так что просто любопытно, является ли SVD сложнее, чем QR / LU для разреженной...

Предположим, у меня есть оригинальная большая, разреженная линейная система: Ax0=b0Ax0=b0A\textbf{x}_0=\textbf{b}_0 . Теперь у меня нет A−1A−1A^{-1} как A слишком велика для разложения или любого вида разложения AAA , но предположим, что у меня есть решение x0x0\textbf{x}_0 найденное с помощью...

Учитывая общую разреженную матрицу A ∈рn × nA∈рN×NA \in \mathbb{R}^{n\times n}с m << n (поправка:м ≪N2м«N2m \ll n^2ненулевые элементы (обычно m ∈ O ( n )м∈О(N)m \in {\cal O}(n)). AAA является общим в том смысле, что не имеет специфических свойств (например, положительной определенности), и не...