tally.h

#ifndef TALLY_H
#define TALLY_H

#include <stdio.h>
#include <stdint.h>

#define COL1 "signal"
#define SEP "\t"
#define PRINTHDR1(col1) printf(col1 "\n")
#define PRINTHDR2(col1, col2) printf(col1 SEP col2 "\n")
#define PRINTHDR3(col1, col2, col3) printf(col1 SEP col2 SEP col3 "\n")
#define PRINTSTDHDR(col2) PRINTHDR2(COL1, col2)
#define TALLY_SZ (65536) /* 2^16 */
#define TALLY_SZ_EXP (8192) /* 2^13 */

static inline void updatetally(uint64_t *tally, int16_t x, int16_t *min,
			       int16_t *max)
{
	tally[(x << 1) ^ (x >> 15)] ++;
	if (x < *min)
		*min = x;
	if (x > *max)
		*max = x;
}
uint64_t *gettally(FILE *fp, int16_t *min, int16_t *max);
uint64_t *gettrally(FILE *fp, uint16_t *min, uint16_t *max);
void printtally(const uint64_t *tally, int16_t min, int16_t max);
void printprob(const uint64_t *tally, uint16_t min, uint16_t max);
void printparity(const uint64_t *tally, uint16_t min, uint16_t max);
void printtrally(const uint64_t *trally, uint16_t min, uint16_t max);

#endif