File Coverage

lib/PDL/Primitive.xs

Criterion	Covered	Total	%
statement	126	163	77.3
branch	430	1258	34.1
condition			n/a
subroutine			n/a
pod			n/a
total	556	1421	39.1

#line 453 "lib/PDL/PP.pm"
/*
 * THIS FILE WAS GENERATED BY PDL::PP from lib/PDL/Primitive.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_Primitive
extern Core* PDL; /* Structure hold core C functions */
#line 23 "lib/PDL/Primitive.xs"
pdl_error pdl_run_inner(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_outer(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_matmult(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_innerwt(pdl  *a,pdl  *b,pdl  *c,pdl  *d);
pdl_error pdl_run_inner2(pdl  *a,pdl  *b,pdl  *c,pdl  *d);
pdl_error pdl_run_inner2d(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_inner2t(pdl  *a,pdl  *b,pdl  *c,pdl  *d);
pdl_error pdl_run_crossp(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_norm(pdl  *vec,pdl  *norm);
pdl_error pdl_run_indadd(pdl  *input,pdl  *ind,pdl  *sum);
pdl_error pdl_run_conv1d(pdl  *a,pdl  *kern,pdl  *b,int  reflect);
pdl_error pdl_run_in(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_hclip(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_lclip(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_clip(pdl  *a,pdl  *l,pdl  *h,pdl  *c);
pdl_error pdl_run_wtstat(pdl  *a,pdl  *wt,pdl  *avg,pdl  *b,int  deg);
pdl_error pdl_run_statsover(pdl  *a,pdl  *w,pdl  *avg,pdl  *prms,pdl  *min,pdl  *max,pdl  *adev,pdl  *rms);
pdl_error pdl_run_histogram(pdl  *in,pdl  *hist,double  step,double  min,IV  msize);
pdl_error pdl_run_whistogram(pdl  *in,pdl  *wt,pdl  *hist,double  step,double  min,IV  msize);
pdl_error pdl_run_histogram2d(pdl  *ina,pdl  *inb,pdl  *hist,double  stepa,double  mina,IV  masize,double  stepb,double  minb,IV  mbsize);
pdl_error pdl_run_whistogram2d(pdl  *ina,pdl  *inb,pdl  *wt,pdl  *hist,double  stepa,double  mina,IV  masize,double  stepb,double  minb,IV  mbsize);
pdl_error pdl_run_fibonacci(pdl  *i,pdl  *x);
pdl_error pdl_run_append(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_axisvalues(pdl  *i,pdl  *a);
pdl_error pdl_run_cmpvec(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_eqvec(pdl  *a,pdl  *b,pdl  *c);
pdl_error pdl_run_enumvec(pdl  *v,pdl  *k);
pdl_error pdl_run_enumvecg(pdl  *v,pdl  *k);
pdl_error pdl_run_vsearchvec(pdl  *find,pdl  *which,pdl  *found);
pdl_error pdl_run_unionvec(pdl  *a,pdl  *b,pdl  *c,pdl  *nc);
pdl_error pdl_run_intersectvec(pdl  *a,pdl  *b,pdl  *c,pdl  *nc);
pdl_error pdl_run_setdiffvec(pdl  *a,pdl  *b,pdl  *c,pdl  *nc);
pdl_error pdl_run_union_sorted(pdl  *a,pdl  *b,pdl  *c,pdl  *nc);
pdl_error pdl_run_intersect_sorted(pdl  *a,pdl  *b,pdl  *c,pdl  *nc);
pdl_error pdl_run_setdiff_sorted(pdl  *a,pdl  *b,pdl  *c,pdl  *nc);
pdl_error pdl_run_vcos(pdl  *a,pdl  *b,pdl  *vcos);
extern int pdl_srand_threads;
extern uint64_t *pdl_rand_state;
void pdl_srand(uint64_t **s, uint64_t seed, int n);
double pdl_drand(uint64_t *s);
#define PDL_MAYBE_SRAND \
  if (pdl_srand_threads < 0) \
    pdl_srand(&pdl_rand_state, PDL->pdl_seed(), PDL->online_cpus());
#define PDL_RAND_SET_OFFSET(v, thr, pdl) \
  if (v < 0) { \
    if (thr.mag_nthr >= 0) { \
      int thr_no = PDL->magic_get_thread(pdl); \
      if (thr_no < 0) return PDL->make_error_simple(PDL_EFATAL, "Invalid pdl_magic_get_thread!"); \
      v = thr_no == 0 ? thr_no : thr_no % PDL->online_cpus(); \
    } else { \
      v = 0; \
    } \
  }
pdl_error pdl_run_srandom(pdl  *a);
pdl_error pdl_run_random(pdl  *a);
pdl_error pdl_run_randsym(pdl  *a);
pdl_error pdl_run_vsearch_sample(pdl  *vals,pdl  *x,pdl  *idx);
pdl_error pdl_run_vsearch_insert_leftmost(pdl  *vals,pdl  *x,pdl  *idx);
pdl_error pdl_run_vsearch_insert_rightmost(pdl  *vals,pdl  *x,pdl  *idx);
pdl_error pdl_run_vsearch_match(pdl  *vals,pdl  *x,pdl  *idx);
pdl_error pdl_run_vsearch_bin_inclusive(pdl  *vals,pdl  *x,pdl  *idx);
pdl_error pdl_run_vsearch_bin_exclusive(pdl  *vals,pdl  *x,pdl  *idx);
pdl_error pdl_run_interpolate(pdl  *xi,pdl  *x,pdl  *y,pdl  *yi,pdl  *err);
pdl_error pdl_run_which(pdl  *mask,pdl  *inds,pdl  *lastout);
pdl_error pdl_run_which_both(pdl  *mask,pdl  *inds,pdl  *notinds,pdl  *lastout,pdl  *lastoutn);
pdl_error pdl_run_whichover(pdl  *a,pdl  *o);
pdl_error pdl_run_approx_artol(pdl  *got,pdl  *expected,pdl  *result,double  atol,double  rtol);
pdl_error pdl_run_pchip_chim(pdl  *x,pdl  *f,pdl  *d,pdl  *ierr);
pdl_error pdl_run_pchip_chic(pdl  *ic,pdl  *vc,pdl  *mflag,pdl  *x,pdl  *f,pdl  *d,pdl  *ierr);
pdl_error pdl_run_pchip_chsp(pdl  *ic,pdl  *vc,pdl  *x,pdl  *f,pdl  *d,pdl  *ierr);
pdl_error pdl_run_pchip_chfd(pdl  *x,pdl  *f,pdl  *d,pdl  *xe,pdl  *fe,pdl  *de,pdl  *ierr,pdl  *skip);
pdl_error pdl_run_pchip_chfe(pdl  *x,pdl  *f,pdl  *d,pdl  *xe,pdl  *fe,pdl  *ierr,pdl  *skip);
pdl_error pdl_run_pchip_chia(pdl  *x,pdl  *f,pdl  *d,pdl  *la,pdl  *lb,pdl  *ans,pdl  *ierr,pdl  *skip);
pdl_error pdl_run_pchip_chid(pdl  *x,pdl  *f,pdl  *d,pdl  *ia,pdl  *ib,pdl  *ans,pdl  *ierr,pdl  *skip);
pdl_error pdl_run_pchip_chbs(pdl  *x,pdl  *f,pdl  *d,pdl  *knotyp,pdl  *t,pdl  *bcoef,pdl  *ierr);
pdl_error pdl_run_pchip_bvalu(pdl  *t,pdl  *a,pdl  *ideriv,pdl  *x,pdl  *ans,pdl  *inbv);

#line 449 "lib/PDL/PP.pm"
#include "pdlperl.h"
#line 103 "lib/PDL/Primitive.xs"

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

MODULE = PDL::Primitive PACKAGE = PDL PREFIX=pdl_run_

PROTOTYPES: DISABLE


void
pdl_run_inner(a, b, c=PDL_XS_PERLINIT_initsv(c_SV))
  SV *c_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *c
 PPCODE:
  PDL->barf_if_error(pdl_run_inner(a,b,c));
  PDL_XS_RETURN(ST(0) = c_SV)

void
pdl_run_outer(a, b, c=PDL_XS_PERLINIT_initsv(c_SV))
  SV *c_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *c
 PPCODE:
  PDL->barf_if_error(pdl_run_outer(a,b,c));
  PDL_XS_RETURN(ST(0) = c_SV)

void
_matmult_int(a,b,c)
  pdl  *a
  pdl  *b
  pdl  *c
 CODE:
  PDL->barf_if_error(pdl_run_matmult(a,b,c));

void
pdl_run_innerwt(a, b, c, d=PDL_XS_PERLINIT_initsv(d_SV))
  SV *d_SV = items > 3 ? ST(3) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 4) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *c
  pdl *d
 PPCODE:
  PDL->barf_if_error(pdl_run_innerwt(a,b,c,d));
  PDL_XS_RETURN(ST(0) = d_SV)

void
pdl_run_inner2(a, b, c, d=PDL_XS_PERLINIT_initsv(d_SV))
  SV *d_SV = items > 3 ? ST(3) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 4) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *c
  pdl *d
 PPCODE:
  PDL->barf_if_error(pdl_run_inner2(a,b,c,d));
  PDL_XS_RETURN(ST(0) = d_SV)

void
pdl_run_inner2d(a, b, c=PDL_XS_PERLINIT_initsv(c_SV))
  SV *c_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *c
 PPCODE:
  PDL->barf_if_error(pdl_run_inner2d(a,b,c));
  PDL_XS_RETURN(ST(0) = c_SV)

void
pdl_run_inner2t(a, b, c, d=PDL_XS_PERLINIT_initsv(d_SV))
  SV *d_SV = items > 3 ? ST(3) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 4) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *c
  pdl *d
 PPCODE:
  PDL->barf_if_error(pdl_run_inner2t(a,b,c,d));
  PDL_XS_RETURN(ST(0) = d_SV)

void
pdl_run_crossp(a, b, c=PDL_XS_PERLINIT_initsv(c_SV))
  SV *c_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *c
 PPCODE:
  PDL->barf_if_error(pdl_run_crossp(a,b,c));
  PDL_XS_RETURN(ST(0) = c_SV)

void
pdl_run_norm(vec, norm=PDL_XS_PERLINIT_initsv(norm_SV))
  SV *norm_SV = items > 1 ? ST(1) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 2) ? 0 : 1);
 INPUT:
  pdl *vec
  pdl *norm
 PPCODE:
  PDL->barf_if_error(pdl_run_norm(vec,norm));
  PDL_XS_RETURN(ST(0) = norm_SV)

void
pdl_run_indadd(input, ind, sum)
  pdl *input
  pdl *ind
  pdl *sum
 PPCODE:
  PDL->barf_if_error(pdl_run_indadd(input,ind,sum));

void
_conv1d_int(a,kern,b,reflect)
  pdl  *a
  pdl  *kern
  pdl  *b
  int  reflect
 CODE:
  PDL->barf_if_error(pdl_run_conv1d(a,kern,b,reflect));

void
pdl_run_in(a, b, c=PDL_XS_PERLINIT_initsv(c_SV))
  SV *c_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *c
 PPCODE:
  PDL->barf_if_error(pdl_run_in(a,b,c));
  PDL_XS_RETURN(ST(0) = c_SV)

void
_hclip_int(a,b,c)
  pdl  *a
  pdl  *b
  pdl  *c
 CODE:
  PDL->barf_if_error(pdl_run_hclip(a,b,c));

void
_lclip_int(a,b,c)
  pdl  *a
  pdl  *b
  pdl  *c
 CODE:
  PDL->barf_if_error(pdl_run_lclip(a,b,c));

void
_clip_int(a,l,h,c)
  pdl  *a
  pdl  *l
  pdl  *h
  pdl  *c
 CODE:
  PDL->barf_if_error(pdl_run_clip(a,l,h,c));

void
pdl_run_wtstat(a, wt, avg, b=b, deg=deg)
  SV *b_SV = (items == 5) ? ST(3) : NULL;
  SV *deg_SV = (items == 5) ? ST(4) : ST(3);
 PREINIT:
  PDL_XS_PREAMBLE((items == 5) ? 0 : 1);
  if (!((items == 4) || (items == 5)))
    croak("Usage: PDL::wtstat(a,wt,avg,[b],deg) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *a
  pdl *wt
  pdl *avg
  pdl *b=NO_INIT
  int  deg; { deg = (int)SvIV(deg_SV); }
 PPCODE:
  if (b_SV) { b = PDL_CORE_(SvPDLV)(b_SV); } else b = PDL_XS_PERLINIT_initsv(b_SV);
  PDL->barf_if_error(pdl_run_wtstat(a,wt,avg,b,deg));
  PDL_XS_RETURN(ST(0) = b_SV)

void
_statsover_int(a,w,avg,prms,min,max,adev,rms)
  pdl  *a
  pdl  *w
  pdl  *avg
  pdl  *prms
  pdl  *min
  pdl  *max
  pdl  *adev
  pdl  *rms
 CODE:
  PDL->barf_if_error(pdl_run_statsover(a,w,avg,prms,min,max,adev,rms));

void
pdl_run_histogram(in, hist=hist, step=step, min=min, msize=msize)
  SV *hist_SV = (items == 5) ? ST(1) : NULL;
  SV *step_SV = (items == 5) ? ST(2) : ST(1);
  SV *min_SV = (items == 5) ? ST(3) : ST(2);
  SV *msize_SV = (items == 5) ? ST(4) : ST(3);
 PREINIT:
  PDL_XS_PREAMBLE((items == 5) ? 0 : 1);
  if (!((items == 4) || (items == 5)))
    croak("Usage: PDL::histogram(in,[hist],step,min,msize) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *in
  pdl *hist=NO_INIT
  double  step; { step = (double)SvNV(step_SV); }
  double  min; { min = (double)SvNV(min_SV); }
  IV  msize; { msize = (IV)SvIV(msize_SV); }
 PPCODE:
  if (hist_SV) { hist = PDL_CORE_(SvPDLV)(hist_SV); } else hist = PDL_XS_PERLINIT_initsv(hist_SV);
  PDL->barf_if_error(pdl_run_histogram(in,hist,step,min,msize));
  PDL_XS_RETURN(ST(0) = hist_SV)

void
pdl_run_whistogram(in, wt, hist=hist, step=step, min=min, msize=msize)
  SV *hist_SV = (items == 6) ? ST(2) : NULL;
  SV *step_SV = (items == 6) ? ST(3) : ST(2);
  SV *min_SV = (items == 6) ? ST(4) : ST(3);
  SV *msize_SV = (items == 6) ? ST(5) : ST(4);
 PREINIT:
  PDL_XS_PREAMBLE((items == 6) ? 0 : 1);
  if (!((items == 5) || (items == 6)))
    croak("Usage: PDL::whistogram(in,wt,[hist],step,min,msize) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *in
  pdl *wt
  pdl *hist=NO_INIT
  double  step; { step = (double)SvNV(step_SV); }
  double  min; { min = (double)SvNV(min_SV); }
  IV  msize; { msize = (IV)SvIV(msize_SV); }
 PPCODE:
  if (hist_SV) { hist = PDL_CORE_(SvPDLV)(hist_SV); } else hist = PDL_XS_PERLINIT_initsv(hist_SV);
  PDL->barf_if_error(pdl_run_whistogram(in,wt,hist,step,min,msize));
  PDL_XS_RETURN(ST(0) = hist_SV)

void
pdl_run_histogram2d(ina, inb, hist=hist, stepa=stepa, mina=mina, masize=masize, stepb=stepb, minb=minb, mbsize=mbsize)
  SV *hist_SV = (items == 9) ? ST(2) : NULL;
  SV *stepa_SV = (items == 9) ? ST(3) : ST(2);
  SV *mina_SV = (items == 9) ? ST(4) : ST(3);
  SV *masize_SV = (items == 9) ? ST(5) : ST(4);
  SV *stepb_SV = (items == 9) ? ST(6) : ST(5);
  SV *minb_SV = (items == 9) ? ST(7) : ST(6);
  SV *mbsize_SV = (items == 9) ? ST(8) : ST(7);
 PREINIT:
  PDL_XS_PREAMBLE((items == 9) ? 0 : 1);
  if (!((items == 8) || (items == 9)))
    croak("Usage: PDL::histogram2d(ina,inb,[hist],stepa,mina,masize,stepb,minb,mbsize) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *ina
  pdl *inb
  pdl *hist=NO_INIT
  double  stepa; { stepa = (double)SvNV(stepa_SV); }
  double  mina; { mina = (double)SvNV(mina_SV); }
  IV  masize; { masize = (IV)SvIV(masize_SV); }
  double  stepb; { stepb = (double)SvNV(stepb_SV); }
  double  minb; { minb = (double)SvNV(minb_SV); }
  IV  mbsize; { mbsize = (IV)SvIV(mbsize_SV); }
 PPCODE:
  if (hist_SV) { hist = PDL_CORE_(SvPDLV)(hist_SV); } else hist = PDL_XS_PERLINIT_initsv(hist_SV);
  PDL->barf_if_error(pdl_run_histogram2d(ina,inb,hist,stepa,mina,masize,stepb,minb,mbsize));
  PDL_XS_RETURN(ST(0) = hist_SV)

void
pdl_run_whistogram2d(ina, inb, wt, hist=hist, stepa=stepa, mina=mina, masize=masize, stepb=stepb, minb=minb, mbsize=mbsize)
  SV *hist_SV = (items == 10) ? ST(3) : NULL;
  SV *stepa_SV = (items == 10) ? ST(4) : ST(3);
  SV *mina_SV = (items == 10) ? ST(5) : ST(4);
  SV *masize_SV = (items == 10) ? ST(6) : ST(5);
  SV *stepb_SV = (items == 10) ? ST(7) : ST(6);
  SV *minb_SV = (items == 10) ? ST(8) : ST(7);
  SV *mbsize_SV = (items == 10) ? ST(9) : ST(8);
 PREINIT:
  PDL_XS_PREAMBLE((items == 10) ? 0 : 1);
  if (!((items == 10) || (items == 9)))
    croak("Usage: PDL::whistogram2d(ina,inb,wt,[hist],stepa,mina,masize,stepb,minb,mbsize) (you may leave [outputs] and values with =defaults out of list)");
 INPUT:
  pdl *ina
  pdl *inb
  pdl *wt
  pdl *hist=NO_INIT
  double  stepa; { stepa = (double)SvNV(stepa_SV); }
  double  mina; { mina = (double)SvNV(mina_SV); }
  IV  masize; { masize = (IV)SvIV(masize_SV); }
  double  stepb; { stepb = (double)SvNV(stepb_SV); }
  double  minb; { minb = (double)SvNV(minb_SV); }
  IV  mbsize; { mbsize = (IV)SvIV(mbsize_SV); }
 PPCODE:
  if (hist_SV) { hist = PDL_CORE_(SvPDLV)(hist_SV); } else hist = PDL_XS_PERLINIT_initsv(hist_SV);
  PDL->barf_if_error(pdl_run_whistogram2d(ina,inb,wt,hist,stepa,mina,masize,stepb,minb,mbsize));
  PDL_XS_RETURN(ST(0) = hist_SV)

void
_fibonacci_int(i,x)
  pdl  *i
  pdl  *x
 CODE:
  PDL->barf_if_error(pdl_run_fibonacci(i,x));

void
_append_int(a,b,c)
  pdl  *a
  pdl  *b
  pdl  *c
 CODE:
  PDL->barf_if_error(pdl_run_append(a,b,c));

void
pdl_run_axisvalues(i, a=a)
  SV *a_SV = items > 1 ? ST(1) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 2) ? 0 : 1);
 INPUT:
  pdl *i
  pdl *a=NO_INIT
 PPCODE:
  PDL_XS_INPLACE(i, a)
  PDL->barf_if_error(pdl_run_axisvalues(i,a));
  PDL_XS_RETURN(ST(0) = a_SV)

void
pdl_run_cmpvec(a, b, c=PDL_XS_PERLINIT_initsv(c_SV))
  SV *c_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *c
 PPCODE:
  PDL->barf_if_error(pdl_run_cmpvec(a,b,c));
  PDL_XS_RETURN(ST(0) = c_SV)

void
pdl_run_eqvec(a, b, c=PDL_XS_PERLINIT_initsv(c_SV))
  SV *c_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *c
 PPCODE:
  PDL->barf_if_error(pdl_run_eqvec(a,b,c));
  PDL_XS_RETURN(ST(0) = c_SV)

void
pdl_run_enumvec(v, k=PDL_XS_PERLINIT_initsv(k_SV))
  SV *k_SV = items > 1 ? ST(1) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 2) ? 0 : 1);
 INPUT:
  pdl *v
  pdl *k
 PPCODE:
  PDL->barf_if_error(pdl_run_enumvec(v,k));
  PDL_XS_RETURN(ST(0) = k_SV)

void
pdl_run_enumvecg(v, k=PDL_XS_PERLINIT_initsv(k_SV))
  SV *k_SV = items > 1 ? ST(1) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 2) ? 0 : 1);
 INPUT:
  pdl *v
  pdl *k
 PPCODE:
  PDL->barf_if_error(pdl_run_enumvecg(v,k));
  PDL_XS_RETURN(ST(0) = k_SV)

void
pdl_run_vsearchvec(find, which, found=PDL_XS_PERLINIT_initsv(found_SV))
  SV *found_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *find
  pdl *which
  pdl *found
 PPCODE:
  PDL->barf_if_error(pdl_run_vsearchvec(find,which,found));
  PDL_XS_RETURN(ST(0) = found_SV)

void
_unionvec_int(a,b,c,nc)
  pdl  *a
  pdl  *b
  pdl  *c
  pdl  *nc
 CODE:
  PDL->barf_if_error(pdl_run_unionvec(a,b,c,nc));

void
_intersectvec_int(a,b,c,nc)
  pdl  *a
  pdl  *b
  pdl  *c
  pdl  *nc
 CODE:
  PDL->barf_if_error(pdl_run_intersectvec(a,b,c,nc));

void
_setdiffvec_int(a,b,c,nc)
  pdl  *a
  pdl  *b
  pdl  *c
  pdl  *nc
 CODE:
  PDL->barf_if_error(pdl_run_setdiffvec(a,b,c,nc));

void
_union_sorted_int(a,b,c,nc)
  pdl  *a
  pdl  *b
  pdl  *c
  pdl  *nc
 CODE:
  PDL->barf_if_error(pdl_run_union_sorted(a,b,c,nc));

void
_intersect_sorted_int(a,b,c,nc)
  pdl  *a
  pdl  *b
  pdl  *c
  pdl  *nc
 CODE:
  PDL->barf_if_error(pdl_run_intersect_sorted(a,b,c,nc));

void
_setdiff_sorted_int(a,b,c,nc)
  pdl  *a
  pdl  *b
  pdl  *c
  pdl  *nc
 CODE:
  PDL->barf_if_error(pdl_run_setdiff_sorted(a,b,c,nc));

void
pdl_run_vcos(a, b, vcos=PDL_XS_PERLINIT_initsv(vcos_SV))
  SV *vcos_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *b
  pdl *vcos
 PPCODE:
  PDL->barf_if_error(pdl_run_vcos(a,b,vcos));
  PDL_XS_RETURN(ST(0) = vcos_SV)

void
_srandom_int(a)
  pdl  *a
 CODE:
  PDL->barf_if_error(pdl_run_srandom(a));

void
_random_int(a)
  pdl  *a
 CODE:
  PDL->barf_if_error(pdl_run_random(a));

void
_randsym_int(a)
  pdl  *a
 CODE:
  PDL->barf_if_error(pdl_run_randsym(a));

void
pdl_run_vsearch_sample(vals, x, idx=PDL_XS_PERLINIT_initsv(idx_SV))
  SV *idx_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *vals
  pdl *x
  pdl *idx
 PPCODE:
  PDL->barf_if_error(pdl_run_vsearch_sample(vals,x,idx));
  PDL_XS_RETURN(ST(0) = idx_SV)

void
pdl_run_vsearch_insert_leftmost(vals, x, idx=PDL_XS_PERLINIT_initsv(idx_SV))
  SV *idx_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *vals
  pdl *x
  pdl *idx
 PPCODE:
  PDL->barf_if_error(pdl_run_vsearch_insert_leftmost(vals,x,idx));
  PDL_XS_RETURN(ST(0) = idx_SV)

void
pdl_run_vsearch_insert_rightmost(vals, x, idx=PDL_XS_PERLINIT_initsv(idx_SV))
  SV *idx_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *vals
  pdl *x
  pdl *idx
 PPCODE:
  PDL->barf_if_error(pdl_run_vsearch_insert_rightmost(vals,x,idx));
  PDL_XS_RETURN(ST(0) = idx_SV)

void
pdl_run_vsearch_match(vals, x, idx=PDL_XS_PERLINIT_initsv(idx_SV))
  SV *idx_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *vals
  pdl *x
  pdl *idx
 PPCODE:
  PDL->barf_if_error(pdl_run_vsearch_match(vals,x,idx));
  PDL_XS_RETURN(ST(0) = idx_SV)

void
pdl_run_vsearch_bin_inclusive(vals, x, idx=PDL_XS_PERLINIT_initsv(idx_SV))
  SV *idx_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *vals
  pdl *x
  pdl *idx
 PPCODE:
  PDL->barf_if_error(pdl_run_vsearch_bin_inclusive(vals,x,idx));
  PDL_XS_RETURN(ST(0) = idx_SV)

void
pdl_run_vsearch_bin_exclusive(vals, x, idx=PDL_XS_PERLINIT_initsv(idx_SV))
  SV *idx_SV = items > 2 ? ST(2) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 3) ? 0 : 1);
 INPUT:
  pdl *vals
  pdl *x
  pdl *idx
 PPCODE:
  PDL->barf_if_error(pdl_run_vsearch_bin_exclusive(vals,x,idx));
  PDL_XS_RETURN(ST(0) = idx_SV)

void
pdl_run_interpolate(xi, x, y, yi=PDL_XS_PERLINIT_initsv(yi_SV), err=PDL_XS_PERLINIT_initsv(err_SV))
  SV *yi_SV = items > 3 ? ST(3) : NULL;
  SV *err_SV = items > 4 ? ST(4) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 5) ? 0 : 2);
 INPUT:
  pdl *xi
  pdl *x
  pdl *y
  pdl *yi
  pdl *err
 PPCODE:
  PDL->barf_if_error(pdl_run_interpolate(xi,x,y,yi,err));
  PDL_XS_RETURN(ST(0) = yi_SV;ST(1) = err_SV)

