Eu tenho usado o GSL como base de muitas de minhas simulações, mas é um pouco exagerado para meus propósitos e define seu próprio tipo complexo por motivos legados. Em vez de codificar meu próprio solucionador Runge-Kutta ODE, o que provavelmente não seria muito eficiente, há algum solucionador ODE de código aberto que use o tipo complexo C99 nativo?
Não sei onde você deseja usá-lo, mas, em geral, o RK é bastante difícil de ser implementado de maneira não eficiente ... Você fez algum benchmark que demonstrasse que tem esse problema?
—
MBq
Nenhum. Não escrevi o meu próprio porque não quero reinventar a roda. Se for preciso, sim, mas encontrar tempo para gastar em algo que não está quebrado não está disponível para mim agora. Se surgir uma resposta, é o que estou procurando, não poderei usá-lo por alguns meses. Além disso, RK nem sempre é o que eu preciso, apenas o que eu sei o algoritmo.
—
qubyte
Aliás, estou fazendo simulações de pequenos sistemas quânticos na maioria das vezes. Não exclusivamente embora.
—
qubyte
Eu recomendaria não implementar a RK de tamanho de etapa variável (exceto para fins educacionais). Existem muitas heurísticas envolvidas na busca do tamanho ideal da etapa.
—
Jitse Niesen 07/12/11
Como eu disse, qualquer um que eu escrevesse rapidamente seria errado ou lento. É particularmente difícil implementar o RK com entrada / saída complexa? Eu sei que você pode dividi-lo em duas partes reais, mas isso é meio chato!
—
qubyte