#ifndef GENLIMIT

typedef struct {
  float e,b;
}EB;

typedef enum {
  flat=10,
  corr=20
} PRIOR;

double cspdf(double s,double norm,
             int nchan,int nens,const int nobs[],const EB* ens,PRIOR prior);

double csint(double s0,int nchan,int nens,const int nobs[],const EB* ens,
	     int* ngl,double xgl[],double lwgl[],
	     PRIOR prior,double* uncertainty);

void csint2(double s1,double s2,
	    int nchan,int nens,const int nobs[],const EB* ens,
	    int* ngl,double xgl[],double lwgl[],PRIOR prior,
	    double* int1,double* int2,
	    double* v11,double* v12,double* v22);

void csint2cut(double s1,double s2,double shi,
	       int nchan,int nens,const int nobs[],const EB* ens,
	       int* ngl,double xgl[],double lwgl[],PRIOR prior,
	       double* int1,double* int2,
	       double* v11,double* v12,double* v22);

double cslimit(double beta,int nchan,int nens,const int nobs[],const EB* ens,
               int* ngl,double xgl[],double lwgl[],
	       PRIOR prior,double* uncertainty);

double cscutlimit(double beta,double smax,
		  int nchan,int nens,const int nobs[],const EB* ens,
		  int* ngl,double xgl[],double lwgl[],
		  PRIOR prior,double* uncertainty);

double galim(double beta,int nchan,int nens,const int nobs[],const EB* ens);

void gausslaguerre(double x[],double lw[],int n,double alpha);

#define GENLIMIT 1
#endif