Pertanyaan yang diberi tag petsc

Saya memiliki beberapa masalah optimisasi global non-cembung yang menantang untuk dipecahkan. Saat ini saya menggunakan MATLAB's Optimization Toolbox (khusus, fmincon()dengan algoritma = 'sqp'), yang cukup efektif . Namun, sebagian besar kode saya menggunakan Python, dan saya ingin melakukan...

Apa yang bisa salah ketika menggunakan metode Krylov yang dikondisikan sebelumnya dari KSP ( paket solver linier PETSc ) untuk menyelesaikan sistem linier yang jarang seperti yang diperoleh dengan mengkritisi dan membuat linierisasi persamaan diferensial parsial? Langkah apa yang bisa saya ambil...

Sejauh yang saya tahu, dua kerangka kerja perangkat lunak ilmu komputasi generik besar Departemen Energi AS adalah PETSc dan Trilinos . Sepintas mereka tampak serupa, di luar perbedaan bahasa (C versus C ++). Apa perbedaan utama antara kedua kerangka kerja, dan faktor-faktor apa yang harus...

Saya menggunakan paket PETSc non -linier SNES untuk menyelesaikan sistem persamaan nonlinier yang diperoleh dengan mendiskritisasi persamaan diferensial parsial. Bagaimana saya bisa menentukan mengapa pemecah tidak konvergen dan apa yang bisa saya lakukan untuk berhasil menyelesaikan persamaan...

Saya berencana menggunakan multigrid untuk menghitung beberapa nilai eigen dan vektor, dan saya perhatikan PETSc memiliki dukungan tingkat tinggi untuk multigrid. Dokumentasi PETSc mengatakan bahwa bagian PETSc ini tidak boleh digunakan, karena akan segera diganti. Perpustakaan lain mana yang...

Saya ingin tahu apakah ada cara cepat untuk menghitung jarak Euclidean dari dua vektor dalam Oktaf. Sepertinya tidak ada fungsi khusus untuk itu, jadi haruskah saya menggunakan rumus saja

Saya memiliki PETSc Matdan ingin memperkirakan jumlah

Saya bermain-main dengan PETSc dan memperhatikan bahwa ketika saya menjalankan program saya dengan lebih dari satu proses melalui MPI tampaknya berjalan lebih lambat ! Bagaimana saya dapat memeriksa untuk melihat apa yang sedang

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...

Biasanya, saya menulis kode seri, dan ketika saya melakukannya, saya menulis unit test dengan beberapa kerangka pengujian gaya xUnit (MATLAB xUnit, PyUnit / nose, atau kerangka pengujian C ++ Google). Berdasarkan pencarian Google sepintas, saya belum melihat banyak tentang bagaimana praktisi...

Pertanyaan: Metode apa yang tersedia untuk secara akurat dan efisien menghitung struktur sparsity dari matriks elemen hingga? Info: Saya sedang mengerjakan pemecah Persamaan Tekanan Poisson, menggunakan metode Galerkin dengan basis Lagrange kuadratik, ditulis dalam C, dan menggunakan PETSc untuk...

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...

Saya perlu memecahkan sistem linear yang jarang sama (300x300 hingga 1000x1000) dengan banyak sisi kanan (300 hingga 1000). Selain masalah pertama ini, saya juga ingin menyelesaikan sistem yang berbeda, tetapi dengan elemen bukan nol yang sama (hanya nilai yang berbeda), yaitu banyak sistem jarang...

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...

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

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 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 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...

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,...

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...