#include "kappa.h"
#include "local_proto.h"


/* function prototypes */
static int same_cats (int a, int b, int nl);


/* within group totals:
    *ns is the first stat
	(updated upon return to point to next stat)
     nl is the layer number (or level) */

long count_sum(int *ns, int nl)
{
  long count;
  int k,n;

  k = n = *ns;
  count = 0;

  if (nl >= 0) {
    while(n < nstats && same_cats(k,n,nl))
      count += Gstats[n++].count;
  }
  else {
    while(n < nstats)
      count += Gstats[n++].count;
  }

  *ns = n;

  return count;
}


static int same_cats (int a, int b, int nl)
{
  long *cat_a,*cat_b;

  cat_a = Gstats[a].cats;
  cat_b = Gstats[b].cats;

  while (nl-- >= 0)
    if (*cat_a++ != *cat_b++) return 0;

  return 1;
}