File Coverage

Levmar.xs

Criterion	Covered	Total	%
statement	57	111	51.3
branch	100	524	19.0
condition			n/a
subroutine			n/a
pod			n/a
total	157	635	24.7

#line 453 "lib/PDL/PP.pm"
/*
 * THIS FILE WAS GENERATED BY PDL::PP from levmar.pd! Do not modify!
 */

#define PDL_FREE_CODE(trans, destroy, comp_free_code, ntpriv_free_code) \
    if (destroy) { \
	comp_free_code \
    } \
    if ((trans)->dims_redone) { \
	ntpriv_free_code \
    }

#include "EXTERN.h"
#include "perl.h"
#include "XSUB.h"
#include "pdl.h"
#include "pdlcore.h"
#define PDL PDL_Fit_Levmar
extern Core* PDL; /* Structure hold core C functions */
#line 23 "Levmar.xs"

#include
#include
#include
#include
#include "pdlperlfunc.h"
#include "levmar.h"
pdl_error pdl_run_levmar_der_(pdl  *p,pdl  *x,pdl  *t,pdl  *iopts,pdl  *opts,pdl  *covar,pdl  *returnval,pdl  *pout,pdl  *info,IV  funcn,IV  sfuncn,IV  jacn,IV  sjacn,IV  indat,int  want_covar);
pdl_error pdl_run_levmar_der_lb_ub(pdl  *p,pdl  *x,pdl  *t,pdl  *lb,pdl  *ub,pdl  *iopts,pdl  *opts,pdl  *covar,pdl  *returnval,pdl  *pout,pdl  *info,IV  funcn,IV  sfuncn,IV  jacn,IV  sjacn,IV  indat,int  want_covar);
pdl_error pdl_run_levmar_der_ub(pdl  *p,pdl  *x,pdl  *t,pdl  *ub,pdl  *iopts,pdl  *opts,pdl  *covar,pdl  *returnval,pdl  *pout,pdl  *info,IV  funcn,IV  sfuncn,IV  jacn,IV  sjacn,IV  indat,int  want_covar);
pdl_error pdl_run_levmar_der_lb(pdl  *p,pdl  *x,pdl  *t,pdl  *lb,pdl  *iopts,pdl  *opts,pdl  *covar,pdl  *returnval,pdl  *pout,pdl  *info,IV  funcn,IV  sfuncn,IV  jacn,IV  sjacn,IV  indat,int  want_covar);
pdl_error pdl_run_levmar_diff_(pdl  *p,pdl  *x,pdl  *t,pdl  *iopts,pdl  *opts,pdl  *covar,pdl  *returnval,pdl  *pout,pdl  *info,IV  funcn,IV  sfuncn,IV  indat,int  want_covar);
pdl_error pdl_run_levmar_diff_lb_ub(pdl  *p,pdl  *x,pdl  *t,pdl  *lb,pdl  *ub,pdl  *iopts,pdl  *opts,pdl  *covar,pdl  *returnval,pdl  *pout,pdl  *info,IV  funcn,IV  sfuncn,IV  indat,int  want_covar);
pdl_error pdl_run_levmar_diff_ub(pdl  *p,pdl  *x,pdl  *t,pdl  *ub,pdl  *iopts,pdl  *opts,pdl  *covar,pdl  *returnval,pdl  *pout,pdl  *info,IV  funcn,IV  sfuncn,IV  indat,int  want_covar);
pdl_error pdl_run_levmar_diff_lb(pdl  *p,pdl  *x,pdl  *t,pdl  *lb,pdl  *iopts,pdl  *opts,pdl  *covar,pdl  *returnval,pdl  *pout,pdl  *info,IV  funcn,IV  sfuncn,IV  indat,int  want_covar);
pdl_error pdl_run__levmar_chkjac(pdl  *p,pdl  *t,pdl  *err,IV  func,IV  sfunc,IV  jac,IV  sjac,IV  indat);
pdl_error pdl_run__levmar_chkjac_no_t(pdl  *p,pdl  *err,IV  func,IV  sfunc,IV  jac,IV  sjac,int  N,IV  indat);

#line 449 "lib/PDL/PP.pm"
#include "pdlperl.h"
#line 44 "Levmar.xs"

Core* PDL = NULL; /* Structure hold core C functions */

MODULE = PDL::Fit::Levmar PACKAGE = PDL PREFIX=pdl_run_

PROTOTYPES: DISABLE


