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...
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...
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
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...
Saya punya proyek di mana saya perlu menggunakan bidang kuadrat Khususnya angka dari bentuk dengansebuah,b∈Q.a + b - 3---√a+b−3a + b \sqrt{-3}a,b∈Qa,b∈Qa,b \in \mathbb{Q} Misalnya di sini adalah bilangan prima dalam bilangan bulat Eisenstein : Saya tidak ingin menggunakan bijak. Saya ingin...
Dalam laporan teknis tentang Galahad [1], penulis menyatakan, dalam konteks masalah pemrograman nonlinier umum, Menurut pendapat kami, tidak pernah ada banyak keraguan bahwa metode SQP [sequential quadratic programming] akan lebih berhasil [daripada metode Augmented Lagrangian] dalam jangka...
Saya melakukan pencarian garis sebagai bagian dari algoritma BFGS kuasi-Newton. Dalam satu langkah pencarian baris saya menggunakan interpolasi kubik untuk bergerak lebih dekat ke minimizer lokal. Biarkan menjadi fungsi yang menarik. Saya ingin mencari x ∗ sehingga f ′ ( x ∗ ) ≈ 0 .f: R → R , f∈...
Saya ingin membuat model pancing (atau tali) dengan bergabung dengan segmen pendek. (Segmen mungkin memiliki panjang yang sama (pendek) tetapi setiap segmen harus diberi massa masing-masing.) Satu segmen akan mempengaruhi berikutnya dengan torsi antar segmen. Untuk saat ini sambungan dapat dianggap...
Saya tertarik dalam menghitung solusi sistem lage ODE menggunakan metode krylov seperti pada [1]. Metode semacam itu melibatkan fungsi-fungsi yang berkaitan dengan fungsi eksponensial (yang disebut fungsi ). Ini pada dasarnya terdiri dari menghitung tindakan fungsi matriks dengan membangun ruang...
Beberapa buku yang lebih tua yang pernah saya lihat mengatakan bahwa jumlah minimum tahapan metode Runge-Kutta eksplisit dari pesanan tertentu tidak diketahui untuk pesanan . Apakah ini masih benar?≥ 9≥9\geq 9 Perpustakaan apa yang ada untuk bekerja secara otomatis dengan metode Runge-Kutta...
Saya bertanya-tanya apakah ada yang punya pengalaman berurusan dengan batas-batas ketika menerapkan diferensiasi chebyshev. Saat ini saya sedang mencoba untuk menerapkan syarat batas tanpa selip untuk menyelesaikan persamaan Navier Stokes yang tidak dapat dikompres dalam 3D, untuk memastikan...
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...
Sebagian besar metode untuk integral osilasi yang saya ketahui berhubungan dengan integral bentuk mana berukuran besar.∫f(x)eiωxdx∫f(x)eiωxdx \int f(x)e^{i\omega x}\,dx ωω\omega Jika saya memiliki integral dari bentuk mana adalah fungsi berosilasi yang akarnya hanya diketahui sekitar, tetapi...
Pertimbangkan situasi di mana Anda ingin menyelesaikan sistem linier menggunakan metode Krylov yang telah dikondisikan sebelumnya, tetapi menerapkan prekondisi itu sendiri melibatkan penyelesaian sistem bantu, yang dilakukan dengan metode Krylov lain yang telah dikondisikan sebelumnya. Pada satu...
Saya telah menemukan bahwa metode garis adalah cara yang sangat alami untuk berpikir tentang diskritisasi PDE. Karena itu saya selalu default ke pola pikir itu ketika disajikan dengan set persamaan baru. Saya belum pernah melihat PDE di mana ini tidak akan berhasil. Yang saya ingin tahu adalah...
Dalam literatur FEM, metode semi-variasional biasanya digunakan dalam solusi PDE bergantung waktu. Saya belum melihat pendekatan yang sepenuhnya variasional yaitu di mana ruang dan waktu didiskritkan oleh FEM, mungkin memungkinkan penggunaan jerat ruang-waktu yang tidak terstruktur. Meskipun metode...
Pertimbangkan, Anda memiliki masalah dalam ruang Hilbert atau Banach dimensi tak terbatas (pikirkan PDE atau masalah optimisasi dalam ruang seperti itu) dan Anda memiliki algoritme yang konvergennya lemah ke suatu solusi. Jika Anda mendiskritisasi masalah dan menerapkan algoritma diskretisasi...
Saya bertanya-tanya bagaimana kondisi batas Dirichlet dalam matriks elemen hingga global jarang benar-benar diterapkan secara efisien. Misalnya, katakanlah matriks elemen hingga global kami adalah: K= ⎡⎣⎢⎢⎢⎢⎢⎢520- 102410001632- 1037000203⎤⎦⎥⎥⎥⎥⎥⎥dan vektor sisi kananb = ⎡⎣⎢⎢⎢⎢⎢⎢b 1b 2b 3b 4b...
Pendahuluan multigrid biasanya menggunakan kotak persegi panjang. Interpolasi nilai kemudian lurus ke depan: Hanya interpolasi secara linear di tepi antara dua node yang berdekatan dari grid kasar untuk menemukan nilai dari node grid halus di tepi itu. Untuk aplikasi FEM saya memiliki kotak yang...
Latar Belakang: Saya hanya membangun satu solusi numerik yang berfungsi untuk 2d Navier-Stokes, untuk kursus. Itu adalah solusi untuk aliran rongga yang digerakkan oleh tutup. Kursus ini, bagaimanapun, membahas sejumlah skema untuk diskritisasi spasial dan diskritisasi waktu. Saya juga telah...