Kürzlich habe ich verschiedene nichtlineare Löser von scipy verglichen und war besonders beeindruckt vom Newton-Krylov-Beispiel im Scipy-Kochbuch, in dem sie eine Differentialgleichungsgleichung zweiter Ordnung mit nichtlinearem Reaktionsterm in etwa 20 Codezeilen lösen. Ich habe den Beispielcode...
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 versuche herauszufinden, ob es einen schnelleren Weg gibt, alle Eigenwerte und Eigenvektoren einer sehr großen und spärlichen Adjazenzmatrix zu berechnen, als mit scipy.sparse.linalg.eigsh Soweit ich weiß, verwendet diese Methode nur die Spärlichkeit und Symmetrieattribute der Matrix. Eine...
Ich habe eine Steuerungssystem-Toolbox von Grund auf neu und rein in Python3 geschrieben (schamloser Plug :) harold. Aus meiner bisherigen Forschung habe ich mich immer care.maus technischen / irrelevanten Gründen über den Riccati-Löser beschwert . Daher habe ich meine eigenen Routinen geschrieben....
Ich suche nach dem minimalen Pfad zwischen den Minima einer potenziellen Oberfläche, die bereits in einem Raster bekannt ist. (Quelle: http://www.math.nus.edu.sg/~matrw/string/ ) Jeder Punkt auf dem Pfad befindet sich in allen Richtungen senkrecht zum Pfad auf einem potenziellen Minimum. Gibt...