#include <R.h>
#include <Rinternals.h>

/* #define TRANSIENT */

void bk_point(double *xp,
	     double *yp,
	     double *zp,
	     double *varp,
	     int *dop,
	      int *np,
	     double *lon,
	     double *lat,
	     double *z,
	     int *extrap,
	     int *n,
	     int *covtype,
	     double *covpar,
	     int *trend,
	     int *ntrend,
	     double *mupr,
	     int *ldmpr,
	     double *phipr,
	     int *ldphpr,
	     double *lonpr,
	     double *latpr,
	     int *npr,
	     int *typpr,
	     double *rsearch,
	     int *nsearch,
	     int *nsmin,
	     int *nsmax,
	     int *lwork,
	     int *mode,
	     double *lambda,
	     double *lambd0,
	     int *searchnb,
	     int *ierr,
	      int *glsmth);


void F77_NAME(bk)(double *xp,
		      double *yp,
		  int *dop,
		  int *inddop,
		  int *np,
		      double *lon,
		      double *lat,
		      double *z,
		      int *n,
		      int *covtype,
		      double *covpar,
		      double *covmat,
		      int *ldcov,
		      double *c0vec,
		  int *ldc0,
		      double *cov0,
		      int *trend,
		      int *ntrend,
		      double *mupr,
		      int *ldmpr,
		      double *phipr,
		      int *ldphpr,
		      double *muwrk,
		      double *phiwrk,
		      int *ldphwk,
		      double *lonpr,
		      double *latpr,
		      double *beta,
		      double *errbta,
		      double *covbta,
		      int *ldcvbt,
		      double *dev,
		      double *errdev,
		      double *cvsrnb,
		      double *zsrnb,
		      int *npr,
		      int *typpr,
		      double *rsearch,
		      int *nsearch,
		      int *nsmin,
		      int *nsmax,
		      double *fwork,
		      double *fwrk2,
		      int *ldfwrk,
		      double *f0work,
		  int *ldf0wk,
		      double *dist,
		      int *indsnb,
		      int *indsnw,
		      int *indsrt,
		      double *kwork,
		      int *ldkwrk,
		      double *rhswork,
		      double *fpwork,
		      double *fpfwork,
		      double *fpf0wrk,
		      double *chlup,
		      int *ldclup,
		      double *cminv,
		      int *ldcinv,
		      double *work,
		      int *lwork,
		      int *ipvt,
		      double *ferr,
		      double *berr,
		      int *ipiv,
		      int *iwork,
		      int *mode,
		      double *mu,
		  int *ldmu,
		  double *zp,
		  int *nap,
		      double *lambda,
		      int *ldlmbd,
		      double *lambd0,
		  double *varp,
		      int *searchnb,
		  int *usesbbt,
		      int *ierr,
		      int *glsmth);