本文整理匯總了Python中amuse.community.hermite0.interface.Hermite.set_dt_param方法的典型用法代碼示例。如果您正苦於以下問題:Python Hermite.set_dt_param方法的具體用法?Python Hermite.set_dt_param怎麽用?Python Hermite.set_dt_param使用的例子?那麽, 這裏精選的方法代碼示例或許可以為您提供幫助。您也可以進一步了解該方法所在類amuse.community.hermite0.interface.Hermite
的用法示例。
在下文中一共展示了Hermite.set_dt_param方法的1個代碼示例,這些例子默認根據受歡迎程度排序。您可以為喜歡或者感覺有用的代碼點讚,您的評價將有助於係統推薦出更棒的Python代碼示例。
示例1: evolve_system
# 需要導入模塊: from amuse.community.hermite0.interface import Hermite [as 別名]
# 或者: from amuse.community.hermite0.interface.Hermite import set_dt_param [as 別名]
def evolve_system(particles):
"""
Evolves the system using the Hermite integrator.
Parameters
----------
particles: amuse.datamodel.particles.Particles instance
"""
times = numpy.linspace(0.0001, args.time, args.steps) |u.yr
intr = Hermite(nbody_to_si(particles.total_mass(), 1 | u.AU))
intr.particles.add_particles(particles)
energy_begin = intr.get_total_energy()
if args.dt is not None:
intr.set_dt_param(args.dt)
for t in times:
intr.evolve_model(t)
energy_error = (intr.get_total_energy() - energy_begin)/energy_begin
print(intr.get_time_step().in_(u.day), intr.get_time().in_(u.yr))
print(energy_error)
print("energy error:{}".format(energy_error))
intr.stop()