Dalam makalahnya http://www.sciencedirect.com/science/article/pii/S0045782509003521 , persamaan elemen-lokal HDG dijelaskan pada halaman 584 persamaan (4), dengan salah satu persamaan mengambil bentuk berikut −(uh,∇q)K=−⟨u^h⋅n,q−q¯⟩∂K−(uh,∇q)K=−⟨u^h⋅n,q−q¯⟩∂K-(u_h,\nabla q)_K =...
Dalam penelitian tentang masalah invers, adalah umum untuk membangun set data sintetik dari set parameter yang diketahui dan kemudian menguji apakah teknik inversi dapat merekonstruksi parameter tersebut. Dengan demikian, penting untuk menambahkan tingkat kebisingan acak yang sesuai ke data...
Biarkan menjadi nyata, kotak, matriks padat. G dan Q simetris. MembiarkanA , G , QA,G,QA, G, QGGGQQQ H= [ A- Q- G- AT]H=[A−G−Q−AT]H = \begin{bmatrix} A & -G \\ -Q &-A^T \end{bmatrix} menjadi matriks Hamilton. Saya ingin menghitung eksponensial matriks . Saya membutuhkan eksponensial matriks...
Setiap matriks nyata dapat direduksi menjadi bentuk Schur nyata T = U T A U menggunakan transformasi U yang mirip ortogonal . Di sini matriks T adalah bentuk kuasi-segitiga dengan 1 oleh 1 atau 2 oleh 2 blok pada diagonal utama. Setiap 1 oleh 1 blok berkorespondensi dengan nilai eigen nyata A dan...
Ada algoritma untuk memecahkan fungsi Odes pegangan , di manay∈Rn. Tetapi dalam banyak sistem fisik, persamaan diferensial adalah otonom, jadidydydt= f( y, t )dydt=f(y,t)\frac{dy}{dt} = f(y, t)y∈Rny∈Rny \in \mathbb R^n,y∈Rn, dengant yangtertinggal. Dengan asumsi yang disederhanakan ini, perbaikan...
Saya menemukan bahwa banyak pemrograman sains komputasi saya memiliki persyaratan pengujian yang tidak dicakup oleh kerangka uji standar: Pengujian waktu komputasi Untuk memastikan bahwa algoritma tidak menjadi lebih lambat. Saya bisa melakukan sesuatu seperti
Saya memecahkan multi-spesies, persamaan Navier-Stokes yang dapat dikompres pada kisi terstruktur 3D. Saya telah mendapatkan solusi pada grid yang diberikan (misalkan yang relatif kasar). Saya ingin sekarang memperbaiki grid saya dan menyisipkan solusi saya sebelumnya pada grid baru saya sebelum...
Apakah ada makalah yang bagus dan atau kode yang dipasangkan pasangan solver elemen hingga galerkin terputus dengan solver Riemann? Saya perlu mengeksplorasi masalah elips dan hiperbolik, tetapi sebagian besar metode pemisahan ad hoc. Karena saya memiliki banyak kode FEniCS, saya hanya ingin...
Ketika memecahkan sistem linear jarang menggunakan metode faktorisasi langsung, strategi pemesanan yang digunakan secara signifikan berdampak pada faktor pengisian elemen yang tidak nol dalam faktor. Salah satu strategi pemesanan tersebut adalah diseksi bersarang. Saya bertanya-tanya apakah mungkin...
Saya mencoba mengintegrasikan ∫10t2 n + 2exp( α r0t) dt∫01t2n+2exp(αr0t)dt\int^1_0 t^{2n+2}\exp\left({\frac{\alpha r_0}{t}}\right)dt yang merupakan transformasi sederhana dari ∫∞1x2 nexp( - α r0x ) dx∫1∞x2nexp(−αr0x)dx\int^{\infty}_1 x^{2n}\exp(-\alpha r_0 x)dx menggunakan karena sulit untuk...
Sebagai contoh, saya ingin numerik menghitung -norm dari u = 1L.2L.2L^2 di beberapa domain yang meliputi nol, saya mencoba Gauss quadrature dan gagal, itu agak jauh dari nyataL2-norm pada unit bola menggunakan koordinat bola untuk mengintegrasikan, adalah adakah cara yang baik untuk melakukan ini?...
Sistem matriks yang tidak terbatas muncul misalnya dalam diskritisasi masalah sadel oleh elemen hingga campuran. Matriks sistem kemudian dapat dimasukkan ke dalam formulir (ABBtC)(ABtBC)\begin{pmatrix} A & B^t \\ B & C\end{pmatrix} di mana negatif (semi) -definite, C positif (semi-) pasti dan B...
Saya telah mendengar bahwa beberapa jurnal mendapat peringkat lebih tinggi daripada yang lain. Apakah ini benar? Dan jika demikian, apa kriteria untuk menilai nilai satu jurnal yang ditinjau oleh sejawat? Bagaimana cara mengetahui peringkatnya? Apakah publikasi saya akan kurang "bernilai" jika...
Saya telah mendengar bahwa transformasi fourier cepat dapat digunakan untuk menyelesaikan masalah poisson ketika kondisi batas semuanya satu jenis ... Seri sinus untuk dirichlet, cosinus untuk neumann, dan keduanya untuk periodik. Mengingat domain persegi panjang 2D, anggap dua sisi yang berlawanan...
Saya punya beberapa pertanyaan mengenai hal berikut: Saya mencoba untuk memecahkan persamaan Schrodinger dalam 1D menggunakan diskritisasi nicolson engkol diikuti dengan membalikkan matriks tridiagonal yang dihasilkan. Masalah saya sekarang telah berkembang menjadi masalah dengan kondisi batas...
Saya ingin memahami lebih detail tentang penerapan Metode Multigrid Aljabar (AMG). Saya telah membaca "Tutorial Multigrid", yang cukup bagus dan menjelaskan semua detail interpolasi, operator grid-kasar, dan pemilihan grid kasar untuk AMG. Namun, saya pikir, tidak ada yang seperti bermain-main dan...
Saya memiliki mata, wajah , tepi , dan simpul , dan saya memiliki daftar kontur set level yang telah ditentukan.FFFEEEVVV Algoritma apa yang dapat saya gunakan untuk membuat kontur dengan cara yang paling efisien? Plot kontur ditampilkan di atas. Garis dengan warna yang sama memiliki nilai sama...
Apakah ada format standar untuk jerat elemen hingga yang banyak digunakan dalam industri ini? Terima
Misalkan sistem linear berikut ini diberikan mana adalah Laplacian tertimbang yang dikenal sebagai semi- pasti positif dengan ruang nol satu dimensi yang dibentang oleh 1_n = (1, \ dots, 1) \ in \ mathbb {R } ^ n , dan variasi terjemahan x \ in \ mathbb {R} ^ {n} , yaitu, x + a1_n tidak mengubah...
Biarkan dan deg f > deg g . Saya mencari algoritma asymptotically cepat dan numerik yang stabil untuk komputasi f mod g . Dalam aplikasi yang dimaksudkan, kedua f , g adalah polinomial padat dengan koefisien floating-point presisi ganda. Tapi, untuk saat ini, saya lebih tertarik pada algoritma...