diff options
Diffstat (limited to 'mpi/make-logtab.txr')
| -rw-r--r-- | mpi/make-logtab.txr | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/mpi/make-logtab.txr b/mpi/make-logtab.txr new file mode 100644 index 00000000..5566dd66 --- /dev/null +++ b/mpi/make-logtab.txr @@ -0,0 +1,26 @@ +@(bind logs @(tuples 4 ^(0.0 0.0 ,*(take 63 (mapcar* [chain log2 /] (range 2)))))) +@(do (set *pprint-flo-format* "~0,9f")) +@(output) +/* + * A table of the logs of 2 for various bases (the 0 and 1 entries of + * this table are meaningless and should not be referenced). + * + * This table is used to compute output lengths for the mp_toradix() + * function. Since a number n in radix r takes up about log_r(n) + * digits, we estimate the output size by taking the least integer + * greater than log_r(n), where: + * + * log_r(n) = log_2(n) * log_r(2) + * + * This table, therefore, is a table of log_r(2) for 2 <= r <= 36, + * which are the output bases supported. + */ + +const double s_logv_2[] = { +@ (repeat) + @(rep)@logs, @(last)@logs,@(end) +@ (last) + @(rep)@logs, @(last)@logs@(end) +@ (end) +}; +@(end) |