/* Copyright (C) 1993 by Thomas Glen Smith.	All Rights Reserved. */
/* take2 APL2 V1.0.0 ***************************************************
* APL2 take.												 *
***********************************************************************/
#define INCLUDES APLCB
#include "includes.h"
Aplcb take2(left, rite, axes)
Aplcb left, rite, axes;
{
	Endoper; Errinit; Errstop; Getcb; Intcopy; Integer; Iscalar;
	Ravel; Reshape; Shape; Take; Vectin;
	extern int aplerr, indxorg;
	Aplcb newleft=NULL,out=NULL;
	int i,*ip,j,k,rank;
	static int izero=0;

	if (axes == NULL) return(take(left,rite));
	for(;;) {
		if (errinit()) break;
		if (left->aplrank > 1) {aplerr = 15; break;}
		if (!(APLINT & left->aplflags))
			if (NULL == (left = integer(left))) break;
		if (!(APLINT & axes->aplflags))
			if (NULL == (axes = integer(axes))) break;
		if (left->aplcount != axes->aplcount ||
			left->aplcount > rite->aplrank) {aplerr = 16; break;}
		newleft = getcb(NULL,rite->aplrank,APLINT+APLTEMP,1,NULL);
		if (newleft == NULL) break;
		ip = intcopy(newleft->aplptr.aplint, /* Default take all. */
			rite->apldim,newleft->aplcount,1);
		for(i = 0; i < axes->aplcount; i++) {
			j = *(axes->aplptr.aplint + i) - indxorg; /* axis, rel 0 */
			if (j < 0 || j >= rite->aplrank) aplerr = 9; /* bad axis */
			else *(newleft->aplptr.aplint + j) =
				*(left->aplptr.aplint + i);
		}
		if (aplerr) break;
		endoper(left);
		endoper(axes);
		return(take(newleft,rite));
	}
	endoper(newleft);
	endoper(axes);
	return(errstop(0,left,rite,NULL));
}