Ilmu Komputasi

Saya mencari solusi eksplisit cepat (berani saya katakan optimal?) 3x3 masalah nyata, , . A x = bSEBUAHx=b\mathbf{A}\mathbf{x} = \mathbf{b}A ∈ R3 × 3, b ∈ R3SEBUAH∈R3×3,b∈R3\mathbf{A} \in \mathbf{R}^{3 \times 3}, \mathbf{b} \in \mathbf{R}^{3} Matriks bersifat umum, tetapi dekat dengan matriks identitas dengan angka kondisi mendekati 1. Karena sebenarnya pengukuran …

9 linear-solver  reference-request 

Banyak pustaka FEM adaptif menggunakan struktur data mesh yang lebih canggih untuk menangani penambahan / penghapusan node, edge, triangles, tetrahedra, dll. Misalnya, p4est library menggunakan struktur data octree untuk perbaikan mesh adaptif; Anda tidak akan sering menemukan oktri yang digunakan untuk perhitungan pada jaring statis. Apa perubahan pada sisi aljabar …

9 linear-algebra  finite-element  sparse-matrix 

Saya selalu mendengar bahwa paralelisasi yang mudah adalah salah satu keuntungan dari metode DG, tetapi saya tidak benar-benar melihat mengapa salah satu dari alasan itu tidak berlaku untuk Galerkin terus menerus.

9 finite-element  parallel-computing  discontinuous-galerkin 

Saya memiliki masalah nilai eigen kubik besar: ( A0+ λ A1+ λ2SEBUAH2+ λ3SEBUAH3) x =0.(SEBUAH0+λSEBUAH1+λ2SEBUAH2+λ3SEBUAH3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Aku bisa mengatasi ini dengan mengkonversi ke masalah nilai eigen linear tetapi akan menghasilkan sistem sebagai besar:32323^2 ⎡⎣⎢- A0000saya000saya⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥= λ ⎡⎣⎢SEBUAH1saya0SEBUAH20sayaSEBUAH300⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥,[-SEBUAH0000saya000saya][xyz]=λ[SEBUAH1SEBUAH2SEBUAH3saya000saya0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 \\ …

9 c++  eigensystem  eigenvalues 

double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Di atas adalah kode C ++ …

9 c++  performance 

Apa perbedaan utama antara FEM dan XFEM? Kapan kita (tidak) menggunakan XFEM intead dari FEM dan sebaliknya? Dengan kata lain, ketika saya menemukan masalah baru, bagaimana saya bisa tahu menggunakan yang mana dari mereka?

9 finite-element  numerical-analysis  adaptive-mesh-refinement  numerics 

Pertanyaan-pertanyaan ini mungkin agak di luar topik di komputer. jika diperlukan tolong sarankan di mana itu cocok dengan. Pertanyaannya adalah tentang bagaimana mengelola semua simulasi berjalan secara efisien. katakanlah, misalnya, simulasi memerlukan perbaikan 2 parameter yang harus ditentukan pada rentang nilai yang disarankan. Untuk menemukan hasil yang lebih baik yang …

9 monte-carlo  data-management 

Pertimbangkan masalah optimisasi cembung yang tidak dibatasi ketatBiarkan x_ \ text {opt} menunjukkan minima uniknya dan x_0 menjadi perkiraan awal yang diberikan kepada x_ \ text {opt}. Kami akan memanggil vektor x solusi \ epsilon- close dari \ mathcal {O} jika \ begin {persamaan} \ frac {|| x - x …

9 convex-optimization  complexity 

Apakah ada aturan sederhana untuk mengatakan apakah itu layak untuk dilakukan SOR daripada Gauss-Seidel? (dan cara yang mungkin bagaimana memperkirakan parameter realxation )ωω\omega Maksud saya hanya dengan melihat pada matriks , atau pengetahuan tentang masalah tertentu yang diwakili oleh matriks? Saya membaca jawaban atas pertanyaan ini: Apakah ada heuristik untuk …

9 linear-solver  iterative-method  heuristics 

Saya menulis paket MATLAB kecil yang akan memecahkan masalah numerik kelas tertentu. Ada 3 tahap algoritma dan pengguna memiliki 5 pilihan untuk setiap tahap. Saya telah menerapkan seluruh masalah menggunakan fungsi dan 3 kasus sakelar (satu untuk setiap tahap algoritma). Ini berfungsi dengan baik tetapi saya mempertimbangkan membuatnya melakukan lebih …

9 matlab 

Saya telah mencoba untuk men-debug kesalahan ini beberapa hari terakhir saya bertanya-tanya apakah ada yang punya saran tentang cara melanjutkan. Saya memecahkan persamaan Poisson untuk distribusi muatan langkah (masalah umum dalam fisika elektrostatik / semikonduktor) pada mesh volume terbatas yang tidak seragam di mana yang tidak diketahui didefinisikan pada pusat …

9 discretization  finite-volume  poisson 

Saya bertanya-tanya apakah ada beberapa milis yang sangat bagus atau grup google untuk Ilmu Komputasi, di mana kita dapat mendiskusikan pertanyaan alih-alih hanya bertanya dan menjawab pertanyaan. Bahkan, saya lebih tertarik dengan komputasi paralel dan solusi numerik dari PDE. Tapi saya tidak tahu apa dan bagaimana orang di daerah ini …

9 parallel-computing  reference-request 

Saya menerapkan metode beda hingga ke sistem 3 persamaan digabungkan. Dua persamaan tidak digabungkan, namun persamaan ketiga berpasangan dengan dua lainnya. Saya perhatikan bahwa dengan mengubah urutan persamaan, katakan dari menjadi ( x , z , y ) bahwa matriks koefisien menjadi simetris.( x , y, z)(x,y,z)(x, y, z)(x,z,y)(x,z,y)(x, z, …

9 finite-difference  symmetry 

Saya bertanya-tanya apa cara terbaik untuk belajar bagaimana memprogram di Fortran? Saya telah mencari buku-buku, contoh-contoh yang tampak oke bagi saya adalah "Fortran 95/2003 Dijelaskan" dan "Fortran Untuk Ilmuwan Dan Insinyur" . Saya tidak tahu apa yang kalian pikirkan tentang ini? Hal lain yang saya ingin tahu adalah kompiler mana …

9 fortran 

Saya telah menerapkan pemecah multigrid V-Cycle menggunakan koreksi cacat linier (LDC) dan skema aproksimasi penuh (FAS). Masalah saya adalah sebagai berikut: Menggunakan LDC residual dikurangi dengan faktor ~ 0,03 per siklus. Implementasi FAS memang konvergen dengan faktor linear juga, tetapi faktornya hanya ~ 0,58. Jadi FAS membutuhkan sekitar 20 kali …

9 multigrid