File Coverage

lib/PDL/Slices-pp-broadcastI.c

Criterion	Covered	Total	%
statement	28	28	100.0
branch	9	16	56.2
condition			n/a
subroutine			n/a
pod			n/a
total	37	44	84.0

#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_broadcastI {
#line 26 "lib/PDL/Slices-pp-broadcastI.c"
  PDL_Indx  id;
  PDL_Indx  *whichdims;
  PDL_Indx  whichdims_count;
  PDL_Indx  nrealwhichdims;
} pdl_params_broadcastI;


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

  #line 1880 "lib/PDL/Slices.pd"
                  PDL_Indx nthc,i,j,flag;
                  PDL_RETERROR(PDL_err, PDL->reallocdims(__it,__privtrans->pdls[0]->ndims));;
                  __privtrans->incs = malloc(sizeof(*__privtrans->incs) * __privtrans->pdls[1]->ndims); /* CType.get_malloc */
  ;
                  __privtrans->offs = 0;
                  nthc=0;
                  for(i=0; i<__privtrans->pdls[0]->ndims; i++) {
                          flag=0;
                          if(__privtrans->pdls[0]->nbroadcastids > __params->id && __params->id >= 0 &&
                             i == __privtrans->pdls[0]->broadcastids[__params->id]) {
                             nthc += __params->whichdims_count;
                          }
                          for(j=0; j<__params->whichdims_count; j++) {
                                  if(__params->whichdims[j] == i) {flag=1; break;}
                          }
                          if(flag) {
                                  continue;
                          }
                          __privtrans->pdls[1]->dims[nthc] = __privtrans->pdls[0]->dims[i];
                          __privtrans->incs[nthc] = __privtrans->pdls[0]->dimincs[i];
                          nthc++;
                  }
                  for(i=0; i<__params->whichdims_count; i++) {
                          PDL_Indx cdim,pdim;
                          cdim = i +
                           (__privtrans->pdls[0]->nbroadcastids > __params->id && __params->id >= 0?
                            __privtrans->pdls[0]->broadcastids[__params->id] : __privtrans->pdls[0]->ndims)
                            - __params->nrealwhichdims;
                          pdim = __params->whichdims[i];
                          if(pdim == -1) {
                                  __privtrans->pdls[1]->dims[cdim] = 1;
                                  __privtrans->incs[cdim] = 0;
                          } else {
                                  __privtrans->pdls[1]->dims[cdim] = __privtrans->pdls[0]->dims[pdim];
                                  __privtrans->incs[cdim] = __privtrans->pdls[0]->dimincs[pdim];
                          }
                  }
                  PDL_RETERROR(PDL_err, PDL->setdims_careful(__it));;
                  PDL_RETERROR(PDL_err, PDL->reallocbroadcastids(__privtrans->pdls[1],
                                  PDLMAX(__params->id+1, __privtrans->pdls[0]->nbroadcastids)));
                  for(i=0; i<__privtrans->pdls[1]->nbroadcastids-1; i++) {
                          __privtrans->pdls[1]->broadcastids[i] =
                           (__privtrans->pdls[0]->nbroadcastids > i ?
                            __privtrans->pdls[0]->broadcastids[i] : __privtrans->pdls[0]->ndims) +
                           (i <= __params->id ? - __params->nrealwhichdims :
                            __params->whichdims_count - __params->nrealwhichdims);
                  }
                  __privtrans->pdls[1]->broadcastids[__privtrans->pdls[1]->nbroadcastids-1] = __privtrans->pdls[1]->ndims;
  #line 91 "lib/PDL/Slices-pp-broadcastI.c"
  return PDL_err;
}


#line 1857 "lib/PDL/PP.pm"
pdl_error pdl_broadcastI_free(pdl_trans *__privtrans, char destroy) {
  pdl_error PDL_err = {0, NULL, 0};
#line 99 "lib/PDL/Slices-pp-broadcastI.c"
  pdl_params_broadcastI *__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_broadcastI_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_broadcastI_vtable_realdims[] = { 0, 0 };
static char *pdl_broadcastI_vtable_parnames[] = { "PARENT","CHILD" };
static short pdl_broadcastI_vtable_parflags[] = {
  0,
  PDL_PARAM_ISCREAT|PDL_PARAM_ISCREATEALWAYS|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
};
static pdl_datatypes pdl_broadcastI_vtable_partypes[] = { -1, -1 };
static PDL_Indx pdl_broadcastI_vtable_realdims_starts[] = { 0, 0 };
static PDL_Indx pdl_broadcastI_vtable_realdims_ind_ids[] = { 0 };
static char *pdl_broadcastI_vtable_indnames[] = { "" };
pdl_transvtable pdl_broadcastI_vtable = {
  0, PDL_ITRANS_ISAFFINE|PDL_ITRANS_TWOWAY|PDL_ITRANS_DO_DATAFLOW_ANY, pdl_broadcastI_vtable_gentypes, 1, 2, NULL /*CORE21*/,
  pdl_broadcastI_vtable_realdims, pdl_broadcastI_vtable_parnames,
  pdl_broadcastI_vtable_parflags, pdl_broadcastI_vtable_partypes,
  pdl_broadcastI_vtable_realdims_starts, pdl_broadcastI_vtable_realdims_ind_ids, 0,
  0, pdl_broadcastI_vtable_indnames,
  pdl_broadcastI_redodims, NULL, NULL,
  pdl_broadcastI_free,
  sizeof(pdl_params_broadcastI),"PDL::Slices::broadcastI"
};


pdl_error pdl_run_broadcastI(pdl  *PARENT,pdl  *CHILD,PDL_Indx  id,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_broadcastI_vtable);
  if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
  pdl_params_broadcastI *__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->id) = (id); /* CType.get_copy */
  __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 1865 "lib/PDL/Slices.pd"
                PDL_Indx i,j;
                __params->nrealwhichdims = 0;
                for(i=0; i<__params->whichdims_count; i++) {
                        for(j=i+1; j<__params->whichdims_count; j++)
                                if(__params->whichdims[i] == __params->whichdims[j] &&
                                   __params->whichdims[i] != -1) {
                                return PDL->make_error(PDL_EUSERERROR, "Error in broadcastI:" "duplicate arg %"IND_FLAG" %"IND_FLAG" %"IND_FLAG"",i,j,__params->whichdims[i]);
                        }
                        if(__params->whichdims[i] != -1) {
                                __params->nrealwhichdims ++;
                        }
                }
#line 165 "lib/PDL/Slices-pp-broadcastI.c"
  PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
  return PDL_err;
}