use Array: all;
use SimplePrint: all;
use SimpleFibre: all;
use Math: all;
use String:{tochar};
import Char:all;
use Bits:all;

/*
% This is the APEX stdlib.sis include file.
% Standard equates and constants for APL compiler
% Also standard coercion functions
*/

#define toB(x) tob((x))
#define toI(x) toi((x))
#define toD(x) tod((x))
#define toC(x) (x)
#define toc(x) ((x))

void APEXERROR(char[.] msg)
{
/* Error function. This needs work, as it should kill
 * the running task, too. 
 */
 /*print(msg); */
return();
}

/* Structural function utility functions */
/* Ravel utility */
#define APEXRavel(y) (reshape([prod(shape(y))],y))

inline int VectorRotateAmount(int x, int y)
{ /* Normalize x rotate for array of shape y on selected axis */
 /* normalize rotation count */
 if (x>0)
        z = _mod_(x,y);
 else
        z = y - _mod_(_abs_(x),y);
 return(z);
}

#define APEXRESHAPE(TYPE)					\
inline TYPE[*] APEXReshape(int[.] x, TYPE[*] y)			\
{ /* APEX vector x reshape, with item reuse */			\
  ry = APEXRavel(y);						\
  zxrho = prod(x); /* THIS NEEDS XRHO FOR CODE SAFETY!! */	\
  yxrho = shape(ry)[0];						\
  if( zxrho <= yxrho) { /* No element resuse case */		\
        z = take([zxrho],ry);					\
 } else {							\
        ncopies = zxrho/yxrho; /* # complete copies of y. */	\
        /* FIXME: y empty case !*/				\
        z = with(. <= [i] <= .)					\
                genarray( [ncopies], y,y);			\
        /* Now append the leftover bits */			\
        z = APEXRavel(z) ++ take([zxrho-(ncopies*yxrho)],ry);	\
 }								\
 return(reshape(x,z));						\
}

APEXRESHAPE(bool)
APEXRESHAPE(int)
APEXRESHAPE(double)
APEXRESHAPE(char)

inline bool SameShape(int[*] x, int[*] y)
{ /* Predicate for two shape vectors having same shape */
 if (dim(x) != dim(y))
	z = false;
 else
	z = with ([0] <= i < [dim(y)])
		fold(&, true,
		(shape(x))[i] == (shape(y))[i]);
 return(z);
}

#define APEXmiota(N) with (. <= [iv] <= .) genarray([N], iv)


#define APEXPI 3.1415926535897932d
#define APEXE  2.718281828d
/*  APEXPINFINITYI largest integer */
#define APEXPINFINITYD  1.7976931348623156D308
/*  APEXMINFINITYI smallest integer */
#define APEXMINFINITYD -1.7976931348623156D308


/*
% Floating-point utilities
% This taken from page 93 of the 1993-01-06 version 
% of Committee Draft 1 of the Extended ISO APL Standard
% NO {quad}ct support yet! 
*/

inline int Dsignum(double y)
{ /* signum double */
 if (0.0 == y)
	z = 0;
 else if (0.0 > y)
	z = -1;
 else 
	z = 1;
 return(z);
}

inline int Isignum(int y)
{ /* signum int */
 if (0 == y)
	z = 0;
 else if (0 > y)
	z = -1;
 else 
	z = 1;
 return(z);
}

inline double Dresidue(double x, double y, double QUADct)
{ /* Double residue double */
  /* See Iresidue for definition */
  if (0.0 == x)
	nx = 1.0;
  else 
	nx = x;
  z = y - x * tod(Dfloor(y/nx, QUADct));
  return(z);
}


inline int Iresidue(int x, int y, double ignoreme)
{ /* Integer residue integer */
  /* Third argument is ignored - QUADct for Dresidue */
  /* This definition is taken from SHARP APL Refman May 1991, p.6-26.
   * It extends the definition of residue to fractional right arguments
   * and to zero, negative and fractional left arguments. 
   * r= y-x times floor y divide x+0=x 
   */
  if (0 == x)
	nx = 1;
  else 
	nx = x;
  z = y - x * (y/nx);
  return(z);
}

inline int DfloorNoFuzz(double y)
{ /* Exact floor (no fuzz) */
 return(toi(floor(y)));
}

inline int Dfloor(double y, double QUADct)
{ /* Fuzzy floor */
  /* Definition taken from SHARP APL Refman May 1991, p.6-23
   * floor:  n <- (signum y) times nofuzzfloor 0.5+abs y)
   *         z <- n-(QUADct times 1 max abs y)<(n-y)
   *
   * If you want a double result,  write: "y - 1| y".
   *
   */
   n = tod(DfloorNoFuzz(0.5+fabs(y)));
   if (y < 0.0)
	n = -n;
   else if (0.0 == y)
	n = 0.0;
   range = fabs(y);
   if (1.0 > range)
	range = 1.0;
   fuzzlim = QUADct*range;
   ny = n-y;
   if (fuzzlim < ny)
	z = n - 1.0;
   else
	z = n;
   return(toi(z)); 
}

inline bool APEXFUZZEQ(double x, double y, double QUADct)
{ /* ISO APL Tolerant equality predicate */
 absx = abs(x);
 absy = abs(y);
 tolerance = QUADct * max(absx,absy);
 z = abs(x-y) <= tolerance;
 return(z);
}

/*
% 1996-12-07 Pearl Harbor Day. Try to fix bug in DBank infdivm 
% benchmark. (4 5 rho 6) + rank 1 (4 1 rho 7) introduces singleton
% vectors into scalar fn calls.
% Hence, we need support for singletons within scalar fns.
% There is an similar, but independent, bug in rank support for
% the extension case.
*/


/* vector-scalar simple search loops 
 * This is origin-0 x1 iota y0
 * cfn is comparator function type suffix to use, e.g, b,i,d,c,
 */

