- fail("input file is not reasonable whole number of doubles\n");
-
- oldsz= stab.st_size / (sizeof(double)*D3);
- for (shift=1;
- shift < XBITS+1 && shift < YBITS+1;
- shift++) {
- oldxbits= XBITS-1;
- oldybits= YBITS-1;
- oldx= 1<<oldxbits;
- oldy= (1<<oldybits)-1;
- fprintf(stderr,"sizeof(double)=%d XYBITS=%d,%d, XY=%d*%d=%d"
- " oldsz=%d shift=%d oldxybits=%d,%d oldxy=%d*%d=%d\n",
- (int)sizeof(double), XBITS,YBITS, X,Y,N,
- oldsz,shift,oldxbits,oldybits,oldx,oldy,oldx*oldy);
- if (oldx*oldy == oldsz) goto found;
- }
- fail("input file size cannot be interpolated to target file size\n");
+ fail("input file is not reasonable whole number of vertices\n");
+
+ fprintf(stderr,"XBITS=%d XY=%d*%d=%d\n", XBITS, X,Y,N);
+ fprintf(stderr,"OXBITS=%d OXY=%d*%d=%d\n", OXBITS, OX,OY,ON);
+ fprintf(stderr,"sizeof(double)=%d st_size=%lu\n",
+ (int)sizeof(double), (unsigned long)stab.st_size);