void
_which_int(mask,inds,lastout)
  pdl  *mask
  pdl  *inds
  pdl  *lastout
 CODE:
  PDL->barf_if_error(pdl_run_which(mask,inds,lastout));

void
_which_both_int(mask,inds,notinds,lastout,lastoutn)
  pdl  *mask
  pdl  *inds
  pdl  *notinds
  pdl  *lastout
  pdl  *lastoutn
 CODE:
  PDL->barf_if_error(pdl_run_which_both(mask,inds,notinds,lastout,lastoutn));

void
pdl_run_whichover(a, o=o)
  SV *o_SV = items > 1 ? ST(1) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 2) ? 0 : 1);
 INPUT:
  pdl *a
  pdl *o=NO_INIT
 PPCODE:
  PDL_XS_INPLACE(a, o)
  PDL->barf_if_error(pdl_run_whichover(a,o));
  PDL_XS_RETURN(ST(0) = o_SV)

void
pdl_run_approx_artol(got, expected, atol=1e-06, rtol=0, result=PDL_XS_PERLINIT_initsv(result_SV))
  SV *result_SV = items > 4 ? ST(4) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE(1);
 INPUT:
  pdl *got
  pdl *expected
  double  atol
  double  rtol
  pdl *result
 PPCODE:
  PDL->barf_if_error(pdl_run_approx_artol(got,expected,result,atol,rtol));
  PDL_XS_RETURN(ST(0) = result_SV)

