本文整理汇总了Python中amuse.community.hermite0.interface.Hermite.get_total_energy方法的典型用法代码示例。如果您正苦于以下问题：Python Hermite.get_total_energy方法的具体用法？Python Hermite.get_total_energy怎么用？Python Hermite.get_total_energy使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类amuse.community.hermite0.interface.Hermite的用法示例。
在下文中一共展示了Hermite.get_total_energy方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: evolve_system
# 需要导入模块: from amuse.community.hermite0.interface import Hermite [as 别名]
# 或者: from amuse.community.hermite0.interface.Hermite import get_total_energy [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()