Merge pull request #120 from hannorein/interpolate_fix

Fixed array size (was too large) used in integrate_or_interpolate. ASSIST also now copies the particles initially, so that the radius, mass, hash, etc are maintained. Fixes #119

Author: hannorein

src/assist.c

@@ -482,8 +482,11 @@ void assist_integrate_or_interpolate(struct assist_extras* ax, double t){
     sim->exact_finish_time = 0;
 
     if (ax->current_state==NULL){
-        ax->current_state = malloc(sizeof(struct reb_particle)*6*sim->N);
-        ax->last_state = malloc(sizeof(struct reb_particle)*6*sim->N);
+        ax->current_state = malloc(sizeof(struct reb_particle)*sim->N);
+        ax->last_state = malloc(sizeof(struct reb_particle)*sim->N);
+        // Initialize new arrays with sim->particles (sets mass, radius, hash, if users use those) 
+        memcpy(ax->current_state, sim->particles, sizeof(struct reb_particle)*sim->N);
+        memcpy(ax->last_state, sim->particles, sizeof(struct reb_particle)*sim->N);
     }else{
         assist_swap_particles(sim);
     }