Оптимизация числодробилки

[al_pas]

Рабочая программа для моделирования эволюции тройной звездной системы с учетом общей теории относительности непрерывно улучшалась и модифицировалась в течение последних десяти лет. В результате ее код превратился в такую кашу, что дальнейшая его поддержка стала невозможной

Comments

[vitus_wagner]

Вообще такие вещи бы оптимизирующий компилятор должен делать.

[al_pas]

Разумеется, но в задачах небесной механики, когда результат расчета на 10000000 шаге виляет хвостом на малейший чих в исходных данных, и весь код заточен не столько на быстродействие, сколько на минимизацию потерь точности, лучше не доверять автоматической оптимизации, поэтому она обычно отключена.

[v1adis1av]

А это проверялось? В смысле, если сравнить два расчёта с одинаковыми начальными условиями, один на программе, компилированной с оптимизацией, второй без, -- действительно появляется разница?

[al_pas]

Да, и очень существенная. У меня в программе используется коррекция, обеспечивающая сохранение полного момента и энергии с учетом потерь на излучение гравитационных волн. Шаг интегрирования тоже постоянно меняется, а про корректность алгоритма оптимизации gcc уже писал Линус Торвальдс.

[dims12]

Странно. Этим должен заниматься компилятор, а не программист.

[al_pas]

См. http://al-pas.livejournal.com/164026.html?thread=1045434#t1045434

[back_in_usa]

Так ведь ещё на первом курсе учат, что деление должно быть последним действием.