void
pdl_run_pchip_chim(x, f, d=PDL_XS_PERLINIT_initsv(d_SV), ierr=PDL_XS_PERLINIT_initsv(ierr_SV))
  SV *d_SV = items > 2 ? ST(2) : NULL;
  SV *ierr_SV = items > 3 ? ST(3) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 4) ? 0 : 2);
 INPUT:
  pdl *x
  pdl *f
  pdl *d
  pdl *ierr
 PPCODE:
  PDL->barf_if_error(pdl_run_pchip_chim(x,f,d,ierr));
  PDL_XS_RETURN(ST(0) = d_SV;ST(1) = ierr_SV)

void
pdl_run_pchip_chic(ic, vc, mflag, x, f, d=PDL_XS_PERLINIT_initsv(d_SV), ierr=PDL_XS_PERLINIT_initsv(ierr_SV))
  SV *d_SV = items > 5 ? ST(5) : NULL;
  SV *ierr_SV = items > 6 ? ST(6) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 7) ? 0 : 2);
 INPUT:
  pdl *ic
  pdl *vc
  pdl *mflag
  pdl *x
  pdl *f
  pdl *d
  pdl *ierr
 PPCODE:
  PDL->barf_if_error(pdl_run_pchip_chic(ic,vc,mflag,x,f,d,ierr));
  PDL_XS_RETURN(ST(0) = d_SV;ST(1) = ierr_SV)

