Di mana pun saya melihat, tutorial / dokumen dll PETSc mengatakan bahwa itu berguna untuk aljabar linier dan biasanya menentukan bahwa sistem jarang akan mendapat manfaat. Bagaimana dengan matriks padat? Saya khawatir tentang menyelesaikan untuk padat .A x = bAx=bAx=bSEBUAHAA Saya telah menulis...
Saya membandingkan dua metode berulang untuk membalikkan matriks kuadrat acak. Karena matriksnya acak, setiap test case mengambil jumlah iterasi dan waktu yang berbeda yang berbeda. Pertanyaan saya adalah, di atas rata-rata waktu CPU, adalah nilai rata-rata dari iterasi yang diambil oleh kedua...
Dalam solusi numerik PDE nilai batas awal, sangat umum untuk menggunakan paralelisme dalam ruang . Adalah jauh lebih jarang untuk menggunakan beberapa bentuk paralelisme dalam diskritisasi waktu , dan bahwa paralelisme biasanya jauh lebih terbatas. Saya menyadari semakin banyak kode dan karya yang...
Saya memiliki masalah ketika saya ingin menggunakan perkiraan perbedaan pusat tingkat tinggi: ( - kamui + 2 , j+ 16 ui + 1 , j- 30 usaya , j+ 16 ui - 1 , j- kamui - 2 ,
Adakah orang di sini yang menggunakan komputasi ilmiah presisi ganda dengan GPU generasi baru (misalnya K20) melalui Python? Saya tahu bahwa teknologi ini berkembang pesat, tetapi apa cara terbaik untuk melakukan ini saat ini? GPU berada di luar ruang lingkup untuk perpustakaan ilmiah Python...
Misalkan saya menggunakan generator nomor semu kongruen-acak (PRNG). Diberikan seed , faktor pengali (a), faktor shift (c) dan faktor modulus (m), bagaimana saya bisa menentukan periode PRNG saya? Apakah saya menentukannya dengan eksperimen / algoritma deteksi pola, atau adakah formula langsung...
Algoritma Remez adalah rutin iteratif yang terkenal untuk memperkirakan suatu fungsi oleh polinomial dalam norma minimum. Tapi, seperti yang dikatakan Nick Trefethen [1] tentang itu: Sebagian besar dari [implementasi] ini telah berlangsung bertahun-tahun dan pada kenyataannya, sebagian besar...
Katakanlah saya memiliki sistem linier Ax=bAx=bA x = b , yang konvergen dengan cepat menggunakan metode Krylov yang sesuai (seperti CG atau GMRES) untuk semua bbb . Jika BBB adalah matriks dengan peringkat rendah rrr, akankah metode Krylov yang sama pada sistem (A+B)x=b(A+B)x=b(A + B) x = b juga...
Saya memiliki masalah pemrograman bilangan bulat campuran. Dan saya saat ini menggunakan GLPK sebagai solver saya. Tetapi saya menemukan bahwa GLPK baik untuk masalah Pemrograman Linier, tetapi untuk pemrograman Mixed Integer, membutuhkan waktu yang lebih lama, oleh karena itu tidak memenuhi...
Saya sedang melakukan beberapa simulasi persamaan Langevin, untuk berbagai kekuatan eksternal. Diberitahu bahwa C's rand()from stdlib.hdapat memperkenalkan bias pada hasil saya, saya menggunakan Twister Mersenne. Namun demikian, saya ingin tahu (dan melihat) persis apa jenis kesalahan yang dapat...
Saya sedang menulis makalah yang dapat direproduksi, dan makalah ini memiliki hasil komputasi yang dihasilkan oleh skrip Python (skrip MATLAB serupa menghasilkan hasil yang hampir identik). Saya merasa bahwa makalah ini akan lebih mudah dipahami oleh pembaca jika mereka dapat mencocokkan...
Saya ingin menerapkan ekspresi berikut dalam Python: mana x dan y adalah array numpy dengan ukuran n , dan k adalah array numpy ukuran n × n . Ukuran n mungkin hingga sekitar 10.000, dan fungsinya adalah bagian dari loop dalam yang akan dievaluasi berkali-kali, jadi kecepatan itu penting.xsaya=...
Saya memparalelkan kode untuk memecahkan model keseimbangan populasi 5 dimensi secara numerik. Saat ini saya memiliki kode paralel MPICH2 yang sangat baik di FORTRAN tetapi ketika kami meningkatkan nilai parameter array menjadi terlalu besar untuk dijalankan dalam mode memori terdistribusi. Saya...
Saya memiliki dataset yang mengalir ke jutaan titik data dalam 3D. Untuk perhitungan yang saya lakukan, saya perlu menghitung tetangga (pencarian rentang) untuk setiap titik data dalam radius, mencoba untuk menyesuaikan suatu fungsi, menghitung kesalahan untuk fit, ulangi ini untuk titik data...
Banyak pustaka C ++ yang lebih terkenal dalam ilmu komputasi seperti Eigen , Trilinos , dan deal.II menggunakan objek pustaka header template C ++ standar std::complex<>,, untuk mewakili angka floating-point yang kompleks. Dalam jawaban Jack Poulson untuk pertanyaan tentang konstruktor...
Saat ini saya sedang menulis tesis PhD saya. Saya menghabiskan sebagian besar PhD saya membersihkan dan memperluas kode ilmiah yang ada, menerapkan praktik terbaik rekayasa perangkat lunak yang sebelumnya tidak digunakan, dan ingin menulis tentang ini dalam tesis saya. Daripada hanya mengatakan...
Saya mencoba mencari beberapa sumber untuk membantu menjelaskan bagaimana memilih kondisi batas saat menggunakan metode beda hingga untuk menyelesaikan PDE. Buku-buku dan catatan yang saat ini saya memiliki akses untuk semua mengatakan hal serupa: Aturan umum yang mengatur stabilitas di hadapan...
Sebanyak saya mencoba menemukan penjelasan singkat di internet, saya sepertinya tidak bisa memahami konsep perbedaan terbatas mimesis, atau bagaimana itu bahkan berhubungan dengan perbedaan hingga standar. Akan sangat membantu untuk melihat beberapa contoh sederhana bagaimana mereka...
Saya telah menghabiskan beberapa bulan terakhir pada pengkodean program Fortran untuk memecahkan sistem PDE tertentu (menggambarkan aliran / pembakaran fluida). Saya mencoba menggunakan Fortran standar terbaru dan kapabilitas baru OOP yang dimiliki Fortran. Saya bekerja sendiri dan tidak memiliki...
Saya baru saja mulai bermain-main dengan FEniCS. Saya memecahkan Poisson dengan elemen urutan ke-3 dan ingin memvisualisasikan hasilnya. Namun, ketika saya menggunakan plot (u), visualisasi hanyalah interpolasi linier dari hasilnya. Saya mendapatkan hal yang sama ketika saya output ke VTK. Dalam...