Baru-baru ini, saya mengalami masalah aneh dengan FORTRAN95. Saya menginisialisasi variabel X dan Y sebagai berikut: X=1.0 Y=0.1 Kemudian saya menambahkannya bersama dan mencetak hasilnya: 1.10000000149012 Setelah memeriksa variabel, sepertinya 0,1 tidak terwakili dalam presisi ganda dengan...
Setelah beberapa langkah pembaruan Bayesian, saya dibiarkan dengan distribusi posterior bentuk campuran distribusi normal, Pr ( θ | data ) = ¬i = 1kwsayaN( μsaya, σ2) .Pr(θ|data)=∑i=1kwiN(μi,σ2).\Pr(\theta| \text{data} ) = \sum_{i=1}^k w_i N(\mu_i, \sigma^2). Yaitu, parameter θθ\theta diambil dari...
Saya memiliki matriks persegi simetris nyata padat. Dimensinya sekitar 1000x1000. Saya perlu menghitung komponen utama pertama dan bertanya-tanya apa algoritma terbaik untuk melakukan ini. Tampaknya MATLAB menggunakan algoritma Arnoldi / Lanczos (untuk eigs). Tetapi dari membaca tentang mereka...
Saya memecahkan banyak PDE secara numerik, tetapi matematika terapan bukan bidang saya. Saya belum mengambil jurnal matematika terapan mana yang harus saya baca untuk mengikuti perkembangan terkini di lapangan. Apa jurnal yang baik untuk dibaca untuk mengikuti perkembangan terakhir dalam...
Saya memiliki ODE: u ( 0 ) = - 1u′=−1000u+sin(t)u′=−1000u+sin(t)u'=-1000u+sin(t) u(0)=−11000001u(0)=−11000001u(0)=-\frac{1}{1000001} Saya tahu bahwa ODE khusus ini kaku, analitis. Saya juga tahu bahwa jika kita menggunakan metode loncatan waktu (maju) yang eksplisit (Euler, Runge-Kutta, Adams,...
Saya memecahkan masalah multiskala menggunakan Metode Multiskala Heterogen (HMM) . Pada dasarnya, prosedur khusus saya menggunakan proses berulang berikut: Memecahkan banyak sistem matriks lokal. Hitung nilai yang menarik dari solusi sistem lokal. Merakit sistem matriks global dari "nilai-nilai...
Dalam metode medan konsisten-sendiri Hartree-Fock dalam menyelesaikan persamaan Schroedinger elektronik yang tidak tergantung waktu, kami berupaya meminimalkan energi keadaan dasar, , dari sistem elektron dalam bidang eksternal berkenaan dengan pilihan putaran orbital, . { χ i
Baru-baru ini saya mengajukan pertanyaan di sepanjang baris yang sama untuk matriks miring-Hermitian. Terinspirasi oleh keberhasilan pertanyaan itu, dan setelah membenturkan kepalaku ke dinding selama beberapa jam, aku melihat matriks eksponensial dari matriks asimetris yang nyata. Rute untuk...
Saya telah mengembangkan solusi kerja dari Metode Elemen Hingga untuk menyelesaikan masalah perpindahan panas menggunakan GPU dan OpenCL menggunakan metode Conjugate Gradient. Kerugian utama dari metode ini adalah tingginya permintaan akan memori. Terlebih lagi, dalam kasus kartu grafis, memori...
Saya mulai belajar OpenFOAM dari tutorial Cavity yang disediakan di situs web . Saat bereksperimen dengan angka Reynolds yang berbeda, di bagian "2.1.8.2 Menjalankan kode", tutorial mengatakan untuk menjalankan kembali solver karena "Masuk akal untuk menambah waktu solusi". Tetapi ketika saya...
Bagaimana pemrograman geometris (umum) berbeda dari pemrograman cembung umum? Program geometrik dapat ditransformasikan menjadi program cembung, dan biasanya diselesaikan dengan metode titik interior. Tetapi apa keuntungan dari merumuskan secara langsung masalah sebagai program cembung dan...
Saya ingin tahu yang mana dari pemecah linear klasik (mis. Gauss-Seidel, Jacobi, SOR) dijamin akan menyatu untuk masalah mana A adalah positif setengah pasti dan tentu saja b ∈ i m ( A )A x = bSEBUAHx=bAx=bSEBUAHSEBUAHAb ∈ i m ( A )b∈sayam(SEBUAH)b \in im(A) (Pemberitahuan adalah semi pasti dan...
Saya memiliki masalah yang mirip dalam perumusan dengan posting ini, dengan beberapa perbedaan penting: Apa metode sederhana yang ada untuk pengambilan sampel fungsi 2D secara adaptif? Seperti di pos itu: Saya memiliki dan evaluasi fungsi ini agak mahal untuk dihitungf( x ,
Sebagai seseorang yang harus mengajar kursus dalam ilmu komputasi, saya dihadapkan dengan pertanyaan kuno: bagaimana saya mengevaluasi kemampuan siswa untuk mempelajari subjek yang tergantung pada aplikasi yang sulit untuk diuji dengan metode pengujian "standar" ( ujian tertulis atau lisan)? Bagian...
Saya tahu bahwa kita dapat menggunakan teknik analisis matematis untuk membuktikan apakah IVP atau BVP memiliki solusi, unik, dan bergantung terus menerus pada batas / nilai awal. Untuk beberapa PDE, terutama PDE non-linear, sangat sulit, jika bukan tidak mungkin, untuk membuktikan posisi yang...
seperti judulnya saya mencoba untuk menghitung integral dari fungsi yang didukung secara kompak (polinomial quintic Wendland) pada segitiga. Perhatikan, bahwa pusat fungsi berada di suatu tempat dalam ruang 3-D. Saya mengintegrasikan fungsi ini pada sembarang segitiga, tetapi kecil ( ). Saat ini...
Apa cara yang disarankan untuk melakukan nonlinier kuadrat, min , dengan kotak kendala l o j < = p j < = h i j ? Menurut saya (orang-orang bodoh masuk) bahwa seseorang dapat membuat batasan kotak kuadratik, dan meminimalkan ∑ i e r r i ( p ) 2 + C ∗ ∑ j t u b ( p j , l...
Misalkan saya tahu bahwa urutan nomor acak dihasilkan oleh generator congruential linier. Itu adalah, xn + 1= ( a Xn+ c ) mod mxn+1=(aXn+c)modmx_{n+1}=(aX_n+c) \bmod m Jika saya diberi seluruh periode (atau paling tidak sebuah urutan besar yang berdekatan), bagaimana saya bisa merekonstruksi...
Misalkan ΩΩ\Omega menjadi domain Lipschitz yang dibatasi secara cembung di R2R2\mathbb R^2 , misalkan f∈L2(Ω)f∈L2(Ω)f \in L^2(\Omega) . Maka solusi dari masalah Dirichlet Δu=fΔu=f\Delta u = f dalam ΩΩ\Omega , traceu=0traceu=0\operatorname{trace} u = 0 pada ∂Ω∂Ω\partial\Omega memiliki solusi unik...
Saya bertanya-tanya bagaimana menemukan nilai eigen dari beberapa matriks jarang dalam interval yang diberikan [a, b] dengan metode iteratif. Menurut pemahaman pribadi saya, lebih jelas menggunakan metode ruang bagian Krylov untuk menemukan nilai eigen ekstrem daripada nilai