Kebanyakan metode numerik untuk quadrature memperlakukan integand sebagai fungsi kotak hitam. Bagaimana jika kita memiliki informasi lebih lanjut? Secara khusus, apa manfaatnya, jika ada, yang dapat kita peroleh dari mengetahui beberapa turunan pertama dari integrand? Informasi apa yang mungkin...
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...
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...
Saya sedang melihat menggunakan OpenFOAM untuk menyelesaikan arus internal dasar dalam CFD. Apa cara terbaik untuk memulai, dan bisakah orang tolong tunjukkan saya ke referensi online yang bagus untuk membahas pertanyaan yang mungkin saya miliki setelah saya menyelam? Saya telah mendengar bahwa itu...
Misalkan saya memiliki persamaan matriks yang didefinisikan secara rekursif sebagai A[n] = inverse([1 - b[n]A[n+1]]) * a[n] Kemudian persamaan untuk A [1] terlihat mirip dengan fraksi lanjutan, di mana ada beberapa metode yang sangat efisien yang menghindari perhitungan ulang yang membosankan...
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...
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...
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
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(:)...
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...
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...
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
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...
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...
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...
Saya mencari perangkat lunak generasi mesh itu gratis dan open source, menyediakan antarmuka scripting yang waras untuk spesifikasi domain, bekerja untuk geometri yang kompleks, dapat menghasilkan jerat 2D dan 3D, Opsi apa yang saya
Ketika perhitungan terbatas bandwidth memori dilakukan dalam lingkungan memori bersama (mis. Berulir melalui OpenMP, Pthreads, atau TBB), ada dilema tentang bagaimana memastikan bahwa memori didistribusikan dengan benar di seluruh memori fisik , sehingga masing-masing thread kebanyakan mengakses...
Saya memiliki fungsi 2-D non-cembung terbatas yang ingin saya temukan minimum. Fungsinya cukup lancar. Mengevaluasi itu mahal. Kesalahan yang dapat diterima adalah sekitar 3% dari domain fungsi di setiap sumbu. Saya mencoba menjalankan implementasi algoritma DIRECT di perpustakaan NLOPT, tetapi...
Pertimbangkan IVP untuk sistem ODE , . Paling umum masalah ini dianggap kaku ketika matriks Jacobi memiliki kedua nilai eigen dengan bagian nyata negatif yang sangat besar dan nilai eigen dengan bagian nyata negatif yang sangat kecil (saya menganggap hanya bagian stabil kasus).y ( x 0 ) = y 0 ∂...
Dalam literatur metode numerik, banyak makalah penelitian terdiri dari deskripsi variasi algoritme baru, diikuti oleh beberapa masalah tes membandingkan metode baru dengan satu atau dua metode yang ada. Ini membuatnya sulit untuk ditentukan Bagaimana metode baru akan dilakukan pada masalah lain...