Preguntas etiquetadas con numpy

Quiero modificar una matriz de transición cuadrada densa en el lugar cambiando el orden de varias de sus filas y columnas, usando la biblioteca numpy de python. Matemáticamente, esto corresponde a la multiplicación previa de la matriz por la matriz de permutación P y la multiplicación posterior por P ^ …

27 linear-algebra  python  numpy 

¿Cómo se escalan las matrices Python / Numpy con el aumento de las dimensiones de la matriz? Esto se basa en un comportamiento que noté al comparar el código Python para esta pregunta: cómo expresar esta expresión complicada usando cortes numpy El problema consistía principalmente en la indexación para llenar …

21 python  performance  numpy 

Estoy desarrollando un código más grande para realizar cálculos de valores propios de enormes matrices dispersas, en el contexto de la física computacional. Pruebo mis rutinas contra el oscilador armónico simple en una dimensión, ya que los valores propios son bien conocidos analíticamente. Al hacerlo y al comparar mis propias …

15 python  sparse-matrix  numpy  eigenvalues 

Deseo implementar la siguiente expresión en Python: donde e son matrices numpy de tamaño , y es una matriz numpy de tamaño . El tamaño puede ser de hasta aproximadamente 10000, y la función es parte de un bucle interno que se evaluará muchas veces, por lo que la velocidad …

14 python  numpy 

He implementado un solucionador de Euler hacia atrás en Python 3 (usando numpy). Para mi propia conveniencia y como ejercicio, también escribí una pequeña función que calcula una aproximación de diferencia finita del gradiente para que no siempre tenga que determinar el jacobiano analíticamente (¡si es posible!). Usando las descripciones …

13 pde  stability  numpy  scipy  newton-method 

Estoy resolviendo ecuaciones diferenciales que requieren invertir matrices cuadradas densas. Esta inversión matricial consume la mayor parte de mi tiempo de cálculo, por lo que me preguntaba si estoy usando el algoritmo más rápido disponible. Mi elección actual es numpy.linalg.inv . De mis números veo que se escala como donde …

11 numpy  dense-matrix  inverse 

Para la reducción del modelo, quiero calcular los vectores singulares izquierdos asociados a los - digamos 20 - valores singulares más grandes de una matriz , donde y . Desafortunadamente, mi matriz será densa sin ninguna estructura. N ≈ 10 6 k ≈ 10 3 AA∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx 10^6k≈103k≈103k\approx …

10 linear-algebra  python  reference-request  numpy  svd 

Soy un físico que intenta modelar las características de voltaje de corriente de una unión superconductor-superconductor. La ecuación para este modelo es: I(V)=1eRn−n∫∞−∞|E|[E2−Δ21]1/2|E+eV|[(E+eV)2−Δ22]1/2[f(E)−f(E+eV)]dEI(V)=1eRn−n∫−∞∞|E|[E2−Δ12]1/2|E+eV|[(E+eV)2−Δ22]1/2[f(E)−f(E+eV)]dE\begin{align} I(V) = \frac{1}{eR_{\mathrm{n-n}}}\int_{-\infty}^{\infty}\frac{|E|}{[E^{2} - \Delta_{1}^{2}]^{1/2}}\frac{|E + eV|}{[(E + eV)^{2} - \Delta_{2}^{2}]^{1/2}}[f(E) - f(E + eV)]\,\mathrm{d}E \end{align} Los valores de corriente ( o en el código) se calculan …

9 python  numpy  integral-equations 

Tengo matrices -dimensional H ( → k ) en función de vector de parámetros → k .nnnH^(k⃗ )H^(k→)\mathrm{\hat{H}}(\vec{k})k⃗ k→\vec{k} Ahora, las rutinas de valores propios devuelven valores propios sin un orden particular (generalmente están ordenados), pero quiero rastrear valores propios como funciones suaves de → k . Debido a que …

8 linear-algebra  matrix  eigensystem  numpy