![[personal profile]](https://www.dreamwidth.org/img/silk/identity/user.png)
Обсчитывать многократные суммы от всяких математических функций.
Типа f(x) = sum(i=-inf..inf)sum(j=-inf..inf) g(i,j,x)
где g(i,j,x) = (besseli(i,i+j,x))/besseli0(i+j,x)
Не смущайтесь бесконечностями, функция имеет более-менее ненулевые значения в ограниченных пределах. Эпсилон для определения пределов на ходу задавать надо.
Требования:
- работать быстро, рассчитывать надо дофига;
- Экономить CPU - кэшировать уже рассчитанные суммы, чтобы не считать повторно;
- Экономить память - кэш не должен есть много;
- Крайне желательно сохранять кэш между запусками;
Спец. пакет, если такой есть, подойдёт. Mathematica, скажем, не тянет.
C - это хорошо, но слишком много их, хороших и разных.
Типа f(x) = sum(i=-inf..inf)sum(j=-inf..inf) g(i,j,x)
где g(i,j,x) = (besseli(i,i+j,x))/besseli0(i+j,x)
Не смущайтесь бесконечностями, функция имеет более-менее ненулевые значения в ограниченных пределах. Эпсилон для определения пределов на ходу задавать надо.
Требования:
- работать быстро, рассчитывать надо дофига;
- Экономить CPU - кэшировать уже рассчитанные суммы, чтобы не считать повторно;
- Экономить память - кэш не должен есть много;
- Крайне желательно сохранять кэш между запусками;
Спец. пакет, если такой есть, подойдёт. Mathematica, скажем, не тянет.
C - это хорошо, но слишком много их, хороших и разных.
Tags:
Курить-1
21/5/08 18:11 (UTC)http://en.wikipedia.org/wiki/Category:Numerical_programming_languages
http://en.wikipedia.org/wiki/J_(programming_language)
http://freemat.sourceforge.net/
(no subject)
21/5/08 18:33 (UTC)C++ в связке с boost,
Посмотри сюда, тут дофига функций: http://www.boost.org/doc/libs/1_35_0/libs/math/doc/sf_and_dist/html/index.html
Даже если там чего-то и нет, C++ элементарно интегрируется с фортраном. Для фортрана есть все.
Посмотри еще в сторону matlab/octave, эта хрень интерпретируемая, плюс не факт, что будет легко интегрироваться. Зато если нужно что-то быстро задрафтить и график нарисовать - вещь идеальная (для дипломного проекта самое оно).
(no subject)
21/5/08 18:51 (UTC)Хотя, что-то для бесконечных сумм там есть...