«sparse-matrix» 태그된 질문

볼록하지 않은 전역 최적화 문제를 해결해야 할 몇 가지 문제가 있습니다. 현재 MATLAB의 Optimization Toolbox (특히 fmincon()algorithm = 사용 'sqp')를 사용하고 있습니다. 그러나 내 코드의 대부분은 Python이며 Python에서도 최적화를하고 싶습니다. 파이썬 바인딩과 경쟁 할 수있는 NLP 솔버가 fmincon()있습니까? 반드시 비선형 평등과 불평등 제약을 다룰 수있다 사용자가 Jacobian을 제공하도록 요구하지 않습니다. …

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 

스파 스 선형 시스템은 응용 분야에서 주파수가 증가함에 따라 나타납니다. 이러한 시스템을 해결하기 위해 선택할 수있는 많은 루틴이 있습니다. 최상위 수준에서는 직접 (예 : 희소 가우시안 제거 또는 Cholesky 분해, 특수 순서 알고리즘 및 다중 정면 방법) 방법과 반복적 (예 : GMRES, (bi) 공액 구배) 방법 사이에 유역이 있습니다. 직접 …

49 linear-algebra  linear-solver  sparse-matrix 

편집 : 고유 값의 크기가 1 이상인지 테스트하고 있습니다. 큰 희소 비대칭 행렬의 최대 절대 고유 값을 찾아야합니다. 나는 R의 eigen()함수를 사용하여 EISPACK 또는 LAPACK의 QR algo를 사용하여 모든 고유 값을 찾은 다음 abs()절대 값을 얻는 데 사용 합니다. 그러나 나는 더 빨리해야합니다. igraphR 패키지 에서 ARPACK 인터페이스를 사용해 보았습니다 …

27 linear-algebra  algorithms  eigensystem  sparse-matrix 

행렬 와 G가 있습니다. A 는 드문 드문하고 n이 매우 큰 n x n 이며 (수백만 정도일 수 있습니다.) G 는 m 이 작은 n x m 높이 행렬 ( 1 < m < 1000 )이며 각 열은 하나만 가질 수 있습니다. 1 명 , 나머지는 존재와 엔트리 0 '이되도록, …

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

객체 지향 프로그래밍을 최대한 활용하도록 가르치는 방법으로 희소 행렬 계산을위한 작은 라이브러리를 작성하고 있습니다. 나는 부품 (연결성 구조를 나타내는 희소 행렬과 그래프)이 매우 느슨하게 결합 된 멋진 객체 모델을 만들기 위해 열심히 노력했습니다. 내 생각에 코드는 훨씬 확장 가능하고 유지 관리가 가능합니다. 그러나 무딘 접근 방식을 사용하는 것보다 다소 느립니다. …

17 linear-algebra  sparse-matrix  oop 

두 행렬이 모두 CSC 또는 CSR 형식이라는 가정하에 희소 행렬 곱셈 연산에 대해 0이 아닌 숫자를 미리 찾는 빠르고 효율적인 방법이 있는지 궁금합니다. 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 

유한 차이로 2D 푸 아송 방정식을 풀려고합니다. 이 과정 에서 각 방정식에 변수 만있는 희소 행렬을 얻습니다 . 예를 들어 변수가 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의 빌트인 솔버와 비교하면서 아래 그림에 표시된 이상한 점을 발견했습니다. 여기에서 처음 100 개의 숫자로 계산 된 고유 …

15 python  sparse-matrix  numpy  eigenvalues 

나는 거의 n ~ 100,000 노드가있는 두 개의 그래프를 가지고 있습니다. 두 그래프에서 각 노드는 정확히 3 개의 다른 노드에 연결되므로 인접 행렬은 대칭이며 매우 드문 드문합니다. 어려운 부분은 인접 행렬의 모든 고유 값이 필요 하지만 고유 벡터 는 필요 하지 않다는 것입니다. 정확히 말하면, 이것은 내 생애에서 한 번 …

13 linear-algebra  sparse-matrix  eigenvalues  matrix 

Thomas 알고리즘이 알고리즘 복잡성 (LAPACK과 같은 구현 패키지를 찾지 않음) 측면에서 대칭 대각선으로 우세한 삼각형 시스템을 해결하는 가장 빠른 방법 (아마도?)인지 궁금합니다. Thomas 알고리즘과 멀티 그리드 둘 다 복잡도이지만 멀티 그리드의 상수 요소가 적다는 것을 알고 있습니까? 멀티 그리드가 더 빠를 수는 있지만 긍정적이지 않습니다.O ( n )영형(엔)O(n) 참고 : …

13 linear-solver  sparse-matrix  complexity  multigrid 

오른쪽 (300 ~ 1000)이 많은 동일한 희소 선형 시스템 (300x300 ~ 1000x1000)을 해결해야합니다. 이 첫 번째 문제 외에도 다른 시스템을 해결하고 싶지만 0이 아닌 동일한 요소 (단지 다른 값)로 일정한 희소성 패턴을 가진 많은 희소 시스템입니다. 내 행렬은 무기한입니다. 분해 및 초기화 성능은 중요하지 않지만 해결 단계의 성능은 중요합니다. 현재 …

12 linear-algebra  petsc  linear-solver  sparse-matrix 

파이썬 패키지 (petsc4py 등)와 지원하는 스파 스 다이렉트 솔버를 나열하십시오. 패키지 당 하나의 (커뮤니티-위키) 답변하십시오.

11 python  linear-solver  sparse-matrix 

예를 들어 내가 사용한 C ++ 희소 행렬 라이브러리-Eigen 및 SuiteSparse는 희소 행렬에 대한 SVD 기능이없는 것 같습니다. 궁금한 점은 SVD가 희소 행렬의 QR / LU보다 어렵다는 것입니다.

10 linear-algebra  sparse-matrix  matrix  svd  eigen 

나는 조건이 좋지 않은 큰 희소 선형 방정식 방정식을 연구하고 있습니다. 더블 더블 산술 또는 쿼드 더블 산술을 사용하여 문제를 해결하고 싶습니다. Nakata, Maho에서 개발 한 MPACK이라는 패키지가 쿼드-이중 산술 하에서 수치 선형 대수 계산을 수행 할 수 있다는 것을 알고 있습니다. 그러나 희소 행렬이 아닌 밀도가 높은 행렬을 위해 …

10 sparse-matrix  matrix  condition-number