File Coverage

lib/PDL/Slices-pp-lags.c

Criterion	Covered	Total	%
statement	21	22	95.4
branch	6	12	50.0
condition			n/a
subroutine			n/a
pod			n/a
total	27	34	79.4

#line 453 "lib/PDL/PP.pm"
/*
 * THIS FILE WAS GENERATED BY PDL::PP from lib/PDL/Slices.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_Slices
extern Core* PDL; /* Structure hold core C functions */

#line 1846 "lib/PDL/PP.pm"
typedef struct pdl_params_lags {
#line 26 "lib/PDL/Slices-pp-lags.c"
  PDL_Indx  nthdim;
  PDL_Indx  step;
  PDL_Indx  nlags;
} pdl_params_lags;


#line 1857 "lib/PDL/PP.pm"
pdl_error pdl_lags_redodims(pdl_trans *__privtrans) {
  pdl_error PDL_err = {0, NULL, 0};
#line 36 "lib/PDL/Slices-pp-lags.c"
  pdl_params_lags *__params = __privtrans->params; (void)__params;
  pdl *__it = __privtrans->pdls[1]; (void) __it;
  PDL->hdr_childcopy(__privtrans); __privtrans->dims_redone = 1;

  #line 1731 "lib/PDL/Slices.pd"
                  PDL_Indx i;
                  if (__params->nthdim < 0)  /* the usual conventions */
                     __params->nthdim += __privtrans->pdls[0]->ndims;
                  if (__params->nthdim < 0 || __params->nthdim >= __privtrans->pdls[0]->ndims)
                     return PDL->make_error(PDL_EUSERERROR, "Error in lags:" "dim out of range");
                  if (__params->nlags < 1)
                     return PDL->make_error(PDL_EUSERERROR, "Error in lags:" "number of lags must be positive");
                  if (__params->step < 1)
                     return PDL->make_error(PDL_EUSERERROR, "Error in lags:" "step must be positive");
                  __privtrans->offs = 0;
                  PDL_RETERROR(PDL_err, PDL->reallocdims(__it,__privtrans->pdls[0]->ndims+1));;
                  __privtrans->incs = malloc(sizeof(*__privtrans->incs) * __privtrans->pdls[1]->ndims); /* CType.get_malloc */
  ;
                  for(i=0; i<__params->nthdim; i++) {
                          __privtrans->pdls[1]->dims[i] = __privtrans->pdls[0]->dims[i];
                          __privtrans->incs[i] = __privtrans->pdls[0]->dimincs[i];
                  }
                  __privtrans->pdls[1]->dims[i] = __privtrans->pdls[0]->dims[i] - __params->step * (__params->nlags-1);
                  if (__privtrans->pdls[1]->dims[i] < 1)
                    return PDL->make_error(PDL_EUSERERROR, "Error in lags:" "product of step size and number of lags too large");
                  __privtrans->pdls[1]->dims[i+1] = __params->nlags;
                  __privtrans->incs[i] = (__privtrans->pdls[0]->dimincs[i]);
                  __privtrans->incs[i+1] = - __privtrans->pdls[0]->dimincs[i] * __params->step;
                  __privtrans->offs += (__privtrans->pdls[1]->dims[i+1] - 1) * (-__privtrans->incs[i+1]);
                  i++;
                  for(; i<__privtrans->pdls[0]->ndims; i++) {
                          __privtrans->pdls[1]->dims[i+1] = __privtrans->pdls[0]->dims[i];
                          __privtrans->incs[i+1] = __privtrans->pdls[0]->dimincs[i];
                  }
                  PDL_RETERROR(PDL_err, PDL->setdims_careful(__it));;
  #line 72 "lib/PDL/Slices-pp-lags.c"
  return PDL_err;
}


#line 1857 "lib/PDL/PP.pm"
pdl_error pdl_lags_free(pdl_trans *__privtrans, char destroy) {
  pdl_error PDL_err = {0, NULL, 0};
#line 80 "lib/PDL/Slices-pp-lags.c"
  pdl_params_lags *__params = __privtrans->params; (void)__params;
  pdl *__it = __privtrans->pdls[1]; (void) __it;
  PDL_FREE_CODE(__privtrans, destroy, , free(__privtrans->incs); /* CType.get_free */
  )  return PDL_err;
}

static pdl_datatypes pdl_lags_vtable_gentypes[] = { PDL_SB, PDL_B, PDL_S, PDL_US, PDL_L, PDL_UL, PDL_IND, PDL_ULL, PDL_LL, PDL_F, PDL_D, PDL_LD, PDL_CF, PDL_CD, PDL_CLD, -1 };
static PDL_Indx pdl_lags_vtable_realdims[] = { 0, 0 };
static char *pdl_lags_vtable_parnames[] = { "PARENT","CHILD" };
static short pdl_lags_vtable_parflags[] = {
  0,
  PDL_PARAM_ISCREAT|PDL_PARAM_ISCREATEALWAYS|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
};
static pdl_datatypes pdl_lags_vtable_partypes[] = { -1, -1 };
static PDL_Indx pdl_lags_vtable_realdims_starts[] = { 0, 0 };
static PDL_Indx pdl_lags_vtable_realdims_ind_ids[] = { 0 };
static char *pdl_lags_vtable_indnames[] = { "" };
pdl_transvtable pdl_lags_vtable = {
  0, PDL_ITRANS_ISAFFINE|PDL_ITRANS_TWOWAY|PDL_ITRANS_DO_DATAFLOW_ANY, pdl_lags_vtable_gentypes, 1, 2, NULL /*CORE21*/,
  pdl_lags_vtable_realdims, pdl_lags_vtable_parnames,
  pdl_lags_vtable_parflags, pdl_lags_vtable_partypes,
  pdl_lags_vtable_realdims_starts, pdl_lags_vtable_realdims_ind_ids, 0,
  0, pdl_lags_vtable_indnames,
  pdl_lags_redodims, NULL, NULL,
  pdl_lags_free,
  sizeof(pdl_params_lags),"PDL::Slices::lags"
};


pdl_error pdl_run_lags(pdl  *PARENT,pdl  *CHILD,PDL_Indx  nthdim,PDL_Indx  step,PDL_Indx  nlags) {
  pdl_error PDL_err = {0, NULL, 0};
  if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
  pdl_trans *__privtrans = PDL->create_trans(&pdl_lags_vtable);
  if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
  pdl_params_lags *__params = __privtrans->params;
  __privtrans->pdls[0] = PARENT;
  __privtrans->pdls[1] = CHILD;
  PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
  (__params->nthdim) = (nthdim); /* CType.get_copy */
  (__params->step) = (step); /* CType.get_copy */
  (__params->nlags) = (nlags); /* CType.get_copy */
  PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
  return PDL_err;
}