Pertanyaan yang diberi tag «sparse-matrix»

Saat ini saya sedang bekerja untuk memecahkan sistem simetris yang sangat besar (tapi tidak pasti positif), yang dihasilkan oleh beberapa algoritma tertentu. Matriks ini memiliki sparsity blok bagus yang dapat digunakan untuk pemecahan paralel. Tetapi saya tidak dapat memutuskan apakah saya harus menggunakan pendekatan langsung (seperti Multi-frontal) atau iteratif (GMRES …

10 parallel-computing  linear-solver  sparse-matrix  linear-system  gmres 

Apakah kompilasi PETSc dengan perpustakaan BLAS / LAPACK eksternal secara signifikan mempengaruhi kinerja pada matriks jarang, atau apakah hanya menggunakan perpustakaan tersebut untuk matematika matriks padat?

10 petsc  sparse-matrix  lapack  blas 

Apakah perkalian matriks (baik Mat * Mat, dan Mat * Vec) skala dengan jumlah non-nol, atau dengan ukuran matriks? Atau kombinasi keduanya. Bagaimana dengan bentuk. Sebagai contoh, saya memiliki 100 x 100 matriks dengan 100 nilai di dalamnya, atau 1000 x 1000 matriks dengan 100 nilai di dalamnya. Ketika mengkuadratkan …

10 linear-algebra  performance  sparse-matrix 

Saya bertanya-tanya bagaimana kondisi batas Dirichlet dalam matriks elemen hingga global jarang benar-benar diterapkan secara efisien. Misalnya, katakanlah matriks elemen hingga global kami adalah: K= ⎡⎣⎢⎢⎢⎢⎢⎢520- 102410001632- 1037000203⎤⎦⎥⎥⎥⎥⎥⎥dan vektor sisi kananb = ⎡⎣⎢⎢⎢⎢⎢⎢b 1b 2b 3b 4b 5⎤⎦⎥⎥⎥⎥⎥⎥K=[520-102410001632-1037000203]dan vektor sisi kananb=[b1b2b3b4b5]K = \begin{bmatrix} 5 & 2 & 0 & -1 …

9 finite-element  sparse-matrix  boundary-conditions 

Banyak pustaka FEM adaptif menggunakan struktur data mesh yang lebih canggih untuk menangani penambahan / penghapusan node, edge, triangles, tetrahedra, dll. Misalnya, p4est library menggunakan struktur data octree untuk perbaikan mesh adaptif; Anda tidak akan sering menemukan oktri yang digunakan untuk perhitungan pada jaring statis. Apa perubahan pada sisi aljabar …

9 linear-algebra  finite-element  sparse-matrix 

Saya bermaksud untuk menyelesaikan Ax = b di mana A adalah matriks persegi atau persegi panjang empat persegi panjang yang rumit, jarang, tidak simetris dan sangat tidak terkondisi (kondisi nomor ~ 1E + 20). Saya telah berhasil menyelesaikan sistem dengan ZGELSS di LAPACK secara akurat. Tetapi seiring dengan meningkatnya derajat …

9 linear-solver  sparse-matrix  lapack  condition-number 

Misalkan saya memiliki sistem linier besar, jarang: . Sekarang, saya tidak memiliki A - 1 karena A terlalu besar untuk faktor atau segala bentuk penguraian A , tetapi berasumsi bahwa saya memiliki solusi x 0 yang ditemukan dengan penyelesaian berulang.Ax0=b0Ax0=b0A\textbf{x}_0=\textbf{b}_0A- 1A−1A^{-1}SEBUAHAAx0x0\textbf{x}_0 Sekarang, saya ingin menerapkan pembaruan peringkat kecil ke diagonal …

9 linear-algebra  iterative-method  sparse-matrix 

Diberikan matriks jarang generik dengan m << n (koreksi: ) elemen tidak nol (biasanya ). adalah generik dalam arti bahwa ia tidak memiliki sifat spesifik (misalnya kepastian positif), dan tidak ada struktur (misalnya pita) yang diasumsikan.A∈Rn×nA∈Rn×nA \in \mathbb{R}^{n\times n}m≪n2m≪n2m \ll n^2m∈O(n)m∈O(n)m \in {\cal O}(n)AAA Apa saja metode numerik yang baik …

9 linear-algebra  algorithms  sparse-matrix