void
pdl_run_pchip_chsp(ic, vc, x, f, d=PDL_XS_PERLINIT_initsv(d_SV), ierr=PDL_XS_PERLINIT_initsv(ierr_SV))
  SV *d_SV = items > 4 ? ST(4) : NULL;
  SV *ierr_SV = items > 5 ? ST(5) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 6) ? 0 : 2);
 INPUT:
  pdl *ic
  pdl *vc
  pdl *x
  pdl *f
  pdl *d
  pdl *ierr
 PPCODE:
  PDL->barf_if_error(pdl_run_pchip_chsp(ic,vc,x,f,d,ierr));
  PDL_XS_RETURN(ST(0) = d_SV;ST(1) = ierr_SV)

void
pdl_run_pchip_chfd(x, f, d, xe, fe=PDL_XS_PERLINIT_initsv(fe_SV), de=PDL_XS_PERLINIT_initsv(de_SV), ierr=PDL_XS_PERLINIT_initsv(ierr_SV), skip=PDL_XS_PERLINIT_initsv(skip_SV))
  SV *fe_SV = items > 4 ? ST(4) : NULL;
  SV *de_SV = items > 5 ? ST(5) : NULL;
  SV *ierr_SV = items > 6 ? ST(6) : NULL;
  SV *skip_SV = items > 7 ? ST(7) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 8) ? 0 : 4);
 INPUT:
  pdl *x
  pdl *f
  pdl *d
  pdl *xe
  pdl *fe
  pdl *de
  pdl *ierr
  pdl *skip
 PPCODE:
  PDL->barf_if_error(pdl_run_pchip_chfd(x,f,d,xe,fe,de,ierr,skip));
  PDL_XS_RETURN(ST(0) = fe_SV;ST(1) = de_SV;ST(2) = ierr_SV;ST(3) = skip_SV)

