File Coverage

lib/PDL/Slices-pp-diagonal.c

Criterion	Covered	Total	%
statement	32	32	100.0
branch	12	22	54.5
condition			n/a
subroutine			n/a
pod			n/a
total	44	54	81.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 23 "lib/PDL/Slices-pp-diagonal.c"
static int cmp_pdll(const void *a_,const void *b_) {
        PDL_Indx *a = (PDL_Indx *)a_; PDL_Indx *b=(PDL_Indx *)b_;
        if(*a>*b) return 1;
        else if(*a==*b) return 0;
        else return -1;
}



#line 1846 "lib/PDL/PP.pm"
typedef struct pdl_params_diagonal {
#line 35 "lib/PDL/Slices-pp-diagonal.c"
  PDL_Indx  *whichdims;
  PDL_Indx  whichdims_count;
} pdl_params_diagonal;


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

  #line 2464 "lib/PDL/Slices.pd"

                  PDL_Indx nthp,nthc,nthd, cd = __params->whichdims[0];
                  PDL_RETERROR(PDL_err, PDL->reallocdims(__it,__privtrans->pdls[0]->ndims-__params->whichdims_count+1));;
                  __privtrans->incs = malloc(sizeof(*__privtrans->incs) * __privtrans->pdls[1]->ndims); /* CType.get_malloc */
  ;
                  __privtrans->offs = 0;
                  if (__params->whichdims[__params->whichdims_count-1] >= __privtrans->pdls[0]->ndims ||
                          __params->whichdims[0] < 0)
                          return PDL->make_error(PDL_EUSERERROR, "Error in diagonal:" "dim out of range");
                  nthd=0; nthc=0;
                  for(nthp=0; nthp<__privtrans->pdls[0]->ndims; nthp++)
                          if (nthd < __params->whichdims_count &&
                              nthp == __params->whichdims[nthd]) {
                                  if (!nthd) {
                                          __privtrans->pdls[1]->dims[cd] = __privtrans->pdls[0]->dims[cd];
                                          nthc++;
                                          __privtrans->incs[cd] = 0;
                                  }
                                  if (nthd && __params->whichdims[nthd] ==
                                      __params->whichdims[nthd-1])
                                         return PDL->make_error(PDL_EUSERERROR, "Error in diagonal:" "dims must be unique");
                                  nthd++; /* advance pointer into whichdims */
                                  if(__privtrans->pdls[1]->dims[cd] !=
                                      __privtrans->pdls[0]->dims[nthp]) {
                                          return PDL->make_error(PDL_EUSERERROR, "Error in diagonal:" "Different dims %"IND_FLAG" and %"IND_FLAG"",__privtrans->pdls[1]->dims[cd],__privtrans->pdls[0]->dims[nthp]);
                                  }
                                  __privtrans->incs[cd] += __privtrans->pdls[0]->dimincs[nthp];
                          } else {
                                  __privtrans->incs[nthc] = __privtrans->pdls[0]->dimincs[nthp];
                                  __privtrans->pdls[1]->dims[nthc] = __privtrans->pdls[0]->dims[nthp];
                                  nthc++;
                          }
                  PDL_RETERROR(PDL_err, PDL->setdims_careful(__it));;

  #line 84 "lib/PDL/Slices-pp-diagonal.c"
  return PDL_err;
}


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

static pdl_datatypes pdl_diagonal_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_diagonal_vtable_realdims[] = { 0, 0 };
static char *pdl_diagonal_vtable_parnames[] = { "PARENT","CHILD" };
static short pdl_diagonal_vtable_parflags[] = {
  0,
  PDL_PARAM_ISCREAT|PDL_PARAM_ISCREATEALWAYS|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
};
static pdl_datatypes pdl_diagonal_vtable_partypes[] = { -1, -1 };
static PDL_Indx pdl_diagonal_vtable_realdims_starts[] = { 0, 0 };
static PDL_Indx pdl_diagonal_vtable_realdims_ind_ids[] = { 0 };
static char *pdl_diagonal_vtable_indnames[] = { "" };
pdl_transvtable pdl_diagonal_vtable = {
  0, PDL_ITRANS_ISAFFINE|PDL_ITRANS_TWOWAY|PDL_ITRANS_DO_DATAFLOW_ANY, pdl_diagonal_vtable_gentypes, 1, 2, NULL /*CORE21*/,
  pdl_diagonal_vtable_realdims, pdl_diagonal_vtable_parnames,
  pdl_diagonal_vtable_parflags, pdl_diagonal_vtable_partypes,
  pdl_diagonal_vtable_realdims_starts, pdl_diagonal_vtable_realdims_ind_ids, 0,
  0, pdl_diagonal_vtable_indnames,
  pdl_diagonal_redodims, NULL, NULL,
  pdl_diagonal_free,
  sizeof(pdl_params_diagonal),"PDL::Slices::diagonal"
};


pdl_error pdl_run_diagonal(pdl  *PARENT,pdl  *CHILD,PDL_Indx  *whichdims,PDL_Indx  whichdims_count) {
  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_diagonal_vtable);
  if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
  pdl_params_diagonal *__params = __privtrans->params;
  __privtrans->pdls[0] = PARENT;
  __privtrans->pdls[1] = CHILD;
  PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
  PARENT = __privtrans->pdls[0];
  CHILD = __privtrans->pdls[1];
  __params->whichdims = malloc(sizeof(*__params->whichdims) * whichdims_count); /* CType.get_malloc */

  			  if(!whichdims) {__params->whichdims=0;} /* CType.get_copy */
  			  else {int __malloc_ind_1;
  				for(__malloc_ind_1 = 0;
  					__malloc_ind_1 < whichdims_count;
  					__malloc_ind_1 ++) { __params->whichdims[__malloc_ind_1] = whichdims[__malloc_ind_1]; }}
  (__params->whichdims_count) = (whichdims_count); /* CType.get_copy */


#line 2450 "lib/PDL/Slices.pd"

                if (__params->whichdims_count < 1)
                  return PDL->make_error(PDL_EUSERERROR, "Error in diagonal:" "must have at least 1 dimension");
                qsort(__params->whichdims, __params->whichdims_count, sizeof(PDL_Indx),
                        cmp_pdll);

#line 151 "lib/PDL/Slices-pp-diagonal.c"
  PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
  return PDL_err;
}