Ilmu Komputasi

Sumber istilah, seperti yang karena batimetri dalam persamaan air dangkal, perlu diintegrasikan dengan cara khusus untuk menjaga keadaan fisik stabil. Apakah ada cara umum untuk membangun metode yang seimbang, atau apakah itu memerlukan teknik khusus untuk setiap persamaan?

13 pde  hyperbolic-pde 

Apa yang akan Anda katakan akan menjadi keuntungan / kerugian dari dua pendekatan untuk pengkodean perpustakaan umum (volume terbatas, fem, dg) untuk Computational Continuum Mechanics? Ini adalah bagaimana saya melihat sesuatu sekarang, jadi tolong berikan pengalaman Anda sendiri dan jangan nyalakan saya untuk saya :): 1) C ++: pemrograman generik, …

13 libraries  programming-paradigms 

Jika kita menggunakan Metode Garis untuk diskritisasi (waktu dan ruang terpisah diskritisasi) PDE hiperbolik yang kita peroleh setelah diskritisasi spasial dengan metode numerik favorit kita (fx. Metode Volume Hingga) dalam praktiknya pemecah ODE mana yang kita gunakan untuk diskretisasi temporal (TVD / SSP / dll)? Beberapa informasi tambahan ditambahkan: Masalah …

13 hyperbolic-pde 

Saya tahu bahwa orang sering mengganti matriks massa yang konsisten dengan matriks diagonal yang disatukan. Di masa lalu, saya juga menerapkan kode di mana vektor beban dirakit dalam mode yang disatukan daripada mode yang konsisten dengan FEM. Tapi saya tidak pernah melihat mengapa kita diizinkan melakukan ini sejak awal. Apa …

13 finite-element 

Saya mencari sumber daya pada pola desain dan prinsip - prinsip untuk perangkat lunak numerik, berpotensi tetapi tidak harus dengan fokus pada pendekatan berorientasi objek untuk kode numerik. Saya mengetahui teksnya, Menulis Perangkat Lunak Ilmiah: Panduan untuk Gaya yang Baik , tetapi ini masih banyak berhubungan dengan hal-hal seperti masalah …

13 design-patterns 

Dalam metode seperti gmres atau bicgstab, mungkin menarik untuk menggunakan metode krylov lain sebagai prasyarat. Bagaimanapun mereka mudah diimplementasikan dalam cara bebas-matriks dan dalam lingkungan paralel. Misalnya, satu dapat menggunakan beberapa (misalkan ~ 5) iterasi bigcstab tanpa syarat sebagai precontioner untuk gmres, atau kombinasi metode krylov lainnya. Saya tidak menemukan …

13 linear-algebra  linear-solver  preconditioning  krylov-method  gmres 

Saya tidak terbiasa dengan format file yang digunakan dalam HDF5, tetapi saya bertanya-tanya apakah file HDF5 cocok untuk kontrol revisi dengan git (atau misalnya Mercurial atau Subversion)? Saya kira apa yang saya maksud adalah: apakah file HDF5 cocok untuk perbedaan berbasis garis atau git harus memperlakukan HDF5 sebagai satu biner …

13 data-management  data-storage 

Persamaan Poisson dengan semua kondisi batas Neumann memiliki ruang nol dimensi konstan tunggal. Ketika menyelesaikan melalui metode Krylov, ruang nol dapat dihapus baik dengan mengurangi rata-rata solusi setiap iterasi atau dengan menyematkan nilai satu titik. Menyematkan satu titik memiliki manfaat kesederhanaan, dan juga menghindari pengurangan global tambahan per proyeksi. Namun, …

13 krylov-method 

Dalam suatu sistem di mana energi secara teoritis harus dilestarikan, simulasi yang paling akurat akan menghemat energi (serta memberikan posisi, kecepatan, dan lain-lain yang akurat). RK4 lebih akurat daripada leapfrog, namun leapfrog menghemat energi dan RK4 tidak. Kenapa ini?

13 computational-physics  time-integration  molecular-dynamics  runge-kutta 

Saya punya pertanyaan tentang quadric fit untuk satu set poin dan normals yang sesuai (atau setara, tangen). Pemasangan permukaan kuadrat ke titik data dieksplorasi dengan baik. Beberapa karya adalah sebagai berikut: Pemasangan Langsung Berbatas Type-Quadric , James Andrews, Carlo H. Sequin Desain & Aplikasi Berbantuan Komputer, 10 (a), 2013, bbb-ccc …

13 computational-geometry  regression  geometry  curve-fitting 

Bentuk PDE yang kuat mensyaratkan bahwa solusi yang tidak diketahui milik . Tetapi bentuk lemah hanya mensyaratkan bahwa solusi yang tidak diketahui milik H 1 .H2H2H^2H1H1H^1 Bagaimana Anda mendamaikan ini?

13 pde  weak-solution 

Saya mencoba memahami persamaan advection dengan koefisien kecepatan variabel sedikit lebih baik. Khususnya saya tidak mengerti bagaimana persamaan bisa konservatif. The persamaan adveksi , ∂u∂t+∂∂x(vu)=0∂u∂t+∂∂x(vu)=0 \frac{\partial u}{\partial t} + \frac{\partial}{\partial x}(\boldsymbol{v}u) = 0 Mari kita menafsirkan sebagai konsentrasi beberapa spesies fisik ( c m - 3 ) atau kuantitas fisik …

13 advection  conservation 

Saya telah membaca bahwa jika, misalnya, saya memiliki forloop ganda yang berjalan di atas indeks matriks, maka menempatkan indeks berjalan kolom di loop luar lebih efisien. Sebagai contoh: a=zeros(1000); for j=1:1000 for i=1:1000 a(i,j)=1; end end Apa cara paling efisien untuk mengkodekannya jika saya memiliki tiga forloop atau lebih ? …

12 matlab  efficiency 

Saya mengerjakan masalah terbalik untuk Ph.D. penelitian, yang demi kesederhanaan kita akan katakan adalah menentukan dalamββ\beta L(β)u≡−∇⋅(k0eβ∇u)=fL(β)u≡−∇⋅(k0eβ∇u)=fL(\beta)u \equiv -\nabla\cdot(k_0e^\beta\nabla u) = f dari beberapa pengamatan ; k 0 adalah konstanta dan f diketahui. Ini biasanya dirumuskan sebagai masalah optimasi untuk ekstremisasiuouou^ok0k0k_0fff J[u,λ;β]=12∫Ω(u(x)−uo(x))2dx+∫Ωλ(L(β)u−f)dxJ[u,λ;β]=12∫Ω(u(x)−uo(x))2dx+∫Ωλ(L(β)u−f)dxJ[u, \lambda; \beta] = \frac{1}{2}\int_\Omega(u(x) - u^o(x))^2dx + \int_\Omega\lambda(L(\beta)u …

12 inverse-problem 

Hanya dua bagian pertama dari pertanyaan panjang ini yang penting. Yang lain hanya untuk ilustrasi. Latar Belakang Kuadratur maju seperti komposit tingkat tinggi Newton-Cotes, Gauß-Legendre, dan Romberg tampaknya terutama ditujukan untuk kasus-kasus di mana seseorang dapat mengambil sampel fungsi secara halus tetapi tidak berintegrasi secara analitis. Namun, untuk fungsi dengan …

12 reference-request  quadrature