def initialize_huayno(bodies, converter, huayno_eta): stars_gravity = Huayno(converter,channel_type="sockets") stars_gravity.particles.add_particles(bodies) # either bodies or just stars stars_gravity.commit_particles() #stars_gravity.set_timestep_parameter(huayno_eta) stars_gravity.parameters.timestep_parameter = huayno_eta stars_gravity.set_inttype_parameter(12) # CC_KEPLER return stars_gravity
def initialize_huayno(bodies, converter, huayno_eta): stars_gravity = Huayno(converter,channel_type="sockets") stars_gravity.particles.add_particles(bodies) # either bodies or just stars stars_gravity.commit_particles() #stars_gravity.set_timestep_parameter(huayno_eta) stars_gravity.parameters.timestep_parameter = huayno_eta stars_gravity.set_inttype_parameter(12) # CC_KEPLER return stars_gravity
def initialize_huayno(bodies, converter, huayno_eta):
""" initialize huayno code """
gravity = Huayno(converter,channel_type="sockets")
gravity.particles.add_particles(bodies) # either bodies or just stars
gravity.commit_particles()
#gravity.set_timestep_parameter(huayno_eta)
gravity.parameters.timestep_parameter = huayno_eta
gravity.set_inttype_parameter(8) # CC_KEPLER
"""
Huayno Options:
(pass / hold / bridge / maybe shared? are described in the Huayno Paper)
SHARED2=1
EXTRAPOLATE=5
PASS_KDK=2
PASS_DKD=7
HOLD_KDK=3
HOLD_DKD=8 <<<<-------
PPASS_DKD=9
BRIDGE_KDK=4
BRIDGE_DKD=10
CC=11
CC_KEPLER=12 <<<<-------
OK=13
KEPLER=14
SHARED4=15
SHARED6=18
SHARED8=19
SHARED10=20
SHAREDBS=21 **BS**
CCC=22
CCC_KEPLER=23
CC_BS=24 **BS**
CCC_BS=25 **BS**
BS_CC_KEPLER=26 **BS**
CC_BSA=27 **BS** ???
CCC_BSA=28 **BS** ???
SHARED2_COLLISIONS=29
SHARED4_COLLISIONS=30
SHARED6_COLLISIONS=31
SHARED8_COLLISIONS=32
SHARED10_COLLISIONS=33
"""
return gravity
