2015-01-21 6 views
2

Я использую версию C LSODA от , что было переведено с FORTRAN.Решатель LSODA ODE в C

Моя проблема заключается в том, что она, похоже, не поддерживает пользователя, поставленного Jacobian. Аргумент jac в lsoda() отсутствует в коде (но не в документации, которая, вероятно, была скопирована).

Есть ли другая версия lsoda, написанная на C, которая позволяет пользователю поставлять Jacobian?

Благодаря

+2

Есть ли важная причина, по которой вы не можете назвать оригинал? –

+0

Да, моя функция f (в ​​y '= f * y) довольно сложна и написана на C, я не хочу писать ее в fortran. – user8476

+0

Это не объясняет, почему вы не можете назвать это. –

ответ

0

Там два основных варианта открыты для вас: Используйте CVODES пакет солнечных часов, или использовать функцию iso C binding FORtran, чтобы написать обертку вокруг вашей функции C.

Смежные вопросы