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