void
pdl_run_pchip_chfe(x, f, d, xe, fe=PDL_XS_PERLINIT_initsv(fe_SV), ierr=PDL_XS_PERLINIT_initsv(ierr_SV), skip=PDL_XS_PERLINIT_initsv(skip_SV))
  SV *fe_SV = items > 4 ? ST(4) : NULL;
  SV *ierr_SV = items > 5 ? ST(5) : NULL;
  SV *skip_SV = items > 6 ? ST(6) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 7) ? 0 : 3);
 INPUT:
  pdl *x
  pdl *f
  pdl *d
  pdl *xe
  pdl *fe
  pdl *ierr
  pdl *skip
 PPCODE:
  PDL->barf_if_error(pdl_run_pchip_chfe(x,f,d,xe,fe,ierr,skip));
  PDL_XS_RETURN(ST(0) = fe_SV;ST(1) = ierr_SV;ST(2) = skip_SV)

void
pdl_run_pchip_chia(x, f, d, la, lb, ans=PDL_XS_PERLINIT_initsv(ans_SV), ierr=PDL_XS_PERLINIT_initsv(ierr_SV), skip=PDL_XS_PERLINIT_initsv(skip_SV))
  SV *ans_SV = items > 5 ? ST(5) : NULL;
  SV *ierr_SV = items > 6 ? ST(6) : NULL;
  SV *skip_SV = items > 7 ? ST(7) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 8) ? 0 : 3);
 INPUT:
  pdl *x
  pdl *f
  pdl *d
  pdl *la
  pdl *lb
  pdl *ans
  pdl *ierr
  pdl *skip
 PPCODE:
  PDL->barf_if_error(pdl_run_pchip_chia(x,f,d,la,lb,ans,ierr,skip));
  PDL_XS_RETURN(ST(0) = ans_SV;ST(1) = ierr_SV;ST(2) = skip_SV)

