Perguntas com a marcação «matrix»

Eu tenho que fazer uma transformação de coordenadas entre dois sistemas de referência (eixos). Para isso, três matrizes ( ) devem ser multiplicadas devido ao uso de alguns eixos intermediários. Pensei em duas abordagens para resolver isso:3 × 33×33\times3 Método # 1 : Fazendo a multiplicação diretamente, ou seja, vf= …

10 matrix  performance  matlab 

Preciso para calcular uma série de 3 × 33×33\times3 matriz inversa (por Newton iteração decomposição polar), com muito pequeno número de casos degenerados ( &lt; 0,1 %&lt;0,1%<0.1\% ). Inverso explícito (via menores de matriz divididos por determinante) parece funcionar e tem cerca de ~ 32 ~ 40 fracassos fundidos (dependendo …

10 matrix  matrix-equations  inverse  matrix-factorization 

Nas aulas do FEM, geralmente é dado como certo que a matriz de rigidez é definida positivamente, mas eu simplesmente não consigo entender o porquê. Alguém poderia dar uma explicação? Por exemplo, podemos considerar o problema de Poisson: −∇2u=f,−∇2u=f, -\nabla^2 u = f, cuja matriz de rigidez é: Kij=∫Ω∇φi⋅∇φjdΩ,Kij=∫Ω∇φi⋅∇φjdΩ,K_{ij} = …

10 finite-element  matrix  stiffness 

Dada uma matriz simétrica definida positiva, qual é o algoritmo mais rápido para calcular a matriz inversa e seu determinante? Para problemas em que estou interessado, a dimensão da matriz é 30 ou menos. Alta precisão e velocidade são realmente necessárias. (milhões de matrizes são realizadas) O determinante é necessário. …

10 algorithms  matrix 

A rotina QR do LAPACK armazena Q como refletores domésticos. Ele escala o vetor de reflexão com , para que o primeiro elemento do resultado se torne , para que não precise ser armazenado. E ele armazena um vetor separado , que contém os fatores de escala necessários. Portanto, uma …

9 linear-algebra  matrix  lapack 

Estou interessado em calcular a solução de um sistema Lage de ODEs usando um método de krylov como em [1]. Esse método envolve funções relacionadas ao exponencial (as chamadas funções ). Consiste essencialmente em calcular a ação da função de matriz, construindo um subespaço de Krylov usando a iteração de …

9 linear-algebra  matrix  numerical-analysis  time-integration  krylov-method 

Precisamos calcular matrizes de covariância com tamanhos que variam de a . Temos acesso a GPUs e clusters, nos perguntamos qual é a melhor abordagem paralela para acelerar esses cálculos.100000 × 10000010000 × 1000010000×1000010000\times10000100000 × 100000100000×100000100000\times100000

9 matrix  parallel-computing  gpu 

A partir da definição do número da condição, parece que é necessária uma inversão da matriz para calculá-la. Será que para uma matriz quadrada genérica (ou melhor, se a definição positiva simétrica positiva) é possível explorar alguma decomposição da matriz para calcular o número da condição em um maneira mais …

9 linear-algebra  matrix  condition-number 

Suponha que o sistema linear a seguir seja dado que é o Laplaciano ponderado conhecido como positivo definido, com um espaço nulo unidimensional abrangido por , e a variação de conversão de , ou seja, não altera o valor da função (cuja derivada é ). As únicas entradas positivas de …

9 linear-algebra  optimization  iterative-method  matrix  linear-solver 

Suponha que eu conheça a esparsidade esperada de uma matriz (ou seja, o número de não-zeros / número total possível de não-zeros). Existe uma regra prática (talvez aproximada) para decidir se o armazenamento de matriz esparsa deve ser usado (especificamente, armazenamento de linhas compactadas) versus armazenado como uma matriz densa? …

9 matrix  sparse-matrix  dense-matrix 

Estou usando o Tatu para fazer multiplicações de matriz muito intensivas com comprimentos laterais , onde pode ser de até 20 ou mais. Estou usando o Armadillo com o OpenBLAS para multiplicação de matrizes, o que parece estar fazendo um trabalho muito bom em núcleos paralelos, exceto que tenho um …

9 c++  matrix  dense-matrix 

Eu tenho a seguinte situação: Eu tenho uma sequência de vetores e para cada eu quero calcular o produto onde é fixado desde o início. Embora não haja informações sobre a estrutura de , normalmente possui um padrão específico em que muitos valores são repetidos e eu gostaria de calcular …

8 matrix  linear-algebra 

Eu tenho a seguinte pergunta: Suponha que eu tenha duas matrizes ,X,YX,YX, Y ambas de tamanho m×pm×pm\times p e uma matriz gaussiana iid aleatória GGG de tamanho m×km×km \times k , m≫p&gt;km≫p&gt;km\gg p>k . Existe uma maneira rápida de calcular exp(−XYT)Gexp⁡(−XYT)G\exp(-XY^T)G ? Talvez usando o fato de que XXX e …

8 linear-algebra  matrix 

dada uma matriz simétrica , e uma matriz arbitrária X ∈ R n × n , e um vetor v ∈ R n × 1 , é possível calcular a seguinte expressão no tempo O ( n 2 ) ?Y∈Rn×nY∈Rn×nY \in \mathbb{R}^{n \times n}X∈Rn×nX∈Rn×nX \in \mathbb{R}^{n \times n}v∈Rn×1v∈Rn×1v \in \mathbb{R}^{n …

8 optimization  algorithms  performance  matrix  complexity 

Eu tenho esse problema no meu livro: Sugira um algoritmo eficiente para resolver sistemas de equações lineares com matriz tridimensional cíclica, que é da forma: sem trocar nenhuma linha e coluna. Estimar a complexidade.⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢uma1c20 0⋮0 0d2b1uma2⋱⋮⋯0 00 0b2⋱cn - 2⋯⋯⋯0 0⋱uman - 2cn - 10 00 0⋮0 0bn - 2uman …

8 matrix  linear-solver