#include <assert.h>
#include <math.h>
#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;
}