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...
¿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...
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...
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 es importante.Xyo= ∑j =...
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...
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...
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...
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
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...