File Coverage

lib/PDL/Stats/Kmeans-pp-_d_p2l.c

Criterion	Covered	Total	%
statement	43	75	57.3
branch	45	380	11.8
condition			n/a
subroutine			n/a
pod			n/a
total	88	455	19.3

#line 453 "lib/PDL/PP.pm"
/*
 * THIS FILE WAS GENERATED BY PDL::PP from lib/PDL/Stats/Kmeans.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_Stats_Kmeans
extern Core* PDL; /* Structure hold core C functions */
#line 23 "lib/PDL/Stats/Kmeans-pp-_d_p2l.c"

#include
#include
#include


#line 1857 "lib/PDL/PP.pm"
pdl_error pdl__d_p2l_readdata(pdl_trans *__privtrans) {
  pdl_error PDL_err = {0, NULL, 0};
#line 33 "lib/PDL/Stats/Kmeans-pp-_d_p2l.c"
  if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in _d_p2l:" "broadcast.incs NULL");
  /* broadcastloop declarations */
  int __brcloopval;
  register PDL_Indx __tind0,__tind1; /* counters along dim */
  register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
  /* dims here are how many steps along those dims */
  register PDL_Indx __tinc0_xc = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
  register PDL_Indx __tinc0_yc = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
  register PDL_Indx __tinc0_xa = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,0);
  register PDL_Indx __tinc0_ya = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,3,0);
  register PDL_Indx __tinc0_xb = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,4,0);
  register PDL_Indx __tinc0_yb = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,5,0);
  register PDL_Indx __tinc0_d = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,6,0);
  register PDL_Indx __tinc1_xc = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
  register PDL_Indx __tinc1_yc = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
  register PDL_Indx __tinc1_xa = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,1);
  register PDL_Indx __tinc1_ya = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,3,1);
  register PDL_Indx __tinc1_xb = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,4,1);
  register PDL_Indx __tinc1_yb = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,5,1);
  register PDL_Indx __tinc1_d = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,6,1);
  #define PDL_BROADCASTLOOP_START__d_p2l_readdata PDL_BROADCASTLOOP_START( \
    readdata, \
    __privtrans->broadcast, \
    __privtrans->vtable, \
    xc_datap += __offsp[0]; \
    yc_datap += __offsp[1]; \
    xa_datap += __offsp[2]; \
    ya_datap += __offsp[3]; \
    xb_datap += __offsp[4]; \
    yb_datap += __offsp[5]; \
    d_datap += __offsp[6]; \
    , \
    (  ,xc_datap += __tinc1_xc - __tinc0_xc * __tdims0 \
    ,yc_datap += __tinc1_yc - __tinc0_yc * __tdims0 \
    ,xa_datap += __tinc1_xa - __tinc0_xa * __tdims0 \
    ,ya_datap += __tinc1_ya - __tinc0_ya * __tdims0 \
    ,xb_datap += __tinc1_xb - __tinc0_xb * __tdims0 \
    ,yb_datap += __tinc1_yb - __tinc0_yb * __tdims0 \
    ,d_datap += __tinc1_d - __tinc0_d * __tdims0 \
    ), \
    (  ,xc_datap += __tinc0_xc \
    ,yc_datap += __tinc0_yc \
    ,xa_datap += __tinc0_xa \
    ,ya_datap += __tinc0_ya \
    ,xb_datap += __tinc0_xb \
    ,yb_datap += __tinc0_yb \
    ,d_datap += __tinc0_d \
    ) \
  )
  #define PDL_BROADCASTLOOP_END__d_p2l_readdata PDL_BROADCASTLOOP_END( \
    __privtrans->broadcast, \
    xc_datap -= __tinc1_xc * __tdims1 + __offsp[0]; \
    yc_datap -= __tinc1_yc * __tdims1 + __offsp[1]; \
    xa_datap -= __tinc1_xa * __tdims1 + __offsp[2]; \
    ya_datap -= __tinc1_ya * __tdims1 + __offsp[3]; \
    xb_datap -= __tinc1_xb * __tdims1 + __offsp[4]; \
    yb_datap -= __tinc1_yb * __tdims1 + __offsp[5]; \
    d_datap -= __tinc1_d * __tdims1 + __offsp[6]; \
  )
  #ifndef PDL_DECLARE_PARAMS__d_p2l_1
  #define PDL_DECLARE_PARAMS__d_p2l_1(PDL_TYPE_OP,PDL_PPSYM_OP) \
    PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, xc, (__privtrans->pdls[0]), 1, PDL_PPSYM_OP) \
    PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, yc, (__privtrans->pdls[1]), 1, PDL_PPSYM_OP) \
    PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, xa, (__privtrans->pdls[2]), 1, PDL_PPSYM_OP) \
    PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, ya, (__privtrans->pdls[3]), 1, PDL_PPSYM_OP) \
    PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, xb, (__privtrans->pdls[4]), 1, PDL_PPSYM_OP) \
    PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, yb, (__privtrans->pdls[5]), 1, PDL_PPSYM_OP) \
    PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, d, (__privtrans->pdls[6]), 1, PDL_PPSYM_OP)
  #endif
  if ( __privtrans->bvalflag ) { /* ** do 'bad' Code ** */
    #define PDL_BAD_CODE
    #define PDL_IF_BAD(t,f) t
    switch (__privtrans->__datatype) { /* Start generic switch */
      case PDL_F: {
        PDL_DECLARE_PARAMS__d_p2l_1(PDL_Float,F)
        PDL_BROADCASTLOOP_START__d_p2l_readdata    {
            PDL_IF_BAD(if (PDL_ISBAD2((xc_datap)[0],xc_badval,F,xc_badval_isnan) || PDL_ISBAD2((yc_datap)[0],yc_badval,F,yc_badval_isnan) || PDL_ISBAD2((xa_datap)[0],xa_badval,F,xa_badval_isnan) || PDL_ISBAD2((ya_datap)[0],ya_badval,F,ya_badval_isnan) || PDL_ISBAD2((xb_datap)[0],xb_badval,F,xb_badval_isnan) || PDL_ISBAD2((yb_datap)[0],yb_badval,F,yb_badval_isnan) ) { (d_datap)[0]=d_badval; continue; },)
            PDL_Float xba, yba;
            xba = (xb_datap)[0] - (xa_datap)[0];
            yba = (yb_datap)[0] - (ya_datap)[0];
            (d_datap)[0] = fabs( xba * ((ya_datap)[0] - (yc_datap)[0]) - ((xa_datap)[0] - (xc_datap)[0]) * yba )
                 / hypot( xba, yba );
              }PDL_BROADCASTLOOP_END__d_p2l_readdata
      } break;
      case PDL_D: {
        PDL_DECLARE_PARAMS__d_p2l_1(PDL_Double,D)
        PDL_BROADCASTLOOP_START__d_p2l_readdata    {
            PDL_IF_BAD(if (PDL_ISBAD2((xc_datap)[0],xc_badval,D,xc_badval_isnan) || PDL_ISBAD2((yc_datap)[0],yc_badval,D,yc_badval_isnan) || PDL_ISBAD2((xa_datap)[0],xa_badval,D,xa_badval_isnan) || PDL_ISBAD2((ya_datap)[0],ya_badval,D,ya_badval_isnan) || PDL_ISBAD2((xb_datap)[0],xb_badval,D,xb_badval_isnan) || PDL_ISBAD2((yb_datap)[0],yb_badval,D,yb_badval_isnan) ) { (d_datap)[0]=d_badval; continue; },)
            PDL_Double xba, yba;
            xba = (xb_datap)[0] - (xa_datap)[0];
            yba = (yb_datap)[0] - (ya_datap)[0];
            (d_datap)[0] = fabs( xba * ((ya_datap)[0] - (yc_datap)[0]) - ((xa_datap)[0] - (xc_datap)[0]) * yba )
                 / hypot( xba, yba );
              }PDL_BROADCASTLOOP_END__d_p2l_readdata
      } break;
      default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in _d_p2l: unhandled datatype(%d), only handles (FD)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
    }
    #undef PDL_BAD_CODE
    #undef PDL_IF_BAD
  } else { /* ** else do 'good' Code ** */
    #define PDL_IF_BAD(t,f) f
    switch (__privtrans->__datatype) { /* Start generic switch */
      case PDL_F: {
        PDL_DECLARE_PARAMS__d_p2l_1(PDL_Float,F)
        PDL_BROADCASTLOOP_START__d_p2l_readdata    {
            PDL_IF_BAD(if (PDL_ISBAD2((xc_datap)[0],xc_badval,F,xc_badval_isnan) || PDL_ISBAD2((yc_datap)[0],yc_badval,F,yc_badval_isnan) || PDL_ISBAD2((xa_datap)[0],xa_badval,F,xa_badval_isnan) || PDL_ISBAD2((ya_datap)[0],ya_badval,F,ya_badval_isnan) || PDL_ISBAD2((xb_datap)[0],xb_badval,F,xb_badval_isnan) || PDL_ISBAD2((yb_datap)[0],yb_badval,F,yb_badval_isnan) ) { (d_datap)[0]=d_badval; continue; },)
            PDL_Float xba, yba;
            xba = (xb_datap)[0] - (xa_datap)[0];
            yba = (yb_datap)[0] - (ya_datap)[0];
            (d_datap)[0] = fabs( xba * ((ya_datap)[0] - (yc_datap)[0]) - ((xa_datap)[0] - (xc_datap)[0]) * yba )
                 / hypot( xba, yba );
              }PDL_BROADCASTLOOP_END__d_p2l_readdata
      } break;
      case PDL_D: {
        PDL_DECLARE_PARAMS__d_p2l_1(PDL_Double,D)
        PDL_BROADCASTLOOP_START__d_p2l_readdata    {
            PDL_IF_BAD(if (PDL_ISBAD2((xc_datap)[0],xc_badval,D,xc_badval_isnan) || PDL_ISBAD2((yc_datap)[0],yc_badval,D,yc_badval_isnan) || PDL_ISBAD2((xa_datap)[0],xa_badval,D,xa_badval_isnan) || PDL_ISBAD2((ya_datap)[0],ya_badval,D,ya_badval_isnan) || PDL_ISBAD2((xb_datap)[0],xb_badval,D,xb_badval_isnan) || PDL_ISBAD2((yb_datap)[0],yb_badval,D,yb_badval_isnan) ) { (d_datap)[0]=d_badval; continue; },)
            PDL_Double xba, yba;
            xba = (xb_datap)[0] - (xa_datap)[0];
            yba = (yb_datap)[0] - (ya_datap)[0];
            (d_datap)[0] = fabs( xba * ((ya_datap)[0] - (yc_datap)[0]) - ((xa_datap)[0] - (xc_datap)[0]) * yba )
                 / hypot( xba, yba );
              }PDL_BROADCASTLOOP_END__d_p2l_readdata
      } break;
      default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in _d_p2l: unhandled datatype(%d), only handles (FD)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
    }
    #undef PDL_IF_BAD
  }
  return PDL_err;
}

