«lapack» 태그된 질문

LAPACK (Linear Algebra PACKage)은 선형 방정식 세트, 선형 최소 제곱, 고유 값 문제 및 특이 값 분해 솔루션을 포함하여 수치 선형 대수 작업에 일반적으로 사용되는 서브 루틴 라이브러리입니다. LAPACK 루틴은 포트란, C 및 친척 및 다양한 기타 언어와 함께 사용할 수 있습니다.

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

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

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


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

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

2
Matlab [Q, R, E] = qr (A) 뒤에 대응하는 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 

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

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



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

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

3
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 …

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

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.