#include #include #include "main.h" #include "force.h" #include "integrator.h" int main(int argc, char **argv) { Data *data; data = (Data *)malloc( sizeof(Data) ); assert(NULL != data); if (parse_input(argc, argv, data)) return 1; if( initialize(data) ) return 1; leapfrog(data); cleanup(data); free(data); return 0; } int initialize(Data *data) { initial_conditions(data); data->sim_out = fopen(data->outfile,"w"); force_poten_eval(data); energy_eval(data); return 0; } int cleanup(Data *data) { free(data->pos); free(data->vel); free(data->f); free(data->mass); fclose(data->sim_out); return 0; }