Pertanyaan yang diberi tag «sparse-matrix»

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 optimasi dengan Python juga. Apakah ada pemecah NLP dengan binding Python yang …

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

Sistem linear yang jarang muncul dengan meningkatnya frekuensi dalam aplikasi. Seseorang memiliki banyak rutinitas untuk dipilih untuk menyelesaikan sistem ini. Pada level tertinggi, ada batas antara metode langsung (mis. Eliminasi Gaussian yang jarang atau dekomposisi Cholesky, dengan algoritma pemesanan khusus, dan metode multirontal) dan iteratif (misalnya metode GMRES, (bi-) konjugat …

49 linear-algebra  linear-solver  sparse-matrix 

EDIT: Saya menguji apakah ada nilai eigen yang besarnya satu atau lebih. Saya perlu menemukan nilai eigen absolut terbesar dari matriks besar non-simetris yang jarang. Saya telah menggunakan eigen()fungsi R , yang menggunakan algo QR baik dari EISPACK atau LAPACK untuk menemukan semua nilai eigen dan kemudian saya gunakan abs()untuk …

27 linear-algebra  algorithms  eigensystem  sparse-matrix 

Saya memiliki matriks dan . jarang dan dengan sangat besar (bisa dalam urutan beberapa juta.) adalah matriks tinggi dengan agak kecil ( ) dan setiap kolom dapat hanya memiliki satu entri dengan sisanya menjadi 's, sehingga . sangat besar, jadi sangat sulit untuk membalikkan, dan saya bisa menyelesaikan sistem linear …

22 linear-algebra  linear-solver  sparse-matrix  krylov-method 

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 

Saya bertanya-tanya apakah ada metode cepat dan efisien untuk menemukan jumlah bukan nol di muka untuk operasi perkalian matriks jarang dengan asumsi kedua matriks dalam format CSC atau CSR. Saya tahu ada satu dalam paket smmp tapi saya butuh sesuatu yang sudah diimplementasikan dalam C atau C ++. Bantuan apa …

17 matrix  sparse-matrix 

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

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

Saya mencoba menyelesaikan persamaan 2D Poisson dengan perbedaan hingga. Dalam prosesnya, saya mendapatkan matriks jarang dengan hanya variabel di setiap persamaan. Misalnya, jika variabelnya adalah , maka diskritisasi akan menghasilkan:555UUU Ui−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jUi−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jU_{i-1,j} + U_{i+1,j} -4U_{i,j} + U_{i,j-1} + U_{i,j+1} = f_{i,j} Saya tahu bahwa saya dapat menyelesaikan sistem ini dengan metode …

15 linear-algebra  pde  finite-difference  sparse-matrix  banded-matrix 

Saya sedang mengembangkan beberapa kode yang lebih besar untuk melakukan perhitungan nilai eigen dari matriks jarang besar, dalam konteks fisika komputasi. Saya menguji rutinitas saya terhadap osilator harmonik sederhana dalam satu dimensi, karena nilai eigen dikenal secara analitis. Melakukannya dan membandingkan rutinitas saya sendiri dengan pemecah bawaan SciPy, saya telah …

15 python  sparse-matrix  numpy  eigenvalues 

Saya memiliki dua grafik dengan masing-masing hampir n ~ 100000 node. Dalam kedua grafik, setiap node terhubung tepat ke 3 node lainnya sehingga matriks adjacency simetris dan sangat jarang. Bagian yang sulit adalah saya membutuhkan semua nilai eigen dari matriks adjacency tetapi bukan vektor eigen. Agar akurat, ini akan menjadi …

13 linear-algebra  sparse-matrix  eigenvalues  matrix 

Saya bertanya-tanya apakah algoritma Thomas adalah cara tercepat (dapat dibuktikan?) Untuk menyelesaikan sistem tridiagonal simetris dominan diagonal dalam hal kompleksitas algoritmik (tidak mencari paket implementasi seperti LAPACK dll). Saya tahu bahwa kedua algoritma Thomas dan multigrid adalah kompleksitas, tetapi mungkin faktor konstan untuk multigrid kurang? Rasanya bagi saya multigrid tidak …

13 linear-solver  sparse-matrix  complexity  multigrid 

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 dengan pola sparsity konstan. Matriks saya …

12 linear-algebra  petsc  linear-solver  sparse-matrix 

Harap cantumkan paket Python (petsc4py, dll ...) dan pemecah langsung yang jarang didukungnya. Tolong, satu (komunitas-wiki) jawaban per paket.

11 python  linear-solver  sparse-matrix 

Misalnya pustaka matriks jarang C ++ yang saya gunakan - Eigen dan SuiteSparse, mereka tampaknya tidak memiliki funcitionality SVD untuk matriks jarang. Jadi hanya ingin tahu, apakah SVD lebih sulit daripada QR / LU untuk matriks jarang?

10 linear-algebra  sparse-matrix  matrix  svd  eigen 

Saya sedang mengerjakan beberapa sistem persamaan linear besar jarang yang tidak dikondisikan. Saya ingin menggunakan aritmatika double-double atau aritmatika quad-double untuk menyelesaikannya. Saya tahu bahwa ada paket bernama MPACK yang dikembangkan oleh Nakata, Maho, yang dapat melakukan perhitungan aljabar linear numerik di bawah aritmatika quad-double. Namun, ini dirancang untuk matriks …

10 sparse-matrix  matrix  condition-number