Я використовую GSL в якості основи багатьох моїх моделювань, але це трохи надмірно для моїх цілей, і він визначає власний складний тип із застарілих причин. Замість того, щоб кодувати мій власний ODE-розв'язувач Runge-Kutta, який, мабуть, не буде дуже ефективним, чи існують якісь відкриті вихідні джерела ODE, які використовують рідний комплексний тип C99?
Я не знаю, де ви хочете його використовувати, але загалом RK досить важко реалізувати неефективним способом ... Ви зробили якісь орієнтири, які б показали, що у вас є ця проблема?
—
mbq
Жоден. Я не написав свого, тому що не хочу винаходити колесо. Якщо мені доведеться, я зроблю, але знайти час, щоб витратити на щось, що не порушено, зараз немає для мене картки. Якщо з’явиться відповідь, це те, що я шукаю, я не зможу реально користуватися протягом декількох місяців. Крім того, РК - це не завжди те, що мені потрібно, а лише те, для чого я знаю алгоритм.
—
квіте
Між іншим, я займаюся симуляцією малих квантових систем більшу частину часу. Не виключно, хоча.
—
квіте
Я радив би самостійно не застосовувати РК із змінним кроком (крім освітніх цілей). Для пошуку оптимального розміру кроку є багато евристики.
—
Jitse Niesen
Як я вже говорив, все, що я швидко напишу, було б або неправильним, або повільним. Чи особливо важко реалізувати RK зі складним входом / виходом? Я знаю, що ви можете просто розділити його на дві реальні частини, але це щось дратує!
—
квіте