J'ai utilisé GSL comme base de plusieurs de mes simulations, mais c'est un peu exagéré pour mes besoins et il définit son propre type complexe pour des raisons héritées. Plutôt que de coder mon propre solveur Runge-Kutta ODE, qui ne serait probablement pas très efficace, existe-t-il des solveurs ODE open source qui utilisent le type complexe natif C99?
12
Je ne sais pas où voulez-vous l'utiliser, mais en général RK est assez difficile à implémenter de manière non efficace ... Avez-vous fait des tests de performance qui ont montré que vous avez ce problème?
mbq
2
Aucun. Je n'ai pas écrit le mien parce que je ne veux pas réinventer la roue. Si je le dois, je le ferai, mais trouver du temps à consacrer à quelque chose qui n'est pas cassé n'est pas prévu pour moi en ce moment. Si une réponse arrive, c'est ce que je cherche, je ne pourrai pas l'utiliser si pendant quelques mois. De plus, RK n'est pas toujours ce dont j'ai besoin, juste ce pour quoi je connais l'algorithme.
qubyte
Soit dit en passant, je fais la plupart du temps des simulations de petits systèmes quantiques. Pas exclusivement cependant.
qubyte
Je déconseille de mettre en œuvre vous - même un RK de taille variable (sauf à des fins éducatives). Il y a beaucoup d'heuristiques impliquées dans la recherche de la taille de pas optimale.
Jitse Niesen
Comme je l'ai dit, tout ce que j'écrirais rapidement serait soit faux, soit lent. Est-il particulièrement difficile d'implémenter RK avec des entrées / sorties complexes? Je sais que vous pouvez simplement le diviser en deux parties réelles, mais c'est un peu ennuyeux!
qubyte