Ilmu Komputasi

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 level referensi BLAS, seperti netlib, di C / C …

11 blas  education 

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 (arbitrer ), dengan mengambil matriks dan:∑i = 0N( A xsaya- bsaya)2→ min∑saya=0N(SEBUAHxsaya-bsaya)2→min\sum\limits_{i=0}^N (Ax_i - …

11 linear-algebra  least-squares  svd 

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 memperkirakan f ( x ) di mana:F( y)F(y)F(y)f( x )f(x)f(x) dan F ( y ) adalah nyata …

11 algorithms  numerical-analysis  approximation  inverse-problem 

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 titik ini. Tapi saya juga bisa melihat mengapa ini …

11 matlab  least-squares 

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 adalah: abs((np.sin((p*np.pi*X)/a))*(np.sin((q*np.pi*Y)/b))*(np.sin((r*np.pi*Z)/c)))**2 Saya menggunakan np.arange()untuk X, Y dan Z. Saya telah membaca bahwa untuk melakukan ini, Anda perlu memplot permukaan plot …

11 probability  plotting  computational-physics 

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, , adalah perbedaan waktu iterasi antara komputer 1 …

11 performance  iterative-method  navier-stokes 

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 saya, pertama-tama saya ingin menunjukkan apa yang sudah saya terapkan …

11 finite-difference  boundary-conditions 

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 …

11 linear-algebra  poisson 

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 menggunakan beberapa array untuk properti yang berbeda, dalam hal kecepatan, …

11 python  finite-volume  data-structures 

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 hampir sama) sederhana dan lebih …

11 finite-difference  precision 

Saya memiliki dataset titik 3-dimensi yang ingin saya buat mesh, menggunakan python. Semua perangkat lunak yang saya lihat mengharuskan Anda memberikan tepi. Apakah ada program dalam python yang mengambil input satu set poin dalam 3D dan output jala? Jika memungkinkan, saya ingin sambungannya seragam.

11 python  mesh-generation 

Di tempat kerja, saya pada dasarnya berfungsi sebagai konsultan independen. Untuk manajemen dan pelanggan, saya perlu memperkirakan jumlah waktu yang diperlukan untuk mengembangkan perangkat lunak sebagai bagian dari penelitian ilmu komputer saya. Namun, perkiraan waktu saya biasanya tidak aktif. Saya tahu bahwa ada metode untuk memperkirakan waktu yang diperlukan untuk …

11 software 

12×1212×1212 \times 12QQQJdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Saat ini saya melakukan ini dengan perpustakaan armadillo tetapi ternyata terlalu lambat. Masalahnya adalah bahwa saya perlu melakukan ini untuk satu triliun matriks dan ternyata menghitung dua penentu adalah hambatan dari program saya. Karena itu saya punya dua pertanyaan Apakah ada trik …

11 linear-algebra  matrix  c++  graph-theory  c 

Saat ini saya sedang mengembangkan metode dekomposisi domain untuk solusi dari masalah hamburan. Pada dasarnya saya memecahkan sistem Helmholtz BVP secara iteratif. Saya discretize persamaan menggunakan metode elemen hingga lebih dari segitiga atau tetrahedral jerat. Saya sedang mengembangkan kode menuju tesis Phd saya. Saya menyadari beberapa pustaka elemen hingga yang …

11 finite-element  parallel-computing  mesh-generation  domain-decomposition  mesh 

Dimulai dengan persamaan advection dalam bentuk konservasi. ut=(a(x)u)xut=(a(x)u)x u_t = (a(x)u)_x di mana adalah kecepatan yang bergantung pada ruang, dan adalah konsentrasi spesies yang dilestarikan.ua(x)a(x)a(x)uuu Mendiskritkan fluks (di mana fluks , didefinisikan pada tepi sel di antara titik-titik mesh) memberi, u t = 1f=a(x)uf=a(x)uf=a(x)uut=1h(fj−12−fj+12)ut=1h(fj−12−fj+12) u_t = \frac{1}{h}\left( f_{j-{\frac{1}{2}}} - f_{j+{\frac{1}{2}}} …

11 discretization  finite-volume  numerical-analysis