Eu quero modificar uma matriz de transição quadrada densa no local, alterando a ordem de várias de suas linhas e colunas, usando a biblioteca numpy do python. Matematicamente, isso corresponde à pré-multiplicação da matriz pela matriz de permutação P e à pós-multiplicação por P ^ -1 = P ^ T, mas...
Como as matrizes Python / Numpy escalam com o aumento das dimensões da matriz? Isso se baseia em algum comportamento que notei ao comparar o código Python para esta pergunta: Como expressar essa expressão complicada usando fatias numpy O problema envolveu principalmente a indexação para preencher...
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...
Desejo implementar o seguinte expressão em Python: onde x e y são matrizes numpy de tamanho n , e k é uma matriz numpy de tamanho n × n . O tamanho n pode ser de até 10000, e a função faz parte de um loop interno que será avaliado muitas vezes, portanto a velocidade é importante.xEu= ∑j = 1i -...
Eu implementei um solucionador de Euler para trás no python 3 (usando numpy). Para minha própria conveniência e como exercício, também escrevi uma pequena função que calcula uma aproximação de diferença finita do gradiente para que nem sempre seja necessário determinar analiticamente o jacobiano...
Estou resolvendo equações diferenciais que exigem inverter matrizes quadradas densas. Essa inversão de matriz consome a maior parte do meu tempo de computação, então eu queria saber se estou usando o algoritmo mais rápido disponível. Minha escolha atual é numpy.linalg.inv . Dos meus números, vejo...
Para redução do modelo, quero calcular os vetores singulares esquerdos associados aos - digamos 20 - maiores valores singulares de uma matriz , em que N ≈ 10 6 e k ≈ 10 3 . Infelizmente, minha matriz A será densa sem nenhuma estrutura.A ∈ RN, kUMA∈RN,kA \in \mathbb R^{N,k}N≈ 106N≈106N\approx 10^6k...
Eu sou um físico que está tentando modelar as características de tensão de corrente de uma junção supercondutor-supercondutor. A equação para este modelo é: Eu( V) = 1e Rn - n∫∞- ∞| E|[ E2- Δ21 1]1 / 2| E+ e V|[ ( E+ e V)2- Δ22]1 / 2[ f( E) - f( E+ e V) ]d
Eu tenho matrizes dimensionais dependendo do parâmetro do vetor .H ( → k ) → knnnH^( k⃗ )H^(k→)\mathrm{\hat{H}}(\vec{k})k⃗ k→\vec{k} Agora, as rotinas de autovalores retornam autovalores em nenhuma ordem específica (geralmente são classificadas), mas quero rastrear os autovalores como funções...