«linear-solver» 태그된 질문

비선형 평등 제약으로 비선형 문제를 해결하고 싶습니다. 그리고 알려진 Lagrangian을 페널티 정규화 용어와 함께 사용하고 있습니다. . 페널티 기간이 클수록 조건 수가 나빠집니다. 누군가가 특정 상황 에서이 나쁜 컨디셔닝을 제거하는 효율적인 방법을 알고 있습니까? 좀 더 구체적으로 말하면, 나는 일반적으로 중복 될 수있는 제약 조건이 많기 때문에 고전적인 증강 래그 …

12 linear-solver  preconditioning  constrained-optimization 

Jacobian-Free Newton-Krylov (JFNK) 방법과 일반적으로 Krylov 방법은 매트릭스의 명시적인 저장이나 구성이 필요하지 않고 행렬-벡터 제품의 결과 만 필요하므로 매우 유용합니다. 실제로 희소 시스템을 구성하면 많은 사전 조건자가 있습니다. 진정한 매트릭스 프리 분석법에는 무엇이 있습니까? 인터넷 검색은 "매트릭스 추정"에 대한 참조와 이것이 가능함을 나타내는 다른 것들을 나타냅니다. 이러한 방법은 일반적으로 어떻게 …

11 fluid-dynamics  linear-solver  preconditioning 

반복적 인 방법으로 다음 선형 시스템을 효율적으로 해결하려는 희망이 있습니까? ∈ Rn × n, x ∈ R엔, b ∈ R엔, 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 와 Δ − 6 6 …

11 linear-solver  dense-matrix  linear-system 

시스템 주어 여기서 경우, I는 경우 코비 반복 솔버로 사용되며, 그 판독 방법은 수렴하지 않을 비 - 제로를 갖는다 의 널 공간의 구성 요소 . 그렇다면 가 의 null 공간에 걸쳐 0이 아닌 구성 요소를 가지고 있다면 Jacobi 방법은 비 수렴 적이 라고 어떻게 공식적으로 말할 수 있습니까? null 공간에 …

11 linear-algebra  optimization  matrix  iterative-method  linear-solver 

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

11 python  linear-solver  sparse-matrix 

matlab에서 linsolve 와 mldivide 는 모두 결정, 과소 결정 및 불확실한 경우 모두 선형 방정식 시스템을 푸는 데 사용됩니다. 그들의 문서를 읽었을 때, 그들 사이에 어떤 차이점이 있는지 궁금했습니다. 그들은 세 가지 경우에 거의 동일한 행렬 분해 및 삼각 알고리즘 알고리즘을 사용하고 있습니까? A에 opts 속성이있는 경우 linsolve는 mldivide보다 빠릅니다. …

10 matlab  linear-solver 

현재 일부 특정 알고리즘에 의해 생성 된 매우 큰 대칭 (그러나 양의 한정이 아닌) 시스템을 해결하기 위해 노력하고 있습니다. 이 행렬에는 병렬 해결에 사용할 수있는 훌륭한 블록 희소성이 있습니다. 그러나 직접 접근 방식 (예 : 다중 정면) 또는 반복적 접근 방식 (사전 조건 GMRES 또는 MINRES)을 사용 해야하는지 결정할 수 …

10 parallel-computing  linear-solver  sparse-matrix  linear-system  gmres 

나는이 문제에 대한 수렴을 보장하는 고전적인 선형 솔버 (예를 들어, 가우스 - 자이 델, 코비, SOR)의 알고 싶은 긍정적 반 명확한 물론A x = bㅏ엑스=비Ax=bㅏㅏAb ∈ i m ( A )비∈나는미디엄(ㅏ)b \in im(A) (공지 는 반 정확한 것이 아니라 정한 것이 아니다)ㅏㅏA

10 linear-algebra  iterative-method  convergence  linear-solver 

Conjugate Gradient 방법을 사용하여 GPU 및 OpenCL을 사용하여 열 전달 문제를 해결하기 위해 유한 요소 방법의 작동 솔루션을 이미 개발했습니다. 이 방법의 주요 단점은 메모리에 대한 수요가 높다는 것입니다. 또한 그래픽 카드의 경우 메모리가 매우 제한적인 경우가 많습니다. 두 가지 옵션이 있습니다. 호스트 메모리로 서브 도메인 생성 및 메시 부분 …

