Какой язык вы бы предложили для решения системы с:
с использованием Рунге Кутта 4-го порядка или подобного.
Скорость имеет много значение, но ею можно пожертвовать ради:
Я в основном между Lisp и OCaml, но любые другие предложения приветствуются.
Спасибо!
Вот реализация RK в Common Lisp:
http://github.com/bld/bld-ode/blob/master/rk.lisp
Преимущество Common Lisp в том, что вы можете начать с простого и элегантного кода, а затем ускорить выполнение критических частей (например, переключившись с преимущественно функциональных вычислений на вычисления с отслеживанием состояния или объявив типы).
SBCL имеет отличный компилятор собственного кода.
RK4 — очень простой метод, и для него уже написано множество отличных реализаций. Воспользуйтесь одним из них, а свои усилия потратьте на другие аспекты проекта.
Я не знаком с Рунге Кутта, но OCaml может обеспечить хорошую скорость и читабельность в целом, по крайней мере, если вы будете немного осторожны. Тогда у вас будет преимущество надежной статической системы типов для остальной части вашего приложения.
помимо всего прочего, вы можете написать привязки ocaml к вашему существующему решателю C runge-kutta.
Трудно сказать, какой язык будет проще всего, есть lisp, C++, C# и т. д. библиотеки для достижения этой цели, так что многое, если это связано с личными предпочтениями. Я бы предположил, что Matlab является наиболее адаптированным и элегантным решением специально для таких типов задач, и в нем есть много встроенная поддержка ODE... Лисп может быть медленным... и я не могу говорить за OCaml.
Я бы предложил использовать python + numpy + scipy, общая математическая и числовая поддержка (превосходные многомерные массивы) превосходна. В любом случае это зависит от конкретных потребностей.
Fortran или C, возможно, стоит изучить подпрограммы NAG. C был бы более гибким и простым для понимания, но Fortran обычно считается лучшим для числовых вычислений.
