Ilmu Komputasi

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 (Lihat "Resep Numerik" untuk beberapa contoh). Namun, saya bertanya-tanya …

18 algorithms 

Saya mengembangkan kode untuk mensimulasikan aliran fluida dengan zat biologis yang ada dalam aliran. Ini melibatkan persamaan standar Navier-Stokes yang digabungkan ke beberapa model biologis tambahan. Ada banyak parameter / konstanta. Saya memiliki fungsi tertulis untuk menangani perhitungan utama, tetapi masalah yang saya miliki adalah banyaknya konstanta / parameter yang …

17 c++ 

Saya seorang mahasiswa Teknik Mesin yang tertarik pada bidang teknik aerospace di mana, saya diberitahu, Fortran masih umum digunakan. Versi Fortran mana yang harus saya investasikan waktu saya untuk belajar?

17 fortran 

Saya hanya ingin tahu mengapa metode tingkat tinggi (yaitu lebih dari 4) metode Runge-Kutta hampir tidak pernah dibahas / dipekerjakan (setidaknya setahu saya). Saya mengerti ini membutuhkan waktu komputasi yang lebih besar per langkah (misalnya RK14 dengan langkah tertanam urutan ke-12 ), tetapi apakah ada kelemahan lain dari menggunakan metode …

17 ode  runge-kutta 

Saya sedang menulis perpustakaan kecil untuk perhitungan matriks jarang sebagai cara untuk belajar sendiri untuk memanfaatkan pemrograman berorientasi objek. Saya telah bekerja sangat keras untuk memiliki model objek yang bagus, di mana bagian-bagiannya (matriks jarang dan grafik yang menggambarkan struktur konektivitas mereka) sangat longgar digabungkan. Dalam pandangan saya sendiri, kode …

17 linear-algebra  sparse-matrix  oop 

Dengan dua implementasi BLAS yang berbeda, dapatkah kita berharap mereka membuat perhitungan floating point yang sama persis dan mengembalikan hasil yang sama? Atau bisa terjadi, misalnya, bahwa salah satu Toedjoe menghitung produk skalar sebagai dan satu sebagai ( x 1 y 1 + x 2 y 2 ) + ( …

17 floating-point  blas 

Ada dua pendekatan umum untuk mewakili solusi dalam metode galerkin diskontinyu: nodal dan modal. Modal : Solusi yang diwakili oleh jumlah koefisien modal dikalikan dengan satu set polinomial, misalnya di mana φ i biasanya ortogonal polinomial, misalnya Legendre . Salah satu keuntungannya adalah polinomial ortogonal menghasilkan matriks massa diagonal.kamu ( …

17 fluid-dynamics  discontinuous-galerkin 

Saya perlu menulis beberapa data dari perhitungan, yang akan dibaca nanti oleh Paraview (file .vtu atau vtk). Ketika datang ke ukuran file, haruskah saya memilih format ASCII atau format Biner?

17 visualization  paraview  vtk 

Mana yang lebih baik: FORTRAN atau Python? Dan saya kira dalam kedua kasus Anda membutuhkan Gnuplot, apakah saya benar? Saya sedang mengerjakan mesin Windows saat ini. Saya ingin menggunakannya untuk mendapatkan solusi numerik untuk masalah fisika, termasuk simulasi Monte-Carlo, integrasi dan diferensiasi numerik, dinamika molekuler, dll. Saya melihat kursus fisika …

17 python  fortran 

Saya telah menjalankan GROMACS kode dinamika molekul (MD) pada cluster Linux Ubuntu yang terdiri dari node yang berisi 24 CPU Intel Xeon. Titik perhatian saya ternyata agak sensitif terhadap presisi aritmatika floating point, jadi saya harus menjalankan GROMACS dalam presisi ganda daripada presisi tunggal - meskipun biaya komputasi presisi ganda …

17 hpc  floating-point  molecular-dynamics  gpu 

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 miliki?

17 mesh-generation 

Banyak pekerjaan saya sendiri berkisar membuat skala algoritma lebih baik, dan salah satu cara yang disukai untuk menunjukkan penskalaan paralel dan / atau efisiensi paralel adalah untuk memplot kinerja algoritma / kode pada jumlah inti, misalnya di mana sumbu mewakili jumlah core dan sumbu beberapa metrik, misalnya pekerjaan yang dilakukan …

17 parallel-computing  hpc  scaling  plotting 

Saya baru saja menginstal kartu grafis Nvidia GT660 di desktop saya dan, setelah beberapa perjuangan, saya berhasil menghubungkannya dengan R. Saya telah bermain dengan beberapa paket R yang menggunakan GPU, terutama gputools, dan saya membandingkan waktu yang dibutuhkan oleh GPU dan CPU saya untuk melakukan beberapa operasi dasar: matriks pembalik …

17 r  gpu 

Misalkan saya memiliki beberapa fungsi fff dan saya ingin mencari sedemikian rupa sehingga . Saya mungkin menggunakan metode Newton-Raphson. Tapi ini mengharuskan saya tahu fungsi turunan . Ekspresi analitik untuk mungkin tidak tersedia. Sebagai contoh, dapat didefinisikan oleh sepotong kode komputer yang rumit yang berkonsultasi dengan basis data nilai eksperimental.f …

17 reference-request  approximation 

Saya tertarik pada kasus di mana konjugasi gradien bekerja jauh lebih baik daripada metode GMRES. Secara umum, CG adalah pilihan yang lebih disukai dalam banyak kasus SPD (symmetric-positive-definite) karena memerlukan lebih sedikit penyimpanan dan ikatan teoritis pada tingkat konvergensi untuk CG dua kali lipat dari GMRES itu. Apakah ada masalah …

17 linear-solver  conjugate-gradient  gmres