/*Header File for numcip.c*/

static double dsqrarg;
#define DSQR(a) ((dsqrarg=(a)) == 0.0 ? 0.0 : dsqrarg*dsqrarg)

static double dmaxarg1,dmaxarg2;
#define DMAX(a,b) (dmaxarg1=(a),dmaxarg2=(b),(dmaxarg1) > (dmaxarg2) ?\
        (dmaxarg1) : (dmaxarg2))

static double dminarg1,dminarg2;
#define DMIN(a,b) (dminarg1=(a),dminarg2=(b),(dminarg1) < (dminarg2) ?\
        (dminarg1) : (dminarg2))

static int imaxarg1,imaxarg2;
#define IMAX(a,b) (imaxarg1=(a),imaxarg2=(b),(imaxarg1) > (imaxarg2) ?\
        (imaxarg1) : (imaxarg2))

static int iminarg1,iminarg2;
#define IMIN(a,b) (iminarg1=(a),iminarg2=(b),(iminarg1) < (iminarg2) ?\
        (iminarg1) : (iminarg2))

#define SIGN(a,b) ((b) >= 0.0 ? fabs(a) : -fabs(a))
void nrerror(char []);

int *ivector(int ,int );
void free_ivector(int *, int , int );

int **imatrix(int,int ,int ,int );
void free_imatrix(int **, int ,int ,int ,int );
double **dmatrix(int,int ,int ,int );
void free_dmatrix(double **, int ,int ,int ,int );
double *dvector(int ,int );
void free_dvector(double *, int , int );
double ***d3matrix(int,int,int,int,int,int);
void free_d3matrix(double ***,int,int,int,int,int,int);
double *bigvector(long ,long );
void free_bigvector(double *, long , long );
double **bmatrix(long,long ,long ,long );
void free_bmatrix(double **, long ,long ,long ,long );

unsigned char **ucmatrix(int ,int ,int,int );
void free_ucmatrix(unsigned char **, int ,int ,int ,int );

float *vector(int,int);
void free_vector(float *,int,int);
float **n2matrix(int,int,int,int);
void free_n2matrix(float **,int,int,int,int);
float ***n3matrix(int,int,int,int,int,int);
void free_n3matrix(float ***,int,int,int,int,int,int);
float ****n4matrix(int,int,int,int,int,int,int,int);
void free_n4matrix(float ****,int,int,int,int,int,int,int,int);