Pertanyaan ini merupakan perpanjangan dari dua diskusi yang muncul baru-baru ini di balasan untuk " C ++ vs Fortran untuk HPC ". Dan ini sedikit lebih menantang daripada sebuah pertanyaan ... Salah satu argumen yang paling sering terdengar dalam mendukung Fortran adalah bahwa penyusunnya lebih...
Saya telah mencari ke dalam perpustakaan aljabar C ++ linier untuk proyek yang telah saya kerjakan. Sesuatu yang saya masih belum pahami adalah koneksi BLAS dan LAPACK ke perpustakaan aljabar linier lainnya. Melihat melalui artikel ini di perpustakaan aljabar linier, saya merasa menarik...
Dengan dua implementasi BLAS yang berbeda, dapatkah kita berharap mereka membuat perhitungan floating point yang sama persis dan mengembalikan hasil yang sama? Atau bisa terjadi, misalnya, bahwa salah satu Toedjoe menghitung produk skalar sebagai dan satu sebagai ( x 1 y 1 + x 2 y 2 ) + ( x 3 y...
Saya ingin tahu apakah ada cara cepat untuk menghitung jarak Euclidean dari dua vektor dalam Oktaf. Sepertinya tidak ada fungsi khusus untuk itu, jadi haruskah saya menggunakan rumus saja
Maaf untuk posting lama tapi saya ingin memasukkan semua yang saya pikir relevan pada langkah pertama. Apa yang saya inginkan Saya menerapkan versi paralel dari Metode Subruang Krylov untuk Matriks Padat.Terutama GMRES, QMR dan CG. Saya menyadari (setelah membuat profil) bahwa rutinitas DGEMV...
Di mana pun saya melihat, tutorial / dokumen dll PETSc mengatakan bahwa itu berguna untuk aljabar linier dan biasanya menentukan bahwa sistem jarang akan mendapat manfaat. Bagaimana dengan matriks padat? Saya khawatir tentang menyelesaikan untuk padat .A x = bAx=bAx=bSEBUAHAA Saya telah menulis...
nVidia, misalnya, memiliki CUBLAS, yang menjanjikan peningkatan 7-14x. Naifnya, ini sama sekali tidak mendekati throughput teoritis kartu GPU nVidia. Apa saja tantangan dalam mempercepat aljabar linier pada GPU, dan apakah sudah ada rute aljabar linier yang lebih
Implementasi netlib BLAS adalah referensi yang sangat baik, sebagian besar tidak dioptimalkan dan didokumentasikan dengan baik (misalnya zgemm ). Namun, itu di Fortran 77, membuatnya agak tidak dapat diakses oleh mereka yang memiliki pendidikan pemrograman yang lebih modern. Apakah ada implementasi...
Produk dalam kompleks memiliki dua definisi berbeda yang diputuskan oleh konvensi: atau . Di BLAS, saya menemukan rutinitas cdotu, zdotu, dan cdotc, zdotc. Dua rutinitas yang pertama sebenarnya menghitung (produk dalam palsu!) Dan dua rutinitas terakhir mengkonjugasikan vektor pertama dalam produk...
Apakah kompilasi PETSc dengan perpustakaan BLAS / LAPACK eksternal secara signifikan mempengaruhi kinerja pada matriks jarang, atau apakah hanya menggunakan perpustakaan tersebut untuk matematika matriks