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

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

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

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

49 linear-algebra  linear-solver  sparse-matrix 

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

27 linear-algebra  algorithms  eigensystem  sparse-matrix 

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

22 linear-algebra  linear-solver  sparse-matrix  krylov-method 

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

17 linear-algebra  sparse-matrix  oop 

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

17 matrix  sparse-matrix 

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

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

Я пытаюсь решить двумерное уравнение Пуассона с помощью конечных разностей. В процессе, я получаю разреженную матрицу только с переменными в каждом уравнении. Например, если переменные были U , то дискретизация даст:555UUU Ui−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jUi−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,JU_{i-1,j} + U_{i+1,j} -4U_{i,j} + U_{i,j-1} + U_{i,j+1} = f_{i,j} Я знаю, что могу решить эту систему итеративным методом, …

15 linear-algebra  pde  finite-difference  sparse-matrix  banded-matrix 

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

15 python  sparse-matrix  numpy  eigenvalues 

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

13 linear-algebra  sparse-matrix  eigenvalues  matrix 

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

13 linear-solver  sparse-matrix  complexity  multigrid 

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

12 linear-algebra  petsc  linear-solver  sparse-matrix 

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

11 python  linear-solver  sparse-matrix 

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

10 linear-algebra  sparse-matrix  matrix  svd  eigen 

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

10 sparse-matrix  matrix  condition-number