/*  cog.h  (C) 1997 Mark A. Livingston */
/* Any non-commercial use of this code is permitted so long as */
/* the above copyright line appears in the file.  Please contact */
/* the author regarding permission for commercial use. */
#define X_OCT 4
#define Y_OCT 2
#define Z_OCT 1
/* first four args are: xyz-quats, weights, sample locations, num samples,
   first three are modified (the arrays are permuted), new value of fourth
   to use for subsequent ops is return value */
int cogReduce( q_xyz_quat_type *q, double *w, q_vec_type *samp_loc, int n,
               q_vec_type cog, /* input and output */
	       q_vec_type curPt );