«lapack» 태그된 질문

내가 작업 한 프로젝트의 C ++ 선형 대수 라이브러리를 살펴 보았습니다. 여전히 이해하지 못하는 것은 BLAS와 LAPACK을 다른 선형 대수 라이브러리에 연결하는 것입니다. 선형 대수 라이브러리에 대한이 기사를 살펴보면 흥미로운 점을 발견했습니다. 일부 라이브러리는 BLAS 및 LAPACK과 독립적입니다. 일부는 BLAS 및 LAPACK이 필요합니다 일부는 BLAS 및 LAPACK에 대한 선택적 인터페이스가 …

26 linear-algebra  software  lapack  blas 

긴 게시물에 대해 죄송하지만 첫 번째 이동과 관련이 있다고 생각되는 모든 것을 포함하고 싶었습니다. 내가 원하는 것 고밀도 행렬에 대한 병렬 버전의 Krylov Subspace Methods를 구현하고 있습니다. 주로 GMRES, QMR 및 CG. 나는 (프로파일 링 후) 내 DGEMV 루틴이 한심하다는 것을 깨달았다. 그래서 나는 그것을 격리함으로써 그것에 집중하기로 결정했습니다. 12 …

15 parallel-computing  matrix  lapack  blas 

내가 작업중 인 소프트웨어 프로젝트에서 밀도가 낮은 하위 행렬의 경우 특정 계산이 훨씬 쉽습니다. 일부 문제 인스턴스에는 밀도가 낮은 하위 행렬이 포함되어 있지만 요인이 아닌 전체가 제공되므로 하위 순위 구조를 활용하려면 순위를 확인하고 행렬을 인수 분해해야합니다. . 문제의 행렬은 일반적으로 완전히 또는 거의 완전히 밀도가 높으며 n은 100에서 수천까지입니다. 행렬의 …

13 linear-algebra  matrix  lapack  rank  matrix-factorization 

ATLAS가 컴파일 된 머신에 대해 자체 최적화 할 수 있으므로 소스에서 컴파일하면 최대의 이점을 얻을 수 있습니다. 소스에서 LAPACK을 컴파일하면 어떤 이점이 있습니까? 사전 빌드 된 패키지를 설치하는 것이 훨씬 쉬울 것입니다.

13 performance  lapack 

행렬의 순위와 방정식의 행렬 시스템의 솔루션을 모두 계산하기 위해 다른 방법을 사용했습니다. linalg.svd 기능을 발견했습니다. 이것을 가우시안 제거 (Gaussian Elimination)로 시스템을 해결하려는 나의 노력과 비교하면 더 빠르고 정확한 것으로 보입니다. 이것이 어떻게 가능한지 이해하려고합니다. 내가 아는 한 linalg.svd 함수는 QR 알고리즘을 사용하여 행렬의 고유 값을 계산합니다. 나는 이것이 수학적으로 어떻게 …

13 linear-algebra  python  lapack 

I는 C ++로 MATLAB에서 기존 코드를 포팅하고 해결하는 선형 시스템이 오전 (오히려 더 일반적인 형태보다 X = B를 )xA=bxA=bxA=bAx=bAx=bAx=b 행렬 는 밀도가 높고 일반적인 형태이지만 1000x1000보다 크지 않습니다. MATLAB에서 솔루션은 함수 또는 슬래시 표기법으로 찾을 수 있습니다.AAAmrdivide(b,A)x = b/A; BLAS 및 LAPACK 루틴을 사용하여 C ++ 코드에서이 문제를 어떻게 해결해야합니까? …

12 linear-solver  c++  lapack 

나는 현재 행렬 대한 좋은 순위 추정치를 저렴하게 계산하려고합니다 . 따라서 나는 사용하여 columnt 피벗 QR 분해를 계산합니다.AAA [Q,R,E]=qr(A) Matlab에서. 사용하고 있는 의 랭크를 추정AAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) 이것은 잘 작동하며 R의 모든 대각선 항목에 대한 플롯은 다음과 같습니다. RRR 입력 행렬은 두 실험에서 모두 동일합니다. 내 …

12 matlab  lapack 

