Sistem saya adalah masalah FE simetris dengan pengganda lagrange (misalnya aliran Stokes yang tidak dapat dikompresi): ( ABBTC)(SEBUAHBTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} di mana adalah kasus khas (saya bahkan telah memastikan bahwa persamaan diberi nomor sehingga pengali Lagrange...
Berikut ini adalah penggunaan khas pointer fungsi di C. Saya ingin melakukan sesuatu yang serupa di Fortran. Saya punya beberapa ide, tetapi saya ingin tahu apakah ada cara kanonik untuk melakukannya. Pointer fungsi dan konteks yang dilewatkan oleh pengguna disimpan, lalu dipanggil...
Komputasi paralel sering dimodelkan menggunakan laju perhitungan lokal deterministik, overhead latensi, dan bandwidth jaringan. Pada kenyataannya, ini adalah variabel spasial dan non-deterministik. Studi seperti Skinner dan Kramer (2005) mengamati distribusi multi-modal, tetapi analisis kinerja...
Seberapa lambatkah petsc4py vs c / c ++ / fortran? Saya menyadari itu akan sangat tergantung pada kode yang dieksekusi, tetapi bagaimana dengan sesuatu yang sederhana seperti produk
12×1212×1212 \times 12QQQJdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Saat ini saya melakukan ini dengan perpustakaan armadillo tetapi ternyata terlalu lambat. Masalahnya adalah bahwa saya perlu melakukan ini untuk satu triliun matriks dan ternyata menghitung...
Apakah kompilasi PETSc dengan perpustakaan BLAS / LAPACK eksternal secara signifikan mempengaruhi kinerja pada matriks jarang, atau apakah hanya menggunakan perpustakaan tersebut untuk matematika matriks
Saya ingin menggunakan metode urutan ke-8 Runge-Kutta (89) dalam aplikasi mekanika / astrodinamika surgawi, yang ditulis dalam C ++, menggunakan mesin Windows. Karena itu saya bertanya-tanya apakah ada yang tahu perpustakaan / implementasi yang baik yang didokumentasikan dan bebas untuk digunakan?...
Manakah yang dihitung lebih cepat, atau atau ? , dan adalah real positif dengan .log a c b √Sebuahbaba^bcatatanSebuahclogac\log_a c abcb>1c√bcb\sqrt[b]{c}Sebuahaabbbcccb > 1b>1b>1 Jenis algoritma apa yang akan Anda gunakan dalam perbandingan? Apa kompleksitasnya? Misalnya, ketika atau c...
Ada banyak metode numerik yang terkenal untuk menyelesaikan persamaan dari tipe misalnya metode pembagian dua bagian, metode Newton, dll.f(x)=0,x∈Rn,f(x)=0,x∈Rn, f(x) = 0, \quad x \in \mathbb{R}^n, Dalam aplikasi saya dihitung dengan metode stokastik (hasilnya rata-rata).f(x)f(x)f(x) Apakah ada...
Saya mencari referensi yang kredibel yang menyatakan berapa banyak sumber daya superkomputer menghabiskan pada koordinasi dibandingkan melakukan pekerjaan yang berhubungan dengan tugas yang sebenarnya. Sumber daya bisa menjadi kekuatan pemrosesan yang tersedia tetapi bahkan Watts tampak seperti...
Saya baru mengenal paket PETSc. Saya memiliki ~ 4000x4000 matriks A dalam format matriks-pasar dan saya ingin mendapatkan PETSc untuk menyelesaikannya dengan menggunakan beberapa prosesor. Saya tahu bagaimana menyelesaikan sistem pada satu prosesor, tetapi saya tidak tahu bagaimana...
Saya memiliki model yang cukup besar (~ 5000 baris) ditulis dalam C. Ini adalah program serial, tanpa generasi nomor acak di mana pun. Itu menggunakan perpustakaan FFTW untuk fungsi menggunakan FFT - Saya tidak tahu detail implementasi FFTW, tapi saya menganggap fungsi di dalamnya juga...
Dalam beberapa tahun terakhir, beberapa proyek perpustakaan / perangkat lunak telah muncul yang menawarkan beberapa bentuk atau paralelisme memori bersama untuk keperluan umum yang didorong oleh data. Gagasan utamanya adalah bahwa alih-alih menulis kode yang diulir secara eksplisit, pemrogram...
Saya sangat baru dalam pemrograman GPGPU jadi tolong maafkan saya jika pertanyaannya tidak tepat. Dari apa yang saya pahami pemrograman GPU adalah bagian yang sangat rumit dari pekerjaan teknik jika dibandingkan dengan pemrograman CPU biasa. Kita harus sangat berhati-hati tentang masalah...
Pada dasarnya, FEM tampaknya menjadi masalah yang cukup banyak "diselesaikan". Ada banyak kerangka kerja yang kuat yang ada, seperti Trilinos, PETSc, FEniCS, Libmesh atau MOOSE. Satu hal yang mereka miliki bersama: Mereka sangat "berat". Pertama, instalasi biasanya sangat menyakitkan. Kedua,...
Saya terjun ke dunia analisis elemen hingga yang menakjubkan dan ingin menyelesaikan masalah termo-mekanis yang besar (hanya termal mekanis, tanpa umpan balik).→→\rightarrow Untuk masalah mekanis, saya sudah memahami jawaban Geoff , bahwa saya harus menggunakan solver iteratif karena ukuran mesh...
Saya perlu menyelesaikan Ax = b, tetapi saya menyadari bahwa meskipun jarang, menyimpan koefisien matriks dari masalah saya akan memakan terlalu banyak memori. Jadi sekarang saya sedang mempertimbangkan menggunakan metode bebas-matriks, karena koefisien yang sama muncul banyak waktu dalam matriks,...
Sebagian besar pemrograman saya adalah kode penelitian satu kali dalam C untuk saya gunakan sendiri. Saya belum pernah mendistribusikan kode apa pun selain dari kolaborator dekat. Saya telah mengembangkan suatu algoritma yang saya terbitkan dalam jurnal ilmiah. Saya ingin memberikan kode sumber dan...
Seperti yang dicontohkan oleh jawaban Jed Brown untuk Biaya pencarian versus perhitungan , menggunakan operasi floating point yang di-vectorized vs non-vectorized menghasilkan kode yang jauh lebih cepat. Banyak kompiler modern mengklaim bahwa mereka dapat melakukan vektorisasi otomatis. Bagaimana...