Ilmu Komputasi

Saya sering menemukan diri saya menulis kode yang sangat mirip untuk versi satu, dua, dan tiga dimensi dari operasi / algoritma yang diberikan. Mempertahankan semua versi ini bisa menjadi membosankan. Pembuatan kode sederhana berfungsi cukup baik, tetapi tampaknya harus ada cara yang lebih baik. Apakah ada cara yang relatif sederhana …

19 algorithms  software  fortran 

Penelitian yang dapat direproduksi dalam perhitungan bertujuan untuk membuat kode yang diperlukan untuk menghasilkan hasil dalam makalah komputasi yang tersedia bagi peneliti lain sehingga mereka dapat menjalankan kode ini untuk mereproduksi hasil dalam makalah itu. Saya ingin membuat semua penelitian saya dapat direproduksi, tetapi saya menghadapi sedikit hambatan: beberapa makalah …

19 software  publications  documentation  reproducibility 

Diberikan keadaan yang diinginkan y0y0y_0 dan parameter regularisasi , pertimbangkan masalah menemukan keadaan dan kontrol untuk meminimalkan fungsional tunduk pada kendala \ begin {persamaan} Ay = u. \ end {persamaan} di mana untuk kesederhanaan kita dapat memikirkan y, y_0, u \ in \ mathbb R ^ n dan A \ …

19 optimization  reference-request  time-integration 

adalahmatriks jarang simetris positif pasti (SPD). adalah matriks diagonal yang jarang. besar (> 10.000) dan jumlah nonzeros dibiasanya 100 ~ 1000.SEBUAHAAn × nn×nn \times nn n GGGGnnnnnnGGG SEBUAHAA telah difaktorkan dalam bentuk Cholesky sebagai .L D LTLDLTLDL^T Bagaimana cara memperbarui dan efisien ketika menjadi ?D A A + GLLLDDDSEBUAHAAA + …

19 linear-algebra 

Untuk solusi sistem linear besar menggunakan metode iteratif, sering kali menarik untuk memperkenalkan prakondisi, misalnya menyelesaikan , di mana sini digunakan untuk prakondisi kiri sistem. . Biasanya, kita harus memiliki M ^ {- 1} \ approx A ^ {- 1} dan memberikan dasar untuk (lebih banyak) solusi efisien atau pengurangan …

19 linear-algebra  iterative-method  preconditioning 

Saat membaca beberapa literatur tentang pemecah PDE saya menemukan istilah pseudo-melangkah hari ini. Tampaknya menjadi istilah umum, namun saya gagal menemukan definisi yang baik atau artikel pengantar untuk itu. Oleh karena itu: Apa itu pseudo-stepping time, dan bagaimana biasanya digunakan?

18 pde  time-integration 

Pada sebagian besar (semua?) Implementasi Metode Multipole Cepat (FMM), oktri digunakan untuk menguraikan domain yang relevan. Secara teoritis, oktri memberikan ikatan volumetrik sederhana, yang berguna untuk membuktikan runtime O (n) FMM. Di luar dasar pemikiran teoretis ini, apakah ada manfaat menggunakan Octree di atas struktur data pohon atau trie lainnya? …

18 algorithms 

Saya mencoba menyelesaikan beberapa masalah optimasi nonlinear pada GPU (CUDA). Fungsi obyektif adalah fungsi nonlinier yang halus, dan gradiennya relatif murah untuk dihitung secara analitis, jadi saya tidak perlu repot dengan perkiraan numerik. Saya ingin menyelesaikan masalah ini dengan sebagian besar ops matematika fp32 (karena berbagai alasan), jadi metode optimasi …

18 optimization  cuda 

Saya mencoba menerapkan fungsi berikut dalam floating point presisi ganda dengan kesalahan relatif rendah : logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log⁡(exp⁡(x)+exp⁡(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Ini digunakan secara luas dalam aplikasi statistik untuk menambahkan probabilitas atau kepadatan probabilitas yang diwakili dalam ruang log. Tentu saja, baik atau dapat dengan mudah overflow atau underflow, yang akan …

18 floating-point  stability  numerics 

Saya ingin menggunakan f2pyFortran modern. Secara khusus saya mencoba untuk mendapatkan contoh dasar berikut ini untuk bekerja. Ini adalah contoh berguna terkecil yang bisa saya hasilkan. ! alloc_test.f90 subroutine f(x, z) implicit none ! Argument Declarations ! real*8, intent(in) :: x(:) real*8, intent(out) :: z(:) ! Variable Declarations ! real*8, …

18 python  fortran 

Apa algoritma (skalabilitas dan efisiensi) terbaik untuk menghasilkan jerat quad tidak terstruktur dalam 2D? Di mana saya dapat menemukan generator quad mesh terstruktur yang baik? (sumber terbuka lebih disukai)

18 mesh-generation  computational-geometry 

Saya ingin belajar bagaimana metode wavelet dapat diterapkan pada PDE, tetapi sayangnya saya tidak tahu sumber yang bagus untuk mempelajari tentang topik ini. Tampaknya banyak pengantar wavelet fokus pada teori interpolasi, misalnya, mengumpulkan sinyal dengan superposisi beberapa wavelet yang disukai. Aplikasi untuk PDE kadang-kadang disebutkan, tanpa masuk lebih dalam ke …

18 pde  wavelet 

Ingin memperbaiki pos ini? Berikan jawaban terinci untuk pertanyaan ini, termasuk kutipan dan penjelasan mengapa jawaban Anda benar. Jawaban tanpa detail yang cukup dapat diedit atau dihapus. Adaptive mesh refinement (AMR) adalah teknik umum untuk menangani masalah berbagai skala spasial dalam solusi numerik PDE. Apa perpustakaan tujuan umum yang ada …

18 pde  libraries  parallel-computing  adaptive-mesh-refinement 

Kuantitas fisik seperti tekanan, kepadatan, energi, suhu, dan konsentrasi harus selalu positif, tetapi metode numerik terkadang menghitung nilai negatif selama proses solusi. Ini tidak oke karena persamaan akan menghitung nilai yang kompleks atau tidak terbatas (biasanya menabrak kode). Metode numerik mana yang dapat digunakan untuk menjamin bahwa jumlah ini tetap …

18 pde  fluid-dynamics  hyperbolic-pde 

Saya ingin mengembangkan perangkat lunak komputasi paralel ilmiah dari awal. Saya ingin beberapa pemikiran tentang bahasa mana untuk memulai. Program ini melibatkan membaca / menulis data ke file txt dan melakukan perhitungan berat secara paralel, dengan banyak faktorisasi LU dan penggunaan pemecah linier yang jarang. Solusi kandidat yang saya pikirkan …

18 python  hpc  c++  fortran  simulation