У меня сейчас очень много работы по расчету эволюции различных конфигураций задачи трёх тел. Для моделирования несколько лет назад была написана программа, которая использует несколько специфических алгоритмических трюков, о которых я здесь не стану распространяться, позволяющих на порядки уменьшить ошибку численного интегрирования и избежать ( Read more... )

Недавно я наблюдал, как студенты подходят к решению одной простой задачи, предлагаемой в курсе астрофизики. Требуется оценить время коллапса межзвездного пылевого облака массы M и радиуса R в предположении, что облако достаточно разреженное, момент количества движения равен нулю, и коллапс происходит только под действием самогравитации и сферически ( Read more... )

Я понимаю, что для ответа на мой вопрос самым правильным было бы посмотреть исходники R, но прежде чем погрузиться в разбор чужого кода, я хотел бы задать вопрос здесь: а вдруг кто-нибудь из моих друзей уже знает ответ? Проблема возникла из-за того, что мне понадобилось заменить скрипт на R скриптом на Perl. Реализация функции, вычисляющей точный ( Read more... )

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

Пишу популярную статью про гравитационные волны и рисую к ней наглядные картинки. Задался вопросом: как наиболее просто и понятно представить искривление пространства-времени (возможно, в ущерб строгости).


( Read more... )

Octave огорчает всё больше и больше. Вот, казалось бы, совсем простая задача: интегрирование системы дифференциальных уравнений методом Эйлера. Проще уже не бывает, и что же мы получаем?

Время работы Matlab - 1.2 c
Время работы Octave - 101.1 c

( Посмотреть код скрипта )

В предыдущем посте я пытался сравнивать различные программы компьютерной алгебры на предмет удобства и быстродействия. В итоге ни одной из них, за исключением разве что Mathematica, но и то только благодаря тому, что я с ней давно работал, не удалось отдать предпочтения ( Read more... )

Подходит к концу моя студенческая лицензия на Mathematica 9 и я задумался о приобретении собственной персональной лицензии на один из программных продуктов. Выбор профессиональных инструментов не особенно велик. Собственно, всего три варианта:

Maple Personal edition - $299.00
Mathematica Home Edition - $295.00Matlab + Image Processing Toolbox + ( Read more... )

( Read more... )

Недавно schegloff инициировал у себя в журнале дискуссию о механизмах ускорения прогресса. И там в комментариях возник вопрос, а можно ли это как-то формализовать и смоделировать.

И вот, как будто в ответ на «Элементах» появляется статья, рассказывающая о том, что американские и британские ученые разработали компьютерную модель культурно-социальной ( Read more... )