+ oN= stab.st_size / (sizeof(double)*D3);
+
+ for (shift=1, inc=2;
+ shift<XBITS-1;
+ shift++, inc<<=1) {
+ printf("\n shift=%d inc=%d ",shift,inc);
+
+ if (X % inc) continue;
+ oX= X/inc; printf("oX=%d ",oX);
+
+ if ((Y-1) % inc) continue;
+ oY= (Y-1)/inc + 1; printf("oY=%d ",oY);
+
+ oN_calc= oX*oY;
+ printf("oN=%d",oN_calc);
+ if (oN_calc == oN) goto found;
+ }
+ fail("\ninput size cannot be reconciled\n");
+
+ found:
+ oXBITS= XBITS-shift;
+ printf("oXBITS=%d\n", oXBITS);