#define FindFirst(x,y,cfn,QUADct)			\
 sx = shape(x)[0];					\
 z = sx; /* if not found */				\
 for(i=0; i<sx; i++)					\
  if (comparatorfn##cfn(to##cfn(x[i]),to##cfn(y),QUADct)) \
   { z = i;						\
     i = sx; 						\
   }

int FindFirstB(bool[.] x, bool y, double QUADct)
{
 /* x iota y (vector-scalar on Booleans */
 FindFirst(x,y,b,QUADct)
 return(z);
}

int FindFirstI(int[.] x, int y, double QUADct)
{
 /* x iota y (vector-scalar on integers */
 FindFirst(x,y,i,QUADct)
 return(z);
}

int FindFirstD(double[.] x, double y,double QUADct)
{
 /* x iota y (vector-scalar on doubles */
 FindFirst(x,y,d,QUADct)
 return(z);
}

int FindFirstC(char[.] x, char y, double QUADct)
{
 /* x iota y (vector-scalar on characters */
 FindFirst(x,y,c,QUADct)
 return(z);
}


inline bool equals(char x, char y,double QUADct)
{ /* Char comparator */
 z = (x == y);
 return(z);
}

inline bool equals(char[+] x, char[+] y, double QUADct)
{ /* Char item comparator */
  /* This definition is sort of flakey - it really compares ravels... */
  /* It oughta also ensure that the arguments shapes match... */
  z = with(0*shape(x) <= iv < shape(x))
	fold(&, true, x[iv] == y[iv]);
 return(z);
}

inline bool lessthan(char[+] x, char[+] y, double QUADct)
{ /* Char item comparator */
  /* It oughta also ensure that the arguments shapes match... */
  shp=shape(x)[0];
  z=false;
  for (i=0;i<shp;i++)
        if (x[i]!=y[i]){
                z= x[i]<y[i];
                i=shp;
        }
 return(z);
}

inline bool greaterthan(char[+] x, char[+] y, double QUADct)
{ /* Char item comparator */
  /* This definition is sort of flakey - it really compares ravels... */
  /* It oughta also ensure that the arguments shapes match... */
  shp=shape(x)[0];
  z=false;
  for (i=0;i<shp;i++)
        if (x[i]!=y[i]){
                z= x[i]>y[i];
                i=shp;
        }
 return(z);
}

inline bool comparatorfnb(bool x, bool y,double QUADct)
{ /* Boolean comparator */
 z = (x == y);
 return(z);
}

inline bool comparatorfni(int x, int y, double QUADct)
{ /* integer comparator */
 z = (x == y);
 return(z);
}

inline bool comparatorfnd(double x, double y, double QUADct)
{ /* double comparator with fuzz */
 L = _abs_(x-y);  /* Tolerant equality */
 R = QUADct*max(_abs_(x),_abs_(y));
 z = (L <= R);
 return(z);
}


inline bool comparatorfnc(char x, char y,double QUADct)
{ /* char comparator */
 z = (x == y);
 return(z);
}

inline bool comparatorfndnf(double x,double y, double QUADct)
{ /* double comparator, no fuzz */
 z = (x == y);
 return(z);
}

/*  Lehmer's random number generator */
#define lehmer(qrl) mod(qrl*16807,2147483647)


/* Transposes */
#define APEXTRANSPOSE(TYPE,OTFILL)			\
inline TYPE[*] APEXTranspose(TYPE[*] y)			\
{							\
  z = with(iv)						\
	( . <= iv <= .) : y[reverse( iv)];		\
        genarray( reverse( shape(y)), OTFILL);		\
  return(z);						\
}

APEXTRANSPOSE(bool,false)
APEXTRANSPOSE(int,0)
APEXTRANSPOSE(double,0.0)
APEXTRANSPOSE(char,' ')
/* End of transposes */

/* End of boilerplate */


inline bool[+] notXBB(bool[+] y)
{ /* Monadic scalar functions on array */
 z = with( . <= iv <= .)
 genarray(shape(y),
 notXBB(toB(y[iv])),false);
 return(z);
}
inline bool notXBB(bool y)
{ return(!tob(y));
}

inline bool eqBIB(bool x, int y)
{ return(toI(x) == toI(y));
}

inline int mpyIII(int x, int y)
{ return(toI(x)*toI(y));
}

inline int plusIII(int x, int y)
{ return(toI(x)+toI(y));
}

inline double divDID(double x, int y)
{ dx = tod(x);
 dy = tod(y);
 if (dx == dy)
 z = 1.0d;
 else
 z = dx/dy;
 return(z);
}

inline double modIDD(int x, double y, double QUADct)
{ return(Dresidue(toD(x),toD(y),QUADct));
}

inline double mpyDID(double x, int y)
{ return(toD(x)*toD(y));
}

inline double barDDD(double x, double y)
{ return(toD(x)-toD(y));
}

inline double divDBD(double x, bool y)
{ dx = tod(x);
 dy = tod(y);
 if (dx == dy)
 z = 1.0d;
 else
 z = dx/dy;
 return(z);
}

inline double modBDD(bool x, double y, double QUADct)
{ return(Dresidue(toD(x),toD(y),QUADct));
}

inline double mpyDBD(double x, bool y)
{ return(toD(x)*toD(y));
}

inline int plusBBI(bool x, bool y)
{ return(toI(x)+toI(y));
}

inline int[.] comaXII(int y)
{
 return([y]);
}

inline int[.] rotrXII(int[.] y)
{ /* Vector reverse */
 n = shape(y);
 cell = 0;
 z = with( . <= iv <= .)
 genarray(n, y[(n[0]-1)-iv],cell);
 return(z);
}

inline int[.] comaXII(int[.] y)
{
 return(y);
}


inline int[.] comaXII(int[+] y)
{ /* Ravel of anything of rank>1 */
 z = reshape([prod(shape(y))],y);
 return(z);
}

inline double[.] comaXDD(double y)
{
 return([y]);
}

inline bool[.] comaXBB(bool[.] y)
{
 return(y);
}


inline bool[.] rotrXBB(bool[.] y)
{ /* Vector reverse */
 n = shape(y);
 cell = false;
 z = with( . <= iv <= .)
 genarray(n, y[(n[0]-1)-iv],cell);
 return(z);
}

inline int[+] rhoIII(int[.] x, int y)
{ /* Vector reshape scalar to matrix) */
 zxrho = prod(toi(x)); /* Result element count */
 z = genarray([zxrho], y); /* allocate result */
 z = reshape(toi(x),z);
 return(z);
}


inline int[.] rhoIII(int x, int y)
{ /* Scalar reshape scalar to vector) */
 zxrho = prod(toi(x)); /* Result element count */
 z = genarray([zxrho], y); /* allocate result */
 return(z);
}

inline double[*] rhoIDD(int[.] x, double[*] y)
{
 z = APEXReshape(toi(x),y);
 return(z);
}


inline int[*] rhoIII(int[.] x, int[*] y)
{
 z = APEXReshape(toi(x),y);
 return(z);
}


inline int[*] iotaCCI(char[.] x, char[*] y,int QUADio)
{ /* Character vector iota character non-scalar */
 table = genarray([256],shape(x)[0]); /* Not found */
 for(i=shape(x)[0]-1; i<=0; i--)
 table[toi(x[i])] = i;
 z = with(. <= iv <= .)
 genarray(shape(y), table[toi(y[iv])],0);
 return(z+QUADio);
}

inline int[.] takeBII(bool x, int[.] y)
/* LOTS of missing code here. (x0<0 for example */
{ /* Scalar take vector */
 z = with( . <= iv <= .)
 genarray([abs(toi(x))], y[iv]);
 return(z);
}

inline int[*] dropBII(bool x, int[*] y)
{ /* Scalar drop non-scalar */
 res = drop([toi(x)], y);
 return( res);
}

inline bool[*] dropBBB(bool x, bool[*] y)
{ /* Scalar drop non-scalar */
 res = drop([toi(x)], y);
 return( res);
}

inline int[0] rhoXII(int y)
{ /* Shape of scalar */
 return(shape(y));
}

inline int[1] rhoXBI(bool[.] y)
{ /* Shape of vector */
 return(shape(y));
}

inline int[1] rhoXII(int[.] y)
{ /* Shape of vector */
 return(shape(y));
}

inline int[.] iotaXII(int shp, int QUADio)
{ /* Index generator on scalar */
/* HELP! Needs domain check for negative shp */
 res = with( . <= [i] <= .)
 genarray( [toI(shp)], i+QUADio);
 return( res);
}
inline int[.] iotaXIINonNegsy(int shp, int QUADio)
{ /* Index generator on ScalarN when N is non-negative integer */
 res = with( . <= [i] <= .)
 genarray( [toI(shp)], i+QUADio);
 return( res);
}

inline int[.] iotaXII(int[1] shp, int QUADio)
{ /* Index generator on 1-element vector */
 /* HELP! Needs length error check */
/* HELP! Needs domain check for negative shp */
 res = with( . <= i <= .)
 genarray( toI(shp), i[0]+QUADio);
 return( res);
}

inline int[.] iotaXIINonNegsy(int[1] shp, int QUADio)
{ /* Index generator on 1-element vector, known to be non-negative integer */
 res = with( . <= [i] <= .)
 genarray( [toI(shp[0])], i+QUADio);
 return( res);
}

inline int[.] rhoXII(int[+] y)
{ /* Shape of rank-3 or higher */
 return(shape(y));
}

inline int[0] rhoXDI(double y)
{ /* Shape of scalar */
 return(shape(y));
}

int quadXBB(bool[*] y, int QUADpp, int QUADpw)
{ /* {quad}{<-} anything */
 r=display(y); /* KLUDGE!!! print doesn't support booleans! */
/* MORE kludges - sac tends to evaporate side effects,
 so we'll help things along abit until that is repaired. */
 return(r);
}

int quadXII(int[*] y, int QUADpp, int QUADpw)
{ /* {quad}{<-} anything */
 r=display(y); /* KLUDGE!!! print doesn't support booleans! */
/* MORE kludges - sac tends to evaporate side effects,
 so we'll help things along abit until that is repaired. */
 return(r);
}

inline int[.] comaIII(int[.] x, int[.] y)
{ /* VxV last-axis catenate */
 return(toI(x)++toI(y));
}

inline int[.] comaIBI(int[.] x, bool y)
{/* VxS catenate first (or last) axis */
 return(toI(x)++[toI(y)]);
}

inline bool[.] comaBBB(bool[.] x, bool y)
{/* VxS catenate first (or last) axis */
 return(toB(x)++[toB(y)]);
}

inline bool[2] comaBBB(bool x, bool y)
{/* SxS catenate first (or last) axis */
 return([toB(x)]++[toB(y)]);
}

inline int[*] dtakIBI(int x, bool[+] y)
{ /* Scalar basevalue rank>1 */
 yt = APEXTranspose(y); /* Dumb, but easy */
 ycols = (shape(y))[dim(y)-1];
 cellshape=take([1],shape(y));
 cell = 0;
 frameshape = drop([1],shape(y));
 weights = genarray(cellshape, toI(1));
 for (i=ycols-2; i>=0; i--)
 weights[i]= weights[i+1]*toI(x);
 z = with(. <= iv <= .)
 genarray(frameshape,
 with([0] <= jv < cellshape)
 fold(+, 0,
 weights[jv] * toI(yt[iv++jv])),cell);
 return(z);
}

inline bool sameIIB(int x, int y)
{ /* Scalar match scalar */
 if (toI(x) == toI(y))
 z = true;
 else
 z = false;
 return(z);
}

inline int[+] utakIII(int[.] x, int[+] y)
{ /* Vector represent non-scalar */
 yt = APEXTranspose(y);
 cell = genarray(shape(x),0);
 z = with (. <= iv <= .)
 genarray(shape(yt),utakIII(x,yt[iv]),cell);
 return(z);
}


inline int[.] utakIII(int[.] x, int y)
{ /* Vector represent scalar */
 /* Taken from ISO Extended APL standard Draft N93.03, page 155 */
 wts = genarray(shape(x),toI(1));
 for(i=shape(x)[0]-2; i>=0; i--)
 wts[i] = wts[i+1] * toI(x[i+1]);

 z = genarray(shape(x),0);
 cy = toI(y);
 for(i=shape(x)[0]-1; i>=0; i--){
 z[i] = Iresidue(toI(x[i]),cy/wts[i],0.0);
 /* Represent is NOT fuzzy (SAPL Ref Man p.6-47, 1991 */
 cy = cy - z[i] * wts[i];
 }
 return(z);
}

inline bool sameIDB(int[+] x, double[+] y,double QUADct)
{ /* Non-scalar match non-scalar */
 if (dim(x) != dim(y))
 z = false;
 else
 if (! SameShape(shape(x),shape(y)))
 z = false;
 else
 z = with(_mul_SxA_(0,shape(y)) <= iv < shape(y))
 fold(&, true,
 eqDDB(toD(x[iv]), toD(y[iv]), QUADct));
 return(z);
}

inline bool eqIDB(int x, double y, double QUADct)
{ return((toD(x) == toD(y)) || APEXFUZZEQ(toD(x),toD(y),QUADct));
}

inline int dtakIII(int x, int y)
{ /* Scalar basevalue Scalar */
 return(y);
}

inline int[*] dtakIII(int x, int[+] y)
{ /* Scalar basevalue rank>1 */
 yt = APEXTranspose(y); /* Dumb, but easy */
 ycols = (shape(y))[dim(y)-1];
 cellshape=take([1],shape(y));
 cell = 0;
 frameshape = drop([1],shape(y));
 weights = genarray(cellshape, toI(1));
 for (i=ycols-2; i>=0; i--)
 weights[i]= weights[i+1]*toI(x);
 z = with(. <= iv <= .)
 genarray(frameshape,
 with([0] <= jv < cellshape)
 fold(+, 0,
 weights[jv] * toI(yt[iv++jv])),cell);
 return(z);
}

inline int dtakIII(int[.] x, int[.] y)
{ /* Vector basevalue vector */
 /* 3 cases - all give 22200:
 * 10 10 10 basevalue 200 200 200
 * 10 10 10 basevalue 200
 * (,10) basevalue 200 200 200
 */
 ycols = (shape(y))[0];
 if (1 == ycols){ /* Maybe extend y */
 ycols = shape(x)[0];
 y = genarray([ycols],y[0]);
 }
 if (1 == shape(x)[0]){ /* Maybe extend x */
 x = genarray([ycols], x[0]);
 }
 weights = genarray([ycols], toI(1));
 for (i=ycols-2; i>=0; i--)
 weights[i]= weights[i+1]*toI(x[i]);
 z = with([0] <= iv < [ycols])
 fold(+, 0,
 weights[iv] * toI(y[iv]));
 return(z);
}

inline int[*] dtakIII(int[.] x, int[+] y)
{ /* Vector basevalue rank>1 */
 yt = APEXTranspose(y); /* Dumb, but easy */
 ycols = (shape(y))[dim(y)-1];
 cellshape=take([1],shape(y));
 cell=genarray(cellshape,0);
 frameshape = drop([1],shape(y));
 weights = genarray(cellshape, toI(1));
 for (i=ycols-2; i>=0; i--)
 weights[i]= weights[i+1]*toI(x);
 z = with(. <= iv <= .)
 genarray(frameshape,
 with([0] <= jv < cellshape)
 fold(+, 0,
 weights[jv] * toI(yt[iv++jv])),cell);
 return(z);
}





inline bool sameIIB(int[+] x, int[+] y)
{ /* Non-scalar match non-scalar */
 if (dim(x) != dim(y))
 z = false;
 else
 if (! SameShape(shape(x),shape(y)))
 z = false;
 else
 z = with(_mul_SxA_(0,shape(y)) <= iv < shape(y))
 fold(&, true,
 eqIIB(toI(x[iv]), toI(y[iv])));
 return(z);
}

inline bool eqIIB(int x, int y)
{ return(toI(x) == toI(y));
}

inline int[.] utakIII(int[.] x, int y)
{ /* Vector represent scalar */
 /* Taken from ISO Extended APL standard Draft N93.03, page 155 */
 wts = genarray(shape(x),toI(1));
 for(i=shape(x)[0]-2; i>=0; i--)
 wts[i] = wts[i+1] * toI(x[i+1]);

 z = genarray(shape(x),0);
 cy = toI(y);
 for(i=shape(x)[0]-1; i>=0; i--){
 z[i] = Iresidue(toI(x[i]),cy/wts[i],0.0);
 /* Represent is NOT fuzzy (SAPL Ref Man p.6-47, 1991 */
 cy = cy - z[i] * wts[i];
 }
 return(z);
}

inline bool sameIDB(int[.] x, double[.] y,double QUADct)
{ /* vector-vector match */
 if (dim(x) != dim(y))
 z = false;
 else
 z = with(_mul_SxA_(0,shape(y)) <= iv < shape(y))
 fold(&, true,
 eqDDB(toD(x[iv]), toD(y[iv]), QUADct));
 return(z);
}

inline bool eqIDB(int x, double y, double QUADct)
{ return((toD(x) == toD(y)) || APEXFUZZEQ(toD(x),toD(y),QUADct));
}

inline double[.] utakIDD(int[.] x, double y)
{ /* Vector represent scalar */
 /* Taken from ISO Extended APL standard Draft N93.03, page 155 */
 wts = genarray(shape(x),toD(1));
 for(i=shape(x)[0]-2; i>=0; i--)
 wts[i] = wts[i+1] * toD(x[i+1]);

 z = genarray(shape(x),0.0d);
 cy = toD(y);
 for(i=shape(x)[0]-1; i>=0; i--){
 z[i] = Dresidue(toD(x[i]),cy/wts[i],0.0);
 /* Represent is NOT fuzzy (SAPL Ref Man p.6-47, 1991 */
 cy = cy - z[i] * wts[i];
 }
 return(z);
}

inline bool sameDDB(double[.] x, double[.] y,double QUADct)
{ /* vector-vector match */
 if (dim(x) != dim(y))
 z = false;
 else
 z = with(_mul_SxA_(0,shape(y)) <= iv < shape(y))
 fold(&, true,
 eqDDB(toD(x[iv]), toD(y[iv]), QUADct));
 return(z);
}

inline bool eqDDB(double x, double y, double QUADct)
{ return((toD(x) == toD(y)) || APEXFUZZEQ(toD(x),toD(y),QUADct));
}

inline double[.] utakBDD(bool[.] x, double y)
{ /* Vector represent scalar */
 /* Taken from ISO Extended APL standard Draft N93.03, page 155 */
 wts = genarray(shape(x),toD(1));
 for(i=shape(x)[0]-2; i>=0; i--)
 wts[i] = wts[i+1] * toD(x[i+1]);

 z = genarray(shape(x),0.0d);
 cy = toD(y);
 for(i=shape(x)[0]-1; i>=0; i--){
 z[i] = Dresidue(toD(x[i]),cy/wts[i],0.0);
 /* Represent is NOT fuzzy (SAPL Ref Man p.6-47, 1991 */
 cy = cy - z[i] * wts[i];
 }
 return(z);
}

inline int[*] indrfr(int fr, int i, int[+] X)
{ /* X[;;;scalari;;;], where i has fr semicolons to its left */
  /* Indexing is origin-0. Caller will correct this */
  /* This could stand some optimization, perhaps, for boolean i,
   * unless SAC avoids building an array-valued temp of toI(i).
   */
 frameshape = take([fr],shape(X));
 cellshape = drop([fr+1], shape(X));
 cell = genarray(cellshape,0);
 z = with (. <= iv <= .)
 genarray(frameshape,(X[iv])[i], cell);
return(z);
}

inline int[*] indrfr(int fr, int[+] i, int[+] X)
{ /* X[;;;i;;;], where i has fr semicolons to its left */
  /* Indexing is origin-0. Caller will correct this */
  /* This could stand some optimization, perhaps, for boolean i,
   * unless SAC avoids building an array-valued temp of toI(i).
   */
 frameshape = take([fr],shape(X));
 cellshape = shape(i)++drop([fr+1],shape(X));
 cell = genarray(cellshape,0); /* not used, but SAC needs help */
 z = with (. <= iv <= .)
 genarray(frameshape,indrfr0(i,X[iv]),cell);
return(z);
}

inline int[*] indrfr0(int i, int[*] X)
{ /* X[i;;;]  i is scalar */
 z = X[[i]];
 return(z);
}

inline int[*] indrfr0(int[+] i, int[*] X)
{ /* X[im;;;]  im is array */
 cellshape = drop([1],shape(X));
 cell = genarray(cellshape, 0);
 z = with (. <= iv <= .)
 genarray(shape(i),
  X[i[iv]],cell);
 return(z);
}

inline double[*] indrfr(int fr, int i, double[+] X)
{ /* X[;;;scalari;;;], where i has fr semicolons to its left */
  /* Indexing is origin-0. Caller will correct this */
  /* This could stand some optimization, perhaps, for boolean i,
   * unless SAC avoids building an array-valued temp of toI(i).
   */
 frameshape = take([fr],shape(X));
 cellshape = drop([fr+1], shape(X));
 cell = genarray(cellshape,0.0d);
 z = with (. <= iv <= .)
 genarray(frameshape,(X[iv])[i], cell);
return(z);
}

inline double[*] indrfr(int fr, int[+] i, double[+] X)
{ /* X[;;;i;;;], where i has fr semicolons to its left */
  /* Indexing is origin-0. Caller will correct this */
  /* This could stand some optimization, perhaps, for boolean i,
   * unless SAC avoids building an array-valued temp of toI(i).
   */
 frameshape = take([fr],shape(X));
 cellshape = shape(i)++drop([fr+1],shape(X));
 cell = genarray(cellshape,0.0d); /* not used, but SAC needs help */
 z = with (. <= iv <= .)
 genarray(frameshape,indrfr0(i,X[iv]),cell);
return(z);
}

inline double[*] indrfr0(int i, double[*] X)
{ /* X[i;;;]  i is scalar */
 z = X[[i]];
 return(z);
}

inline double[*] indrfr0(int[+] i, double[*] X)
{ /* X[im;;;]  im is array */
 cellshape = drop([1],shape(X));
 cell = genarray(cellshape, 0.0d);
 z = with (. <= iv <= .)
 genarray(shape(i),
  X[i[iv]],cell);
 return(z);
}

inline bool[*] indrfr(int fr, int i, bool[+] X)
{ /* X[;;;scalari;;;], where i has fr semicolons to its left */
  /* Indexing is origin-0. Caller will correct this */
  /* This could stand some optimization, perhaps, for boolean i,
   * unless SAC avoids building an array-valued temp of toI(i).
   */
 frameshape = take([fr],shape(X));
 cellshape = drop([fr+1], shape(X));
 cell = genarray(cellshape,false);
 z = with (. <= iv <= .)
 genarray(frameshape,(X[iv])[i], cell);
return(z);
}

inline bool[*] indrfr(int fr, int[+] i, bool[+] X)
{ /* X[;;;i;;;], where i has fr semicolons to its left */
  /* Indexing is origin-0. Caller will correct this */
  /* This could stand some optimization, perhaps, for boolean i,
   * unless SAC avoids building an array-valued temp of toI(i).
   */
 frameshape = take([fr],shape(X));
 cellshape = shape(i)++drop([fr+1],shape(X));
 cell = genarray(cellshape,false); /* not used, but SAC needs help */
 z = with (. <= iv <= .)
 genarray(frameshape,indrfr0(i,X[iv]),cell);
return(z);
}

inline bool[*] indrfr0(int i, bool[*] X)
{ /* X[i;;;]  i is scalar */
 z = X[[i]];
 return(z);
}

inline bool[*] indrfr0(int[+] i, bool[*] X)
{ /* X[im;;;]  im is array */
 cellshape = drop([1],shape(X));
 cell = genarray(cellshape, false);
 z = with (. <= iv <= .)
 genarray(shape(i),
  X[i[iv]],cell);
 return(z);
}

inline double[+] indsfr(int fr, int[*] i, double[+] X, double[+] Y)
{ /* X[;;;i;;;]<- nonscalar Y, where i has fr semicolons to its left */
 cellshape = shape(i)++drop([fr],shape(X));
 cell = genarray(cellshape,0.0d); /* not used, but SAC needs help */
 frameshape = take([fr],shape(X));
 z = with (. <= iv <= .)
 genarray(frameshape,indsfr0(i,X[iv], Y[iv]),cell);
 zshape = frameshape++cellshape;
 return(reshape(zshape,z));
}


inline double[+] indsfr(int fr, int[+] i, double[+] X, double Y)
{ /* X[;;;i;;;]<- scalar Y, where i has fr semicolons to its left */
 cellshape = drop([fr+1],shape(X));
 cell = genarray(cellshape,0.0d); /* not used, but SAC needs help */
 frameshape = take([fr],shape(X));
 z = with (. <= iv <= .)
 genarray(frameshape,indsfr0(i,X[iv], Y),cell);
 return(z);
}

inline double[+] indsfr(int fr, int i, double[+] X, double Y)
{ /* X[;;;i;;;]<- scalar Y, where i has fr semicolons to its left */
 cellshape = drop([fr+1],shape(X));
 cell = genarray(cellshape,0.0d);
 frameshape = take([fr],shape(X));
 z = with (. <= iv <= .)
 genarray(frameshape,indsfr0(i,X[iv], Y),cell);
 zshape = frameshape++cellshape;
 return(z);
}

inline double[+] indsfr0(int i, double[+] X, double Y)
{ /* Case 1. X[scalarI;;]<- scalarY  NB. Leading axis  */
 cell = genarray(drop([1],shape(X)),tod(Y));
 z = tod(X);
 z[[i]] = cell;
 return(z);
}

inline double[+] indsfr0(int i, double[+] X, double[+] Y)
{ /* Case 1. X[scalarI;;]<- non-scalarY  NB. Leading axis  */
  /* This has to match on shape. FIXME! */
 z = tod(X);
 z[[i]] = tod(Y);
 return(z);
}

inline double[+] indsfr0(int[.] iv, double[+] X, double Y)
{ /* 2. X[non-scalarIV;;]<- scalarY      NB. Leading axis */
  /* This would almost work under a with-loop, but the potential
   * for duplicates in iv scuppers that. Ergo, FOR loop.
   */
 z = tod(X);
 cellshape = drop([1],shape(X));
 cell = genarray (cellshape, tod(Y));
 raveli = APEXRavel(iv);
 for(i=0; i<shape(raveli)[0]; i++)
   z[raveli[i]] = cell;
 return(z);
}

inline double[+] indsfr0(int[+] i, double[+] X, double[+] Y)
{ /* 4. x[non-scalar;;]<- non-scalar  NB. Leading axis */
 z = tod(X);
 raveli = APEXRavel(i);
 for(k=0; k<shape(raveli)[0]; k++)
   z[raveli[k]] = tod(Y[k]);
 return(z);
}

inline int plusdotmpyIBI(int[.] x, bool[.] y)
{ /* Vector-Vector inner product */
 z = plusslIII(mpyIII(toi(x),toi(y)));
 return(z);
}

inline int[+] mpyIBI(int[*] x, bool y)
{ /* AxS scalar function */
 yel = toI(y);
 z = with( . <= iv <= .) {
 xel = toI(x[iv]);
 }
 genarray( shape(x),
 mpyIII(xel,yel),0);
 return(z);
}

inline int plusslIII(int[.] y)
{ /* Last axis reduction of vector */
 z = with (_mul_SxA_(0,shape(y)) <= iv < shape(y))
 fold( plusIII, toI(0), toI(y[iv]));
 return(z);
}

inline int plusdotmpyIII(int[.] x, int y)
{ /* Vector-Scalar inner product */
 z = plusslIII(mpyIII(toi(x),toi(y)));
 return(z);
}

inline int[+] mpyIII(int[*] x, int y)
{ /* AxS scalar function */
 yel = toI(y);
 z = with( . <= iv <= .) {
 xel = toI(x[iv]);
 }
 genarray( shape(x),
 mpyIII(xel,yel),0);
 return(z);
}

inline int plusdotmpyIII(int[.] x, int[.] y)
{ /* Vector-Vector inner product */
 z = plusslIII(mpyIII(toi(x),toi(y)));
 return(z);
}

inline int[*] plusdotmpyIII(int[.] x, int[*] y)
{ /* TRANSPOSE case of inner product z = vector f.g y */
 yt = APEXTranspose(toi(y));
 xct = toi(x);
 /* if (1 != shape(xct)[0]) FIXME; length error check */

 shp = drop([-1],shape(xct)) ++ drop([1], shape(y));
 z = with(iv)
 (. <= iv <= .) {
 vx = xct[take([dim(x)-1], iv)];
 vy = yt[ reverse(take([1-dim(y)], iv))];
 } : plusslIII(mpyIII(vx,vy));
 genarray(shp, 0);
 return(z);
}
inline int[+] mpyIII(int[+] x, int[+] y)
{ /* AxA Dyadic scalar fn, shapes unknown, shapes may or may not match */
 /* Insert length error checking code here!! */
 z = with( . <= iv <= .) {
 xel = toI(x[iv]);
 yel = toI(y[iv]);
 }
 genarray(shape(x),
 mpyIII(xel,yel),0);
 return(z);
}

inline int[.] slBII(bool x, int[.] y)
{ /* Scalar replicate vector */
 cell = genarray([toi(x)], 0);
 z = with(. <= iv <= .)
 genarray(shape(y),
 genarray([toi(x)], y[iv]),cell);
 return(APEXRavel(z));
}

inline int mpyslXII(int[.] y)
{ /* Last axis reduction of vector */
 z = with (_mul_SxA_(0,shape(y)) <= iv < shape(y))
 fold( mpyIII, toI(1), toI(y[iv]));
 return(z);
}

inline int[.] mpybslXII(int[.] y)
{
/* This does the scan in the wrong direction, but since
 * we assume associative functions only, it should be ok.
 */
 size = shape(y);
 z = genarray(size,toI(1));
 if (0 != size[0]) {
 /* real work to do */
 z[[0]] = toI(y[0]); /* Not sure about the coercion... */
 for ( i=1; i<size[0]; i++) {
 z[i] = mpyIII(toI(z[i-1]),toI(y[i]));
 }
 }
 return(z);
}

inline bool[.] mpybslXBB(bool[.] y)
{
/* This does the scan in the wrong direction, but since
 * we assume associative functions only, it should be ok.
 */
 size = shape(y);
 z = genarray(size,toB(1));
 if (0 != size[0]) {
 /* real work to do */
 z[[0]] = toB(y[0]); /* Not sure about the coercion... */
 for ( i=1; i<size[0]; i++) {
 z[i] = mpyBBB(toB(z[i-1]),toB(y[i]));
 }
 }
 return(z);
}

inline bool andslXBB(bool[.] y)
{ /* Last axis reduction of vector */
 z = with (_mul_SxA_(0,shape(y)) <= iv < shape(y))
 fold( andBBB, toB(1), toB(y[iv]));
 return(z);
}

inline bool mpyBBB(bool x, bool y)
{ return(tob(x) & tob(y));
}

inline bool andBBB(bool x, bool y)
{ return(tob(x)&tob(y));
}

inline int BaseValueIBI(int x, bool[.] y,double QUADct)
{ /* Start of function BaseValue 2006-01-08 23:32:20.000 */

/*
r{<-}x BaseValue y;nx;ny;wts
@ APL model of Base Value
@@:if 1={times}/{rho}x @ Extend singleton x
@@  nx{<-}(1{take}{rho}y){rho}x
@@:else
@@  nx{<-},x
nx{<-}1={times}/{rho}x
nx{<-}(nx/(1{take}{rho}y){rho}x),(~nx)/,x
@@:endif
@@if 1={times}/{rho}y @ Extend singleton y
@@ ny{<-}({rho}nx){rho}y
@@lse
@@ny{<-}y
@@:endif
wts{<-}{reverse}{times}\{reverse}(1{drop}nx),1
r{<-}wts+.{times}y
*/
TMP_29=rhoXII( x) ;
 TMP_30=mpyslXII( TMP_29 ) ;
 nx_0=eqBIB(true,TMP_30 );
 TMP_36=comaXII( x) ;
 TMP_37=notXBB( nx_0) ;
 TMP_38=slBII(TMP_37 ,TMP_36 ) ;
 TMP_39=rhoXBI( y) ;
 TMP_40=takeBII(true,TMP_39 ) ;
 TMP_41=rhoIII(TMP_40 ,x) ;
 TMP_42=slBII(nx_0,TMP_41 ) ;
 nx_1=comaIII(TMP_42 ,TMP_38 ) ;
 TMP_44=dropBII(true,nx_1) ;
 TMP_45=comaIBI(TMP_44 ,true) ;
 TMP_46=rotrXII( TMP_45 ) ;
 TMP_47=mpybslXII( TMP_46 ) ;
 wts_0=rotrXII( TMP_47 ) ;
 r_0=plusdotmpyIBI(wts_0,y) ;
 r=r_0;
return(r_0);
}

inline double[.,.] RepresentIID(int[.] x, int[.] y,int QUADio)
{ /* Start of function Represent 2006-01-08 23:32:21.000 */

/*
z{<-}x Represent y;weights;i;cy;shpy;shpz;#ct;yi;zi
@ APL model of Represent.
#ct{<-}0
shpz{<-}({rho}x),{rho}y
x{<-},x
weights{<-}1{drop}{reverse}{times}\{reverse}x,1
shpy{<-}{rho}y
y{<-},y
z{<-}({times}/shpz){rho}{neg}1
:for yi :in {iota}{times}/shpy
 cy{<-}y[yi]
 :for i :in {reverse}{iota}{rho}x
  zi{<-}(i{times}{times}/shpy)+yi
  z[zi]{<-}x[i]|cy{divide}weights[i]
  cy{<-}cy-z[zi]{times}weights[i]
 :endfor
:endfor
z{<-}shpz{rho}z
*/
QUADct_0=tod(( false) );
 TMP_59=rhoXII( y) ;
 TMP_60=rhoXII( x) ;
 shpz_0=comaIII(TMP_60 ,TMP_59 ) ;
 x_0=comaXII( x) ;
 TMP_63=comaIBI(x_0,true) ;
 TMP_64=rotrXII( TMP_63 ) ;
 TMP_65=mpybslXII( TMP_64 ) ;
 TMP_69=rotrXII( TMP_65 ) ;
 weights_0=dropBII(true,TMP_69 ) ;
 shpy_0=rhoXII( y) ;
 y_0=comaXII( y) ;
 TMP_73=mpyslXII( shpz_0) ;
 z_0=rhoIII(TMP_73 ,-1) ;
 TMP_78=mpyslXII( shpy_0) ;
 TMP_83=iotaXII( TMP_78 ,QUADio) ;
 CTR084_= 0;
CTR084z_ = (_shape_(TMP_83 )[0])-1;
z_3=tod(z_0);
for(; CTR084_ <= CTR084z_; CTR084_++){
yi_0 = TMP_83 [CTR084_];
 cy_1=y_0[[toi(yi_0)-QUADio]];
 TMP_90=rhoXII( x_0) ;
 TMP_92=iotaXII( TMP_90 ,QUADio) ;
 TMP_93=rotrXII( TMP_92 ) ;
 CTR094_= 0;
CTR094z_ = (_shape_(TMP_93 )[0])-1;
z_3=tod(z_3);
cy_3=tod(cy_1);
for(; CTR094_ <= CTR094z_; CTR094_++){
i_0 = TMP_93 [CTR094_];
 TMP_97=mpyslXII( shpy_0) ;
 TMP_101=mpyIII(i_0,TMP_97 ) ;
 zi_0=plusIII(TMP_101 ,yi_0) ;
 TMP_105=weights_0[[toi(i_0)-QUADio]];
 TMP_106=divDID(cy_3,TMP_105 ) ;
 TMP_109=x_0[[toi(i_0)-QUADio]];
 TMP_111=modIDD(TMP_109 ,TMP_106 ,QUADct_0) ;
 TMP_114=modarray(z_3,[toi(zi_0)-QUADio],TMP_111 ) ;
 z_3=( TMP_114 ) ;
 TMP_118=weights_0[[toi(i_0)-QUADio]];
 TMP_121=z_3[[toi(zi_0)-QUADio]];
 TMP_122=mpyDID(TMP_121 ,TMP_118 ) ;
 cy_3=barDDD(cy_3,TMP_122 ) ;
 }
 }
 z_4=rhoIDD(shpz_0,z_3) ;
 z=z_4;
return(z_4);
}

inline int BaseValueCLONE0III(int x, int y,double QUADct)
{ /* Start of function BaseValueCLONE0 2006-01-08 23:32:21.000 */

/*
r{<-}x BaseValue y;nx;ny;wts
@ APL model of Base Value
@@:if 1={times}/{rho}x @ Extend singleton x
@@  nx{<-}(1{take}{rho}y){rho}x
@@:else
@@  nx{<-},x
nx{<-}1={times}/{rho}x
nx{<-}(nx/(1{take}{rho}y){rho}x),(~nx)/,x
@@:endif
@@if 1={times}/{rho}y @ Extend singleton y
@@ ny{<-}({rho}nx){rho}y
@@lse
@@ny{<-}y
@@:endif
wts{<-}{reverse}{times}\{reverse}(1{drop}nx),1
r{<-}wts+.{times}y
*/
TMP_29=rhoXII( x) ;
 TMP_30=mpyslXII( TMP_29 ) ;
 nx_0=eqBIB(true,TMP_30 );
 TMP_36=comaXII( x) ;
 TMP_37=notXBB( nx_0) ;
 TMP_38=slBII(TMP_37 ,TMP_36 ) ;
 TMP_39=rhoXII( y) ;
 TMP_40=takeBII(true,TMP_39 ) ;
 TMP_41=rhoIII(TMP_40 ,x) ;
 TMP_42=slBII(nx_0,TMP_41 ) ;
 nx_1=comaIII(TMP_42 ,TMP_38 ) ;
 TMP_44=dropBII(true,nx_1) ;
 TMP_45=comaIBI(TMP_44 ,true) ;
 TMP_46=rotrXII( TMP_45 ) ;
 TMP_47=mpybslXII( TMP_46 ) ;
 wts_0=rotrXII( TMP_47 ) ;
 r_0=plusdotmpyIII(wts_0,y) ;
 r=r_0;
return(r_0);
}

inline int BaseValueCLONE1III(int x, int[.] y,double QUADct)
{ /* Start of function BaseValueCLONE1 2006-01-08 23:32:22.000 */

/*
r{<-}x BaseValue y;nx;ny;wts
@ APL model of Base Value
@@:if 1={times}/{rho}x @ Extend singleton x
@@  nx{<-}(1{take}{rho}y){rho}x
@@:else
@@  nx{<-},x
nx{<-}1={times}/{rho}x
nx{<-}(nx/(1{take}{rho}y){rho}x),(~nx)/,x
@@:endif
@@if 1={times}/{rho}y @ Extend singleton y
@@ ny{<-}({rho}nx){rho}y
@@lse
@@ny{<-}y
@@:endif
wts{<-}{reverse}{times}\{reverse}(1{drop}nx),1
r{<-}wts+.{times}y
*/
TMP_29=rhoXII( x) ;
 TMP_30=mpyslXII( TMP_29 ) ;
 nx_0=eqBIB(true,TMP_30 );
 TMP_36=comaXII( x) ;
 TMP_37=notXBB( nx_0) ;
 TMP_38=slBII(TMP_37 ,TMP_36 ) ;
 TMP_39=rhoXII( y) ;
 TMP_40=takeBII(true,TMP_39 ) ;
 TMP_41=rhoIII(TMP_40 ,x) ;
 TMP_42=slBII(nx_0,TMP_41 ) ;
 nx_1=comaIII(TMP_42 ,TMP_38 ) ;
 TMP_44=dropBII(true,nx_1) ;
 TMP_45=comaIBI(TMP_44 ,true) ;
 TMP_46=rotrXII( TMP_45 ) ;
 TMP_47=mpybslXII( TMP_46 ) ;
 wts_0=rotrXII( TMP_47 ) ;
 r_0=plusdotmpyIII(wts_0,y) ;
 r=r_0;
return(r_0);
}

inline int BaseValueCLONE2III(int[.] x, int[.] y,double QUADct)
{ /* Start of function BaseValueCLONE2 2006-01-08 23:32:22.000 */

/*
r{<-}x BaseValue y;nx;ny;wts
@ APL model of Base Value
@@:if 1={times}/{rho}x @ Extend singleton x
@@  nx{<-}(1{take}{rho}y){rho}x
@@:else
@@  nx{<-},x
nx{<-}1={times}/{rho}x
nx{<-}(nx/(1{take}{rho}y){rho}x),(~nx)/,x
@@:endif
@@if 1={times}/{rho}y @ Extend singleton y
@@ ny{<-}({rho}nx){rho}y
@@lse
@@ny{<-}y
@@:endif
wts{<-}{reverse}{times}\{reverse}(1{drop}nx),1
r{<-}wts+.{times}y
*/
TMP_29=rhoXII( x) ;
 TMP_30=mpyslXII( TMP_29 ) ;
 nx_0=eqBIB(true,TMP_30 );
 TMP_36=comaXII( x) ;
 TMP_37=notXBB( nx_0) ;
 TMP_38=slBII(TMP_37 ,TMP_36 ) ;
 TMP_39=rhoXII( y) ;
 TMP_40=takeBII(true,TMP_39 ) ;
 TMP_41=rhoIII(TMP_40 ,x) ;
 TMP_42=slBII(nx_0,TMP_41 ) ;
 nx_1=comaIII(TMP_42 ,TMP_38 ) ;
 TMP_44=dropBII(true,nx_1) ;
 TMP_45=comaIBI(TMP_44 ,true) ;
 TMP_46=rotrXII( TMP_45 ) ;
 TMP_47=mpybslXII( TMP_46 ) ;
 wts_0=rotrXII( TMP_47 ) ;
 r_0=plusdotmpyIII(wts_0,y) ;
 r=r_0;
return(r_0);
}

inline int[.,.,.] BaseValueCLONE3III(int[.] x, int[.,.,.,.] y,double QUADct)
{ /* Start of function BaseValueCLONE3 2006-01-08 23:32:23.000 */

/*
r{<-}x BaseValue y;nx;ny;wts
@ APL model of Base Value
@@:if 1={times}/{rho}x @ Extend singleton x
@@  nx{<-}(1{take}{rho}y){rho}x
@@:else
@@  nx{<-},x
nx{<-}1={times}/{rho}x
nx{<-}(nx/(1{take}{rho}y){rho}x),(~nx)/,x
@@:endif
@@if 1={times}/{rho}y @ Extend singleton y
@@ ny{<-}({rho}nx){rho}y
@@lse
@@ny{<-}y
@@:endif
wts{<-}{reverse}{times}\{reverse}(1{drop}nx),1
r{<-}wts+.{times}y
*/
TMP_29=rhoXII( x) ;
 TMP_30=mpyslXII( TMP_29 ) ;
 nx_0=eqBIB(true,TMP_30 );
 TMP_36=comaXII( x) ;
 TMP_37=notXBB( nx_0) ;
 TMP_38=slBII(TMP_37 ,TMP_36 ) ;
 TMP_39=rhoXII( y) ;
 TMP_40=takeBII(true,TMP_39 ) ;
 TMP_41=rhoIII(TMP_40 ,x) ;
 TMP_42=slBII(nx_0,TMP_41 ) ;
 nx_1=comaIII(TMP_42 ,TMP_38 ) ;
 TMP_44=dropBII(true,nx_1) ;
 TMP_45=comaIBI(TMP_44 ,true) ;
 TMP_46=rotrXII( TMP_45 ) ;
 TMP_47=mpybslXII( TMP_46 ) ;
 wts_0=rotrXII( TMP_47 ) ;
 r_0=plusdotmpyIII(wts_0,y) ;
 r=r_0;
return(r_0);
}

inline double[.,.] RepresentCLONE0IID(int[.] x, int[.] y,int QUADio)
{ /* Start of function RepresentCLONE0 2006-01-08 23:32:23.000 */

/*
z{<-}x Represent y;weights;i;cy;shpy;shpz;#ct;yi;zi
@ APL model of Represent.
#ct{<-}0
shpz{<-}({rho}x),{rho}y
x{<-},x
weights{<-}1{drop}{reverse}{times}\{reverse}x,1
shpy{<-}{rho}y
y{<-},y
z{<-}({times}/shpz){rho}{neg}1
:for yi :in {iota}{times}/shpy
 cy{<-}y[yi]
 :for i :in {reverse}{iota}{rho}x
  zi{<-}(i{times}{times}/shpy)+yi
  z[zi]{<-}x[i]|cy{divide}weights[i]
  cy{<-}cy-z[zi]{times}weights[i]
 :endfor
:endfor
z{<-}shpz{rho}z
*/
QUADct_0=tod(( false) );
 TMP_59=rhoXII( y) ;
 TMP_60=rhoXII( x) ;
 shpz_0=comaIII(TMP_60 ,TMP_59 ) ;
 x_0=comaXII( x) ;
 TMP_63=comaIBI(x_0,true) ;
 TMP_64=rotrXII( TMP_63 ) ;
 TMP_65=mpybslXII( TMP_64 ) ;
 TMP_69=rotrXII( TMP_65 ) ;
 weights_0=dropBII(true,TMP_69 ) ;
 shpy_0=rhoXII( y) ;
 y_0=comaXII( y) ;
 TMP_73=mpyslXII( shpz_0) ;
 z_0=rhoIII(TMP_73 ,-1) ;
 TMP_78=mpyslXII( shpy_0) ;
 TMP_83=iotaXII( TMP_78 ,QUADio) ;
 CTR084_= 0;
CTR084z_ = (_shape_(TMP_83 )[0])-1;
z_3=tod(z_0);
for(; CTR084_ <= CTR084z_; CTR084_++){
yi_0 = TMP_83 [CTR084_];
 cy_1=y_0[[toi(yi_0)-QUADio]];
 TMP_90=rhoXII( x_0) ;
 TMP_92=iotaXII( TMP_90 ,QUADio) ;
 TMP_93=rotrXII( TMP_92 ) ;
 CTR094_= 0;
CTR094z_ = (_shape_(TMP_93 )[0])-1;
z_3=tod(z_3);
cy_3=tod(cy_1);
for(; CTR094_ <= CTR094z_; CTR094_++){
i_0 = TMP_93 [CTR094_];
 TMP_97=mpyslXII( shpy_0) ;
 TMP_101=mpyIII(i_0,TMP_97 ) ;
 zi_0=plusIII(TMP_101 ,yi_0) ;
 TMP_105=weights_0[[toi(i_0)-QUADio]];
 TMP_106=divDID(cy_3,TMP_105 ) ;
 TMP_109=x_0[[toi(i_0)-QUADio]];
 TMP_111=modIDD(TMP_109 ,TMP_106 ,QUADct_0) ;
 TMP_114=modarray(z_3,[toi(zi_0)-QUADio],TMP_111 ) ;
 z_3=( TMP_114 ) ;
 TMP_118=weights_0[[toi(i_0)-QUADio]];
 TMP_121=z_3[[toi(zi_0)-QUADio]];
 TMP_122=mpyDID(TMP_121 ,TMP_118 ) ;
 cy_3=barDDD(cy_3,TMP_122 ) ;
 }
 }
 z_4=rhoIDD(shpz_0,z_3) ;
 z=z_4;
return(z_4);
}

inline double[.,.,.,.,.] RepresentCLONE1IID(int[.] x, int[.,.,.,.] y,int QUADio)
{ /* Start of function RepresentCLONE1 2006-01-08 23:32:24.000 */

/*
z{<-}x Represent y;weights;i;cy;shpy;shpz;#ct;yi;zi
@ APL model of Represent.
#ct{<-}0
shpz{<-}({rho}x),{rho}y
x{<-},x
weights{<-}1{drop}{reverse}{times}\{reverse}x,1
shpy{<-}{rho}y
y{<-},y
z{<-}({times}/shpz){rho}{neg}1
:for yi :in {iota}{times}/shpy
 cy{<-}y[yi]
 :for i :in {reverse}{iota}{rho}x
  zi{<-}(i{times}{times}/shpy)+yi
  z[zi]{<-}x[i]|cy{divide}weights[i]
  cy{<-}cy-z[zi]{times}weights[i]
 :endfor
:endfor
z{<-}shpz{rho}z
*/
QUADct_0=tod(( false) );
 TMP_59=rhoXII( y) ;
 TMP_60=rhoXII( x) ;
 shpz_0=comaIII(TMP_60 ,TMP_59 ) ;
 x_0=comaXII( x) ;
 TMP_63=comaIBI(x_0,true) ;
 TMP_64=rotrXII( TMP_63 ) ;
 TMP_65=mpybslXII( TMP_64 ) ;
 TMP_69=rotrXII( TMP_65 ) ;
 weights_0=dropBII(true,TMP_69 ) ;
 shpy_0=rhoXII( y) ;
 y_0=comaXII( y) ;
 TMP_73=mpyslXII( shpz_0) ;
 z_0=rhoIII(TMP_73 ,-1) ;
 TMP_78=mpyslXII( shpy_0) ;
 TMP_83=iotaXII( TMP_78 ,QUADio) ;
 CTR084_= 0;
CTR084z_ = (_shape_(TMP_83 )[0])-1;
z_3=tod(z_0);
for(; CTR084_ <= CTR084z_; CTR084_++){
yi_0 = TMP_83 [CTR084_];
 cy_1=y_0[[toi(yi_0)-QUADio]];
 TMP_90=rhoXII( x_0) ;
 TMP_92=iotaXII( TMP_90 ,QUADio) ;
 TMP_93=rotrXII( TMP_92 ) ;
 CTR094_= 0;
CTR094z_ = (_shape_(TMP_93 )[0])-1;
z_3=tod(z_3);
cy_3=tod(cy_1);
for(; CTR094_ <= CTR094z_; CTR094_++){
i_0 = TMP_93 [CTR094_];
 TMP_97=mpyslXII( shpy_0) ;
 TMP_101=mpyIII(i_0,TMP_97 ) ;
 zi_0=plusIII(TMP_101 ,yi_0) ;
 TMP_105=weights_0[[toi(i_0)-QUADio]];
 TMP_106=divDID(cy_3,TMP_105 ) ;
 TMP_109=x_0[[toi(i_0)-QUADio]];
 TMP_111=modIDD(TMP_109 ,TMP_106 ,QUADct_0) ;
 TMP_114=modarray(z_3,[toi(zi_0)-QUADio],TMP_111 ) ;
 z_3=( TMP_114 ) ;
 TMP_118=weights_0[[toi(i_0)-QUADio]];
 TMP_121=z_3[[toi(zi_0)-QUADio]];
 TMP_122=mpyDID(TMP_121 ,TMP_118 ) ;
 cy_3=barDDD(cy_3,TMP_122 ) ;
 }
 }
 z_4=rhoIDD(shpz_0,z_3) ;
 z=z_4;
return(z_4);
}

inline double[.] RepresentCLONE2IID(int[.] x, int y,int QUADio)
{ /* Start of function RepresentCLONE2 2006-01-08 23:32:25.000 */

/*
z{<-}x Represent y;weights;i;cy;shpy;shpz;#ct;yi;zi
@ APL model of Represent.
#ct{<-}0
shpz{<-}({rho}x),{rho}y
x{<-},x
weights{<-}1{drop}{reverse}{times}\{reverse}x,1
shpy{<-}{rho}y
y{<-},y
z{<-}({times}/shpz){rho}{neg}1
:for yi :in {iota}{times}/shpy
 cy{<-}y[yi]
 :for i :in {reverse}{iota}{rho}x
  zi{<-}(i{times}{times}/shpy)+yi
  z[zi]{<-}x[i]|cy{divide}weights[i]
  cy{<-}cy-z[zi]{times}weights[i]
 :endfor
:endfor
z{<-}shpz{rho}z
*/
QUADct_0=tod(( false) );
 TMP_59=rhoXII( y) ;
 TMP_60=rhoXII( x) ;
 shpz_0=comaIII(TMP_60 ,TMP_59 ) ;
 x_0=comaXII( x) ;
 TMP_63=comaIBI(x_0,true) ;
 TMP_64=rotrXII( TMP_63 ) ;
 TMP_65=mpybslXII( TMP_64 ) ;
 TMP_69=rotrXII( TMP_65 ) ;
 weights_0=dropBII(true,TMP_69 ) ;
 shpy_0=rhoXII( y) ;
 y_0=comaXII( y) ;
 TMP_73=mpyslXII( shpz_0) ;
 z_0=rhoIII(TMP_73 ,-1) ;
 TMP_78=mpyslXII( shpy_0) ;
 TMP_83=iotaXII( TMP_78 ,QUADio) ;
 CTR084_= 0;
CTR084z_ = (_shape_(TMP_83 )[0])-1;
z_3=tod(z_0);
for(; CTR084_ <= CTR084z_; CTR084_++){
yi_0 = TMP_83 [CTR084_];
 cy_1=y_0[[toi(yi_0)-QUADio]];
 TMP_90=rhoXII( x_0) ;
 TMP_92=iotaXII( TMP_90 ,QUADio) ;
 TMP_93=rotrXII( TMP_92 ) ;
 CTR094_= 0;
CTR094z_ = (_shape_(TMP_93 )[0])-1;
z_3=tod(z_3);
cy_3=tod(cy_1);
for(; CTR094_ <= CTR094z_; CTR094_++){
i_0 = TMP_93 [CTR094_];
 TMP_97=mpyslXII( shpy_0) ;
 TMP_101=mpyIII(i_0,TMP_97 ) ;
 zi_0=plusIII(TMP_101 ,yi_0) ;
 TMP_105=weights_0[[toi(i_0)-QUADio]];
 TMP_106=divDID(cy_3,TMP_105 ) ;
 TMP_109=x_0[[toi(i_0)-QUADio]];
 TMP_111=modIDD(TMP_109 ,TMP_106 ,QUADct_0) ;
 TMP_114=modarray(z_3,[toi(zi_0)-QUADio],TMP_111 ) ;
 z_3=( TMP_114 ) ;
 TMP_118=weights_0[[toi(i_0)-QUADio]];
 TMP_121=z_3[[toi(zi_0)-QUADio]];
 TMP_122=mpyDID(TMP_121 ,TMP_118 ) ;
 cy_3=barDDD(cy_3,TMP_122 ) ;
 }
 }
 z_4=rhoIDD(shpz_0,z_3) ;
 z=z_4;
return(z_4);
}

inline double[.] RepresentCLONE3IID(int[.] x, int y,int QUADio)
{ /* Start of function RepresentCLONE3 2006-01-08 23:32:25.000 */

/*
z{<-}x Represent y;weights;i;cy;shpy;shpz;#ct;yi;zi
@ APL model of Represent.
#ct{<-}0
shpz{<-}({rho}x),{rho}y
x{<-},x
weights{<-}1{drop}{reverse}{times}\{reverse}x,1
shpy{<-}{rho}y
y{<-},y
z{<-}({times}/shpz){rho}{neg}1
:for yi :in {iota}{times}/shpy
 cy{<-}y[yi]
 :for i :in {reverse}{iota}{rho}x
  zi{<-}(i{times}{times}/shpy)+yi
  z[zi]{<-}x[i]|cy{divide}weights[i]
  cy{<-}cy-z[zi]{times}weights[i]
 :endfor
:endfor
z{<-}shpz{rho}z
*/
QUADct_0=tod(( false) );
 TMP_59=rhoXII( y) ;
 TMP_60=rhoXII( x) ;
 shpz_0=comaIII(TMP_60 ,TMP_59 ) ;
 x_0=comaXII( x) ;
 TMP_63=comaIBI(x_0,true) ;
 TMP_64=rotrXII( TMP_63 ) ;
 TMP_65=mpybslXII( TMP_64 ) ;
 TMP_69=rotrXII( TMP_65 ) ;
 weights_0=dropBII(true,TMP_69 ) ;
 shpy_0=rhoXII( y) ;
 y_0=comaXII( y) ;
 TMP_73=mpyslXII( shpz_0) ;
 z_0=rhoIII(TMP_73 ,-1) ;
 TMP_78=mpyslXII( shpy_0) ;
 TMP_83=iotaXII( TMP_78 ,QUADio) ;
 CTR084_= 0;
CTR084z_ = (_shape_(TMP_83 )[0])-1;
z_3=tod(z_0);
for(; CTR084_ <= CTR084z_; CTR084_++){
yi_0 = TMP_83 [CTR084_];
 cy_1=y_0[[toi(yi_0)-QUADio]];
 TMP_90=rhoXII( x_0) ;
 TMP_92=iotaXII( TMP_90 ,QUADio) ;
 TMP_93=rotrXII( TMP_92 ) ;
 CTR094_= 0;
CTR094z_ = (_shape_(TMP_93 )[0])-1;
z_3=tod(z_3);
cy_3=tod(cy_1);
for(; CTR094_ <= CTR094z_; CTR094_++){
i_0 = TMP_93 [CTR094_];
 TMP_97=mpyslXII( shpy_0) ;
 TMP_101=mpyIII(i_0,TMP_97 ) ;
 zi_0=plusIII(TMP_101 ,yi_0) ;
 TMP_105=weights_0[[toi(i_0)-QUADio]];
 TMP_106=divDID(cy_3,TMP_105 ) ;
 TMP_109=x_0[[toi(i_0)-QUADio]];
 TMP_111=modIDD(TMP_109 ,TMP_106 ,QUADct_0) ;
 TMP_114=modarray(z_3,[toi(zi_0)-QUADio],TMP_111 ) ;
 z_3=( TMP_114 ) ;
 TMP_118=weights_0[[toi(i_0)-QUADio]];
 TMP_121=z_3[[toi(zi_0)-QUADio]];
 TMP_122=mpyDID(TMP_121 ,TMP_118 ) ;
 cy_3=barDDD(cy_3,TMP_122 ) ;
 }
 }
 z_4=rhoIDD(shpz_0,z_3) ;
 z=z_4;
return(z_4);
}

inline double[.] RepresentCLONE4IDD(int[.] x, double y,int QUADio)
{ /* Start of function RepresentCLONE4 2006-01-08 23:32:26.000 */

/*
z{<-}x Represent y;weights;i;cy;shpy;shpz;#ct;yi;zi
@ APL model of Represent.
#ct{<-}0
shpz{<-}({rho}x),{rho}y
x{<-},x
weights{<-}1{drop}{reverse}{times}\{reverse}x,1
shpy{<-}{rho}y
y{<-},y
z{<-}({times}/shpz){rho}{neg}1
:for yi :in {iota}{times}/shpy
 cy{<-}y[yi]
 :for i :in {reverse}{iota}{rho}x
  zi{<-}(i{times}{times}/shpy)+yi
  z[zi]{<-}x[i]|cy{divide}weights[i]
  cy{<-}cy-z[zi]{times}weights[i]
 :endfor
:endfor
z{<-}shpz{rho}z
*/
QUADct_0=tod(( false) );
 TMP_59=rhoXDI( y) ;
 TMP_60=rhoXII( x) ;
 shpz_0=comaIII(TMP_60 ,TMP_59 ) ;
 x_0=comaXII( x) ;
 TMP_63=comaIBI(x_0,true) ;
 TMP_64=rotrXII( TMP_63 ) ;
 TMP_65=mpybslXII( TMP_64 ) ;
 TMP_69=rotrXII( TMP_65 ) ;
 weights_0=dropBII(true,TMP_69 ) ;
 shpy_0=rhoXDI( y) ;
 y_0=comaXDD( y) ;
 TMP_73=mpyslXII( shpz_0) ;
 z_0=rhoIII(TMP_73 ,-1) ;
 TMP_78=mpyslXII( shpy_0) ;
 TMP_83=iotaXII( TMP_78 ,QUADio) ;
 CTR084_= 0;
CTR084z_ = (_shape_(TMP_83 )[0])-1;
z_3=tod(z_0);
for(; CTR084_ <= CTR084z_; CTR084_++){
yi_0 = TMP_83 [CTR084_];
 cy_1=y_0[[toi(yi_0)-QUADio]];
 TMP_90=rhoXII( x_0) ;
 TMP_92=iotaXII( TMP_90 ,QUADio) ;
 TMP_93=rotrXII( TMP_92 ) ;
 CTR094_= 0;
CTR094z_ = (_shape_(TMP_93 )[0])-1;
z_3=tod(z_3);
cy_3=tod(cy_1);
for(; CTR094_ <= CTR094z_; CTR094_++){
i_0 = TMP_93 [CTR094_];
 TMP_97=mpyslXII( shpy_0) ;
 TMP_101=mpyIII(i_0,TMP_97 ) ;
 zi_0=plusIII(TMP_101 ,yi_0) ;
 TMP_105=weights_0[[toi(i_0)-QUADio]];
 TMP_106=divDID(cy_3,TMP_105 ) ;
 TMP_109=x_0[[toi(i_0)-QUADio]];
 TMP_111=modIDD(TMP_109 ,TMP_106 ,QUADct_0) ;
 TMP_114=modarray(z_3,[toi(zi_0)-QUADio],TMP_111 ) ;
 z_3=( TMP_114 ) ;
 TMP_118=weights_0[[toi(i_0)-QUADio]];
 TMP_121=z_3[[toi(zi_0)-QUADio]];
 TMP_122=mpyDID(TMP_121 ,TMP_118 ) ;
 cy_3=barDDD(cy_3,TMP_122 ) ;
 }
 }
 z_4=rhoIDD(shpz_0,z_3) ;
 z=z_4;
return(z_4);
}

inline double[.] RepresentCLONE5BDD(bool[.] x, double y,int QUADio)
{ /* Start of function RepresentCLONE5 2006-01-08 23:32:26.000 */

/*
z{<-}x Represent y;weights;i;cy;shpy;shpz;#ct;yi;zi
@ APL model of Represent.
#ct{<-}0
shpz{<-}({rho}x),{rho}y
x{<-},x
weights{<-}1{drop}{reverse}{times}\{reverse}x,1
shpy{<-}{rho}y
y{<-},y
z{<-}({times}/shpz){rho}{neg}1
:for yi :in {iota}{times}/shpy
 cy{<-}y[yi]
 :for i :in {reverse}{iota}{rho}x
  zi{<-}(i{times}{times}/shpy)+yi
  z[zi]{<-}x[i]|cy{divide}weights[i]
  cy{<-}cy-z[zi]{times}weights[i]
 :endfor
:endfor
z{<-}shpz{rho}z
*/
QUADct_0=tod(( false) );
 TMP_59=rhoXDI( y) ;
 TMP_60=rhoXBI( x) ;
 shpz_0=comaIII(TMP_60 ,TMP_59 ) ;
 x_0=comaXBB( x) ;
 TMP_63=comaBBB(x_0,true) ;
 TMP_64=rotrXBB( TMP_63 ) ;
 TMP_65=mpybslXBB( TMP_64 ) ;
 TMP_69=rotrXBB( TMP_65 ) ;
 weights_0=dropBBB(true,TMP_69 ) ;
 shpy_0=rhoXDI( y) ;
 y_0=comaXDD( y) ;
 TMP_73=mpyslXII( shpz_0) ;
 z_0=rhoIII(TMP_73 ,-1) ;
 TMP_78=mpyslXII( shpy_0) ;
 TMP_83=iotaXII( TMP_78 ,QUADio) ;
 CTR084_= 0;
CTR084z_ = (_shape_(TMP_83 )[0])-1;
z_3=tod(z_0);
for(; CTR084_ <= CTR084z_; CTR084_++){
yi_0 = TMP_83 [CTR084_];
 cy_1=y_0[[toi(yi_0)-QUADio]];
 TMP_90=rhoXBI( x_0) ;
 TMP_92=iotaXII( TMP_90 ,QUADio) ;
 TMP_93=rotrXII( TMP_92 ) ;
 CTR094_= 0;
CTR094z_ = (_shape_(TMP_93 )[0])-1;
z_3=tod(z_3);
cy_3=tod(cy_1);
for(; CTR094_ <= CTR094z_; CTR094_++){
i_0 = TMP_93 [CTR094_];
 TMP_97=mpyslXII( shpy_0) ;
 TMP_101=mpyIII(i_0,TMP_97 ) ;
 zi_0=plusIII(TMP_101 ,yi_0) ;
 TMP_105=weights_0[[toi(i_0)-QUADio]];
 TMP_106=divDBD(cy_3,TMP_105 ) ;
 TMP_109=x_0[[toi(i_0)-QUADio]];
 TMP_111=modBDD(TMP_109 ,TMP_106 ,QUADct_0) ;
 TMP_114=modarray(z_3,[toi(zi_0)-QUADio],TMP_111 ) ;
 z_3=( TMP_114 ) ;
 TMP_118=weights_0[[toi(i_0)-QUADio]];
 TMP_121=z_3[[toi(zi_0)-QUADio]];
 TMP_122=mpyDBD(TMP_121 ,TMP_118 ) ;
 cy_3=barDDD(cy_3,TMP_122 ) ;
 }
 }
 z_4=rhoIDD(shpz_0,z_3) ;
 z=z_4;
return(z_4);
}

inline double[.] RepresentCLONE6IID(int[.] x, int y,int QUADio)
{ /* Start of function RepresentCLONE6 2006-01-08 23:32:27.000 */

/*
z{<-}x Represent y;weights;i;cy;shpy;shpz;#ct;yi;zi
@ APL model of Represent.
#ct{<-}0
shpz{<-}({rho}x),{rho}y
x{<-},x
weights{<-}1{drop}{reverse}{times}\{reverse}x,1
shpy{<-}{rho}y
y{<-},y
z{<-}({times}/shpz){rho}{neg}1
:for yi :in {iota}{times}/shpy
 cy{<-}y[yi]
 :for i :in {reverse}{iota}{rho}x
  zi{<-}(i{times}{times}/shpy)+yi
  z[zi]{<-}x[i]|cy{divide}weights[i]
  cy{<-}cy-z[zi]{times}weights[i]
 :endfor
:endfor
z{<-}shpz{rho}z
*/
QUADct_0=tod(( false) );
 TMP_59=rhoXII( y) ;
 TMP_60=rhoXII( x) ;
 shpz_0=comaIII(TMP_60 ,TMP_59 ) ;
 x_0=comaXII( x) ;
 TMP_63=comaIBI(x_0,true) ;
 TMP_64=rotrXII( TMP_63 ) ;
 TMP_65=mpybslXII( TMP_64 ) ;
 TMP_69=rotrXII( TMP_65 ) ;
 weights_0=dropBII(true,TMP_69 ) ;
 shpy_0=rhoXII( y) ;
 y_0=comaXII( y) ;
 TMP_73=mpyslXII( shpz_0) ;
 z_0=rhoIII(TMP_73 ,-1) ;
 TMP_78=mpyslXII( shpy_0) ;
 TMP_83=iotaXII( TMP_78 ,QUADio) ;
 CTR084_= 0;
CTR084z_ = (_shape_(TMP_83 )[0])-1;
z_3=tod(z_0);
for(; CTR084_ <= CTR084z_; CTR084_++){
yi_0 = TMP_83 [CTR084_];
 cy_1=y_0[[toi(yi_0)-QUADio]];
 TMP_90=rhoXII( x_0) ;
 TMP_92=iotaXII( TMP_90 ,QUADio) ;
 TMP_93=rotrXII( TMP_92 ) ;
 CTR094_= 0;
CTR094z_ = (_shape_(TMP_93 )[0])-1;
z_3=tod(z_3);
cy_3=tod(cy_1);
for(; CTR094_ <= CTR094z_; CTR094_++){
i_0 = TMP_93 [CTR094_];
 TMP_97=mpyslXII( shpy_0) ;
 TMP_101=mpyIII(i_0,TMP_97 ) ;
 zi_0=plusIII(TMP_101 ,yi_0) ;
 TMP_105=weights_0[[toi(i_0)-QUADio]];
 TMP_106=divDID(cy_3,TMP_105 ) ;
 TMP_109=x_0[[toi(i_0)-QUADio]];
 TMP_111=modIDD(TMP_109 ,TMP_106 ,QUADct_0) ;
 TMP_114=modarray(z_3,[toi(zi_0)-QUADio],TMP_111 ) ;
 z_3=( TMP_114 ) ;
 TMP_118=weights_0[[toi(i_0)-QUADio]];
 TMP_121=z_3[[toi(zi_0)-QUADio]];
 TMP_122=mpyDID(TMP_121 ,TMP_118 ) ;
 cy_3=barDDD(cy_3,TMP_122 ) ;
 }
 }
 z_4=rhoIDD(shpz_0,z_3) ;
 z=z_4;
return(z_4);
}

inline double[.] RepresentCLONE7BDD(bool[.] x, double y,int QUADio)
{ /* Start of function RepresentCLONE7 2006-01-08 23:32:27.000 */

/*
z{<-}x Represent y;weights;i;cy;shpy;shpz;#ct;yi;zi
@ APL model of Represent.
#ct{<-}0
shpz{<-}({rho}x),{rho}y
x{<-},x
weights{<-}1{drop}{reverse}{times}\{reverse}x,1
shpy{<-}{rho}y
y{<-},y
z{<-}({times}/shpz){rho}{neg}1
:for yi :in {iota}{times}/shpy
 cy{<-}y[yi]
 :for i :in {reverse}{iota}{rho}x
  zi{<-}(i{times}{times}/shpy)+yi
  z[zi]{<-}x[i]|cy{divide}weights[i]
  cy{<-}cy-z[zi]{times}weights[i]
 :endfor
:endfor
z{<-}shpz{rho}z
*/
QUADct_0=tod(( false) );
 TMP_59=rhoXDI( y) ;
 TMP_60=rhoXBI( x) ;
 shpz_0=comaIII(TMP_60 ,TMP_59 ) ;
 x_0=comaXBB( x) ;
 TMP_63=comaBBB(x_0,true) ;
 TMP_64=rotrXBB( TMP_63 ) ;
 TMP_65=mpybslXBB( TMP_64 ) ;
 TMP_69=rotrXBB( TMP_65 ) ;
 weights_0=dropBBB(true,TMP_69 ) ;
 shpy_0=rhoXDI( y) ;
 y_0=comaXDD( y) ;
 TMP_73=mpyslXII( shpz_0) ;
 z_0=rhoIII(TMP_73 ,-1) ;
 TMP_78=mpyslXII( shpy_0) ;
 TMP_83=iotaXII( TMP_78 ,QUADio) ;
 CTR084_= 0;
CTR084z_ = (_shape_(TMP_83 )[0])-1;
z_3=tod(z_0);
for(; CTR084_ <= CTR084z_; CTR084_++){
yi_0 = TMP_83 [CTR084_];
 cy_1=y_0[[toi(yi_0)-QUADio]];
 TMP_90=rhoXBI( x_0) ;
 TMP_92=iotaXII( TMP_90 ,QUADio) ;
 TMP_93=rotrXII( TMP_92 ) ;
 CTR094_= 0;
CTR094z_ = (_shape_(TMP_93 )[0])-1;
z_3=tod(z_3);
cy_3=tod(cy_1);
for(; CTR094_ <= CTR094z_; CTR094_++){
i_0 = TMP_93 [CTR094_];
 TMP_97=mpyslXII( shpy_0) ;
 TMP_101=mpyIII(i_0,TMP_97 ) ;
 zi_0=plusIII(TMP_101 ,yi_0) ;
 TMP_105=weights_0[[toi(i_0)-QUADio]];
 TMP_106=divDBD(cy_3,TMP_105 ) ;
 TMP_109=x_0[[toi(i_0)-QUADio]];
 TMP_111=modBDD(TMP_109 ,TMP_106 ,QUADct_0) ;
 TMP_114=modarray(z_3,[toi(zi_0)-QUADio],TMP_111 ) ;
 z_3=( TMP_114 ) ;
 TMP_118=weights_0[[toi(i_0)-QUADio]];
 TMP_121=z_3[[toi(zi_0)-QUADio]];
 TMP_122=mpyDBD(TMP_121 ,TMP_118 ) ;
 cy_3=barDDD(cy_3,TMP_122 ) ;
 }
 }
 z_4=rhoIDD(shpz_0,z_3) ;
 z=z_4;
return(z_4);
}

inline double[.,.] RepresentCLONE8IID(int[.] x, int[.] y,int QUADio)
{ /* Start of function RepresentCLONE8 2006-01-08 23:32:28.000 */

/*
z{<-}x Represent y;weights;i;cy;shpy;shpz;#ct;yi;zi
@ APL model of Represent.
#ct{<-}0
shpz{<-}({rho}x),{rho}y
x{<-},x
weights{<-}1{drop}{reverse}{times}\{reverse}x,1
shpy{<-}{rho}y
y{<-},y
z{<-}({times}/shpz){rho}{neg}1
:for yi :in {iota}{times}/shpy
 cy{<-}y[yi]
 :for i :in {reverse}{iota}{rho}x
  zi{<-}(i{times}{times}/shpy)+yi
  z[zi]{<-}x[i]|cy{divide}weights[i]
  cy{<-}cy-z[zi]{times}weights[i]
 :endfor
:endfor
z{<-}shpz{rho}z
*/
QUADct_0=tod(( false) );
 TMP_59=rhoXII( y) ;
 TMP_60=rhoXII( x) ;
 shpz_0=comaIII(TMP_60 ,TMP_59 ) ;
 x_0=comaXII( x) ;
 TMP_63=comaIBI(x_0,true) ;
 TMP_64=rotrXII( TMP_63 ) ;
 TMP_65=mpybslXII( TMP_64 ) ;
 TMP_69=rotrXII( TMP_65 ) ;
 weights_0=dropBII(true,TMP_69 ) ;
 shpy_0=rhoXII( y) ;
 y_0=comaXII( y) ;
 TMP_73=mpyslXII( shpz_0) ;
 z_0=rhoIII(TMP_73 ,-1) ;
 TMP_78=mpyslXII( shpy_0) ;
 TMP_83=iotaXII( TMP_78 ,QUADio) ;
 CTR084_= 0;
CTR084z_ = (_shape_(TMP_83 )[0])-1;
z_3=tod(z_0);
for(; CTR084_ <= CTR084z_; CTR084_++){
yi_0 = TMP_83 [CTR084_];
 cy_1=y_0[[toi(yi_0)-QUADio]];
 TMP_90=rhoXII( x_0) ;
 TMP_92=iotaXII( TMP_90 ,QUADio) ;
 TMP_93=rotrXII( TMP_92 ) ;
 CTR094_= 0;
CTR094z_ = (_shape_(TMP_93 )[0])-1;
z_3=tod(z_3);
cy_3=tod(cy_1);
for(; CTR094_ <= CTR094z_; CTR094_++){
i_0 = TMP_93 [CTR094_];
 TMP_97=mpyslXII( shpy_0) ;
 TMP_101=mpyIII(i_0,TMP_97 ) ;
 zi_0=plusIII(TMP_101 ,yi_0) ;
 TMP_105=weights_0[[toi(i_0)-QUADio]];
 TMP_106=divDID(cy_3,TMP_105 ) ;
 TMP_109=x_0[[toi(i_0)-QUADio]];
 TMP_111=modIDD(TMP_109 ,TMP_106 ,QUADct_0) ;
 TMP_114=modarray(z_3,[toi(zi_0)-QUADio],TMP_111 ) ;
 z_3=( TMP_114 ) ;
 TMP_118=weights_0[[toi(i_0)-QUADio]];
 TMP_121=z_3[[toi(zi_0)-QUADio]];
 TMP_122=mpyDID(TMP_121 ,TMP_118 ) ;
 cy_3=barDDD(cy_3,TMP_122 ) ;
 }
 }
 z_4=rhoIDD(shpz_0,z_3) ;
 z=z_4;
return(z_4);
}

int main()
{ /* Start of function main 2006-01-08 23:32:29.000 */

/*
r{<-}main;y;S;BV;M4;Hex;Byte
S{<-}2
BV{<-}1 1 0 1 1 0 1 1
M4{<-}2 3 4 5{rho}{iota}999
Hex{<-}16 16
Byte{<-}256 256 256 256
r{<-}  (2{basevalue}BV){match}S BaseValue BV                               @ S{basevalue}BV
r{<-}r,(Hex{represent}#av{iota}#av){match}Hex Represent #av{iota}#av                 @ V{represent}IV
r{<-}r,((8{rho}2){represent}#av{iota}#av){match}(8{rho}2)Represent #av{iota}#av              @ IV{represent}IV
r{<-}r,(2{basevalue}1000){match}2 BaseValue 1000                           @ IV{basevalue}IS
r{<-}r,(2{basevalue}{iota}10){match}2 BaseValue {iota}10                             @ IS{basevalue}IV
r{<-}r,(0 100 100{basevalue}10 20 30){match}0 100 100 BaseValue 10 20 30   @ IV{basevalue}IV
r{<-}r,(2 3{basevalue}M4){match}2 3 BaseValue M4                           @ IV{basevalue}IM
r{<-}r,(2 3{represent}M4){match}2 3 Represent M4                           @ IV{represent}IM
r{<-}r,(10 10 10{represent}3247){match}10 10 10 Represent 3247
r{<-}r,(0  10 10{represent}3247){match} 0 10 10 Represent 3247
r{<-}r,(10 10 10{represent}247.33){match}10 10 10 Represent 247.33
r{<-}r,(0 1{represent} 247.33){match}0 1 Represent 247.33
r{<-}r,((4{rho}10){represent}{neg}1){match}(4{rho}10)Represent {neg}1
r{<-}r,(0 1{represent}{neg}247.33){match}0 1 Represent {neg}247.33
r{<-}r,((4{rho}10){represent}1984 2001 {neg}44){match}(4{rho}10)Represent 1984 2001 {neg}44
#{<-}r
r{<-}#{<-}1+^/r
*/
QUADio_0=toi(( false) );
 QUADct_0=( 1.0e-13) ;
 QUADpp_0=( 10) ;
 QUADpw_0=( 80) ;
 QUADrl_0=( 16807) ;
 S_0=( 2) ;
 BV_0=( [true,true,false,true,true,false,true,true]) ;
 TMP_119=iotaXII( 999,QUADio_0) ;
 M4_0=rhoIII([2, 3, 4, 5],TMP_119 ) ;
 Hex_0=( [16, 16]) ;
 Byte_0=( [256, 256, 256, 256]) ;
 TMP_128=BaseValueIBI(S_0,BV_0,QUADct_0) ;
 TMP_129=dtakIBI(2,BV_0) ;
 r_0=sameIIB(TMP_129 ,TMP_128 ) ;
 TMP_134=iotaCCI(['\x00','\x01','\x02','\x03','\x04','\x05','\x06','\x07',
'\x08','\x09','\x0a','\x0b','\x0c','\x0d','\x0e','\x0f',
'\x10','\x11','\x12','\x13','\x14','\x15','\x16','\x17',
'\x18','\x19','\x1a','\x1b','\x1c','\x1d','\x1e','\x1f',
'\x20','\x21','\x22','\x23','\x24','\x25','\x26','\x27',
'\x28','\x29','\x2a','\x2b','\x2c','\x2d','\x2e','\x2f',
'\x30','\x31','\x32','\x33','\x34','\x35','\x36','\x37',
'\x38','\x39','\x3a','\x3b','\x3c','\x3d','\x3e','\x3f',
'\x40','\x41','\x42','\x43','\x44','\x45','\x46','\x47',
'\x48','\x49','\x4a','\x4b','\x4c','\x4d','\x4e','\x4f',
'\x50','\x51','\x52','\x53','\x54','\x55','\x56','\x57',
'\x58','\x59','\x5a','\x5b','\x5c','\x5d','\x5e','\x5f',
'\x60','\x61','\x62','\x63','\x64','\x65','\x66','\x67',
'\x68','\x69','\x6a','\x6b','\x6c','\x6d','\x6e','\x6f',
'\x70','\x71','\x72','\x73','\x74','\x75','\x76','\x77',
'\x78','\x79','\x7a','\x7b','\x7c','\x7d','\x7e','\x7f',
'\x80','\x81','\x82','\x83','\x84','\x85','\x86','\x87',
'\x88','\x89','\x8a','\x8b','\x8c','\x8d','\x8e','\x8f',
'\x90','\x91','\x92','\x93','\x94','\x95','\x96','\x97',
'\x98','\x99','\x9a','\x9b','\x9c','\x9d','\x9e','\x9f',
'\xa0','\xa1','\xa2','\xa3','\xa4','\xa5','\xa6','\xa7',
'\xa8','\xa9','\xaa','\xab','\xac','\xad','\xae','\xaf',
'\xb0','\xb1','\xb2','\xb3','\xb4','\xb5','\xb6','\xb7',
'\xb8','\xb9','\xba','\xbb','\xbc','\xbd','\xbe','\xbf',
'\xc0','\xc1','\xc2','\xc3','\xc4','\xc5','\xc6','\xc7',
'\xc8','\xc9','\xca','\xcb','\xcc','\xcd','\xce','\xcf',
'\xd0','\xd1','\xd2','\xd3','\xd4','\xd5','\xd6','\xd7',
'\xd8','\xd9','\xda','\xdb','\xdc','\xdd','\xde','\xdf',
'\xe0','\xe1','\xe2','\xe3','\xe4','\xe5','\xe6','\xe7',
'\xe8','\xe9','\xea','\xeb','\xec','\xed','\xee','\xef',
'\xf0','\xf1','\xf2','\xf3','\xf4','\xf5','\xf6','\xf7',
'\xf8','\xf9','\xfa','\xfb','\xfc','\xfd','\xfe','\xff'
],['\x00','\x01','\x02','\x03','\x04','\x05','\x06','\x07',
'\x08','\x09','\x0a','\x0b','\x0c','\x0d','\x0e','\x0f',
'\x10','\x11','\x12','\x13','\x14','\x15','\x16','\x17',
'\x18','\x19','\x1a','\x1b','\x1c','\x1d','\x1e','\x1f',
'\x20','\x21','\x22','\x23','\x24','\x25','\x26','\x27',
'\x28','\x29','\x2a','\x2b','\x2c','\x2d','\x2e','\x2f',
'\x30','\x31','\x32','\x33','\x34','\x35','\x36','\x37',
'\x38','\x39','\x3a','\x3b','\x3c','\x3d','\x3e','\x3f',
'\x40','\x41','\x42','\x43','\x44','\x45','\x46','\x47',
'\x48','\x49','\x4a','\x4b','\x4c','\x4d','\x4e','\x4f',
'\x50','\x51','\x52','\x53','\x54','\x55','\x56','\x57',
'\x58','\x59','\x5a','\x5b','\x5c','\x5d','\x5e','\x5f',
'\x60','\x61','\x62','\x63','\x64','\x65','\x66','\x67',
'\x68','\x69','\x6a','\x6b','\x6c','\x6d','\x6e','\x6f',
'\x70','\x71','\x72','\x73','\x74','\x75','\x76','\x77',
'\x78','\x79','\x7a','\x7b','\x7c','\x7d','\x7e','\x7f',
'\x80','\x81','\x82','\x83','\x84','\x85','\x86','\x87',
'\x88','\x89','\x8a','\x8b','\x8c','\x8d','\x8e','\x8f',
'\x90','\x91','\x92','\x93','\x94','\x95','\x96','\x97',
'\x98','\x99','\x9a','\x9b','\x9c','\x9d','\x9e','\x9f',
'\xa0','\xa1','\xa2','\xa3','\xa4','\xa5','\xa6','\xa7',
'\xa8','\xa9','\xaa','\xab','\xac','\xad','\xae','\xaf',
'\xb0','\xb1','\xb2','\xb3','\xb4','\xb5','\xb6','\xb7',
'\xb8','\xb9','\xba','\xbb','\xbc','\xbd','\xbe','\xbf',
'\xc0','\xc1','\xc2','\xc3','\xc4','\xc5','\xc6','\xc7',
'\xc8','\xc9','\xca','\xcb','\xcc','\xcd','\xce','\xcf',
'\xd0','\xd1','\xd2','\xd3','\xd4','\xd5','\xd6','\xd7',
'\xd8','\xd9','\xda','\xdb','\xdc','\xdd','\xde','\xdf',
'\xe0','\xe1','\xe2','\xe3','\xe4','\xe5','\xe6','\xe7',
'\xe8','\xe9','\xea','\xeb','\xec','\xed','\xee','\xef',
'\xf0','\xf1','\xf2','\xf3','\xf4','\xf5','\xf6','\xf7',
'\xf8','\xf9','\xfa','\xfb','\xfc','\xfd','\xfe','\xff'
],QUADio_0) ;
 TMP_139=RepresentIID(Hex_0,TMP_134 ,QUADio_0) ;
 TMP_142=iotaCCI(['\x00','\x01','\x02','\x03','\x04','\x05','\x06','\x07',
'\x08','\x09','\x0a','\x0b','\x0c','\x0d','\x0e','\x0f',
'\x10','\x11','\x12','\x13','\x14','\x15','\x16','\x17',
'\x18','\x19','\x1a','\x1b','\x1c','\x1d','\x1e','\x1f',
'\x20','\x21','\x22','\x23','\x24','\x25','\x26','\x27',
'\x28','\x29','\x2a','\x2b','\x2c','\x2d','\x2e','\x2f',
'\x30','\x31','\x32','\x33','\x34','\x35','\x36','\x37',
'\x38','\x39','\x3a','\x3b','\x3c','\x3d','\x3e','\x3f',
'\x40','\x41','\x42','\x43','\x44','\x45','\x46','\x47',
'\x48','\x49','\x4a','\x4b','\x4c','\x4d','\x4e','\x4f',
'\x50','\x51','\x52','\x53','\x54','\x55','\x56','\x57',
'\x58','\x59','\x5a','\x5b','\x5c','\x5d','\x5e','\x5f',
'\x60','\x61','\x62','\x63','\x64','\x65','\x66','\x67',
'\x68','\x69','\x6a','\x6b','\x6c','\x6d','\x6e','\x6f',
'\x70','\x71','\x72','\x73','\x74','\x75','\x76','\x77',
'\x78','\x79','\x7a','\x7b','\x7c','\x7d','\x7e','\x7f',
'\x80','\x81','\x82','\x83','\x84','\x85','\x86','\x87',
'\x88','\x89','\x8a','\x8b','\x8c','\x8d','\x8e','\x8f',
'\x90','\x91','\x92','\x93','\x94','\x95','\x96','\x97',
'\x98','\x99','\x9a','\x9b','\x9c','\x9d','\x9e','\x9f',
'\xa0','\xa1','\xa2','\xa3','\xa4','\xa5','\xa6','\xa7',
'\xa8','\xa9','\xaa','\xab','\xac','\xad','\xae','\xaf',
'\xb0','\xb1','\xb2','\xb3','\xb4','\xb5','\xb6','\xb7',
'\xb8','\xb9','\xba','\xbb','\xbc','\xbd','\xbe','\xbf',
'\xc0','\xc1','\xc2','\xc3','\xc4','\xc5','\xc6','\xc7',
'\xc8','\xc9','\xca','\xcb','\xcc','\xcd','\xce','\xcf',
'\xd0','\xd1','\xd2','\xd3','\xd4','\xd5','\xd6','\xd7',
'\xd8','\xd9','\xda','\xdb','\xdc','\xdd','\xde','\xdf',
'\xe0','\xe1','\xe2','\xe3','\xe4','\xe5','\xe6','\xe7',
'\xe8','\xe9','\xea','\xeb','\xec','\xed','\xee','\xef',
'\xf0','\xf1','\xf2','\xf3','\xf4','\xf5','\xf6','\xf7',
'\xf8','\xf9','\xfa','\xfb','\xfc','\xfd','\xfe','\xff'
],['\x00','\x01','\x02','\x03','\x04','\x05','\x06','\x07',
'\x08','\x09','\x0a','\x0b','\x0c','\x0d','\x0e','\x0f',
'\x10','\x11','\x12','\x13','\x14','\x15','\x16','\x17',
'\x18','\x19','\x1a','\x1b','\x1c','\x1d','\x1e','\x1f',
'\x20','\x21','\x22','\x23','\x24','\x25','\x26','\x27',
'\x28','\x29','\x2a','\x2b','\x2c','\x2d','\x2e','\x2f',
'\x30','\x31','\x32','\x33','\x34','\x35','\x36','\x37',
'\x38','\x39','\x3a','\x3b','\x3c','\x3d','\x3e','\x3f',
'\x40','\x41','\x42','\x43','\x44','\x45','\x46','\x47',
'\x48','\x49','\x4a','\x4b','\x4c','\x4d','\x4e','\x4f',
'\x50','\x51','\x52','\x53','\x54','\x55','\x56','\x57',
'\x58','\x59','\x5a','\x5b','\x5c','\x5d','\x5e','\x5f',
'\x60','\x61','\x62','\x63','\x64','\x65','\x66','\x67',
'\x68','\x69','\x6a','\x6b','\x6c','\x6d','\x6e','\x6f',
'\x70','\x71','\x72','\x73','\x74','\x75','\x76','\x77',
'\x78','\x79','\x7a','\x7b','\x7c','\x7d','\x7e','\x7f',
'\x80','\x81','\x82','\x83','\x84','\x85','\x86','\x87',
'\x88','\x89','\x8a','\x8b','\x8c','\x8d','\x8e','\x8f',
'\x90','\x91','\x92','\x93','\x94','\x95','\x96','\x97',
'\x98','\x99','\x9a','\x9b','\x9c','\x9d','\x9e','\x9f',
'\xa0','\xa1','\xa2','\xa3','\xa4','\xa5','\xa6','\xa7',
'\xa8','\xa9','\xaa','\xab','\xac','\xad','\xae','\xaf',
'\xb0','\xb1','\xb2','\xb3','\xb4','\xb5','\xb6','\xb7',
'\xb8','\xb9','\xba','\xbb','\xbc','\xbd','\xbe','\xbf',
'\xc0','\xc1','\xc2','\xc3','\xc4','\xc5','\xc6','\xc7',
'\xc8','\xc9','\xca','\xcb','\xcc','\xcd','\xce','\xcf',
'\xd0','\xd1','\xd2','\xd3','\xd4','\xd5','\xd6','\xd7',
'\xd8','\xd9','\xda','\xdb','\xdc','\xdd','\xde','\xdf',
'\xe0','\xe1','\xe2','\xe3','\xe4','\xe5','\xe6','\xe7',
'\xe8','\xe9','\xea','\xeb','\xec','\xed','\xee','\xef',
'\xf0','\xf1','\xf2','\xf3','\xf4','\xf5','\xf6','\xf7',
'\xf8','\xf9','\xfa','\xfb','\xfc','\xfd','\xfe','\xff'
],QUADio_0) ;
 TMP_143=utakIII(Hex_0,TMP_142 ) ;
 TMP_145=sameIDB(TMP_143 ,TMP_139 ,QUADct_0) ;
 r_1=comaBBB(r_0,TMP_145 ) ;
 TMP_149=iotaCCI(['\x00','\x01','\x02','\x03','\x04','\x05','\x06','\x07',
'\x08','\x09','\x0a','\x0b','\x0c','\x0d','\x0e','\x0f',
'\x10','\x11','\x12','\x13','\x14','\x15','\x16','\x17',
'\x18','\x19','\x1a','\x1b','\x1c','\x1d','\x1e','\x1f',
'\x20','\x21','\x22','\x23','\x24','\x25','\x26','\x27',
'\x28','\x29','\x2a','\x2b','\x2c','\x2d','\x2e','\x2f',
'\x30','\x31','\x32','\x33','\x34','\x35','\x36','\x37',
'\x38','\x39','\x3a','\x3b','\x3c','\x3d','\x3e','\x3f',
'\x40','\x41','\x42','\x43','\x44','\x45','\x46','\x47',
'\x48','\x49','\x4a','\x4b','\x4c','\x4d','\x4e','\x4f',
'\x50','\x51','\x52','\x53','\x54','\x55','\x56','\x57',
'\x58','\x59','\x5a','\x5b','\x5c','\x5d','\x5e','\x5f',
'\x60','\x61','\x62','\x63','\x64','\x65','\x66','\x67',
'\x68','\x69','\x6a','\x6b','\x6c','\x6d','\x6e','\x6f',
'\x70','\x71','\x72','\x73','\x74','\x75','\x76','\x77',
'\x78','\x79','\x7a','\x7b','\x7c','\x7d','\x7e','\x7f',
'\x80','\x81','\x82','\x83','\x84','\x85','\x86','\x87',
'\x88','\x89','\x8a','\x8b','\x8c','\x8d','\x8e','\x8f',
'\x90','\x91','\x92','\x93','\x94','\x95','\x96','\x97',
'\x98','\x99','\x9a','\x9b','\x9c','\x9d','\x9e','\x9f',
'\xa0','\xa1','\xa2','\xa3','\xa4','\xa5','\xa6','\xa7',
'\xa8','\xa9','\xaa','\xab','\xac','\xad','\xae','\xaf',
'\xb0','\xb1','\xb2','\xb3','\xb4','\xb5','\xb6','\xb7',
'\xb8','\xb9','\xba','\xbb','\xbc','\xbd','\xbe','\xbf',
'\xc0','\xc1','\xc2','\xc3','\xc4','\xc5','\xc6','\xc7',
'\xc8','\xc9','\xca','\xcb','\xcc','\xcd','\xce','\xcf',
'\xd0','\xd1','\xd2','\xd3','\xd4','\xd5','\xd6','\xd7',
'\xd8','\xd9','\xda','\xdb','\xdc','\xdd','\xde','\xdf',
'\xe0','\xe1','\xe2','\xe3','\xe4','\xe5','\xe6','\xe7',
'\xe8','\xe9','\xea','\xeb','\xec','\xed','\xee','\xef',
'\xf0','\xf1','\xf2','\xf3','\xf4','\xf5','\xf6','\xf7',
'\xf8','\xf9','\xfa','\xfb','\xfc','\xfd','\xfe','\xff'
],['\x00','\x01','\x02','\x03','\x04','\x05','\x06','\x07',
'\x08','\x09','\x0a','\x0b','\x0c','\x0d','\x0e','\x0f',
'\x10','\x11','\x12','\x13','\x14','\x15','\x16','\x17',
'\x18','\x19','\x1a','\x1b','\x1c','\x1d','\x1e','\x1f',
'\x20','\x21','\x22','\x23','\x24','\x25','\x26','\x27',
'\x28','\x29','\x2a','\x2b','\x2c','\x2d','\x2e','\x2f',
'\x30','\x31','\x32','\x33','\x34','\x35','\x36','\x37',
'\x38','\x39','\x3a','\x3b','\x3c','\x3d','\x3e','\x3f',
'\x40','\x41','\x42','\x43','\x44','\x45','\x46','\x47',
'\x48','\x49','\x4a','\x4b','\x4c','\x4d','\x4e','\x4f',
'\x50','\x51','\x52','\x53','\x54','\x55','\x56','\x57',
'\x58','\x59','\x5a','\x5b','\x5c','\x5d','\x5e','\x5f',
'\x60','\x61','\x62','\x63','\x64','\x65','\x66','\x67',
'\x68','\x69','\x6a','\x6b','\x6c','\x6d','\x6e','\x6f',
'\x70','\x71','\x72','\x73','\x74','\x75','\x76','\x77',
'\x78','\x79','\x7a','\x7b','\x7c','\x7d','\x7e','\x7f',
'\x80','\x81','\x82','\x83','\x84','\x85','\x86','\x87',
'\x88','\x89','\x8a','\x8b','\x8c','\x8d','\x8e','\x8f',
'\x90','\x91','\x92','\x93','\x94','\x95','\x96','\x97',
'\x98','\x99','\x9a','\x9b','\x9c','\x9d','\x9e','\x9f',
'\xa0','\xa1','\xa2','\xa3','\xa4','\xa5','\xa6','\xa7',
'\xa8','\xa9','\xaa','\xab','\xac','\xad','\xae','\xaf',
'\xb0','\xb1','\xb2','\xb3','\xb4','\xb5','\xb6','\xb7',
'\xb8','\xb9','\xba','\xbb','\xbc','\xbd','\xbe','\xbf',
'\xc0','\xc1','\xc2','\xc3','\xc4','\xc5','\xc6','\xc7',
'\xc8','\xc9','\xca','\xcb','\xcc','\xcd','\xce','\xcf',
'\xd0','\xd1','\xd2','\xd3','\xd4','\xd5','\xd6','\xd7',
'\xd8','\xd9','\xda','\xdb','\xdc','\xdd','\xde','\xdf',
'\xe0','\xe1','\xe2','\xe3','\xe4','\xe5','\xe6','\xe7',
'\xe8','\xe9','\xea','\xeb','\xec','\xed','\xee','\xef',
'\xf0','\xf1','\xf2','\xf3','\xf4','\xf5','\xf6','\xf7',
'\xf8','\xf9','\xfa','\xfb','\xfc','\xfd','\xfe','\xff'
],QUADio_0) ;
 TMP_150=rhoIII(8,2) ;
 TMP_155=RepresentCLONE0IID(TMP_150 ,TMP_149 ,QUADio_0) ;
 TMP_158=iotaCCI(['\x00','\x01','\x02','\x03','\x04','\x05','\x06','\x07',
'\x08','\x09','\x0a','\x0b','\x0c','\x0d','\x0e','\x0f',
'\x10','\x11','\x12','\x13','\x14','\x15','\x16','\x17',
'\x18','\x19','\x1a','\x1b','\x1c','\x1d','\x1e','\x1f',
'\x20','\x21','\x22','\x23','\x24','\x25','\x26','\x27',
'\x28','\x29','\x2a','\x2b','\x2c','\x2d','\x2e','\x2f',
'\x30','\x31','\x32','\x33','\x34','\x35','\x36','\x37',
'\x38','\x39','\x3a','\x3b','\x3c','\x3d','\x3e','\x3f',
'\x40','\x41','\x42','\x43','\x44','\x45','\x46','\x47',
'\x48','\x49','\x4a','\x4b','\x4c','\x4d','\x4e','\x4f',
'\x50','\x51','\x52','\x53','\x54','\x55','\x56','\x57',
'\x58','\x59','\x5a','\x5b','\x5c','\x5d','\x5e','\x5f',
'\x60','\x61','\x62','\x63','\x64','\x65','\x66','\x67',
'\x68','\x69','\x6a','\x6b','\x6c','\x6d','\x6e','\x6f',
'\x70','\x71','\x72','\x73','\x74','\x75','\x76','\x77',
'\x78','\x79','\x7a','\x7b','\x7c','\x7d','\x7e','\x7f',
'\x80','\x81','\x82','\x83','\x84','\x85','\x86','\x87',
'\x88','\x89','\x8a','\x8b','\x8c','\x8d','\x8e','\x8f',
'\x90','\x91','\x92','\x93','\x94','\x95','\x96','\x97',
'\x98','\x99','\x9a','\x9b','\x9c','\x9d','\x9e','\x9f',
'\xa0','\xa1','\xa2','\xa3','\xa4','\xa5','\xa6','\xa7',
'\xa8','\xa9','\xaa','\xab','\xac','\xad','\xae','\xaf',
'\xb0','\xb1','\xb2','\xb3','\xb4','\xb5','\xb6','\xb7',
'\xb8','\xb9','\xba','\xbb','\xbc','\xbd','\xbe','\xbf',
'\xc0','\xc1','\xc2','\xc3','\xc4','\xc5','\xc6','\xc7',
'\xc8','\xc9','\xca','\xcb','\xcc','\xcd','\xce','\xcf',
'\xd0','\xd1','\xd2','\xd3','\xd4','\xd5','\xd6','\xd7',
'\xd8','\xd9','\xda','\xdb','\xdc','\xdd','\xde','\xdf',
'\xe0','\xe1','\xe2','\xe3','\xe4','\xe5','\xe6','\xe7',
'\xe8','\xe9','\xea','\xeb','\xec','\xed','\xee','\xef',
'\xf0','\xf1','\xf2','\xf3','\xf4','\xf5','\xf6','\xf7',
'\xf8','\xf9','\xfa','\xfb','\xfc','\xfd','\xfe','\xff'
],['\x00','\x01','\x02','\x03','\x04','\x05','\x06','\x07',
'\x08','\x09','\x0a','\x0b','\x0c','\x0d','\x0e','\x0f',
'\x10','\x11','\x12','\x13','\x14','\x15','\x16','\x17',
'\x18','\x19','\x1a','\x1b','\x1c','\x1d','\x1e','\x1f',
'\x20','\x21','\x22','\x23','\x24','\x25','\x26','\x27',
'\x28','\x29','\x2a','\x2b','\x2c','\x2d','\x2e','\x2f',
'\x30','\x31','\x32','\x33','\x34','\x35','\x36','\x37',
'\x38','\x39','\x3a','\x3b','\x3c','\x3d','\x3e','\x3f',
'\x40','\x41','\x42','\x43','\x44','\x45','\x46','\x47',
'\x48','\x49','\x4a','\x4b','\x4c','\x4d','\x4e','\x4f',
'\x50','\x51','\x52','\x53','\x54','\x55','\x56','\x57',
'\x58','\x59','\x5a','\x5b','\x5c','\x5d','\x5e','\x5f',
'\x60','\x61','\x62','\x63','\x64','\x65','\x66','\x67',
'\x68','\x69','\x6a','\x6b','\x6c','\x6d','\x6e','\x6f',
'\x70','\x71','\x72','\x73','\x74','\x75','\x76','\x77',
'\x78','\x79','\x7a','\x7b','\x7c','\x7d','\x7e','\x7f',
'\x80','\x81','\x82','\x83','\x84','\x85','\x86','\x87',
'\x88','\x89','\x8a','\x8b','\x8c','\x8d','\x8e','\x8f',
'\x90','\x91','\x92','\x93','\x94','\x95','\x96','\x97',
'\x98','\x99','\x9a','\x9b','\x9c','\x9d','\x9e','\x9f',
'\xa0','\xa1','\xa2','\xa3','\xa4','\xa5','\xa6','\xa7',
'\xa8','\xa9','\xaa','\xab','\xac','\xad','\xae','\xaf',
'\xb0','\xb1','\xb2','\xb3','\xb4','\xb5','\xb6','\xb7',
'\xb8','\xb9','\xba','\xbb','\xbc','\xbd','\xbe','\xbf',
'\xc0','\xc1','\xc2','\xc3','\xc4','\xc5','\xc6','\xc7',
'\xc8','\xc9','\xca','\xcb','\xcc','\xcd','\xce','\xcf',
'\xd0','\xd1','\xd2','\xd3','\xd4','\xd5','\xd6','\xd7',
'\xd8','\xd9','\xda','\xdb','\xdc','\xdd','\xde','\xdf',
'\xe0','\xe1','\xe2','\xe3','\xe4','\xe5','\xe6','\xe7',
'\xe8','\xe9','\xea','\xeb','\xec','\xed','\xee','\xef',
'\xf0','\xf1','\xf2','\xf3','\xf4','\xf5','\xf6','\xf7',
'\xf8','\xf9','\xfa','\xfb','\xfc','\xfd','\xfe','\xff'
],QUADio_0) ;
 TMP_159=rhoIII(8,2) ;
 TMP_160=utakIII(TMP_159 ,TMP_158 ) ;
 TMP_162=sameIDB(TMP_160 ,TMP_155 ,QUADct_0) ;
 r_2=comaBBB(r_1,TMP_162 ) ;
 TMP_169=BaseValueCLONE0III(2,1000,QUADct_0) ;
 TMP_170=dtakIII(2,1000) ;
 TMP_172=sameIIB(TMP_170 ,TMP_169 ) ;
 r_3=comaBBB(r_2,TMP_172 ) ;
 TMP_175=iotaXII( 10,QUADio_0) ;
 TMP_181=BaseValueCLONE1III(2,TMP_175 ,QUADct_0) ;
 TMP_183=iotaXII( 10,QUADio_0) ;
 TMP_184=dtakIII(2,TMP_183 ) ;
 TMP_186=sameIIB(TMP_184 ,TMP_181 ) ;
 r_4=comaBBB(r_3,TMP_186 ) ;
 TMP_193=BaseValueCLONE2III([0, 100, 100],[10, 20, 30],QUADct_0) ;
 TMP_194=dtakIII([0, 100, 100],[10, 20, 30]) ;
 TMP_196=sameIIB(TMP_194 ,TMP_193 ) ;
 r_5=comaBBB(r_4,TMP_196 ) ;
 TMP_203=BaseValueCLONE3III([2, 3],M4_0,QUADct_0) ;
 TMP_204=dtakIII([2, 3],M4_0) ;
 TMP_206=sameIIB(TMP_204 ,TMP_203 ) ;
 r_6=comaBBB(r_5,TMP_206 ) ;
 TMP_212=RepresentCLONE1IID([2, 3],M4_0,QUADio_0) ;
 TMP_213=utakIII([2, 3],M4_0) ;
 TMP_215=sameIDB(TMP_213 ,TMP_212 ,QUADct_0) ;
 r_7=comaBBB(r_6,TMP_215 ) ;
 TMP_221=RepresentCLONE2IID([10, 10, 10],3247,QUADio_0) ;
 TMP_222=utakIII([10, 10, 10],3247) ;
 TMP_224=sameIDB(TMP_222 ,TMP_221 ,QUADct_0) ;
 r_8=comaBBB(r_7,TMP_224 ) ;
 TMP_230=RepresentCLONE3IID([0, 10, 10],3247,QUADio_0) ;
 TMP_231=utakIII([0, 10, 10],3247) ;
 TMP_233=sameIDB(TMP_231 ,TMP_230 ,QUADct_0) ;
 r_9=comaBBB(r_8,TMP_233 ) ;
 TMP_239=RepresentCLONE4IDD([10, 10, 10],247.33,QUADio_0) ;
 TMP_240=utakIDD([10, 10, 10],247.33) ;
 TMP_242=sameDDB(TMP_240 ,TMP_239 ,QUADct_0) ;
 r_10=comaBBB(r_9,TMP_242 ) ;
 TMP_248=RepresentCLONE5BDD([false,true],247.33,QUADio_0) ;
 TMP_249=utakBDD([false,true],247.33) ;
 TMP_251=sameDDB(TMP_249 ,TMP_248 ,QUADct_0) ;
 r_11=comaBBB(r_10,TMP_251 ) ;
 TMP_253=rhoIII(4,10) ;
 TMP_258=RepresentCLONE6IID(TMP_253 ,-1,QUADio_0) ;
 TMP_259=rhoIII(4,10) ;
 TMP_260=utakIII(TMP_259 ,-1) ;
 TMP_262=sameIDB(TMP_260 ,TMP_258 ,QUADct_0) ;
 r_12=comaBBB(r_11,TMP_262 ) ;
 TMP_268=RepresentCLONE7BDD([false,true],-247.33,QUADio_0) ;
 TMP_269=utakBDD([false,true],-247.33) ;
 TMP_271=sameDDB(TMP_269 ,TMP_268 ,QUADct_0) ;
 r_13=comaBBB(r_12,TMP_271 ) ;
 TMP_273=rhoIII(4,10) ;
 TMP_278=RepresentCLONE8IID(TMP_273 ,[1984, 2001, -44],QUADio_0) ;
 TMP_279=rhoIII(4,10) ;
 TMP_280=utakIII(TMP_279 ,[1984, 2001, -44]) ;
 TMP_282=sameIDB(TMP_280 ,TMP_278 ,QUADct_0) ;
 r_14=comaBBB(r_13,TMP_282 ) ;
 TMP_286=quadXII(toi( r_14),QUADpp_0,QUADpw_0) ;
 TMP_287=andslXBB( r_14) ;
 TMP_291=plusBBI(true,TMP_287 );
 r_15=quadXII( TMP_291 ,QUADpp_0,QUADpw_0) ;
 r=r_15;
return(r_15+TMP_286);
}