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

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

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

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, por isso é realmente difícil de inverter, e eu posso resolver um...

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

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

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

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

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

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

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 dos autovetores. Para ser...

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

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

A multiplicação da matriz (Mat * Mat e Mat * Vec) é dimensionada com o número de não-zeros ou com o tamanho da matriz? Ou alguma combinação dos dois. Que tal com forma. Por exemplo, eu tenho uma matriz 100 x 100 com 100 valores ou uma matriz 1000 x 1000 com 100 valores. Ao esquadrar essas...

A compilação do PETSc com uma biblioteca externa BLAS / LAPACK afeta significativamente o desempenho em matrizes esparsas ou apenas usa essas bibliotecas para matemática densa de

Atualmente, estou trabalhando na solução de sistemas simétricos muito grandes (mas não positivos definidos), gerados por alguns algoritmos. Essas matrizes têm uma boa esparsidade de bloco que pode ser usada para solução paralela. Mas não consigo decidir se devo usar uma abordagem direta (como...

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

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

Dada uma matriz esparsa genérica com m << n (correção: m ≪ n 2 ) elementos diferentes de zero (normalmente m ∈ O ( n ) ). A é genérico no sentido de que não possui propriedades específicas (por exemplo, definição positiva), e nenhuma estrutura (por exemplo, bandagem) é assumida.A∈Rn×nA∈Rn×nA...

Suponha que eu tenha o sistema linear grande e esparso original: . Agora, eu não tenho pois A é muito grande para fatorar ou qualquer tipo de decomposição de , mas suponha que eu tenha a solução encontrada com uma solução iterativa.A - 1 A x 0A x0 0= b0 0Ax0=b0A\textbf{x}_0=\textbf{b}_0UMA-...