void
pdl_run_levmar_der_(p, x, t, iopts, opts, covar=covar, returnval=returnval, pout=pout, info=info, funcn=funcn, sfuncn=sfuncn, jacn=jacn, sjacn=sjacn, indat=indat, want_covar=want_covar)
  SV *covar_SV = (items == 15) ? ST(5) : NULL;
  SV *returnval_SV = (items == 15) ? ST(6) : NULL;
  SV *pout_SV = (items == 15) ? ST(7) : NULL;
  SV *info_SV = (items == 15) ? ST(8) : NULL;
  SV *funcn_SV = (items == 15) ? ST(9) : ST(5);
  SV *sfuncn_SV = (items == 15) ? ST(10) : ST(6);
  SV *jacn_SV = (items == 15) ? ST(11) : ST(7);
  SV *sjacn_SV = (items == 15) ? ST(12) : ST(8);
  SV *indat_SV = (items == 15) ? ST(13) : ST(9);
  SV *want_covar_SV = (items == 15) ? ST(14) : ST(10);
 PREINIT:
  PDL_XS_PREAMBLE((items == 15) ? 0 : 4);
  if (!((items == 11) || (items == 15)))
    croak("Usage: PDL::levmar_der_(p,x,t,iopts,opts,[covar],[returnval],[pout],[info],funcn,sfuncn,jacn,sjacn,indat,want_covar) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *p
  pdl *x
  pdl *t
  pdl *iopts
  pdl *opts
  pdl *covar=NO_INIT
  pdl *returnval=NO_INIT
  pdl *pout=NO_INIT
  pdl *info=NO_INIT
  IV  funcn; { funcn = (IV)SvIV(funcn_SV); }
  IV  sfuncn; { sfuncn = (IV)SvIV(sfuncn_SV); }
  IV  jacn; { jacn = (IV)SvIV(jacn_SV); }
  IV  sjacn; { sjacn = (IV)SvIV(sjacn_SV); }
  IV  indat; { indat = (IV)SvIV(indat_SV); }
  int  want_covar; { want_covar = (int)SvIV(want_covar_SV); }
 PPCODE:
  if (covar_SV) { covar = PDL_CORE_(SvPDLV)(covar_SV); } else covar = PDL_XS_PERLINIT_initsv(covar_SV);
  if (returnval_SV) { returnval = PDL_CORE_(SvPDLV)(returnval_SV); } else returnval = PDL_XS_PERLINIT_initsv(returnval_SV);
  if (pout_SV) { pout = PDL_CORE_(SvPDLV)(pout_SV); } else pout = PDL_XS_PERLINIT_initsv(pout_SV);
  if (info_SV) { info = PDL_CORE_(SvPDLV)(info_SV); } else info = PDL_XS_PERLINIT_initsv(info_SV);
  PDL->barf_if_error(pdl_run_levmar_der_(p,x,t,iopts,opts,covar,returnval,pout,info,funcn,sfuncn,jacn,sjacn,indat,want_covar));
  PDL_XS_RETURN(ST(0) = covar_SV;ST(1) = returnval_SV;ST(2) = pout_SV;ST(3) = info_SV)

void
pdl_run_levmar_der_lb_ub(p, x, t, lb, ub, iopts, opts, covar=covar, returnval=returnval, pout=pout, info=info, funcn=funcn, sfuncn=sfuncn, jacn=jacn, sjacn=sjacn, indat=indat, want_covar=want_covar)
  SV *covar_SV = (items == 17) ? ST(7) : NULL;
  SV *returnval_SV = (items == 17) ? ST(8) : NULL;
  SV *pout_SV = (items == 17) ? ST(9) : NULL;
  SV *info_SV = (items == 17) ? ST(10) : NULL;
  SV *funcn_SV = (items == 17) ? ST(11) : ST(7);
  SV *sfuncn_SV = (items == 17) ? ST(12) : ST(8);
  SV *jacn_SV = (items == 17) ? ST(13) : ST(9);
  SV *sjacn_SV = (items == 17) ? ST(14) : ST(10);
  SV *indat_SV = (items == 17) ? ST(15) : ST(11);
  SV *want_covar_SV = (items == 17) ? ST(16) : ST(12);
 PREINIT:
  PDL_XS_PREAMBLE((items == 17) ? 0 : 4);
  if (!((items == 13) || (items == 17)))
    croak("Usage: PDL::levmar_der_lb_ub(p,x,t,lb,ub,iopts,opts,[covar],[returnval],[pout],[info],funcn,sfuncn,jacn,sjacn,indat,want_covar) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *p
  pdl *x
  pdl *t
  pdl *lb
  pdl *ub
  pdl *iopts
  pdl *opts
  pdl *covar=NO_INIT
  pdl *returnval=NO_INIT
  pdl *pout=NO_INIT
  pdl *info=NO_INIT
  IV  funcn; { funcn = (IV)SvIV(funcn_SV); }
  IV  sfuncn; { sfuncn = (IV)SvIV(sfuncn_SV); }
  IV  jacn; { jacn = (IV)SvIV(jacn_SV); }
  IV  sjacn; { sjacn = (IV)SvIV(sjacn_SV); }
  IV  indat; { indat = (IV)SvIV(indat_SV); }
  int  want_covar; { want_covar = (int)SvIV(want_covar_SV); }
 PPCODE:
  if (covar_SV) { covar = PDL_CORE_(SvPDLV)(covar_SV); } else covar = PDL_XS_PERLINIT_initsv(covar_SV);
  if (returnval_SV) { returnval = PDL_CORE_(SvPDLV)(returnval_SV); } else returnval = PDL_XS_PERLINIT_initsv(returnval_SV);
  if (pout_SV) { pout = PDL_CORE_(SvPDLV)(pout_SV); } else pout = PDL_XS_PERLINIT_initsv(pout_SV);
  if (info_SV) { info = PDL_CORE_(SvPDLV)(info_SV); } else info = PDL_XS_PERLINIT_initsv(info_SV);
  PDL->barf_if_error(pdl_run_levmar_der_lb_ub(p,x,t,lb,ub,iopts,opts,covar,returnval,pout,info,funcn,sfuncn,jacn,sjacn,indat,want_covar));
  PDL_XS_RETURN(ST(0) = covar_SV;ST(1) = returnval_SV;ST(2) = pout_SV;ST(3) = info_SV)

void
pdl_run_levmar_der_ub(p, x, t, ub, iopts, opts, covar=covar, returnval=returnval, pout=pout, info=info, funcn=funcn, sfuncn=sfuncn, jacn=jacn, sjacn=sjacn, indat=indat, want_covar=want_covar)
  SV *covar_SV = (items == 16) ? ST(6) : NULL;
  SV *returnval_SV = (items == 16) ? ST(7) : NULL;
  SV *pout_SV = (items == 16) ? ST(8) : NULL;
  SV *info_SV = (items == 16) ? ST(9) : NULL;
  SV *funcn_SV = (items == 16) ? ST(10) : ST(6);
  SV *sfuncn_SV = (items == 16) ? ST(11) : ST(7);
  SV *jacn_SV = (items == 16) ? ST(12) : ST(8);
  SV *sjacn_SV = (items == 16) ? ST(13) : ST(9);
  SV *indat_SV = (items == 16) ? ST(14) : ST(10);
  SV *want_covar_SV = (items == 16) ? ST(15) : ST(11);
 PREINIT:
  PDL_XS_PREAMBLE((items == 16) ? 0 : 4);
  if (!((items == 12) || (items == 16)))
    croak("Usage: PDL::levmar_der_ub(p,x,t,ub,iopts,opts,[covar],[returnval],[pout],[info],funcn,sfuncn,jacn,sjacn,indat,want_covar) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *p
  pdl *x
  pdl *t
  pdl *ub
  pdl *iopts
  pdl *opts
  pdl *covar=NO_INIT
  pdl *returnval=NO_INIT
  pdl *pout=NO_INIT
  pdl *info=NO_INIT
  IV  funcn; { funcn = (IV)SvIV(funcn_SV); }
  IV  sfuncn; { sfuncn = (IV)SvIV(sfuncn_SV); }
  IV  jacn; { jacn = (IV)SvIV(jacn_SV); }
  IV  sjacn; { sjacn = (IV)SvIV(sjacn_SV); }
  IV  indat; { indat = (IV)SvIV(indat_SV); }
  int  want_covar; { want_covar = (int)SvIV(want_covar_SV); }
 PPCODE:
  if (covar_SV) { covar = PDL_CORE_(SvPDLV)(covar_SV); } else covar = PDL_XS_PERLINIT_initsv(covar_SV);
  if (returnval_SV) { returnval = PDL_CORE_(SvPDLV)(returnval_SV); } else returnval = PDL_XS_PERLINIT_initsv(returnval_SV);
  if (pout_SV) { pout = PDL_CORE_(SvPDLV)(pout_SV); } else pout = PDL_XS_PERLINIT_initsv(pout_SV);
  if (info_SV) { info = PDL_CORE_(SvPDLV)(info_SV); } else info = PDL_XS_PERLINIT_initsv(info_SV);
  PDL->barf_if_error(pdl_run_levmar_der_ub(p,x,t,ub,iopts,opts,covar,returnval,pout,info,funcn,sfuncn,jacn,sjacn,indat,want_covar));
  PDL_XS_RETURN(ST(0) = covar_SV;ST(1) = returnval_SV;ST(2) = pout_SV;ST(3) = info_SV)

void
pdl_run_levmar_der_lb(p, x, t, lb, iopts, opts, covar=covar, returnval=returnval, pout=pout, info=info, funcn=funcn, sfuncn=sfuncn, jacn=jacn, sjacn=sjacn, indat=indat, want_covar=want_covar)
  SV *covar_SV = (items == 16) ? ST(6) : NULL;
  SV *returnval_SV = (items == 16) ? ST(7) : NULL;
  SV *pout_SV = (items == 16) ? ST(8) : NULL;
  SV *info_SV = (items == 16) ? ST(9) : NULL;
  SV *funcn_SV = (items == 16) ? ST(10) : ST(6);
  SV *sfuncn_SV = (items == 16) ? ST(11) : ST(7);
  SV *jacn_SV = (items == 16) ? ST(12) : ST(8);
  SV *sjacn_SV = (items == 16) ? ST(13) : ST(9);
  SV *indat_SV = (items == 16) ? ST(14) : ST(10);
  SV *want_covar_SV = (items == 16) ? ST(15) : ST(11);
 PREINIT:
  PDL_XS_PREAMBLE((items == 16) ? 0 : 4);
  if (!((items == 12) || (items == 16)))
    croak("Usage: PDL::levmar_der_lb(p,x,t,lb,iopts,opts,[covar],[returnval],[pout],[info],funcn,sfuncn,jacn,sjacn,indat,want_covar) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *p
  pdl *x
  pdl *t
  pdl *lb
  pdl *iopts
  pdl *opts
  pdl *covar=NO_INIT
  pdl *returnval=NO_INIT
  pdl *pout=NO_INIT
  pdl *info=NO_INIT
  IV  funcn; { funcn = (IV)SvIV(funcn_SV); }
  IV  sfuncn; { sfuncn = (IV)SvIV(sfuncn_SV); }
  IV  jacn; { jacn = (IV)SvIV(jacn_SV); }
  IV  sjacn; { sjacn = (IV)SvIV(sjacn_SV); }
  IV  indat; { indat = (IV)SvIV(indat_SV); }
  int  want_covar; { want_covar = (int)SvIV(want_covar_SV); }
 PPCODE:
  if (covar_SV) { covar = PDL_CORE_(SvPDLV)(covar_SV); } else covar = PDL_XS_PERLINIT_initsv(covar_SV);
  if (returnval_SV) { returnval = PDL_CORE_(SvPDLV)(returnval_SV); } else returnval = PDL_XS_PERLINIT_initsv(returnval_SV);
  if (pout_SV) { pout = PDL_CORE_(SvPDLV)(pout_SV); } else pout = PDL_XS_PERLINIT_initsv(pout_SV);
  if (info_SV) { info = PDL_CORE_(SvPDLV)(info_SV); } else info = PDL_XS_PERLINIT_initsv(info_SV);
  PDL->barf_if_error(pdl_run_levmar_der_lb(p,x,t,lb,iopts,opts,covar,returnval,pout,info,funcn,sfuncn,jacn,sjacn,indat,want_covar));
  PDL_XS_RETURN(ST(0) = covar_SV;ST(1) = returnval_SV;ST(2) = pout_SV;ST(3) = info_SV)

void
pdl_run_levmar_diff_(p, x, t, iopts, opts, covar=covar, returnval=returnval, pout=pout, info=info, funcn=funcn, sfuncn=sfuncn, indat=indat, want_covar=want_covar)
  SV *covar_SV = (items == 13) ? ST(5) : NULL;
  SV *returnval_SV = (items == 13) ? ST(6) : NULL;
  SV *pout_SV = (items == 13) ? ST(7) : NULL;
  SV *info_SV = (items == 13) ? ST(8) : NULL;
  SV *funcn_SV = (items == 13) ? ST(9) : ST(5);
  SV *sfuncn_SV = (items == 13) ? ST(10) : ST(6);
  SV *indat_SV = (items == 13) ? ST(11) : ST(7);
  SV *want_covar_SV = (items == 13) ? ST(12) : ST(8);
 PREINIT:
  PDL_XS_PREAMBLE((items == 13) ? 0 : 4);
  if (!((items == 13) || (items == 9)))
    croak("Usage: PDL::levmar_diff_(p,x,t,iopts,opts,[covar],[returnval],[pout],[info],funcn,sfuncn,indat,want_covar) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *p
  pdl *x
  pdl *t
  pdl *iopts
  pdl *opts
  pdl *covar=NO_INIT
  pdl *returnval=NO_INIT
  pdl *pout=NO_INIT
  pdl *info=NO_INIT
  IV  funcn; { funcn = (IV)SvIV(funcn_SV); }
  IV  sfuncn; { sfuncn = (IV)SvIV(sfuncn_SV); }
  IV  indat; { indat = (IV)SvIV(indat_SV); }
  int  want_covar; { want_covar = (int)SvIV(want_covar_SV); }
 PPCODE:
  if (covar_SV) { covar = PDL_CORE_(SvPDLV)(covar_SV); } else covar = PDL_XS_PERLINIT_initsv(covar_SV);
  if (returnval_SV) { returnval = PDL_CORE_(SvPDLV)(returnval_SV); } else returnval = PDL_XS_PERLINIT_initsv(returnval_SV);
  if (pout_SV) { pout = PDL_CORE_(SvPDLV)(pout_SV); } else pout = PDL_XS_PERLINIT_initsv(pout_SV);
  if (info_SV) { info = PDL_CORE_(SvPDLV)(info_SV); } else info = PDL_XS_PERLINIT_initsv(info_SV);
  PDL->barf_if_error(pdl_run_levmar_diff_(p,x,t,iopts,opts,covar,returnval,pout,info,funcn,sfuncn,indat,want_covar));
  PDL_XS_RETURN(ST(0) = covar_SV;ST(1) = returnval_SV;ST(2) = pout_SV;ST(3) = info_SV)

void
pdl_run_levmar_diff_lb_ub(p, x, t, lb, ub, iopts, opts, covar=covar, returnval=returnval, pout=pout, info=info, funcn=funcn, sfuncn=sfuncn, indat=indat, want_covar=want_covar)
  SV *covar_SV = (items == 15) ? ST(7) : NULL;
  SV *returnval_SV = (items == 15) ? ST(8) : NULL;
  SV *pout_SV = (items == 15) ? ST(9) : NULL;
  SV *info_SV = (items == 15) ? ST(10) : NULL;
  SV *funcn_SV = (items == 15) ? ST(11) : ST(7);
  SV *sfuncn_SV = (items == 15) ? ST(12) : ST(8);
  SV *indat_SV = (items == 15) ? ST(13) : ST(9);
  SV *want_covar_SV = (items == 15) ? ST(14) : ST(10);
 PREINIT:
  PDL_XS_PREAMBLE((items == 15) ? 0 : 4);
  if (!((items == 11) || (items == 15)))
    croak("Usage: PDL::levmar_diff_lb_ub(p,x,t,lb,ub,iopts,opts,[covar],[returnval],[pout],[info],funcn,sfuncn,indat,want_covar) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *p
  pdl *x
  pdl *t
  pdl *lb
  pdl *ub
  pdl *iopts
  pdl *opts
  pdl *covar=NO_INIT
  pdl *returnval=NO_INIT
  pdl *pout=NO_INIT
  pdl *info=NO_INIT
  IV  funcn; { funcn = (IV)SvIV(funcn_SV); }
  IV  sfuncn; { sfuncn = (IV)SvIV(sfuncn_SV); }
  IV  indat; { indat = (IV)SvIV(indat_SV); }
  int  want_covar; { want_covar = (int)SvIV(want_covar_SV); }
 PPCODE:
  if (covar_SV) { covar = PDL_CORE_(SvPDLV)(covar_SV); } else covar = PDL_XS_PERLINIT_initsv(covar_SV);
  if (returnval_SV) { returnval = PDL_CORE_(SvPDLV)(returnval_SV); } else returnval = PDL_XS_PERLINIT_initsv(returnval_SV);
  if (pout_SV) { pout = PDL_CORE_(SvPDLV)(pout_SV); } else pout = PDL_XS_PERLINIT_initsv(pout_SV);
  if (info_SV) { info = PDL_CORE_(SvPDLV)(info_SV); } else info = PDL_XS_PERLINIT_initsv(info_SV);
  PDL->barf_if_error(pdl_run_levmar_diff_lb_ub(p,x,t,lb,ub,iopts,opts,covar,returnval,pout,info,funcn,sfuncn,indat,want_covar));
  PDL_XS_RETURN(ST(0) = covar_SV;ST(1) = returnval_SV;ST(2) = pout_SV;ST(3) = info_SV)

void
pdl_run_levmar_diff_ub(p, x, t, ub, iopts, opts, covar=covar, returnval=returnval, pout=pout, info=info, funcn=funcn, sfuncn=sfuncn, indat=indat, want_covar=want_covar)
  SV *covar_SV = (items == 14) ? ST(6) : NULL;
  SV *returnval_SV = (items == 14) ? ST(7) : NULL;
  SV *pout_SV = (items == 14) ? ST(8) : NULL;
  SV *info_SV = (items == 14) ? ST(9) : NULL;
  SV *funcn_SV = (items == 14) ? ST(10) : ST(6);
  SV *sfuncn_SV = (items == 14) ? ST(11) : ST(7);
  SV *indat_SV = (items == 14) ? ST(12) : ST(8);
  SV *want_covar_SV = (items == 14) ? ST(13) : ST(9);
 PREINIT:
  PDL_XS_PREAMBLE((items == 14) ? 0 : 4);
  if (!((items == 10) || (items == 14)))
    croak("Usage: PDL::levmar_diff_ub(p,x,t,ub,iopts,opts,[covar],[returnval],[pout],[info],funcn,sfuncn,indat,want_covar) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *p
  pdl *x
  pdl *t
  pdl *ub
  pdl *iopts
  pdl *opts
  pdl *covar=NO_INIT
  pdl *returnval=NO_INIT
  pdl *pout=NO_INIT
  pdl *info=NO_INIT
  IV  funcn; { funcn = (IV)SvIV(funcn_SV); }
  IV  sfuncn; { sfuncn = (IV)SvIV(sfuncn_SV); }
  IV  indat; { indat = (IV)SvIV(indat_SV); }
  int  want_covar; { want_covar = (int)SvIV(want_covar_SV); }
 PPCODE:
  if (covar_SV) { covar = PDL_CORE_(SvPDLV)(covar_SV); } else covar = PDL_XS_PERLINIT_initsv(covar_SV);
  if (returnval_SV) { returnval = PDL_CORE_(SvPDLV)(returnval_SV); } else returnval = PDL_XS_PERLINIT_initsv(returnval_SV);
  if (pout_SV) { pout = PDL_CORE_(SvPDLV)(pout_SV); } else pout = PDL_XS_PERLINIT_initsv(pout_SV);
  if (info_SV) { info = PDL_CORE_(SvPDLV)(info_SV); } else info = PDL_XS_PERLINIT_initsv(info_SV);
  PDL->barf_if_error(pdl_run_levmar_diff_ub(p,x,t,ub,iopts,opts,covar,returnval,pout,info,funcn,sfuncn,indat,want_covar));
  PDL_XS_RETURN(ST(0) = covar_SV;ST(1) = returnval_SV;ST(2) = pout_SV;ST(3) = info_SV)

void
pdl_run_levmar_diff_lb(p, x, t, lb, iopts, opts, covar=covar, returnval=returnval, pout=pout, info=info, funcn=funcn, sfuncn=sfuncn, indat=indat, want_covar=want_covar)
  SV *covar_SV = (items == 14) ? ST(6) : NULL;
  SV *returnval_SV = (items == 14) ? ST(7) : NULL;
  SV *pout_SV = (items == 14) ? ST(8) : NULL;
  SV *info_SV = (items == 14) ? ST(9) : NULL;
  SV *funcn_SV = (items == 14) ? ST(10) : ST(6);
  SV *sfuncn_SV = (items == 14) ? ST(11) : ST(7);
  SV *indat_SV = (items == 14) ? ST(12) : ST(8);
  SV *want_covar_SV = (items == 14) ? ST(13) : ST(9);
 PREINIT:
  PDL_XS_PREAMBLE((items == 14) ? 0 : 4);
  if (!((items == 10) || (items == 14)))
    croak("Usage: PDL::levmar_diff_lb(p,x,t,lb,iopts,opts,[covar],[returnval],[pout],[info],funcn,sfuncn,indat,want_covar) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *p
  pdl *x
  pdl *t
  pdl *lb
  pdl *iopts
  pdl *opts
  pdl *covar=NO_INIT
  pdl *returnval=NO_INIT
  pdl *pout=NO_INIT
  pdl *info=NO_INIT
  IV  funcn; { funcn = (IV)SvIV(funcn_SV); }
  IV  sfuncn; { sfuncn = (IV)SvIV(sfuncn_SV); }
  IV  indat; { indat = (IV)SvIV(indat_SV); }
  int  want_covar; { want_covar = (int)SvIV(want_covar_SV); }
 PPCODE:
  if (covar_SV) { covar = PDL_CORE_(SvPDLV)(covar_SV); } else covar = PDL_XS_PERLINIT_initsv(covar_SV);
  if (returnval_SV) { returnval = PDL_CORE_(SvPDLV)(returnval_SV); } else returnval = PDL_XS_PERLINIT_initsv(returnval_SV);
  if (pout_SV) { pout = PDL_CORE_(SvPDLV)(pout_SV); } else pout = PDL_XS_PERLINIT_initsv(pout_SV);
  if (info_SV) { info = PDL_CORE_(SvPDLV)(info_SV); } else info = PDL_XS_PERLINIT_initsv(info_SV);
  PDL->barf_if_error(pdl_run_levmar_diff_lb(p,x,t,lb,iopts,opts,covar,returnval,pout,info,funcn,sfuncn,indat,want_covar));
  PDL_XS_RETURN(ST(0) = covar_SV;ST(1) = returnval_SV;ST(2) = pout_SV;ST(3) = info_SV)

void
pdl_run__levmar_chkjac(p, t, err=err, func=func, sfunc=sfunc, jac=jac, sjac=sjac, indat=indat)
  SV *err_SV = (items == 8) ? ST(2) : NULL;
  SV *func_SV = (items == 8) ? ST(3) : ST(2);
  SV *sfunc_SV = (items == 8) ? ST(4) : ST(3);
  SV *jac_SV = (items == 8) ? ST(5) : ST(4);
  SV *sjac_SV = (items == 8) ? ST(6) : ST(5);
  SV *indat_SV = (items == 8) ? ST(7) : ST(6);
 PREINIT:
  PDL_XS_PREAMBLE((items == 8) ? 0 : 1);
  if (!((items == 7) || (items == 8)))
    croak("Usage: PDL::_levmar_chkjac(p,t,[err],func,sfunc,jac,sjac,indat) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *p
  pdl *t
  pdl *err=NO_INIT
  IV  func; { func = (IV)SvIV(func_SV); }
  IV  sfunc; { sfunc = (IV)SvIV(sfunc_SV); }
  IV  jac; { jac = (IV)SvIV(jac_SV); }
  IV  sjac; { sjac = (IV)SvIV(sjac_SV); }
  IV  indat; { indat = (IV)SvIV(indat_SV); }
 PPCODE:
  if (err_SV) { err = PDL_CORE_(SvPDLV)(err_SV); } else err = PDL_XS_PERLINIT_initsv(err_SV);
  PDL->barf_if_error(pdl_run__levmar_chkjac(p,t,err,func,sfunc,jac,sjac,indat));
  PDL_XS_RETURN(ST(0) = err_SV)

void
pdl_run__levmar_chkjac_no_t(p, err=err, func=func, sfunc=sfunc, jac=jac, sjac=sjac, N=N, indat=indat)
  SV *err_SV = (items == 8) ? ST(1) : NULL;
  SV *func_SV = (items == 8) ? ST(2) : ST(1);
  SV *sfunc_SV = (items == 8) ? ST(3) : ST(2);
  SV *jac_SV = (items == 8) ? ST(4) : ST(3);
  SV *sjac_SV = (items == 8) ? ST(5) : ST(4);
  SV *N_SV = (items == 8) ? ST(6) : ST(5);
  SV *indat_SV = (items == 8) ? ST(7) : ST(6);
 PREINIT:
  PDL_XS_PREAMBLE((items == 8) ? 0 : 1);
  if (!((items == 7) || (items == 8)))
    croak("Usage: PDL::_levmar_chkjac_no_t(p,[err],func,sfunc,jac,sjac,N,indat) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *p
  pdl *err=NO_INIT
  IV  func; { func = (IV)SvIV(func_SV); }
  IV  sfunc; { sfunc = (IV)SvIV(sfunc_SV); }
  IV  jac; { jac = (IV)SvIV(jac_SV); }
  IV  sjac; { sjac = (IV)SvIV(sjac_SV); }
  int  N; { N = (int)SvIV(N_SV); }
  IV  indat; { indat = (IV)SvIV(indat_SV); }
 PPCODE:
  if (err_SV) { err = PDL_CORE_(SvPDLV)(err_SV); } else err = PDL_XS_PERLINIT_initsv(err_SV);
  PDL->barf_if_error(pdl_run__levmar_chkjac_no_t(p,err,func,sfunc,jac,sjac,N,indat));
  PDL_XS_RETURN(ST(0) = err_SV)

MODULE = PDL::Fit::Levmar PACKAGE = PDL::Fit::Levmar




void *
DFP_create()

     CODE:
     DFP * dat;
     dat = (DFP *)malloc(sizeof( DFP ));
     if ( NULL == dat ) {
       fprintf(stderr, "Can't allocate storage for dat in DFP_create\n");
       exit(1);
     }
     RETVAL = dat;
     OUTPUT:
     RETVAL

void
DFP_free( dat )
    IV  dat

    CODE:
    if ( (DFP *) dat) free( (DFP *) dat);


void
DFP_set_perl_funcs ( data,  perl_fit_func,  perl_jac_func )
     IV data
     SV* perl_fit_func
     SV* perl_jac_func

     CODE:
     DFP *dat = (DFP *) data;
     if ( dat == NULL ) {
        fprintf(stderr, "DFP_set_perl_funcs got null struct\n");
        exit(1);
     }
     dat->perl_fit_func = perl_fit_func;
     dat->perl_jac_func = perl_jac_func;


double
get_dbl_max(  )
     CODE:
     RETVAL = DBL_MAX;
     OUTPUT:
     RETVAL

double
LM_INIT_MU(  )
     CODE:
     RETVAL = LM_INIT_MU;
     OUTPUT:
     RETVAL

double
LM_STOP_THRESH(  )
     CODE:
     RETVAL = LM_STOP_THRESH;
     OUTPUT:
     RETVAL

double
LM_DIFF_DELTA(  )
     CODE:
     RETVAL = LM_DIFF_DELTA;
     OUTPUT:
     RETVAL


=pod


=begin comment

  The next two routines just get pointers to functions from
  the C code and return them as SVs to perl. The functions
  are the C wrappers to the users _perl_ fit and jacobian
  functions.  LEVFUNC and JLEVFUNC should be re-entrant so
  the following two routines can called once when the
  Levmar::Func module is loaded and can be shared by all
  object instances.

=end comment

=cut

void *
get_perl_func_wrapper( )

    CODE:
    RETVAL = &LEVFUNC;
    OUTPUT:
    RETVAL

void *
get_perl_jac_wrapper( )

    CODE:
    RETVAL = &JLEVFUNC;
    OUTPUT:
    RETVAL




MODULE = PDL::Fit::Levmar PACKAGE = PDL PREFIX=pdl_run_



#line 483 "lib/PDL/PP.pm"
BOOT:
   /* Get pointer to structure of core shared C routines */
   /* make sure PDL::Core is loaded */
#line 534 "Levmar.xs"
   perl_require_pv ("PDL/Core.pm"); /* make sure PDL::Core is loaded */
#ifndef aTHX_
#define aTHX_
#endif
   if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
   SV* CoreSV = perl_get_sv("PDL::SHARE",FALSE); /* var with core structure */
   if (!CoreSV)
     Perl_croak(aTHX_ "We require the PDL::Core module, which was not found");
   if (!(PDL = INT2PTR(Core*,SvIV( CoreSV )))) /* Core* value */
     Perl_croak(aTHX_ "Got NULL pointer for PDL");
   if (PDL->Version != PDL_CORE_VERSION)
     Perl_croak(aTHX_ "[PDL->Version: %ld PDL_CORE_VERSION: %ld XS_VERSION: %s] PDL::Fit::Levmar needs to be recompiled against the newly installed PDL", (long int)PDL->Version, (long int)PDL_CORE_VERSION, XS_VERSION);