예를 들어, nVidia에는 ​​CUBLAS가있어 속도가 7-14 배 향상됩니다. 기본적으로 이것은 nVidia의 GPU 카드의 이론적 인 처리량에 거의 근접하지 않습니다. GPU에서 선형 대수의 속도를 높이는 데 어떤 어려움이 있으며, 이미 사용 가능한 빠른 선형 대수 라우팅이 있습니까?

11 linear-algebra  lapack  blas  gpu 

몇 가지 양자 역학 시뮬레이션을 위해 포트란 95에 집중하고 있습니다. 솔직히 옥타브에 의해 망쳐 졌으므로 당연히 행렬 지수를 취했습니다. 크기가 인 (작고 ) skew- Hermitian 행렬이 주어지면 LAPACK을 사용하여이 문제를 해결하는 가장 효율적인 방법은 무엇입니까? LAPACK95 래퍼를 사용하지 않고 LAPACK을 직접 호출합니다.n ≤ 36n≤36n\leq 36n × nn×nn\times n

11 fortran  lapack 

모든 실수 행렬 는 직교 치열 변환 U를 사용하여 실수 형태 T = U T A U 로 축소 될 수 있습니다 . 여기서 행렬 T 는 주 대각선에 1 x 1 또는 2 x 2 블록의 준 삼각형입니다. 각 1 x 1 블록은 A 의 실수 고유 값에 해당 …

10 eigenvalues  eigensystem  lapack  numerics  scalapack 

외부 BLAS / LAPACK 라이브러리를 사용하여 PETSc를 컴파일하면 희소 행렬의 성능에 큰 영향을 미치거나 밀도가 높은 행렬 수학에 해당 라이브러리 만 사용합니까?

10 petsc  sparse-matrix  lapack  blas 

최근 에 skew-Hermitian 행렬에 대해 같은 줄을 따라 질문했습니다 . 그 질문의 성공에 영감을 받아 벽에 머리를 두 시간 동안 두드린 후 실제 비대칭 행렬의 행렬 지수를보고 있습니다. 고유 값과 고유 벡터를 찾는 경로는 다소 복잡해 보이므로 잃어버린 것 같습니다. 배경 : 얼마 전에 나는 이론 물리 물리학 SE에 대해이 …

10 linear-algebra  fortran  lapack 

저는 컴퓨터 과학에 익숙하지 않고 이미 C ++에서 RK4, Numerov 등과 같은 통합, 보간법, 방법에 대한 기본 방법을 배웠지 만 최근 교수님은 매트릭스와 관련된 문제를 해결하기 위해 LAPACK을 사용하는 방법을 배우도록 요청했습니다. 예를 들어 복소수 행렬의 고유 값을 찾는 것과 같습니다. 나는 타사 라이브러리를 사용한 적이 없으며 거의 ​​항상 내 …

10 lapack 

LAPACK의 QR 루틴은 Q를 가정용 반사기로 저장합니다. 반사 벡터의 크기를 조정합니다VVv 와 1 /V11/V11/v_1결과의 첫 번째 요소는 111저장하지 않아도됩니다. 그리고 별도 저장ττ\tau필요한 스케일 팩터를 포함하는 벡터. 리플렉터 매트릭스는 다음과 같습니다.H= 나는− τVV티,H=나는−τVV티,H=I-\tau v v^T, 어디 VVv정규화되지 않았습니다. 교과서에서 반사기 매트릭스는 H= 나는− 2 vV티,H=나는−2VV티,H = I-2vv^T, 여기서 는 정규화됩니다.VVv LAPACK …

9 linear-algebra  matrix  lapack 

제어 시스템 도구 상자를 처음부터 순수하게 Python3 (shameless plug :)으로 작성했습니다 harold. 과거의 연구에서 나는 care.m기술적 / 관련이없는 이유로 Riccati 솔버 에 대해 항상 불만을 제기했습니다 . 따라서 나는 내 자신의 루틴 세트를 작성해 왔습니다. 내가 찾을 수없는 한 가지 방법은 적어도만큼 좋은 고성능 밸런싱 알고리즘을 얻는 것 balance.m입니다. 언급하기 …

9 linear-algebra  matlab  lapack  scipy  matrix-equations