File Coverage

lib/PDL/Slices-pp-splitdim.c

Criterion	Covered	Total	%
statement	21	21	100.0
branch	6	12	50.0
condition			n/a
subroutine			n/a
pod			n/a
total	27	33	81.8

#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_splitdim {
#line 26 "lib/PDL/Slices-pp-splitdim.c"
  PDL_Indx  nthdim;
  PDL_Indx  nsp;
} pdl_params_splitdim;


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

  #line 1789 "lib/PDL/Slices.pd"
                  PDL_Indx i = __params->nthdim;
                  PDL_Indx nsp = __params->nsp;
                  if(nsp == 0) {return PDL->make_error(PDL_EUSERERROR, "Error in splitdim:" "Cannot split to 0\n");}
                  if (i < 0)
                          i = __params->nthdim += __privtrans->pdls[0]->ndims;
                  if (i < 0 || i >= __privtrans->pdls[0]->ndims)
                          return PDL->make_error(PDL_EUSERERROR, "Error in splitdim:" "nthdim %"IND_FLAG" after adjusting for negative must not be negative or greater or equal to number of dims %"IND_FLAG"\n",i,__privtrans->pdls[0]->ndims);
                  if (nsp > __privtrans->pdls[0]->dims[i])
                          return PDL->make_error(PDL_EUSERERROR, "Error in splitdim:" "nsp %"IND_FLAG" cannot be greater than dim %"IND_FLAG"\n",nsp,__privtrans->pdls[0]->dims[i]);
                  if ((__privtrans->pdls[0]->dims[i] % nsp) != 0)
                          return PDL->make_error(PDL_EUSERERROR, "Error in splitdim:" "nsp %"IND_FLAG" non-divisible into dim %"IND_FLAG"\n",nsp,__privtrans->pdls[0]->dims[i]);
                  __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] = __params->nsp;
                  __privtrans->pdls[1]->dims[i+1] = __privtrans->pdls[0]->dims[i] / __params->nsp;
                  __privtrans->incs[i] = __privtrans->pdls[0]->dimincs[i];
                  __privtrans->incs[i+1] = __privtrans->pdls[0]->dimincs[i] * __params->nsp;
                  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 70 "lib/PDL/Slices-pp-splitdim.c"
  return PDL_err;
}


#line 1857 "lib/PDL/PP.pm"
pdl_error pdl_splitdim_free(pdl_trans *__privtrans, char destroy) {
  pdl_error PDL_err = {0, NULL, 0};
#line 78 "lib/PDL/Slices-pp-splitdim.c"
  pdl_params_splitdim *__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_splitdim_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_splitdim_vtable_realdims[] = { 0, 0 };
static char *pdl_splitdim_vtable_parnames[] = { "PARENT","CHILD" };
static short pdl_splitdim_vtable_parflags[] = {
  0,
  PDL_PARAM_ISCREAT|PDL_PARAM_ISCREATEALWAYS|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
};
static pdl_datatypes pdl_splitdim_vtable_partypes[] = { -1, -1 };
static PDL_Indx pdl_splitdim_vtable_realdims_starts[] = { 0, 0 };
static PDL_Indx pdl_splitdim_vtable_realdims_ind_ids[] = { 0 };
static char *pdl_splitdim_vtable_indnames[] = { "" };
pdl_transvtable pdl_splitdim_vtable = {
  0, PDL_ITRANS_ISAFFINE|PDL_ITRANS_TWOWAY|PDL_ITRANS_DO_DATAFLOW_ANY, pdl_splitdim_vtable_gentypes, 1, 2, NULL /*CORE21*/,
  pdl_splitdim_vtable_realdims, pdl_splitdim_vtable_parnames,
  pdl_splitdim_vtable_parflags, pdl_splitdim_vtable_partypes,
  pdl_splitdim_vtable_realdims_starts, pdl_splitdim_vtable_realdims_ind_ids, 0,
  0, pdl_splitdim_vtable_indnames,
  pdl_splitdim_redodims, NULL, NULL,
  pdl_splitdim_free,
  sizeof(pdl_params_splitdim),"PDL::Slices::splitdim"
};


pdl_error pdl_run_splitdim(pdl  *PARENT,pdl  *CHILD,PDL_Indx  nthdim,PDL_Indx  nsp) {
  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_splitdim_vtable);
  if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
  pdl_params_splitdim *__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->nsp) = (nsp); /* CType.get_copy */
  PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
  return PDL_err;
}