Perguntas com a marcação «sparse-matrix»

Eu tenho vários problemas desafiadores de otimização global não convexa para resolver. Atualmente, uso o Optimization Toolbox do MATLAB (especificamente, fmincon()com o algoritmo = 'sqp'), o que é bastante eficaz . No entanto, a maior parte do meu código está em Python, e eu adoraria fazer a otimização em Python …

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 

Sistemas lineares esparsos aparecem com crescente frequência nas aplicações. Um deles tem muitas rotinas para escolher para resolver esses sistemas. No nível mais alto, existe um divisor de águas entre os métodos diretos (por exemplo, eliminação Gaussiana esparsa ou decomposição de Cholesky, com algoritmos especiais de ordenação e métodos multifrontais) …

49 linear-algebra  linear-solver  sparse-matrix 

EDIT: Estou testando se algum autovalor tem uma magnitude de um ou mais. Preciso encontrar o maior autovalor absoluto de uma grande matriz esparsa e não simétrica. Eu tenho usado a eigen()função de R , que usa o QR algo do EISPACK ou LAPACK para encontrar todos os autovalores e, …

27 linear-algebra  algorithms  eigensystem  sparse-matrix 

Eu tenho matrizes e . é escasso e é com muito grande (pode ser da ordem de vários milhões.) é uma matriz de altura com bastante pequeno ( ) e cada coluna pode só tem uma única entrada com o resto sendo 's, de tal modo que . é enorme, …

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

Estou escrevendo uma pequena biblioteca para cálculos matriciais esparsos como uma maneira de me ensinar a fazer o melhor uso da programação orientada a objetos. Eu trabalhei muito duro para ter um bom modelo de objeto, onde as partes (matrizes esparsas e os gráficos que descrevem sua estrutura de conectividade) …

17 linear-algebra  sparse-matrix  oop 

Eu queria saber se existe um método rápido e eficiente para encontrar antecipadamente o número de não-zeros para operações de multiplicação de matrizes esparsas, assumindo que ambas as matrizes estejam no formato CSC ou CSR. Eu sei que existe um no pacote smmp, mas preciso de algo que já esteja …

17 matrix  sparse-matrix 

Gostaria de saber se existe uma maneira rápida de calcular a distância euclidiana de dois vetores na oitava. Parece que não há função especial para isso, então devo usar apenas a fórmula com 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 

Estou tentando resolver uma equação de Poisson 2D por diferenças finitas. No processo, obtenho uma matriz esparsa com apenas variáveis ​​em cada equação. Por exemplo, se as variáveis ​​fossem U , a discretização produziria: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} Eu sei que posso resolver esse …

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

Estou desenvolvendo um código maior para realizar cálculos de autovalor de enormes matrizes esparsas, no contexto da física computacional. Testo minhas rotinas contra o oscilador harmônico simples em uma dimensão, uma vez que os autovalores são bem conhecidos analiticamente. Fazendo isso e comparando minhas próprias rotinas aos solucionadores embutidos do …

15 python  sparse-matrix  numpy  eigenvalues 

Eu tenho dois gráficos com quase n ~ 100000 nós cada. Nos dois gráficos, cada nó é conectado a exatamente 3 outros nós, portanto a matriz de adjacência é simétrica e muito esparsa. A parte difícil é que eu preciso de todos os autovalores da matriz de adjacência, mas não …

13 linear-algebra  sparse-matrix  eigenvalues  matrix 

Eu estou querendo saber se o algoritmo Thomas é a maneira mais rápida (provàvelmente?) De resolver um sistema tridiagonal esparso com dominação diagonal simétrica em termos de complexidade algorítmica (sem procurar pacotes de implementação como LAPACK etc.). Eu sei que o algoritmo Thomas e o multigrid são de complexidade , …

13 linear-solver  sparse-matrix  complexity  multigrid 

Preciso resolver o mesmo sistema linear esparso (300x300 a 1000x1000) com muitos lados do lado direito (300 a 1000). Além deste primeiro problema, eu também gostaria de resolver sistemas diferentes, mas com os mesmos elementos diferentes de zero (apenas valores diferentes), ou seja, muitos sistemas esparsos com padrão de esparsidade …

12 linear-algebra  petsc  linear-solver  sparse-matrix 

Por favor, liste o pacote Python (petsc4py, etc ...) e os solucionadores diretos esparsos que ele suporta. Uma resposta (wiki da comunidade) por pacote, por favor.

11 python  linear-solver  sparse-matrix 

Por exemplo, as bibliotecas de matriz esparsa C ++ que eu usei - Eigen e SuiteSparse, elas parecem não ter nenhuma funcionalidade SVD para matriz esparsa. Então, apenas curioso, SVD é mais difícil que QR / LU para matriz esparsa?

10 linear-algebra  sparse-matrix  matrix  svd  eigen 

Estou trabalhando em um sistema de equações linear esparso grande e mal condicionado. Quero usar aritmética duplo-duplo ou aritmético quad-duplo para resolvê-los. Eu sei que existe um pacote chamado MPACK desenvolvido por Nakata, Maho, que pode executar cálculos algébricos lineares numéricos sob aritmética quádrupla. No entanto, ele é projetado para …

10 sparse-matrix  matrix  condition-number