@@ -161,8 +161,13 @@ int get_total_radius(double * radius){
161161 return 0 ;
162162}
163163
164- int new_particle (int * index_of_the_particle, double mass, double x,
165- double y, double z, double vx, double vy, double vz, double radius, int code_index){
164+ int new_particle (
165+ int * index_of_the_particle,
166+ double mass,
167+ double x, double y, double z,
168+ double vx, double vy, double vz,
169+ double radius, int code_index
170+ ){
166171 if (code_index < 0 || code_index >= (signed ) codes.size ()){
167172 *index_of_the_particle=0 ;
168173 return -10 ;
@@ -859,20 +864,20 @@ int _set_integrator(int value, int code_index){
859864 code->integrator = reb_simulation::REB_INTEGRATOR_JANUS;
860865 break ;
861866 case 9 :
862- code->integrator = reb_simulation::REB_INTEGRATOR_WHFAST512 ;
867+ code->integrator = reb_simulation::REB_INTEGRATOR_MERCURIUS ;
863868 break ;
864869 case 10 :
865870 code->integrator = reb_simulation::REB_INTEGRATOR_SABA;
866871 break ;
867872 case 11 :
868- code->integrator = reb_simulation::REB_INTEGRATOR_MERCURIUS;
869- break ;
870- case 12 :
871873 code->integrator = reb_simulation::REB_INTEGRATOR_EOS;
872874 break ;
873- case 13 :
875+ case 12 :
874876 code->integrator = reb_simulation::REB_INTEGRATOR_BS;
875877 break ;
878+ case 21 :
879+ code->integrator = reb_simulation::REB_INTEGRATOR_WHFAST512;
880+ break ;
876881 default :
877882 code->integrator = reb_simulation::REB_INTEGRATOR_NONE;
878883 return -1 ;
0 commit comments