Insinyur sering bersikeras menggunakan metode konservatif lokal seperti volume hingga, perbedaan hingga konservatif, atau metode Galerkin diskontinyu untuk menyelesaikan PDE. Apa yang salah ketika menggunakan metode yang tidak konservatif secara lokal? Oke, jadi konservasi lokal itu penting untuk...
Adakah yang bisa merekomendasikan saya sistem manajemen alur kerja yang baik (WMS), lebih disukai di Python? Sejauh ini saya telah menggunakan GNU Make, tetapi memperkenalkan lapisan kompleksitas yang ingin saya hindari. WMS yang baik harus memiliki fitur-fitur berikut: mengintegrasikan dengan...
Saya punya pertanyaan sederhana yang benar-benar sulit bagi Google (di samping kanonik Apa Yang Harus Diketahui Setiap Ilmuwan Komputer tentang kertas Aritmatika Titik Apung ). Kapan seharusnya fungsi seperti log1patau expm1digunakan, bukan logdan exp? Kapan mereka tidak digunakan? Bagaimana...
Baru-baru ini saya menemukan formulasi dari fenomena-meta : " dua itu mudah, tiga itu sulit " (diucapkan seperti ini oleh Federico Poloni), yang dapat dijelaskan, sebagai berikut: Ketika masalah tertentu dirumuskan untuk dua entitas, itu relatif mudah untuk diselesaikan; Namun, algoritma untuk...
Jurnal Fisika Komputasi telah menjadi outlet penting untuk ilmu komputasi di masa lalu, dan saya telah menerbitkan di sana sebelumnya. Demi kepentingan orang-orang (seperti saya) yang telah menandatangani boikot Elsevier , jurnal non-Elsevier apa yang akan menjadi tempat yang tepat untuk...
det(A)≈0det(A)≈0\det(A) \approx 0 Apakah kebalikannya juga benar? Apakah matriks yang dikondisikan buruk memiliki determinan yang hampir nol? Ini adalah sesuatu yang saya coba di Octave: a = rand(4,4); det(a) %0.008 cond(a)%125 a(:,4) = 1*a(:,1) + 2*a(:,2) =
Sebagai seseorang yang memegang gelar BA dalam fisika saya agak tersinggung ketika saya mulai bekerja dengan simulasi molekuler. Agak mengejutkan mengetahui bahwa bahkan simulasi yang paling terperinci dan mahal secara komputasi tidak dapat secara kuantitatif mereproduksi perilaku penuh air dari...
Latar Belakang: Saya pikir saya mungkin ingin mem-porting beberapa kode yang menghitung produk matriks eksponensial-vektor menggunakan metode subruang Krylov dari MATLAB ke Python. (Secara khusus, fungsi expmvp Jitse Niesen , yang menggunakan algoritma yang dijelaskan dalam makalah ini .) Namun,...
Masalah Saat ini saya sedang mengerjakan simulasi Elemen Hingga Navier Stokes dan saya ingin menyelidiki efek dari berbagai parameter. Beberapa parameter ditentukan dalam file input atau melalui opsi baris perintah; parameter lain disediakan sebagai flag dalam Makefile sehingga kode saya harus...
Versi 3.0 dari standar MPI secara formal menghapus antarmuka C ++ (sebelumnya tidak digunakan lagi). Sementara implementasi mungkin masih mendukungnya, fitur-fitur yang baru di MPI-3 tidak memiliki antarmuka C ++ yang didefinisikan dalam standar MPI. Lihat
Dalam statistik dan berbagai aplikasinya, kita sering menghitung matriks kovarians , yang pasti positif (dalam kasus yang dipertimbangkan) dan simetris, untuk berbagai kegunaan. Terkadang, kita membutuhkan invers dari matriks ini untuk berbagai perhitungan (bentuk kuadratik dengan invers ini...
Pertanyaan ini dimigrasikan dari Stack Overflow karena dapat dijawab di Computational Science Stack Exchange. Bermigrasi 7 tahun yang lalu . Adakah masalah / algoritme terkenal dalam komputasi ilmiah yang tidak dapat dipercepat dengan paralisis? Sepertinya saya ketika
EDIT: Saya menguji apakah ada nilai eigen yang besarnya satu atau lebih. Saya perlu menemukan nilai eigen absolut terbesar dari matriks besar non-simetris yang jarang. Saya telah menggunakan eigen()fungsi R , yang menggunakan algo QR baik dari EISPACK atau LAPACK untuk menemukan semua nilai...
Saya tahu, dan telah menggunakan f2py2e untuk membungkus beberapa kode Fortran 77 lama, tetapi pemahaman saya adalah bahwa itu tidak bekerja dengan kode Fortran 95 yang lebih baru. Saya telah meneliti apa yang harus saya gunakan, dan telah menemukan fwrap dan G3 f2py, yang keduanya tidak memberikan...
Saya ingin memodifikasi matriks transisi persegi yang padat di tempat dengan mengubah urutan beberapa baris dan kolomnya, menggunakan perpustakaan numpy python. Secara matematis ini berhubungan dengan pra-mengalikan matriks dengan matriks permutasi P dan mengalikannya dengan P ^ -1 = P ^ T, tetapi...
Apa yang bisa salah ketika menggunakan metode Krylov yang dikondisikan sebelumnya dari KSP ( paket solver linier PETSc ) untuk menyelesaikan sistem linier yang jarang seperti yang diperoleh dengan mengkritisi dan membuat linierisasi persamaan diferensial parsial? Langkah apa yang bisa saya ambil...
Apakah ada alat yang dirancang khusus untuk mengompresi data ilmiah floating point? Jika suatu fungsi halus, jelas ada banyak korelasi antara angka-angka yang mewakili fungsi itu, sehingga data harus dikompres dengan baik. Zipping / gzipping data floating point biner tidak kompres dengan baik....
Salah satu masalah utama yang harus kita hadapi dalam simulasi molekuler adalah perhitungan gaya yang bergantung pada jarak. Jika kita dapat membatasi fungsi gaya dan jarak agar memiliki kekuatan genap dari jarak pemisahan , maka kita bisa menghitung kuadrat jarak r 2 = r ⋅ r dan tidak perlu...
Saya tidak begitu akrab dengan skema diskritisasi umum untuk PDE. Saya tahu bahwa Crank-Nicolson adalah skema populer untuk mendiskritisasi persamaan difusi. Apakah juga pilihan yang baik untuk masa adveksi? Saya menarik dalam menyelesaikan persamaan Reaksi-Difusi-Adveksi , ∂kamu∂t+ ∇ ⋅ ( v u - D...
Pengembang perangkat lunak memiliki pilihan untuk memilih lisensi yang sesuai sesuai dengan tujuan pekerjaan. Adakah yang bisa memberikan rekomendasi / pengalaman tentang lisensi mana yang harus dipilih untuk perangkat lunak? Apa pro / kontra dari "memberikan" semua pekerjaan yang dikodekan...