Pytania otagowane jako scipy

1
Kiedy Newton-Krylov nie jest odpowiednim rozwiązaniem?
Ostatnio porównywałem różne nieliniowe solwery z scipy i byłem pod szczególnym wrażeniem przykładu Newtona-Kryłowa w książce kucharskiej Scipy, w którym rozwiązują równanie różniczkowe drugiego rzędu z nieliniowym terminem reakcji w około 20 wierszach kodu. Zmodyfikowałem przykładowy kod, aby rozwiązać nieliniowe równanie Poissona ( zwane również równaniem Poissona-Boltzmanna , patrz strona …

1
Czy przybliżony jakobian ze skończonymi różnicami może powodować niestabilność w metodzie Newtona?
Zaimplementowałem solver z Eulerem wstecznym w Pythonie 3 (używając numpy). Dla własnej wygody i jako ćwiczenie napisałem również małą funkcję, która oblicza przybliżoną różnicę skończoną różnicy gradientu, aby nie zawsze musiałem określać analitycznie jakobian (jeśli to w ogóle możliwe!). Korzystając z opisów podanych w Ascher i Petzold 1998 , napisałem …

2
Jaki jest najszybszy sposób obliczenia wszystkich wartości własnych bardzo dużej i rzadkiej macierzy przylegania w pythonie?
Próbuję dowiedzieć się, czy istnieje szybszy sposób obliczenia wszystkich wartości własnych i wektorów własnych bardzo dużej i rzadkiej macierzy przylegania niż przy użyciu scipy.sparse.linalg.eigsh O ile mi wiadomo, metody te wykorzystują tylko rzadkość i atrybuty symetrii macierzy. Macierz przylegania jest również binarna, co sprawia, że ​​myślę, że istnieje szybszy sposób, …

1
Algorytm równoważenia macierzy
Pisałem przybornik systemu sterowania od zera i wyłącznie w Python3 (bezwstydna wtyczka:) harold. Na podstawie moich wcześniejszych badań zawsze narzekałem na solver Riccati care.mz powodów technicznych / nieistotnych. Dlatego piszę własny zestaw procedur. Jednej rzeczy, której nie mogę znaleźć, jest uzyskanie algorytmu równoważenia o wysokiej wydajności, przynajmniej tak dobrego balance.m. …
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.