Откуда берет свое начало метод адаптивного размера шага SciPy для конечных разностей? ⇐ Python

[Anonymous]

Внутри класса KrylovJacobian из SciPy есть этот метод:

Код: Выделить всё

def _update_diff_step(self):
mx = abs(self.x0).max()
mf = abs(self.f0).max()
self.omega = self.rdiff * max(1, mx) / max(1, mf)

что будет то же самое, что:
[img]https://i.sstatic. net/FyfKhk1V.png[/img]

Это изменяет размер шага, который использует метод конечных разностей, однако я не могу найти происхождение этого выражения или то, почему оно работает.
Кто-нибудь знает источник этого метода или причина этого?

Подробнее здесь: https://stackoverflow.com/questions/792 ... -originate