Ada banyak skema FD untuk persamaan advection diskusikan di web. Misalnya di sini: http://farside.ph.utexas.edu/teaching/329/lectures/node89.html∂T∂t+u∂T∂x=0∂T∂t+u∂T∂x=0\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}=0 Tapi saya belum melihat ada yang mengusulkan skema angin...
apa metode canggih untuk solusi numerik ODE dengan sisi kanan terputus-putus? Saya sebagian besar tertarik sedikit demi sedikit fungsi sisi kanan, misalnya tanda. Saya mencoba memecahkan persamaan dari tipe berikut:
Saya tertarik untuk menyelesaikan persamaan Poisson menggunakan pendekatan finite-difference. Saya ingin lebih memahami bagaimana menulis persamaan matriks dengan kondisi batas Neumann. Akankah seseorang meninjau yang berikut, apakah itu benar? Matriks beda hingga Persamaan
Saya menjalankan simulasi dinamika molekul (MD) menggunakan beberapa paket perangkat lunak, seperti Gromacs dan DL_POLY. Gromacs sekarang mendukung algoritma dekomposisi partikel dan domain dekomposisi. Secara default, simulasi Gromacs menggunakan dekomposisi domain, meskipun selama...
Apakah ada algoritma (efisien) untuk memilih himpunan bagian dari titik dari satu set titik ( ) sedemikian rupa sehingga mereka "menutupi" sebagian besar wilayah (di atas semua himpunan bagian ukuran )?MM.MNNNM< NM<NM < NM.M.M Saya berasumsi poinnya ada di bidang 2D. Algoritme naif itu...
Biarkan x, ymenjadi dua angka floating-point. Apa cara yang benar untuk menghitung artinya? Cara naif (x+y)/2dapat menghasilkan luapan saat xdan yterlalu besar. Saya pikir 0.5 * x + 0.5 * ymungkin lebih baik, tetapi melibatkan dua perkalian (yang mungkin tidak efisien), dan saya tidak yakin apakah...
Apakah ada kasus di mana Hartree-Fock bukanlah pendekatan yang baik untuk menghitung geometri keseimbangan ketika molekul berada dalam kondisi tidak mengikat?
Saya baru dalam simulasi dinamika molekul (MD). Apa kompleksitas simulasi dinamika molekul dalam hal waktu simulasi? Dengan kata lain, jika saya ingin meningkatkan waktu simulasi dari 10 nanodetik ke 20 nanodetik, apa yang bisa saya harapkan dalam hal peningkatan
Metode multigrid biasanya menyelesaikan masalah Dirichlet pada level (mis. Titik Jacobi atau Gauss-Seidel). Ketika menggunakan metode elemen hingga terus menerus, jauh lebih murah untuk merakit masalah Neumann kecil daripada merakit masalah Dirichlet kecil. Metode dekomposisi domain yang tidak...
Saya bekerja cukup banyak dalam biologi matematika / epidemiologi, di mana sebagian besar karya ilmu pemodelan / komputasi masih didominasi oleh set ODE, diakui kadang-kadang set cukup rumit dari mereka. Salah satu kelebihan model-model ini adalah model-model itu agak mudah digambarkan dan ditiru....
Simulasi dinamika plasma kami sering menghasilkan terlalu banyak informasi. Selama simulasi, kami merekam berbagai properti fisik pada kisi (x, y, z, t) yang sebesar (8192x1024x1024x1500), untuk setidaknya 10 properti. Informasi ini diproses setelah simulasi selesai. Dengan itu kita membuat film...
Saya membandingkan dua metode berulang untuk membalikkan matriks kuadrat acak. Karena matriksnya acak, setiap test case mengambil jumlah iterasi dan waktu yang berbeda yang berbeda. Pertanyaan saya adalah, di atas rata-rata waktu CPU, adalah nilai rata-rata dari iterasi yang diambil oleh kedua...
Saya ingin menghitung spektrum ( semua nilai eigen) dari matriks jarang besar (ratusan ribu baris). Ini sulit. Saya bersedia menerima perkiraan. Apakah ada metode perkiraan untuk melakukan ini? Sementara saya berharap jawaban umum untuk pertanyaan ini, saya juga akan puas dengan jawaban dalam...
Dalam simulasi semikonduktor, adalah umum bahwa persamaan diskalakan sehingga memiliki nilai yang dinormalisasi. Sebagai contoh, dalam kasus ekstrim kerapatan elektron dalam semikonduktor dapat bervariasi di atas 18 orde magnitudo, dan medan listrik dapat berubah bentuk, lebih dari 6 (atau lebih)...
Seperti yang saya pahami, karena solusi untuk program linier selalu terjadi pada titik set yang layak polyhedral (jika ada solusi dan nilai fungsi objektif yang optimal dibatasi dari bawah, dengan asumsi masalah minimisasi), bagaimana bisa pencarian melalui interior wilayah yang layak menjadi lebih...
Saya baru saja mulai bermain-main dengan FEniCS. Saya memecahkan Poisson dengan elemen urutan ke-3 dan ingin memvisualisasikan hasilnya. Namun, ketika saya menggunakan plot (u), visualisasi hanyalah interpolasi linier dari hasilnya. Saya mendapatkan hal yang sama ketika saya output ke VTK. Dalam...
Dalam praktiknya, runtime penyelesaian IVP ˙ x ( t ) = f ( t , x ( t ) ) secara numerik x˙(t)=f(t,x(t)) for t∈[t0,t1]x˙(t)=f(t,x(t)) for t∈[t0,t1] \dot{x}(t) = f(t, x(t)) \quad \text{ for } t \in [t_0, t_1] x(t0)=x0x(t0)=x0 x(t_0) = x_0 sering didominasi oleh lamanya mengevaluasi sisi kanan...
Saya memecahkan sistem dua digabungkan PDE dalam dua dimensi spasial dan waktu komputasi. Karena evaluasi fungsi mahal, saya ingin menggunakan metode multistep (diinisialisasi menggunakan Runge-Kutta 4-5). Metode Adams-Bashforth menggunakan lima evaluasi fungsi sebelumnya memiliki kesalahan global...
Banyak pustaka C ++ yang lebih terkenal dalam ilmu komputasi seperti Eigen , Trilinos , dan deal.II menggunakan objek pustaka header template C ++ standar std::complex<>,, untuk mewakili angka floating-point yang kompleks. Dalam jawaban Jack Poulson untuk pertanyaan tentang konstruktor...
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...