PETSc có bao giờ sử dụng các thư viện LAPACK cho toán học ma trận thưa thớt không?

Việc biên dịch PETSc với thư viện BLAS / LAPACK bên ngoài có ảnh hưởng đáng kể đến hiệu suất trên các ma trận thưa thớt không, hay nó chỉ sử dụng các thư viện đó cho toán học ma trận dày đặc?

PETSc sử dụng BLAS cho một vài nguyên hàm vector, nhưng chúng thường bị giới hạn bởi băng thông bộ nhớ và không có nhiều sự khác biệt trong "tối ưu hóa", do đó nó có xu hướng không tạo ra nhiều khác biệt về hiệu năng.

Nó cũng sử dụng Lapack cho một số phân tích như ước tính của Lanczos hoặc Arnoldi về giá trị bản địa và giá trị số ít, nhưng chúng thường không nhạy cảm với hiệu suất.

Các hoạt động "cấp 3" dày đặc thường chỉ xuất hiện trong bối cảnh nhạy cảm về hiệu năng khi sử dụng các bộ giải trực tiếp thưa thớt từ các thư viện của bên thứ ba (ví dụ: MUMPS, SuperLU, UMFPACK), trong trường hợp đó, việc lấp đầy cuối cùng dẫn đến các vấn đề dày đặc đủ lớn để hưởng lợi từ gọi BLAS.

Nếu bạn dựa vào các bộ giải trực tiếp thưa thớt này được áp dụng cho các vấn đề lớn, thì đáng để xây dựng với triển khai BLAS được điều chỉnh, nếu không, nó sẽ tạo ra rất ít sự khác biệt.