LINUX.ORG.RU

История изменений

Исправление AIv, (текущая версия) :

Насколько я понял вики, d omega = 0 это как раз и есть интеграл движения.

В общем, ошибка округления приводит к нарушению этой самой «симплектичности» только если ее сохраняет числ. схема (схема консервативна). Для многих задач схемы неконсервативны (напр. для задачи многих теЛ полностью консервативная схема вообще не строится), и энергия нарушается за счет схемного источника.

Если же схема консервативна, то увидев что на тестах та же энергия сохраняется с машинной точностью, все кричат «алиллуйя» и идут квасить, поскольку для борьбы с ошибкой округления существует double;-)

Исходная версия AIv, :

Насколько я понял вики, d omega = 0 это как раз и есть интеграл движения.

В общем, ошибка округления приводит к нарушению этой самой «симплектичности» только если ее сохраняет числ. схема (схема консервативна). Для многих задач схемы неконсервативны (напр. для задачи многих тех полностью консервативная схема вообще не строится), и энергия нарушается за счет схемного источника.

Если же схема консервативна, то увидев что на тестах та же энергия сохраняется с машинной точностью, все кричат «алиллуйя» и идут квасить, поскольку для борьбы с ошибкой округления существует double;-)