Ich möchte eine dichte quadratische Übergangsmatrix direkt ändern, indem ich die Reihenfolge mehrerer Zeilen und Spalten mithilfe der Numpy-Bibliothek von Python ändere. Mathematisch entspricht dies einer Vormultiplikation der Matrix mit der Permutationsmatrix P und einer Nachmultiplikation mit P ^...
Wie skalieren Python / Numpy-Arrays mit zunehmenden Array-Dimensionen? Dies basiert auf einem Verhalten, das ich beim Benchmarking von Python-Code für diese Frage festgestellt habe: Wie kann man diesen komplizierten Ausdruck mit numpy-Slices ausdrücken? Das Problem bestand hauptsächlich in der...
λn u mλnum\lambda_{num}λa n aλeinnein\lambda_{ana} Um den Eigenwert 40 beginnen die numerischen Ergebnisse von den analytischen zu abweichen. Das überrascht mich nicht (ich gehe hier nicht auf das Warum ein, es sei denn, es wird in der Diskussion erwähnt). Was mich jedoch überrascht , ist, dass...
Ich möchte den folgenden Ausdruck in Python implementieren: wobei x und y numpy Arrays der Größe n sind und k ein numpy Array ist der Größe n × n . Die Größe n kann bis zu 10000 betragen, und die Funktion ist Teil einer inneren Schleife, die viele Male ausgewertet wird. Daher ist die...
Ich habe einen Rückwärts-Euler-Löser in Python 3 implementiert (mit Numpy). Zu meiner eigenen Bequemlichkeit und als Übung habe ich auch eine kleine Funktion geschrieben, die eine endliche Differenzapproximation des Gradienten berechnet, so dass ich den Jacobian nicht immer analytisch bestimmen...
Ich löse Differentialgleichungen, die das Invertieren dichter quadratischer Matrizen erfordern. Diese Matrixinversion nimmt den größten Teil meiner Rechenzeit in Anspruch, daher habe ich mich gefragt, ob ich den schnellsten verfügbaren Algorithmus verwende. Meine aktuelle Wahl ist numpy.linalg.inv...
Zur Modellreduktion möchte ich die linken Singularvektoren berechnen, die den - sagen wir 20 - größten Singularwerten einer Matrix , wobei und . Leider wird meine Matrix ohne Struktur dicht sein. N ≤ 10 6 k ≤ 10 3 A.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx 10^6k≈103k≈103k\approx 10^3AAA...
Ich bin ein Physiker, der versucht, die Strom-Spannungs-Eigenschaften eines Supraleiter-Supraleiter-Übergangs zu modellieren. Die Gleichung für dieses Modell lautet: ich( V.) = 1e R.n - n∫∞- ∞| E.|[ E.2- Δ21]]1 / 2| E.+ e V.|[ ( E.+ e V.)2- Δ22]]1 / 2[ f( E.) - f( E.+ e V.) ]d E.ich(V.)=1eR.n-...
Ich habe -dimensionalen Matrizen H ( → K ) in Abhängigkeit von Vektorparametern → k .nnnH^(k⃗ )H^(k→)\mathrm{\hat{H}}(\vec{k})k⃗ k→\vec{k} Nun geben Eigenwertroutinen Eigenwerte in keiner bestimmten Reihenfolge zurück (sie werden normalerweise sortiert), aber ich möchte Eigenwerte als glatte...