3 각 행렬의 고유 시스템에 대한 병렬 알고리즘


11

나는 큰 희소 행렬 (~ 2 백만 요소)의 Lanczos 대각선 화를하고 있습니다. Lanzcos 알고리즘의 거의 모든 단계는 Lanczos 행렬을 대각선으로 모아 수렴을 확인하는 것을 제외하고 GPU에서 병렬로 수행됩니다. 이를 위해 Numerical Recipes의 TQLI 알고리즘을 사용하고 있습니다. 평행하거나 쉽게 병렬화 할 수있는 3 각형 행렬의 고유 시스템을 찾는 방법이 있습니까? 병렬 버전의 TQLI가 있습니까?

답변:


4

SLEPc 과 같은 라이브러리를 사용하는 것이 좋습니다.이 라이브러리 에는 고유 시스템을 직렬 또는 병렬로 해결하기위한 다양한 방법에 대한 인터페이스가 포함되어 있습니다. 사용자 설명서는 고유치 문제를 해결하기위한 여러 가지 방법에 대한 참조를 포함합니다.


실제로, 기존 희소 고유 해석기는 Rayleigh 몫에 대해 평행 선형 대수를 사용하지 않습니다. 나는 이번 여름에 그러한 고유 해석기를 썼지 만 불행히도 비공개 소스입니다.
잭 폴슨

9

TQL은 병렬화 할 수 없습니다.

표준 병렬 알고리즘은 Cuppen의 표준 알고리즘입니다.

JJM Cuppen, 대칭 삼각형 고유 문제의 분할 및 정복 방법, 1980.
http://www.springerlink.com/content/t21365q2gh702714/

또한보십시오:

F. Tisseur, 분산 메모리 아키텍처에서 대칭 고유 값 문제에 대한 병렬 분할 및 정복 알고리즘, 1999
http://eprints.ma.man.ac.uk/981/01/covered/MIMS_ep2007_225.pdf

http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.75.4109&rep=rep1&type=pdf

http://www14.in.tum.de/konferenzen/Jass09/courses/2/Kleine_Albers_paper.pdf


Arvo 링크는 이제 슬프게도 깨졌습니다. :(
Geoffrey Irving

@ GeoffreyIrving : 모든 사람에게 무료는 아니지만 작동하는 것으로 교체했습니다. 그리고 나는 Tisseur의 논문에 대한 새로운 참조를 추가했습니다.
Arnold Neumaier

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