Saya mencoba mencari tahu apakah ada cara yang lebih cepat untuk menghitung semua nilai eigen dan vektor eigen dari matriks adjacency yang sangat besar dan jarang daripada menggunakan scipy.sparse.linalg.eigsh Sejauh yang saya tahu, metode ini hanya menggunakan sparseness dan atribut simetri dari...
Dalam kedua metode dekomposisi domain (DD) dan multigrid (MG), seseorang dapat menyusun aplikasi pembaruan blok atau koreksi kasar baik sebagai aditif atau multiplikatif . Untuk pemecah poin point, ini adalah perbedaan antara iterasi Jacobi dan Gauss-Seidel. Multiplicative smoother untuk bertindak...
Misalkan saya memiliki dua matriks Nx2, Mx2 masing-masing mewakili N, M 2d vektor. Apakah ada cara sederhana dan bagus untuk menghitung jarak antara masing-masing pasangan vektor (n, m)? Cara yang mudah namun tidak efisien tentu saja: d = zeros(N, M); for i = 1:N, for j = 1:M, d(i,j) =...
Komputasi paralel sering dimodelkan menggunakan laju perhitungan lokal deterministik, overhead latensi, dan bandwidth jaringan. Pada kenyataannya, ini adalah variabel spasial dan non-deterministik. Studi seperti Skinner dan Kramer (2005) mengamati distribusi multi-modal, tetapi analisis kinerja...
Ini adalah masalah enumerasi simetri sederhana. Saya memberikan latar belakang penuh di sini, tetapi tidak ada pengetahuan tentang kimia kuantum diperlukan. Kedua partikel terpisahkan adalah: ⟨ i j | k l ⟩ = ∫ ψ * i ( x ) ψ * j ( x ' ) ψ k ( x ) ψ l ( x ' )⟨ I j | k l ⟩⟨ij|kl⟩\langle ij|kl\rangle...
Pertimbangkan dua komputer dengan konfigurasi perangkat keras dan perangkat lunak yang berbeda. Saat menjalankan kode Navier-Stokes serial yang sama persis pada setiap platform, diperlukan waktu x dan y untuk menjalankan satu iterasi untuk komputer 1 dan 2, secara berurutan. Dalam kasus ini, ,...
Saya mencari literatur paling penting dan referensi slide untuk memodelkan biaya energi saat ini dan masa depan dari operasi floating-point dan transfer data di seluruh CPU, memori, jaringan, dan penyimpanan. Saya telah menandai pertanyaan ini sebagai wiki komunitas, dan saya lebih suka Anda...
Apakah perkalian matriks (baik Mat * Mat, dan Mat * Vec) skala dengan jumlah non-nol, atau dengan ukuran matriks? Atau kombinasi keduanya. Bagaimana dengan bentuk. Sebagai contoh, saya memiliki 100 x 100 matriks dengan 100 nilai di dalamnya, atau 1000 x 1000 matriks dengan 100 nilai di...
C ++ 11 memperkenalkan langkah semantik yang dapat, misalnya, meningkatkan kinerja kode dalam situasi di mana C ++ 03 perlu melakukan konstruksi salin atau tugas penyalinan. Artikel ini melaporkan bahwa kode berikut mengalami kecepatan 5x saat dikompilasi dengan C + 11: vector<vector<int>...
Saya harus membuat transformasi koordinat antara dua sistem referensi (sumbu). Untuk itu, tiga matriks ( ) harus dikalikan karena beberapa sumbu perantara yang digunakan. Saya telah memikirkan dua pendekatan untuk menyelesaikan ini:3×33×33\times3 Metode # 1 : Membuat perkalian secara langsung,...
Saya mencoba memahami perbedaan antara kedua kartu grafis ini untuk komputasi akademik, khususnya untuk komponen DGEMM. Jika kita melihat statistik mentah, keduanya memiliki chip GK110 yang sama, memiliki statistik yang sebanding di hampir setiap kategori, dan, saya percaya, memiliki arsitektur...
Saya memiliki masalah di mana saya harus menemukan semua positif (seperti dalam nilai eigen positif) pasang eigen dari matriks nonsimetrik kecil (biasanya lebih kecil dari 60x60). Saya bisa berhenti menghitung ketika nilai eigen lebih kecil dari ambang tertentu. Saya tahu bahwa nilai eigennya...
double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Di atas...