Ir para conteúdo
Fórum Script Brasil
  • 0

Como obter informações adicionais dos solvers do “scipy.integrate.ode”?


Ayanami

Pergunta

Oi pessoal,

Eu estou comparando diferentes solvers ( dopri5, dop853, vode, etc) do “scipy.integrate.ode” e gostaria de saber se há  algum modo de obter informações como: numero de iterações, numero de passos utilizados, tempo de simulação, ...

Eu sei como conseguir esse tipo de informação no "scipy.integrate.odeint" usando, por exemplo, "infodict". No entanto, não descobri nenhuma maneira de fazer isso para os outros solvers que estou testando. A proposito, para descobrir o tempo necessário, estou usando "time.clock()" como abaixo (não sei se existe um modo melhor de fazer isso) :

time_start_RK4 = time.clock()                       

solver = ode(reactor_ode)
solver.set_integrator('dopri5')
solver.set_initial_value(Initial_conditions, t0)


solution_RK4 = np.empty((Nt, 8*N))
solution_RK4[0, :] = Initial_conditions
k = 1


while solver.successful() and solver.t < tf:
   solver.integrate(t[k])
   solution_RK4[k] = solver.y
   k += 1

time_elapsed_RK4 = (time.clock() - time_start_RK4)  

print(time_elapsed_RK4, 's')

Muito obrigado desde já! 

Link para o comentário
Compartilhar em outros sites

0 respostass a esta questão

Posts Recomendados

Até agora não há respostas para essa pergunta

Participe da discussão

Você pode postar agora e se registrar depois. Se você já tem uma conta, acesse agora para postar com sua conta.

Visitante
Responder esta pergunta...

×   Você colou conteúdo com formatação.   Remover formatação

  Apenas 75 emoticons são permitidos.

×   Seu link foi incorporado automaticamente.   Exibir como um link em vez disso

×   Seu conteúdo anterior foi restaurado.   Limpar Editor

×   Você não pode colar imagens diretamente. Carregar ou inserir imagens do URL.



  • Estatísticas dos Fóruns

    • Tópicos
      152,3k
    • Posts
      652,5k
×
×
  • Criar Novo...