Модифицированную функцию Бесселя, чем лучше посчитать? Опции

[Pathfinder]

Ищу библиотеку, которая позволяет вычислить модифицированную функцию Бесселя с экспоненциальным масштабным множителем. Пока что использую GSL из GNUWin32, но она не считает функцию с порядком выше 32-го.
Насколько я понял, в Матлабе используется вот эта библиотека: http://www.netlib.org/amos, но она на фортране ... боюсь, проблематично её будет с MSVC-шным проектом использовать.

[Pathfinder]

mdmitry
Насколько я понял, в Scientific Linux используется GSL, и как раз её я использую сейчас. Правда приходится пользоваться версией 1.8, поскольку свежие не понял, как под винду собрать...
За наводку про F2C спасибо, нашёл и попробовал. Исходник конечно не фонтан, но лучше чем ничего. Самое неприятное, что в результате сборки бесселя получается dll-ка размером 132КБ, поскольку приходится фортрановский рантайм прилинковывать. При таких размерах она скорее всего не влезет в страницу кэша и будет тормозить.

=GM=
Согласен, решение логичное, и в моём случае удобное, поскольку мне нужно считать I_0, I_n и I_2n для одного аргумента. Но смущает то, что его не используют (как мне показалось) ни в netlib.org/amos ни в GSL. Может при рекурсивном вычислении ошибка накапливается? Там ведь дробно-полиномиальный множитель будет...
А не помните, какая точность получалась в сравнении с таблицами? Мне очень желательно 12 знаков получить.

[=GM=]

=GM=
А не помните, какая точность получалась в сравнении с таблицами? Мне очень желательно 12 знаков получить

Давно было, лет 25 назад, считал на фортране ряды Пуассона, плохо сходящиеся, помню, что точность переменных была удвоенная, то ли 64 бита, то ли 80, с таблицами совпадало до последнего знака.