static pdl_datatypes pdl__d_p2l_vtable_gentypes[] = { PDL_F, PDL_D, -1 };
static PDL_Indx pdl__d_p2l_vtable_realdims[] = { 0, 0, 0, 0, 0, 0, 0 };
static char *pdl__d_p2l_vtable_parnames[] = { "xc","yc","xa","ya","xb","yb","d" };
static short pdl__d_p2l_vtable_parflags[] = {
  0,
  0,
  0,
  0,
  0,
  0,
  PDL_PARAM_ISCREAT|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
};
static pdl_datatypes pdl__d_p2l_vtable_partypes[] = { -1, -1, -1, -1, -1, -1, -1 };
static PDL_Indx pdl__d_p2l_vtable_realdims_starts[] = { 0, 0, 0, 0, 0, 0, 0 };
static PDL_Indx pdl__d_p2l_vtable_realdims_ind_ids[] = { 0 };
static char *pdl__d_p2l_vtable_indnames[] = { "" };
pdl_transvtable pdl__d_p2l_vtable = {
  PDL_TRANS_DO_BROADCAST|PDL_TRANS_BADPROCESS, 0, pdl__d_p2l_vtable_gentypes, 6, 7, NULL /*CORE21*/,
  pdl__d_p2l_vtable_realdims, pdl__d_p2l_vtable_parnames,
  pdl__d_p2l_vtable_parflags, pdl__d_p2l_vtable_partypes,
  pdl__d_p2l_vtable_realdims_starts, pdl__d_p2l_vtable_realdims_ind_ids, 0,
  0, pdl__d_p2l_vtable_indnames,
  NULL, pdl__d_p2l_readdata, NULL,
  NULL,
  0,"PDL::Stats::Kmeans::_d_p2l"
};


pdl_error pdl_run__d_p2l(pdl  *xc,pdl  *yc,pdl  *xa,pdl  *ya,pdl  *xb,pdl  *yb,pdl  *d) {
  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__d_p2l_vtable);
  if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
  __privtrans->pdls[0] = xc;
  __privtrans->pdls[1] = yc;
  __privtrans->pdls[2] = xa;
  __privtrans->pdls[3] = ya;
  __privtrans->pdls[4] = xb;
  __privtrans->pdls[5] = yb;
  __privtrans->pdls[6] = d;
  PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
  PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
  return PDL_err;
}