Adakah orang di sini yang menggunakan komputasi ilmiah presisi ganda dengan GPU generasi baru (misalnya K20) melalui Python? Saya tahu bahwa teknologi ini berkembang pesat, tetapi apa cara terbaik untuk melakukan ini saat ini? GPU berada di luar ruang lingkup untuk perpustakaan ilmiah Python...
Dalam solusi numerik PDE nilai batas awal, sangat umum untuk menggunakan paralelisme dalam ruang . Adalah jauh lebih jarang untuk menggunakan beberapa bentuk paralelisme dalam diskritisasi waktu , dan bahwa paralelisme biasanya jauh lebih terbatas. Saya menyadari semakin banyak kode dan karya yang...
Sumber istilah, seperti yang karena batimetri dalam persamaan air dangkal, perlu diintegrasikan dengan cara khusus untuk menjaga keadaan fisik stabil. Apakah ada cara umum untuk membangun metode yang seimbang, atau apakah itu memerlukan teknik khusus untuk setiap
Saya ingin menulis solver saya sendiri untuk persamaan Euler yang dapat dikompres, dan yang paling penting saya ingin itu bekerja dengan kuat dalam semua situasi. Saya ingin ini berbasis FE (DG tidak apa-apa). Apa metode yang mungkin? Saya sadar melakukan urutan ke-0 DG (volume terbatas) dan itu...
Baru-baru ini, saya telah mendorong grup saya untuk memasukkan lebih banyak pengujian saat menulis kode mereka. Ada beberapa bug utama yang membutuhkan waktu lebih lama untuk ditangkap daripada yang mungkin diperlukan, karena kami tidak memiliki rezim pengujian yang baik. Namun, saya menduga bahwa...
Saya mengalami masalah dalam mengimplementasikan fungsi secara numerik. Itu menderita dari kenyataan bahwa pada nilai input besar hasilnya adalah jumlah yang sangat besar kali jumlah yang sangat kecil. Saya tidak yakin apakah pembatalan katastropik adalah istilah yang tepat, jadi mohon perbaiki...
Saya mengalami beberapa frustrasi atas cara matlab menangani integrasi numerik vs Scipy. Saya mengamati perbedaan berikut dalam kode pengujian saya di bawah ini: Versi Matlab berjalan rata-rata 24 kali lebih cepat dari yang setara dengan python saya! Versi Matlab mampu menghitung integral tanpa...
Ketika mengevaluasi jumlah FLOPs dalam fungsi yang sederhana, kita seringkali hanya bisa menurunkan ekspresi penghitungan operator aritmatika dasar. Namun, dalam kasus pernyataan matematis yang melibatkan pembagian yang merata, seseorang tidak dapat melakukan ini dan berharap dapat membandingkan...
Saya mencoba belajar tentang menyelesaikan sendiri PDE secara numerik. Saya sudah mulai dengan metode beda hingga (FDM) untuk beberapa waktu karena saya mendengar bahwa FDM adalah dasar dari banyak metode numerik untuk PDE. Sejauh ini saya punya beberapa pemahaman dasar untuk FDM dan bisa menulis...
Saya telah membaca buku Computer Simulation of Liquids oleh Allen dan Tildesley. Mulai dari halaman 71, penulis membahas berbagai algoritma yang digunakan untuk mengintegrasikan persamaan gerak Newton dalam simulasi dinamika molekul (MD). Mulai dari halaman 78, penulis membahas algoritma Verlet,...
Saya ingin memodelkan aliran air laminar dari akar ke batang tanaman. Di ujung akar, ukuran dan panjang tabung bervariasi dari milimeter hingga sentimeter. Saat kami semakin dekat ke batang, akarnya semakin besar panjang dan diameternya. Saya ingin membuat domain 3D acak yang mewakili jaringan akar...
Saya bertanya-tanya apakah algoritma Thomas adalah cara tercepat (dapat dibuktikan?) Untuk menyelesaikan sistem tridiagonal simetris dominan diagonal dalam hal kompleksitas algoritmik (tidak mencari paket implementasi seperti LAPACK dll). Saya tahu bahwa kedua algoritma Thomas dan multigrid adalah...
Saya telah menggunakan metode yang berbeda untuk menghitung peringkat matriks dan solusi sistem persamaan matriks. Saya menemukan fungsi linalg.svd. Membandingkan ini dengan upaya saya sendiri untuk memecahkan sistem dengan Gaussian Elimination, tampaknya lebih cepat dan lebih tepat. Saya mencoba...
Katakanlah saya menjalankan komputasi superkomputer pada 100k core selama 4 jam di http://www.nersc.gov/users/computational-systems/edison/configuration , bertukar sekitar 4 PB data melalui jaringan dan melakukan sekitar 4 TB I / HAI. Penghitungan semuanya bilangan bulat, sehingga hasilnya benar...
Saya mencoba untuk menyelesaikan persamaan MHD Ideal menggunakan metode semi diskrit, rekonstruksi spasial ENO dan loncatan waktu TVD RK. Saya mendapatkan solusi kondisi tunak yang berbeda dengan urutan waktu yang berbeda. Apakah itu
Saya mencoba memahami persamaan advection dengan koefisien kecepatan variabel sedikit lebih baik. Khususnya saya tidak mengerti bagaimana persamaan bisa konservatif. The persamaan adveksi , ∂u∂t+∂∂x(vu)=0∂u∂t+∂∂x(vu)=0 \frac{\partial u}{\partial t} + \frac{\partial}{\partial x}(\boldsymbol{v}u) =...
Dalam metode seperti gmres atau bicgstab, mungkin menarik untuk menggunakan metode krylov lain sebagai prasyarat. Bagaimanapun mereka mudah diimplementasikan dalam cara bebas-matriks dan dalam lingkungan paralel. Misalnya, satu dapat menggunakan beberapa (misalkan ~ 5) iterasi bigcstab tanpa syarat...
Menurut Wikipedia , laju konvergensi dinyatakan sebagai rasio spesifik dari norma-norma vektor. Saya mencoba memahami perbedaan antara tingkat "linear" dan "kuadrat", pada titik waktu yang berbeda (pada dasarnya, "di awal" dari iterasi, dan "di akhir"). Mungkinkah dinyatakan
Judulnya adalah pertanyaan. Teknik ini melibatkan penggunaan "matriks kofaktor", atau "matriks adjugate", dan memberikan formula eksplisit untuk komponen kebalikan dari matriks kuadrat. Tidak mudah dilakukan dengan tangan untuk sebuah matriks yang lebih besar dari, katakanlah, 3×33×33\times 3 ....
Kesimpulannya di sini: Seberapa jauh sebenarnya penyusun Fortran? adalah gfortran dan gcc yang cepat untuk kode sederhana. Jadi saya ingin mencoba sesuatu yang lebih rumit. Saya mengambil contoh baku tembak norma spektral. Saya pertama-tama menghitung ulang matriks 2D A (:, :), dan kemudian...