Saya sedang mengerjakan beberapa sistem persamaan linear besar jarang yang tidak dikondisikan. Saya ingin menggunakan aritmatika double-double atau aritmatika quad-double untuk menyelesaikannya. Saya tahu bahwa ada paket bernama MPACK yang dikembangkan oleh Nakata, Maho, yang dapat melakukan perhitungan aljabar linear numerik di bawah aritmatika quad-double. Namun, ini dirancang untuk matriks padat, bukan matriks jarang. Apakah Anda tahu apakah ada paket matriks aritmatika quad-double?
sparse-matrix
matrix
condition-number
Hanyu Ye
sumber
sumber
Jawaban:
Pada versi 3.2 , PETSc mendukung perhitungan presisi quad jarang di gcc / gfortran 4.6 dan yang lebih baru.
Anda akan membutuhkan BLAS dan LAPACK quad-precision, yang dapat disediakan PETSc kepada Anda (bersama dengan dukungan quad) dengan perintah konfigurasi (parsial) berikut:
Lihat FAQ untuk informasi lebih lanjut.
Juga, saya setuju dengan nOOb , jika memungkinkan, cobalah untuk mengatur sistem sebelum beralih ke quad-precision.
sumber
Sistem yang dikondisikan buruk lebih baik diselesaikan dengan regularisasi daripada dengan meningkatkan presisi numerik. Cari "regularisasi ill-posed" untuk detail berdarah.
sumber
Saya mungkin akan mencoba perpustakaan Trilinos. Mereka telah mem-templar perpustakaan matriks jarang di bawah Tpetra (yang seharusnya menggantikan Epetra, perpustakaan matriks jarang asli mereka). Anda dapat templat ganda, kompleks, quad, dll, dan mereka mungkin memiliki pilihan solver terbesar (baik langsung dan berulang) di sebelah PETSc.
Sunting: setelah membaca komentar, kegunaan langsung Tpetra tampaknya agak dipertanyakan presisi wrt quad ...
sumber
__float128
didukung oleh Tpetra, apalagi quad-double. Tpetra tidak berdiri sendiri dan tidak semua header, dan bahkan jika itu, hal-hal sepertistd::complex
hanya bekerja denganfloat
dandouble
.Kotak alat komputasi multiprecision untuk MATLAB memiliki dukungan untuk matriks jarang dan secara khusus dioptimalkan untuk komputasi dengan presisi empat kali lipat.
Berikut ini rincian waktu untuk pemecah jarang presisi quad: Pemecah Langsung untuk Matriks Jarang
(Saya penulis toolbox).
sumber