| 1282 | |
| 1283 | |
| 1284 | static double ratio(Bigint *a, Bigint *b) |
| 1285 | { |
| 1286 | U da, db; |
| 1287 | int k, ka, kb; |
| 1288 | |
| 1289 | dval(&da)= b2d(a, &ka); |
| 1290 | dval(&db)= b2d(b, &kb); |
| 1291 | k= ka - kb + 32*(a->wds - b->wds); |
| 1292 | if (k > 0) |
| 1293 | word0(&da)+= k*Exp_msk1; |
| 1294 | else |
| 1295 | { |
| 1296 | k= -k; |
| 1297 | word0(&db)+= k*Exp_msk1; |
| 1298 | } |
| 1299 | return dval(&da) / dval(&db); |
| 1300 | } |
| 1301 | |
| 1302 | static const double tens[] = |
| 1303 | { |