Geografic-Information-System/raster/r.kappa/sum.c

#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;
}