#include <float.h>
#include <fstream.h>
#include <iostream.h>
#include <math.h>
#include <unistd.h>
#include <stdio.h>
#include <stdlib.h>
#include <GL/glut.h>

#define LENGTH_NAME 80
#define size_W 300
#define size_H 300
#define F77NAME(x) x## _
#define F77_NAME(name) name## __

extern "C" {
  void readinput(int, char**, int*, int*, int*, int*, double*, double*, int*);
  void writeout(int, int, double*, char*);
  void reinitiate(int, int, int, double*);
  void Plot_Swarming(int, int, int, double*, double*, double*, double*);
  void Plot_RadialDensity(int, double, double*);
  void AddParticle(void);
  void gtDrawAxis(double, double, double, double);
  void gtWrite(const char*, float, float);
  void F77NAME(initsl)(const int &ncells,
      const int &fx,const int &lx,const int &ifirst,const int &ilast,
    double *xold);
  void F77NAME(method)(const double &t, const double &dt,
      const int &fm,const int &lm, const int &ifirst,const int &ilast,
      const double *xold, double *xnew,
      double *xm, double *vm, double *km);
  void F77NAME(method2)(const double &t, const double &dt,
      const int &fm,const int &lm, const int &ifirst,const int &ilast,
      const double *xold, double *xnew,
      double *xm, double *vm, double *km);
  void F77NAME(fourstep)(const double &t, double &dt, const int &fm,
      const int &lm, const int &ifirst, const int&ilast, const double *xold,
      double *xnew, double *km, const double *km1, const double *km2, 
      const double *km3, double *Gforce);
  void F77NAME(group)(const int &fm,const int &lm, const int &ifirst,
      const int &ilast, const double *xnew, double &polar, double &momen,
      const int &rhoPoints, double *rho, double &dr, int &cw, int &ccw,
      double *Gforce);
}