Saya harus menulis kode volume terbatas untuk Magnetohydrodynamics (MHD). Saya telah menulis kode numerik sebelumnya tetapi tidak pada skala ini. Saya hanya ingin bertanya yang akan menjadi pilihan yang baik, menggunakan struktur data (pendekatan berorientasi objek) dengan kelas atau hanya...
Untuk sebuah proyek, saya harus mengimplementasikan kedua metode ini dan membandingkan kinerjanya pada fungsi yang berbeda. Sepertinya metode gradien konjugasi dimaksudkan untuk menyelesaikan sistem persamaan linear untuk Ax=bAx=b A\mathbf{x} = \mathbf{b} Di mana adalah matriks n-by-n yang...
Saya punya pertanyaan yang mirip dengan yang diajukan sebelumnya kecuali dalam 3D, dan saya hanya perlu volume, bukan bentuk lambung yang sebenarnya. Lebih tepatnya, saya diberi satu set kecil poin (katakanlah, 10-15) dalam 3D, yang semuanya diketahui terletak pada cembung lambung set point...
Pertimbangkan sistem linear tridiagonal pasti positif simetris mana A ∈ R n × n dan b ∈ R n . Diberikan tiga indeks 0 ≤ i < j < k < n , jika kita mengasumsikan hanya baris persamaan antara i dan k tahan, kita dapat menghilangkan variabel perantara untuk mendapatkan persamaan bentuk u x...
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 punya pertanyaan mengenai kondisi batas pengkodean untuk mekanika solid (elastisitas linier). Dalam kasus khusus saya harus menggunakan perbedaan hingga (3D). Saya sangat baru dalam topik ini, jadi mungkin beberapa pertanyaan berikut mungkin sangat mendasar. Untuk mengarah ke masalah spesifik...
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...
Saya mencoba membalikkan transformasi integral berikut secara numerik: F(y) = ∫∞0yexp[ - 12( y2+ x2) ] Saya0( X y) f( x )d xF(y)=∫0∞yexp[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Jadi untuk saya perlu...
Saya mencoba merencanakan fungsi gelombang untuk sebuah partikel dalam kotak 3D. Ini mengharuskan saya untuk memplot 4 variabel: sumbu x, y, z dan fungsi kepadatan probabilitas. Fungsi kepadatan probabilitas
Apakah ada algoritma SVD terpotong yang menghitung nilai singular satu per satu? Masalah saya: Saya ingin menghitung nilai singular pertama (dan vektor singular) dari matriks padat , tapi saya tidak tahu apa nilai tepat. itu besar, jadi untuk alasan efisiensi, saya lebih suka tidak mengevaluasi...
Adakah yang bisa merekomendasikan metode untuk masalah kuadrat-terkecil berikut: R ∈ R3 × 3R∈R3×3R \in \mathbb{R}^{3 \times 3}∑i = 0N( R xsaya- bsaya)2→ min∑saya=0N(Rxsaya-bsaya)2→min\sum\limits_{i=0}^N (Rx_i - b_i)^2 \rightarrow \minRRR Saya bisa mendapatkan solusi perkiraan dengan meminimalkan...
Saya membaca di suatu tempat bahwa Pandas pertama kali dikembangkan untuk dunia keuangan, setidaknya tidak terutama untuk ilmu alam (fisika, biologi, dan sebagainya), jadi apakah ada analisis data serupa paket Python yang lebih “berorientasi ilmu alam”? Saya baru saja mulai menggunakan Pandas dan...
Saya mengambil kursus perhitungan ilmiah, dan kami baru saja mendekati pendekatan kuadrat. Pertanyaan saya secara khusus tentang perkiraan menggunakan polinomial. Saya mengerti bahwa jika Anda memiliki titik data n +1, Anda dapat menemukan polinomial tingkat n yang unik yang menggambarkan semua...
Ketika seseorang ingin menghitung turunan numerik, metode yang disajikan oleh Bengt Fornberg di sini (dan dilaporkan di sini ) sangat mudah (baik tepat dan mudah diimplementasikan). Sebagai kertas asli tanggal dari 1988, saya ingin tahu apakah ada alternatif yang lebih baik hari ini (sebagai (atau...
Metode Jacobian-Free Newton-Krylov (JFNK), dan metode Krylov secara umum, dapat sangat berguna karena mereka tidak memerlukan penyimpanan eksplisit atau konstruksi matriks, hanya hasil dari produk vektor-matriks. Jika Anda benar-benar membentuk sistem sparse, ada banyak prasyarat di luar sana untuk...
Saya memiliki daftar dari matriks simetris yang saya perlu periksa untuk semi-definiteness positif (yaitu nilai eigennya non-negatif.)LL{\cal L} Komentar di atas menyiratkan bahwa seseorang dapat melakukannya dengan menghitung nilai eigen masing-masing dan memeriksa apakah nilainya tidak negatif...
Bagaimana keadaan terkini dari seni untuk memecahkan PDE parabola dimensi tinggi (3-10) dalam domain kompleks dengan kutub sederhana (dari bentuk ) dan menyerap kondisi batas?1|r⃗ 1−r⃗ 2|1|r→1−r→2| \frac{1}{|\vec{r}_1 - \vec{r}_2|} Secara khusus, saya tertarik untuk menyelesaikan persamaan...
Saya mencoba memahami bagaimana metode pengoptimalan berbasis adjoint bekerja untuk optimisasi terbatas PDE. Khususnya, saya mencoba memahami mengapa metode adjoint lebih efisien untuk masalah di mana jumlah variabel desain besar, tetapi "jumlah persamaan kecil". Apa yang saya
pengantar Ijinkan saya menyatakan beberapa pernyataan yang bertentangan tentang masalah ini untuk menggambarkan apa masalahnya. Secara pribadi saya ingin kode saya terbuka di setiap tahap pengembangan, sejak itu yang lain akan melihat dan memanfaatkan apa yang saya lakukan Saya juga suka...
Saya ingin menyelesaikan beberapa PDE pada manifold, katakanlah misalnya persamaan elips pada bola. Di mana saya memulai? Saya ingin menemukan sesuatu yang digunakan sudah ada sebelumnya kode / perpustakaan di 2d, tidak ada yang begitu mewah (untuk saat ini) Ditambahkan Kemudian: Artikel dan...