void
pdl_run_pchip_chid(x, f, d, ia, ib, ans=PDL_XS_PERLINIT_initsv(ans_SV), ierr=PDL_XS_PERLINIT_initsv(ierr_SV), skip=PDL_XS_PERLINIT_initsv(skip_SV))
  SV *ans_SV = items > 5 ? ST(5) : NULL;
  SV *ierr_SV = items > 6 ? ST(6) : NULL;
  SV *skip_SV = items > 7 ? ST(7) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 8) ? 0 : 3);
 INPUT:
  pdl *x
  pdl *f
  pdl *d
  pdl *ia
  pdl *ib
  pdl *ans
  pdl *ierr
  pdl *skip
 PPCODE:
  PDL->barf_if_error(pdl_run_pchip_chid(x,f,d,ia,ib,ans,ierr,skip));
  PDL_XS_RETURN(ST(0) = ans_SV;ST(1) = ierr_SV;ST(2) = skip_SV)

void
pdl_run_pchip_chbs(x, f, d, knotyp, t=PDL_XS_PERLINIT_initsv(t_SV), bcoef=PDL_XS_PERLINIT_initsv(bcoef_SV), ierr=PDL_XS_PERLINIT_initsv(ierr_SV))
  SV *t_SV = items > 4 ? ST(4) : NULL;
  SV *bcoef_SV = items > 5 ? ST(5) : NULL;
  SV *ierr_SV = items > 6 ? ST(6) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 7) ? 0 : 3);
 INPUT:
  pdl *x
  pdl *f
  pdl *d
  pdl *knotyp
  pdl *t
  pdl *bcoef
  pdl *ierr
 PPCODE:
  PDL->barf_if_error(pdl_run_pchip_chbs(x,f,d,knotyp,t,bcoef,ierr));
  PDL_XS_RETURN(ST(0) = t_SV;ST(1) = bcoef_SV;ST(2) = ierr_SV)

void
pdl_run_pchip_bvalu(t, a, ideriv, x, ans=PDL_XS_PERLINIT_initsv(ans_SV), inbv=PDL_XS_PERLINIT_initsv(inbv_SV))
  SV *ans_SV = items > 4 ? ST(4) : NULL;
  SV *inbv_SV = items > 5 ? ST(5) : NULL;
 PREINIT:
  PDL_XS_PREAMBLE((items == 6) ? 0 : 2);
 INPUT:
  pdl *t
  pdl *a
  pdl *ideriv
  pdl *x
  pdl *ans
  pdl *inbv
 PPCODE:
  PDL->barf_if_error(pdl_run_pchip_bvalu(t,a,ideriv,x,ans,inbv));
  PDL_XS_RETURN(ST(0) = ans_SV;ST(1) = inbv_SV)


#line 483 "lib/PDL/PP.pm"
BOOT:
   /* Get pointer to structure of core shared C routines */
   /* make sure PDL::Core is loaded */
#line 900 "lib/PDL/Primitive.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::Primitive needs to be recompiled against the newly installed PDL", (long int)PDL->Version, (long int)PDL_CORE_VERSION, XS_VERSION);