행렬 라이브러리는 행렬 어떤 식으로 표현, 관리하냐?

[익명(175.223)]

그냥 2차원 배열?

[커헠(msca8h)]

1차원 배열 ㅇㅇ. 좋은 라이브러리는 패딩까지 하기도 하고. 어쨋든 1차원 배열에 인덱스 차원변환해서 쓴다 - return 0;

[커헠(msca8h)]

보통은 포트란 과의 호환성 때문에 칼럼 메이저로 많이 쓴다. numpy 는 특이하게 row major 긴 한데 그 때문에 성능을 조금 희생하는편 - return 0;

[익명(118.40)]

대각 행렬이나 희소 행렬 같은 것도 마찬가지로 처리함?

[커헠(msca8h)]

대각 행렬은 라이브러리마다 다른데 특수취급하는 경우에는 대각만 1차원 배열로 저장하고 대각행렬에 대한 연산도 전부 따로 만들어놓음 - return 0;

[커헠(msca8h)]

희소의 경우 여러가지 방법들이 공존하는데 CPU 에서 연산 돌릴 경우 CRS, COO 포맷을 주로 쓴다. 전통적으로는 CRS 가 강세 - return 0;