Jadi saya mengerti (atau setidaknya saya percaya saya tahu) bagaimana siklus-V berjalan. Saya telah menulis dalam Matlab versi 1-D, versi rekursif dari siklus-V. Namun, ketika saya menjalankan kode saya untuk FMG, solusi saya tidak konvergen. Saya percaya masalah saya terletak pada pemahaman saya...
Saya perlu memecahkan sistem linear yang jarang sama (300x300 hingga 1000x1000) dengan banyak sisi kanan (300 hingga 1000). Selain masalah pertama ini, saya juga ingin menyelesaikan sistem yang berbeda, tetapi dengan elemen bukan nol yang sama (hanya nilai yang berbeda), yaitu banyak sistem jarang...
Saya menggunakan MATLAB untuk memecahkan masalah yang melibatkan penyelesaian di setiap timestep, di mana berubah seiring waktu. Saat ini, saya menyelesaikan ini menggunakan MATLAB :bAx=bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}bb\mathbf{b}mldivide x = A\b Saya memiliki fleksibilitas untuk membuat...
Saya memiliki koleksi model komputasi yang dapat digambarkan sebagai automata seluler asinkron. Model-model ini menyerupai model Ising, tetapi sedikit lebih rumit. Sepertinya model seperti itu akan mendapat manfaat dari dijalankan pada GPU daripada CPU. Sayangnya itu tidak cukup mudah untuk...
Saya porting kode yang ada dari MATLAB ke C ++ dan memiliki sistem linier untuk menyelesaikan (daripada bentuk yang lebih khas A x = b )xA=bxA=bxA=bAx=bAx=bAx=b Matriks padat, dan berbentuk umum, tetapi tidak lebih besar dari 1000x1000. Jadi dalam MATLAB, solusinya ditemukan oleh fungsi atau...
Saya mensimulasikan aliran yang tidak dapat dimampatkan di atas silinder pada bilangan Reynold 500. Saya memecahkan persamaan nokes stokes menggunakan metode koreksi tekanan. Solusi saya menjadi tidak stabil setelah waktu tertentu (sekitar 5 detik). Saya telah mencoba memperbaiki mesh saya,...
Ada beberapa keluarga metode untuk quadrature numerik. Jika saya memiliki kelas khusus integrand, bagaimana cara memilih metode yang ideal? Apa pertanyaan yang relevan untuk ditanyakan baik tentang integrand (misalnya apakah lancar? Apakah memiliki singularitas?) Dan masalah komputasi (misalnya...
Misalkan saya memiliki masalah nilai batas: d2udx2+dvdx=f in Ωd2udx2+dvdx=f in Ω\frac{d^2u}{dx^2} + \frac{dv}{dx}=f \text{ in } \Omega dudx+d2vdx2=g in Ωdudx+d2vdx2=g in Ω\frac{du}{dx} +\frac{d^2v}{dx^2} =g \text{ in } \Omega u=h in ∂Ωu=h in ∂Ωu=h \text{ in } \partial\Omega Tujuan saya adalah...
Ketika diskritisasi FEM dan menyelesaikan masalah difusi reaksi, misalnya, dengan (perturbasi singular), solusi dari masalah diskrit biasanya akan menunjukkan lapisan berosilasi dekat dengan batas. Dengan , dan elemen hingga linier, solusi terlihat seperti0 < ε ≪ 1 Ω = ( 0 , 1 ) ε = 10 - 5 u...
Sistem saya adalah masalah FE simetris dengan pengganda lagrange (misalnya aliran Stokes yang tidak dapat dikompresi): ( ABBTC)(SEBUAHBTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} di mana adalah kasus khas (saya bahkan telah memastikan bahwa persamaan diberi nomor sehingga pengali Lagrange...
Hai dan terima kasih telah meluangkan waktu untuk melihat pertanyaan saya. Ini adalah versi terbaru dari pertanyaan saya yang saya posting sebelumnya di physics.stackexchange.com Saat ini saya sedang mempelajari spin exciton 2D Bose-Einstein Condensate dan saya ingin tahu tentang keadaan dasar...
Saya mencoba mencari tahu apakah ada cara yang lebih cepat untuk menghitung semua nilai eigen dan vektor eigen dari matriks adjacency yang sangat besar dan jarang daripada menggunakan scipy.sparse.linalg.eigsh Sejauh yang saya tahu, metode ini hanya menggunakan sparseness dan atribut simetri dari...
Dalam persamaan Navier-Stokes yang tidak dapat dimampatkan, istilah tekanan sering disebut sebagai pengganda Lagrange yang menegakkan kondisi tidak tertahankan.ρ(ut+(u⋅∇)u)∇⋅u=−∇p+μΔu+f=0ρ(ut+(u⋅∇)u)=−∇p+μΔu+f∇⋅u=0\begin{align*} \rho\left(\mathbf{u}_t + (\mathbf{u} \cdot \nabla)\mathbf{u}\right)...
Saya tertarik secara global untuk memaksimalkan fungsi dari banyak ( ) parameter nyata (hasil dari simulasi yang kompleks). Namun, fungsi yang dimaksud relatif mahal untuk dievaluasi, membutuhkan sekitar 2 hari untuk setiap set parameter. Saya membandingkan opsi yang berbeda, dan bertanya-tanya...
Kami sedang mengerjakan model Bayesian untuk proses ruang-waktu, dan menggunakan No-U-Turn sampler (NUTS) yang membutuhkan model untuk probabilitas log dan gradiennya sehubungan dengan parameter model. Lebih ringkasnya, kami memiliki fungsi probabilitas log yang cukup rumit , yang melibatkan...
Biarkan dan f ( → x ) : [ 0 , 1 ] n → C menjadi fungsi dalam variabel-variabel ini.x⃗ =(x1,x2,…,xn)∈[0,1]nx→=(x1,x2,…,xn)∈[0,1]n\vec{x} = (x_1, x_2, \dots, x_n) \in [0,1]^nf(x⃗ ):[0,1]n→Cf(x→):[0,1]n→Cf(\vec{x}): [0,1]^n \to \mathbb{C} Apakah ada skema rekursif untuk integral iterasi...
Katakanlah saya memiliki versi kode yang lambat dan cepat, dan ingin melaporkan nomor percepatan yang membandingkan keduanya. Saya menjalankan versi lambat kali dan versi cepat m kali, menghasilkan kali ( s 1 , ... , s n ) dan ( f 1 , ... , f m ) . Cara paling sederhana untuk menghasilkan speedup...
Ini mungkin pertanyaan yang lebih cocok untuk sisi Rekomendasi Perangkat Lunak dari SE, namun saya percaya bahwa orang yang sering menggunakan bagian SE ini cenderung menjawab pertanyaan ini. Saya mencari alternatif gratis (tidak hanya dalam kebebasan) untuk Comsol Multiphysics . Berikut ini...
Jika saya memecahkan masalah yang relatif kecil, yaitu, masalah yang dapat ditangani dengan metode langsung seperti LU, lalu apakah nomor kondisi operator linier memengaruhi keakuratan solusi? Salah satu masalah penelitian yang saya kerjakan berfokus pada pengembangan teknik optimisasi untuk...
Saya baru saja membaca Bab 3 dalam "A Multigrid Tutorial" oleh Briggs / Henson / McCormick, tautan . Teksnya tentang siklus Multigrid seperti siklus-V, siklus-mu, FMG. Apa yang menarik perhatian saya: Dalam kebanyakan prosedur iteratif seseorang memeriksa apakah ia telah menyatu dengan toleransi /...