Có bất kỳ gói ma trận thưa số học bốn nhân đôi?

Tôi đang làm việc trên một số hệ phương trình tuyến tính thưa thớt lớn điều hòa. Tôi muốn sử dụng số học hai lần hoặc số học gấp đôi để giải quyết chúng. Tôi biết rằng có một gói có tên MPACK được phát triển bởi Nakata, Maho, có thể thực hiện các phép tính đại số tuyến tính theo số học bốn nhân đôi. Tuy nhiên, nó được thiết kế cho ma trận dày đặc, không phải ma trận thưa thớt. Bạn có biết liệu có bất kỳ gói ma trận thưa số học bốn nhân đôi nào không?

Kể từ phiên bản 3.2 , PETSc hỗ trợ tính toán độ chính xác tứ giác thưa thớt trên gcc / gfortran 4.6 và mới hơn.

Bạn sẽ cần một BLAS và LAPACK bốn độ chính xác, mà PETSc có thể cung cấp cho bạn (cùng với hỗ trợ quad) với lệnh cấu hình (một phần) sau:

./configure --with-precision=__float128 --download-f2cblaslapack

Xem FAQ để biết thêm thông tin.

Ngoài ra, tôi đồng ý với nOOb , nếu có thể, hãy cố gắng thường xuyên hóa hệ thống trước khi chuyển sang độ chính xác bốn.

Các hệ thống điều hòa không khí được giải quyết tốt hơn bằng cách chính quy hơn là tăng độ chính xác bằng số. Tìm kiếm "chính quy hóa sai" để biết chi tiết về tin đồn.

Tôi có thể thử thư viện Trilinos. Họ có các thư viện ma trận thưa thớt theo Tpetra (được cho là thay thế Epetra, thư viện ma trận thưa thớt ban đầu của họ). Bạn có thể tạo mẫu gấp đôi, phức tạp, quad, v.v. và chúng có thể có nhiều lựa chọn bộ giải nhất (cả trực tiếp và lặp) bên cạnh PETSc.

Chỉnh sửa: sau khi đọc các bình luận, tính hữu dụng ngay lập tức của Tpetra dường như là một chút nghi vấn chính xác ...

Hộp công cụ tính toán đa điểm cho MATLAB có hỗ trợ cho ma trận thưa thớt và được tối ưu hóa đặc biệt cho các tính toán với độ chính xác gấp bốn lần.

Dưới đây là chi tiết thời gian cho bộ giải thưa thưa chính xác bốn phần: Bộ giải trực tiếp cho ma trận thưa

(Tôi là tác giả của hộp công cụ).