# include "external.h"
# include "shell.h"
# include "dvomath.h"
# include "convert.h"
# include "display.h"

# ifndef DATA_H
# define DATA_H

/*** typedef structs used by math functions ***/
typedef struct {
  int   NLINES;
  int   Nlines;
  char **lines;
  char  *name;
} Queue;

typedef struct {
  char *name;
  int NWORDS;
  int Nwords;
  char **words;
  char **value;
} Page;

typedef struct {
  char *name;
  int NPAGES;
  int Npages;
  Page **pages;
  int *index;
  char **pageIDs;
} Book;

void InitData ();

/* in book.c */
void InitBooks ();
void InitBook (Book *book, char *name);
void FreeBook (Book *book);
Book *FindBook (char *name);
Book *GetBook (int where);
Book *CreateBook (char *name);
int DeleteBook (Book *book);
void ListBooks ();

/* in page.c */
void InitPage (Page *page, char *name);
void FreePage (Page *page);
Page *FindPage (Book *book, char *name);
Page *GetPage (Book *book, int where);
Page *GetPageRestricted (Book *book, int where, char *keyName, char *keyValue);
Page *CreatePage (Book *book, char *name);
int DeletePage (Book *book, Page *page);
void ListPages (Book *book);
void ListWords (Page *page);
int BookSetWord (Page *page, char *word, char *value);
char *BookGetWord (Page *page, char *word);

/* in queues.c */
void InitQueues ();
void ListQueues ();
Queue *FindQueue (char *name);
Queue *CreateQueue (char *name);
void PushQueue (Queue *queue, char *line);
void PushNamedQueue (char *name, char *line);
char *PopQueue (Queue *queue);
char *PopQueueMatch (Queue *queue, char *Key, char *value);
void PushQueueUnique (Queue *queue, char *line, char *Key);
void PushQueueReplace (Queue *queue, char *line, char *Key);
int InitQueue (Queue *queue);
int DeleteQueue (Queue *queue);
int PrintQueue (Queue *queue);

/* in sort.c */
void sort (double *value, int N);
void fsort (float *value, int N);
void sortpair (double *value1, double *value2, int N);
void fsortpair (float *X, float *Y, int N);
void sortthree (float *X, float *Y, float *Z, int N);
void sort_lists (float *X, float *Y, int *S, int N);
void dsort_lists (double *X, double *Y, int *S, int N);
void isort_pair (int *X, int *Y, int N);

/* in fft.c */
void fft (float *Data, int N, int isign);
void fftold (float *Data, int N, int isign);
void fftN (float *data, int *nn, int ndim, int isign);
int IsBinary (int N);
void fourn (float *data, int *nn, int ndim, int isign);

/* in gaussj.c */
int gaussj (double **a, int n, double **b, int m);
int fgaussj (float **a, int n, float **b, int m);

/* in svdcmp.c */
int svdcmp (float *a, float *w, float *v, int Nx, int Ny);

/* in spline.c */
void spline (float *x, float *y, int N, float *y2);
float splint (float *x, float *y, float *y2, int N, float X);

/* mrqmin.c */
float mrqcof (float *x, float *y, float *dy, int Npts, 
	      float *par, int Npar, float **ta, float **tb, 
	      float (funcs)(float, float *, int, float *));

float mrqmin (float *x, float *y, float *dy, int Npts, 
	      float *par, int Npar, 
	      float (funcs)(float, float *, int, float *), int VERBOSE);

float mrqinit (float *x, float *y, float *dy, int Npts, 
	       float *par, int Npar, 
	       float (funcs)(float, float *, int, float *), int VERBOSE);

float **mrqcovar (int Npar);

void mrqfree (int Npar);

/* mrq2dmin.c */
float mrq2dcof (float *x, float *t, float *y, float *dy, int Npts, 
		float *par, int Npar, float **ta, float **tb, 
		float (funcs)(float, float, float *, int, float *));

float mrq2dmin (float *x, float *t, float *y, float *dy, int Npts, 
		float *par, int Npar, 
		float (funcs)(float, float, float *, int, float *), int VERBOSE);

float mrq2dinit (float *x, float *t, float *y, float *dy, int Npts, 
		 float *par, int Npar, 
		 float (funcs)(float, float, float *, int, float *), int VERBOSE);

float mrq2dchi (float *x, float *t, float *y, float *dy, int Npts, 
		float *par, int Npar, 
		float (funcs)(float, float, float *, int, float *));

int mrq2dlimits (float *pmin, float *pmax, int Npar);

float **mrq2dcovar (int Npar);

void mrq2dfree (int Npar);

/* powell.c */
void powell (float *p, int Npar, float (func)() );

/* gaussian.c */
double gaussian (double x, double mean, double sigma);
void gauss_init (int Nbin);
double rnd_gauss (double mean, double sigma);

/* starfuncs.c */
double get_aperture_stats (Matrix *matrix, int X, int Y, int Npix, int Nborder, double max);
int set_rough_radii (double Ra, double Ri, double Ro);
int get_rough_star (float *data, int Nx, int Ny, int x, int y, float *xc, float *yc, float *sx, float *sy, float *sxy, float *zs, float *zp, float *sk);

/* precess.c */
double BtoJ (double in_epoch);
double get_epoch (char *in_epoch, char mode);

/* graphtools.c */
void          SetLimits             PROTO((Vector *xvec, Vector *yvec, Graphdata *graphmode));
void          SetLimitsRaw          PROTO((float *xvec, float *yvec, int Npts, Graphdata *graphmode));
void          ApplyLimits           PROTO((int Xgraph, Graphdata *graphmode, int apply));
int           style_args            PROTO((Graphdata *graphmode, int *argc, char **argv, int Nforce));

int read_table_vectors (int argc, char **argv, char *extname);
# endif