10 parallel-computing  opencl  linear-solver 

작은 매트릭스 (10x10), 때로는 작은 매트릭스 라고도하는 선형 시스템 해결에서 지옥을 최적화하는 데 매우 관심이 있습니다. 이에 대한 준비된 해결책이 있습니까? 행렬은 단수로 가정 할 수 있습니다. 이 솔버는 Intel CPU에서 10 만 회 이상 마이크로 초 단위로 실행됩니다. 컴퓨터 게임에 사용되는 최적화 수준에 대해 이야기하고 있습니다. 어셈블리 및 아키텍처 …

9 linear-algebra  linear-solver  compiling  open-source  assembly 

3x3 선형 실제 문제인 빠른 (최적이라고 말합니까?) 명시 적 솔루션을 찾고 있습니다. A x = bAx=b\mathbf{A}\mathbf{x} = \mathbf{b}, ∈아르 자형3 × 3, b ∈아르 자형삼A∈R3×3,b∈R3\mathbf{A} \in \mathbf{R}^{3 \times 3}, \mathbf{b} \in \mathbf{R}^{3}. 매트릭스 ㅏA\mathbf{A} 일반적이지만 조건 번호가 1에 가까운 항등 행렬에 가깝습니다. 비b\mathbf{b} 실제로 약 5 자리의 정밀도로 센서를 측정하는 …

9 linear-solver  reference-request 

Gauss-Seidel 대신 SOR을 수행 할 가치 가 있는지에 대한 간단한 경험 법칙이 있습니까? (및 Realxation 매개 변수를 추정하는 가능한 방법ωω\omega) 나는 매트릭스 를 보거나 매트릭스가 나타내는 특정 문제에 대한 지식 만으로 의미 합니까? 나는이 질문들에 대한 답을 읽고 있었다 : 연속적인 이완법 (SOR) 방법을 최적화하기위한 휴리스틱이 있습니까? 하지만 너무 정교합니다. …

9 linear-solver  iterative-method  heuristics 

A가 복잡하고 희박하며 비대칭이며 조건이 좋지 않은 (조건 번호 ~ 1E + 20) 정사각형 또는 직사각형 행렬 인 Ax = b를 풀려고합니다. LAPACK의 ZGELSS로 시스템을 정확하게 해결할 수있었습니다. 그러나 시스템의 자유도가 커짐에 따라 희소성이 악용되지 않기 때문에 ZGELSS가 설치된 PC에서 시스템을 해결하는 데 시간이 오래 걸립니다. 최근에 동일한 시스템에 대해 …

9 linear-solver  sparse-matrix  lapack  condition-number 

Ax = b를 풀어야하지만 희소하더라도 문제의 행렬 계수를 저장하면 너무 많은 메모리가 필요하다는 것을 알고 있습니다. 그래서 이제는 매트릭스에 동일한 계수가 많은 시간을 나타 내기 때문에 매트릭스없는 방법을 사용하려고합니다. 그래서 개인 스토리지 체계를 사용하고 캐시 효율성을 높일 수 있습니다. 나는 매트릭스가없는 선형 연산자에 대한 인터페이스를 제공하는 petsc를보고 있지만 실제로 이해하지 …

9 linear-solver  petsc  sparse-operator 

내 프로그래밍의 대부분은 내가 사용하는 C의 일회용 연구 코드입니다. 가까운 공동 작업자 이외의 다른 코드를 배포 한 적이 없습니다. 과학 저널에 게시하는 알고리즘을 개발했습니다. 기사의 온라인 보충 자료에 소스 코드와 실행 코드를 제공하고 싶습니다. 동료가 알고리즘을 일반화하여 C ++ (ack!)로 작성해야하고 작은 조밀 한 선형 시스템을 풀어야한다고 요구했습니다. 알고리즘에 대한 …

9 linear-solver  libraries  c++  c