Pull to refresh

Comments 5

Странный выбор. Вообще-то, один из самых популярных методов в современных коммерческих пакетах — метод Ньютона-Крылова. Очень быстро сходится.

Вот пример использования scipy.optimize.newton_krylov (даже без preconditioner'а) для решения аж интегро-дифференциального уравнения:

https://docs.scipy.org/doc/scipy/reference/optimize.nonlin.html

Ну, справедливости ради, бывает так, что релаксация может оказаться не хуже, а то и лучше прочих продвинутых методов (в нелинейных задачах правда). Я, во всяком случае, такое встречал.
Ну вот пусть автор публикации и сравнит. Хоть на своем примере, хоть на том, что из моей ссылки.
Нехорошо как-то граничные условия Неймана аппроксимированы в одномерной задаче, в двумерной, я так полагаю, кругом Дирихле? Второе, что такое «квадратичная матрице {a}»? Ну и центральные разности для конвекции в сочетании с релаксацией… Не очень понял, кому предназначена эта публикация? Школьникам старших классов? Студентам первого курса?
У Петра Николаевича Вабищевича защищался как раз по этой теме много лет назад. Эх, ностальгия
Sign up to leave a comment.

Articles