static void normalize( int root, int n, int64 * a )
{
int e1,e2;
while( 2 * root + 1 < n )
{
e1 = 2 * root + 1; // first
e2 = e1 + 1; // second
/*
check if second ancessor exists and greater than first
*/
#ifndef SORT_REVERSE
if( e2 < n && cmp( a[e1], a[e2] ) )
#else
if( e2 < n && rcmp( a[e1], a[e2] ) )
#endif
e1 = e2; // swap
// if a[root] < max( a[e1],a[e2] ), do a swap
#ifndef SORT_REVERSE
if( cmp( a[root], a[e1] ) )
#else
if( rcmp( a[root], a[e1] ) )
#endif
{
swap(&a[root],&a[e1]);
root = e1;
}
else
break;
}
}
int main(){
int count = 0;
char sysout[1000];
int sysrc;
int vdrc;
double value;
int w,p;
static char *format[] = { "%*.*g","%-*.*g","%#*.*g","%+*.*g", "% *.*G",
"%*.*e", "%#*.*E", "%*.*f", "%#*.*f", "%*%%*g"};
int fcnt;
for(fcnt=0; fcnt<sizeof(format)/sizeof(char*); fcnt++){
for(w=0; w<=12; w += 3){
for(p=0; p<=12; p += 3){
for(value=1.234567890123456789e-32; value<1.3e+32; value *= 1e8 ){
outidx = 0;
sprintf(sysout,format[fcnt],w,p,value);
vdrc = testprintf(format[fcnt],w,p,value);
sysrc = strlen(sysout);
outbuf[outidx] = 0;
count++;
if( vdrc!=sysrc || !rcmp(outbuf,sysout) ){
printf("Count=%d Format=%s Width=%d Precision=%d Value=%.16g\n",
count,format[fcnt],w,p,value);
printf(" system: [%s] returned %d\n",sysout,sysrc);
printf(" vxprintf: [%s] returned %d\n",outbuf,vdrc);
}
}
}
}
}
printf("Checked %d combinations.\n",count);
return 0;
}
static void rPsort2(double *x, int lo, int hi, int k)
{
double v, w;
int L, R, i, j;
for (L = lo, R = hi; L < R; ) {
v = x[k];
for(i = L, j = R; i <= j;) {
while (rcmp(x[i], v) < 0) i++;
while (rcmp(v, x[j]) < 0) j--;
if (i <= j) { w = x[i]; x[i++] = x[j]; x[j--] = w; }
}
if (j < k) L = i;
if (k < i) R = j;
}
}
void rsort_with_index(double *x, unsigned *indx, int n) {
double v;
int i, j, h, iv;
for (h = 1; h <= n / 9; h = 3 * h + 1)
;
for (; h > 0; h /= 3)
for (i = h; i < n; i++) {
v = x[i];
iv = indx[i];
j = i;
while (j >= h && rcmp(x[j - h], v) > 0) {
x[j] = x[j - h];
indx[j] = indx[j - h];
j -= h;
}
x[j] = v;
indx[j] = iv;
}
}
int rcmp(struct vertex *a, struct vertex *b)
{
#ifdef HASH
if (a->hash < b->hash)
return -1;
if (a->hash > b->hash)
return 1;
#endif
if (a->deg < b->deg)
return -1;
if (a->deg > b->deg)
return 1;
if (a->depth < b->depth)
return -1;
if (a->depth > b->depth)
return 1;
for (int i=0; i<a->deg; i++)
{
int res = rcmp(a->son[i], b->son[i]);
if (res)
return res;
}
return 0;
}
int cmp(const void *A, const void *B)
{
struct vertex * const *AA = A, * const *BB = B;
struct vertex *a = *AA, *b = *BB;
return rcmp(a, b);
}
main(int argc, char *argv[])
#endif
{
struct stat sb1, sb2;
off_t skip1, skip2;
int ch, fd1, fd2, special;
char *file1, *file2;
#ifdef __STACK_CHECK__
atexit(cleanup);
_beginStackCheck();
#endif
#ifdef __GNO__
while ((ch = getopt(argc, argv, "-lsr")) != -1)
#else
while ((ch = getopt(argc, argv, "-ls")) != -1)
#endif
switch (ch) {
case 'l': /* print all differences */
lflag = 1;
break;
case 's': /* silent run */
sflag = 1;
break;
#ifdef __GNO__
case 'r':
rflag = 1; /* ignore resource fork */
break;
#endif
case '-': /* stdin (must be after options) */
--optind;
goto endargs;
case '?':
default:
usage();
}
endargs:
#ifndef __ORCAC__
argv += optind;
#else
argv = argv + optind;
#endif
argc -= optind;
if (lflag && sflag)
errx(ERR_EXIT, "only one of -l and -s may be specified");
if (argc < 2 || argc > 4)
usage();
/* Backward compatibility -- handle "-" meaning stdin. */
special = 0;
if (strcmp(file1 = argv[0], "-") == 0) {
special = 1;
fd1 = 0;
file1 = "stdin";
}
#ifndef __GNO__
else if ((fd1 = open(file1, O_RDONLY, 0)) < 0) {
#else
/* open() for GNO/ME requires 3rd arg iff creating file */
else if ((fd1 = open(file1, O_RDONLY)) < 0) {
#endif
if (!sflag)
err(ERR_EXIT, "%s", file1);
else
exit(1);
}
if (strcmp(file2 = argv[1], "-") == 0) {
if (special)
errx(ERR_EXIT,
"standard input may only be specified once");
special = 1;
fd2 = 0;
file2 = "stdin";
}
#ifndef __GNO__
else if ((fd2 = open(file2, O_RDONLY, 0)) < 0) {
#else
else if ((fd2 = open(file2, O_RDONLY)) < 0) {
#endif
if (!sflag)
err(ERR_EXIT, "%s", file2);
else
exit(1);
}
skip1 = argc > 2 ? strtol(argv[2], NULL, 10) : 0;
skip2 = argc == 4 ? strtol(argv[3], NULL, 10) : 0;
#ifndef __GNO__
if (!special) {
if (fstat(fd1, &sb1)) {
if (!sflag)
err(ERR_EXIT, "%s", file1);
else
exit(1);
}
if (!S_ISREG(sb1.st_mode))
special = 1;
else {
if (fstat(fd2, &sb2)) {
if (!sflag)
err(ERR_EXIT, "%s", file2);
else
exit(1);
}
if (!S_ISREG(sb2.st_mode))
special = 1;
}
}
#else
special = 1; /* GNO doesn't have mmap.h, so treat every file as
a special file and process byte by byte */
#endif
if (special)
#ifdef __GNO__
{
#endif
c_special(fd1, file1, skip1, fd2, file2, skip2);
#ifdef __GNO__
close(fd1);
close(fd2);
}
if (rflag != 1) {
rflag = 2;
rcmp(&fd1, &fd2, file1, file2);
close(fd1);
close(fd2);
}
#else
else
c_regular(fd1, file1, skip1, sb1.st_size,
fd2, file2, skip2, sb2.st_size);
#endif
/* exit(0); */
}
static void
#ifndef __STDC__
usage()
#else
usage(void)
#endif
{
(void)fprintf(stderr,
#ifndef __GNO__
"usage: cmp [-l | -s] file1 file2 [skip1 [skip2]]\n");
#else
"usage: cmp [-l | -s] [-r] file1 file2 [skip1 [skip2]]\n");
#endif
exit(ERR_EXIT);
}
