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')
Pergunta
Ayanami
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
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.