30 return(**((
long **)a)-**((
long **)b));
33 int main(
int argc,
char *argv[]){
45 fprintf(stderr,
"Usage: distribution {data.vqd [bins]| book.vqh} \n\n");
51 in=fopen(argv[1],
"r");
53 fprintf(stderr,
"Could not open input file %s\n",argv[1]);
57 if(strrchr(argv[1],
'.') && strcmp(strrchr(argv[1],
'.'),
".vqh")==0){
68 printf(
"entropy codebook only; no mappings\n");
87 printf(
"Minimum scalar value: %f\n",min);
88 printf(
"Maximum scalar value: %f\n",max);
96 long **sort=
calloc(bins,
sizeof(
long *));
98 countarray=
calloc(bins,
sizeof(
long));
100 for(i=0;i<bins;i++)sort[i]=c->
quantlist+i;
111 for(k=0;k<b->
dim;k++){
112 int index= (j/indexdiv)%bins;
126 for(k=0;k<b->
dim;k++){
128 int index= (j/indexdiv)%bins;
140 if(countarray[i]>maxcount)maxcount=countarray[
i];
144 int stars=rint(50./maxcount*countarray[ptr]);
145 printf(
"%+08f (%8ld) |",c->
quantlist[ptr]*delta+min,countarray[ptr]);
146 for(j=0;j<stars;j++)printf(
"*");
155 printf(
"Can't do probability dump of a trained [type 2] codebook (yet)\n");
169 sprintf(buf,
"getting min/max (%.2f::%.2f). lines...",min,max);
170 if(!(lines&0xff))
spinnit(buf,lines);
172 while(!flag && sscanf(line,
"%f",&code)==1){
173 line=strchr(line,
',');
178 while(line && sscanf(line,
"%f",&code)==1){
179 line=strchr(line,
',');
181 if(code<min)min=code;
182 if(code>max)max=code;
189 if((
int)(max-min)==min-max){
197 printf(
"Minimum scalar value: %f\n",min);
198 printf(
"Maximum scalar value: %f\n",max);
202 printf(
"\n counting hits into %ld bins...\n",bins+1);
203 countarray=
calloc(bins+1,
sizeof(
long));
210 if(!(lines&0xff))
spinnit(
"counting distribution. lines so far...",lines);
212 while(line && sscanf(line,
"%f",&code)==1){
213 line=strchr(line,
',');
219 countarray[(
int)rint(code)]++;
229 for(
i=0;
i<bins+1;
i++)
230 if(countarray[
i]>maxcount)maxcount=countarray[
i];
233 printf(
"Total scalars: %ld\n",total);
234 for(
i=0;
i<bins+1;
i++){
235 int stars=rint(50./maxcount*countarray[
i]);
236 printf(
"%08f (%8ld) |",(max-min)/bins*i+min,countarray[i]);
237 for(j=0;j<stars;j++)printf(
"*");
GLuint const GLfloat * val
int main(int argc, char **argv)
const static_codebook * c
GLboolean GLboolean GLboolean GLboolean a
int ascend(const void *a, const void *b)
float _float32_unpack(long val)
void spinnit(char *s, int n)
char * setup_line(FILE *in)
GLenum GLuint GLsizei const GLchar * buf
long _book_maptype1_quantvals(const static_codebook *b)
GLdouble GLdouble GLdouble b
codebook * codebook_load(char *filename)
void qsort(void *base, size_t nmemb, size_t size, int(*compare)(const void *, const void *))