«lapack» 태그된 질문

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

9 linear-solver  sparse-matrix  lapack  condition-number 

내 프로젝트에서 매 단계마다 두 개의 삼각 행렬을 풀어야하므로 좋은 솔버를 갖는 것이 중요합니다. 나는 Wikipedia에 설명 된 고전적인 방법 인 내 자신의 구현을했습니다. 그런 다음 Lapack을 대신 사용해 보았습니다. 놀랍게도 느려졌습니다! 이제 Lapack 내부에서는 LU 인수 분해로 해결하는 것처럼 보이며 왜 그렇게 복잡하지 않은지 궁금합니다. 또한 nr.com의 "Numerical Recipes"책에서 …

9 banded-matrix  lapack