Add a new mp_todec_fast...

that uses Barrett reduction to speed up stringifying large integers.
libtom · Aug 23, 2019 · 78be8b8 · 78be8b8
1 parent 584405f
commit 78be8b8
Show file tree

Hide file tree

Showing 8 changed files with 261 additions and 35 deletions.
diff --git a/bn_mp_todecimal_fast.c b/bn_mp_todecimal_fast.c
@@ -0,0 +1,199 @@
+#include "tommath_private.h"
+#include <string.h>
+#ifdef BN_MP_TODECIMAL_FAST_C
+/* LibTomMath, multiple-precision integer library -- Tom St Denis */
+/* SPDX-License-Identifier: Unlicense */
+
+/* store a bignum as a decimal ASCII string */
+mp_err mp_todecimal_fast_rec(mp_int *number, mp_int *nL, mp_int *shiftL, mp_int *mL, int index, int left, char **result) {
+    mp_int q, nLq, r;
+    mp_err err;
+
+    if (index < 0) {
+        char *next_piece = calloc(4, sizeof(char));
+        int s_s = snprintf(next_piece, 4, left ? "%u" : "%03u", mp_get_u32(number));
+        int r_s = strlen(*result);
+        (*result) = realloc(*result, r_s + s_s + 2);
+        strcat(*result, next_piece);
+        return MP_OKAY;
+    }
+
+    if ((err = mp_init_multi(&q, &nLq, &r, NULL)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+    if ((err = mp_mul(number, &mL[index], &q)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+    if ((err = mp_div_2d(&q, mp_get_u32(&shiftL[index]), &q, NULL)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    if ((err = mp_mul(&nL[index], &q, &nLq)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    if ((err = mp_sub(number, &nLq, &r)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    if (mp_isneg(&r)) {
+        if ((err = mp_sub_d(&q, 1, &q)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+        if ((err = mp_add(&r, &nL[index], &r)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+    }
+
+    --index;
+    if (left && mp_iszero(&q)) {
+        if ((err = mp_todecimal_fast_rec(&r, nL, shiftL, mL, index, 1, result)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+    } else {
+        if ((err = mp_todecimal_fast_rec(&q, nL, shiftL, mL, index, left, result)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+        if ((err = mp_todecimal_fast_rec(&r, nL, shiftL, mL, index, 0, result)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+    }
+
+    err = MP_OKAY;
+
+LBL_ERR:mp_clear_multi (&q, &nLq, &r, NULL);
+    return err;
+}
+
+mp_err mp_todecimal_fast(mp_int *number, char **result) {
+    mp_int n, shift, M, M2, M22, M4, M44;
+    mp_int *nL, *shiftL, *mL;
+    mp_err err;
+    int index = 1;
+
+    if ((err = mp_init_multi(&M2, &M22, &M4, &M44, NULL)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    if (mp_isneg(number)) {
+        if ((err = mp_neg(number, number)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+        *result[0] = '-';
+    }
+    if ((err = mp_init_set(&n, 1000)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    nL = malloc(20 * sizeof(mp_int));
+    if ((err = mp_init_copy(&nL[0], &n)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    if ((err = mp_init_set(&shift, 20)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    shiftL = malloc(20 * sizeof(mp_int));
+    if ((err = mp_init_copy(&shiftL[0], &shift)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    /* (8 * 2**$shift) / $n rounded up */
+    if ((err = mp_init_set(&M, 8389)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    /* $M / 8, rounded up */
+    mL = malloc(20 * sizeof(mp_int));
+    if ((err = mp_init_set(&mL[0], 1049)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    while (1) {
+        if ((err = mp_sqr(&n, &n)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+        if (mp_cmp(&n, number) == MP_GT) {
+            break;
+        }
+
+        if ((err = mp_mul_2(&shift, &shift)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+
+        /* The following is a Newton-Raphson step, to restore the invariant
+         * that $M is (8 * 2**$shift) / $n, rounded up. */
+        {
+            if ((err = mp_sqr(&M, &M2)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_sqr(&M2, &M4)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+
+            if ((err = mp_mul(&M4, &n, &M4)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_div_2d(&M4, mp_get_ul(&shift) + 6, &M4, NULL)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_mul_2(&M2, &M2)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_sub(&M4, &M2, &M4)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_add_d(&M4, 1, &M4)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_div_2d(&M4, 3, &M4, NULL)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_sub_d(&M4, 1, &M4)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_neg(&M4, &M)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+        }
+
+        if ((err = mp_init_copy(&nL[index], &n)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+        if ((err = mp_init_copy(&shiftL[index], &shift)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+
+        /* Divide by 8, round up */
+        {
+            if ((err = mp_add_d(&M4, 1, &M4)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_div_2d(&M4, 3, &M4, NULL)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_sub_d(&M4, 1, &M4)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+            if ((err = mp_neg(&M4, &M4)) != MP_OKAY) {
+                goto LBL_ERR;
+            }
+        }
+        if ((err = mp_init_copy(&mL[index], &M4)) != MP_OKAY) {
+            goto LBL_ERR;
+        }
+        index++;
+    }
+
+    if ((err = mp_todecimal_fast_rec(number, nL, shiftL, mL, index - 1, 1, result)) != MP_OKAY) {
+        goto LBL_ERR;
+    }
+
+    err = MP_OKAY;
+
+LBL_ERR:mp_clear_multi (&n, &shift, &M, &M2, &M22, &M4, &M44, NULL);
+    return err;
+}
+
+#endif
diff --git a/libtommath_VS2008.vcproj b/libtommath_VS2008.vcproj
@@ -832,6 +832,10 @@
 			RelativePath="bn_mp_to_unsigned_bin_n.c"
 			>
 		</File>
+		<File
+			RelativePath="bn_mp_todecimal_fast.c"
+			>
+		</File>
 		<File
 			RelativePath="bn_mp_toradix.c"
 			>

diff --git a/makefile b/makefile
@@ -49,13 +49,13 @@ bn_mp_set.o bn_mp_set_double.o bn_mp_set_i32.o bn_mp_set_i64.o bn_mp_set_l.o bn_
 bn_mp_set_u32.o bn_mp_set_u64.o bn_mp_set_ul.o bn_mp_set_ull.o bn_mp_shrink.o bn_mp_signed_bin_size.o \
 bn_mp_signed_rsh.o bn_mp_sqr.o bn_mp_sqrmod.o bn_mp_sqrt.o bn_mp_sqrtmod_prime.o bn_mp_sub.o bn_mp_sub_d.o \
 bn_mp_submod.o bn_mp_to_signed_bin.o bn_mp_to_signed_bin_n.o bn_mp_to_unsigned_bin.o \
-bn_mp_to_unsigned_bin_n.o bn_mp_toradix.o bn_mp_toradix_n.o bn_mp_unsigned_bin_size.o bn_mp_xor.o \
-bn_mp_zero.o bn_prime_tab.o bn_s_mp_add.o bn_s_mp_balance_mul.o bn_s_mp_exptmod.o bn_s_mp_exptmod_fast.o \
-bn_s_mp_get_bit.o bn_s_mp_invmod_fast.o bn_s_mp_invmod_slow.o bn_s_mp_karatsuba_mul.o \
-bn_s_mp_karatsuba_sqr.o bn_s_mp_montgomery_reduce_fast.o bn_s_mp_mul_digs.o bn_s_mp_mul_digs_fast.o \
-bn_s_mp_mul_high_digs.o bn_s_mp_mul_high_digs_fast.o bn_s_mp_prime_is_divisible.o \
-bn_s_mp_rand_jenkins.o bn_s_mp_rand_platform.o bn_s_mp_reverse.o bn_s_mp_sqr.o bn_s_mp_sqr_fast.o \
-bn_s_mp_sub.o bn_s_mp_toom_mul.o bn_s_mp_toom_sqr.o
+bn_mp_to_unsigned_bin_n.o bn_mp_todecimal_fast.o bn_mp_toradix.o bn_mp_toradix_n.o \
+bn_mp_unsigned_bin_size.o bn_mp_xor.o bn_mp_zero.o bn_prime_tab.o bn_s_mp_add.o bn_s_mp_balance_mul.o \
+bn_s_mp_exptmod.o bn_s_mp_exptmod_fast.o bn_s_mp_get_bit.o bn_s_mp_invmod_fast.o bn_s_mp_invmod_slow.o \
+bn_s_mp_karatsuba_mul.o bn_s_mp_karatsuba_sqr.o bn_s_mp_montgomery_reduce_fast.o bn_s_mp_mul_digs.o \
+bn_s_mp_mul_digs_fast.o bn_s_mp_mul_high_digs.o bn_s_mp_mul_high_digs_fast.o \
+bn_s_mp_prime_is_divisible.o bn_s_mp_rand_jenkins.o bn_s_mp_rand_platform.o bn_s_mp_reverse.o \
+bn_s_mp_sqr.o bn_s_mp_sqr_fast.o bn_s_mp_sub.o bn_s_mp_toom_mul.o bn_s_mp_toom_sqr.o
 
 #END_INS
 

diff --git a/makefile.mingw b/makefile.mingw
@@ -52,13 +52,13 @@ bn_mp_set.o bn_mp_set_double.o bn_mp_set_i32.o bn_mp_set_i64.o bn_mp_set_l.o bn_
 bn_mp_set_u32.o bn_mp_set_u64.o bn_mp_set_ul.o bn_mp_set_ull.o bn_mp_shrink.o bn_mp_signed_bin_size.o \
 bn_mp_signed_rsh.o bn_mp_sqr.o bn_mp_sqrmod.o bn_mp_sqrt.o bn_mp_sqrtmod_prime.o bn_mp_sub.o bn_mp_sub_d.o \
 bn_mp_submod.o bn_mp_to_signed_bin.o bn_mp_to_signed_bin_n.o bn_mp_to_unsigned_bin.o \
-bn_mp_to_unsigned_bin_n.o bn_mp_toradix.o bn_mp_toradix_n.o bn_mp_unsigned_bin_size.o bn_mp_xor.o \
-bn_mp_zero.o bn_prime_tab.o bn_s_mp_add.o bn_s_mp_balance_mul.o bn_s_mp_exptmod.o bn_s_mp_exptmod_fast.o \
-bn_s_mp_get_bit.o bn_s_mp_invmod_fast.o bn_s_mp_invmod_slow.o bn_s_mp_karatsuba_mul.o \
-bn_s_mp_karatsuba_sqr.o bn_s_mp_montgomery_reduce_fast.o bn_s_mp_mul_digs.o bn_s_mp_mul_digs_fast.o \
-bn_s_mp_mul_high_digs.o bn_s_mp_mul_high_digs_fast.o bn_s_mp_prime_is_divisible.o \
-bn_s_mp_rand_jenkins.o bn_s_mp_rand_platform.o bn_s_mp_reverse.o bn_s_mp_sqr.o bn_s_mp_sqr_fast.o \
-bn_s_mp_sub.o bn_s_mp_toom_mul.o bn_s_mp_toom_sqr.o
+bn_mp_to_unsigned_bin_n.o bn_mp_todecimal_fast.o bn_mp_toradix.o bn_mp_toradix_n.o \
+bn_mp_unsigned_bin_size.o bn_mp_xor.o bn_mp_zero.o bn_prime_tab.o bn_s_mp_add.o bn_s_mp_balance_mul.o \
+bn_s_mp_exptmod.o bn_s_mp_exptmod_fast.o bn_s_mp_get_bit.o bn_s_mp_invmod_fast.o bn_s_mp_invmod_slow.o \
+bn_s_mp_karatsuba_mul.o bn_s_mp_karatsuba_sqr.o bn_s_mp_montgomery_reduce_fast.o bn_s_mp_mul_digs.o \
+bn_s_mp_mul_digs_fast.o bn_s_mp_mul_high_digs.o bn_s_mp_mul_high_digs_fast.o \
+bn_s_mp_prime_is_divisible.o bn_s_mp_rand_jenkins.o bn_s_mp_rand_platform.o bn_s_mp_reverse.o \
+bn_s_mp_sqr.o bn_s_mp_sqr_fast.o bn_s_mp_sub.o bn_s_mp_toom_mul.o bn_s_mp_toom_sqr.o
 
 HEADERS_PUB=tommath.h
 HEADERS=tommath_private.h tommath_class.h tommath_superclass.h $(HEADERS_PUB)

diff --git a/makefile.msvc b/makefile.msvc
@@ -44,13 +44,13 @@ bn_mp_set.obj bn_mp_set_double.obj bn_mp_set_i32.obj bn_mp_set_i64.obj bn_mp_set
 bn_mp_set_u32.obj bn_mp_set_u64.obj bn_mp_set_ul.obj bn_mp_set_ull.obj bn_mp_shrink.obj bn_mp_signed_bin_size.obj \
 bn_mp_signed_rsh.obj bn_mp_sqr.obj bn_mp_sqrmod.obj bn_mp_sqrt.obj bn_mp_sqrtmod_prime.obj bn_mp_sub.obj bn_mp_sub_d.obj \
 bn_mp_submod.obj bn_mp_to_signed_bin.obj bn_mp_to_signed_bin_n.obj bn_mp_to_unsigned_bin.obj \
-bn_mp_to_unsigned_bin_n.obj bn_mp_toradix.obj bn_mp_toradix_n.obj bn_mp_unsigned_bin_size.obj bn_mp_xor.obj \
-bn_mp_zero.obj bn_prime_tab.obj bn_s_mp_add.obj bn_s_mp_balance_mul.obj bn_s_mp_exptmod.obj bn_s_mp_exptmod_fast.obj \
-bn_s_mp_get_bit.obj bn_s_mp_invmod_fast.obj bn_s_mp_invmod_slow.obj bn_s_mp_karatsuba_mul.obj \
-bn_s_mp_karatsuba_sqr.obj bn_s_mp_montgomery_reduce_fast.obj bn_s_mp_mul_digs.obj bn_s_mp_mul_digs_fast.obj \
-bn_s_mp_mul_high_digs.obj bn_s_mp_mul_high_digs_fast.obj bn_s_mp_prime_is_divisible.obj \
-bn_s_mp_rand_jenkins.obj bn_s_mp_rand_platform.obj bn_s_mp_reverse.obj bn_s_mp_sqr.obj bn_s_mp_sqr_fast.obj \
-bn_s_mp_sub.obj bn_s_mp_toom_mul.obj bn_s_mp_toom_sqr.obj
+bn_mp_to_unsigned_bin_n.obj bn_mp_todecimal_fast.obj bn_mp_toradix.obj bn_mp_toradix_n.obj \
+bn_mp_unsigned_bin_size.obj bn_mp_xor.obj bn_mp_zero.obj bn_prime_tab.obj bn_s_mp_add.obj bn_s_mp_balance_mul.obj \
+bn_s_mp_exptmod.obj bn_s_mp_exptmod_fast.obj bn_s_mp_get_bit.obj bn_s_mp_invmod_fast.obj bn_s_mp_invmod_slow.obj \
+bn_s_mp_karatsuba_mul.obj bn_s_mp_karatsuba_sqr.obj bn_s_mp_montgomery_reduce_fast.obj bn_s_mp_mul_digs.obj \
+bn_s_mp_mul_digs_fast.obj bn_s_mp_mul_high_digs.obj bn_s_mp_mul_high_digs_fast.obj \
+bn_s_mp_prime_is_divisible.obj bn_s_mp_rand_jenkins.obj bn_s_mp_rand_platform.obj bn_s_mp_reverse.obj \
+bn_s_mp_sqr.obj bn_s_mp_sqr_fast.obj bn_s_mp_sub.obj bn_s_mp_toom_mul.obj bn_s_mp_toom_sqr.obj
 
 HEADERS_PUB=tommath.h
 HEADERS=tommath_private.h tommath_class.h tommath_superclass.h $(HEADERS_PUB)

diff --git a/makefile.shared b/makefile.shared
@@ -46,13 +46,13 @@ bn_mp_set.o bn_mp_set_double.o bn_mp_set_i32.o bn_mp_set_i64.o bn_mp_set_l.o bn_
 bn_mp_set_u32.o bn_mp_set_u64.o bn_mp_set_ul.o bn_mp_set_ull.o bn_mp_shrink.o bn_mp_signed_bin_size.o \
 bn_mp_signed_rsh.o bn_mp_sqr.o bn_mp_sqrmod.o bn_mp_sqrt.o bn_mp_sqrtmod_prime.o bn_mp_sub.o bn_mp_sub_d.o \
 bn_mp_submod.o bn_mp_to_signed_bin.o bn_mp_to_signed_bin_n.o bn_mp_to_unsigned_bin.o \
-bn_mp_to_unsigned_bin_n.o bn_mp_toradix.o bn_mp_toradix_n.o bn_mp_unsigned_bin_size.o bn_mp_xor.o \
-bn_mp_zero.o bn_prime_tab.o bn_s_mp_add.o bn_s_mp_balance_mul.o bn_s_mp_exptmod.o bn_s_mp_exptmod_fast.o \
-bn_s_mp_get_bit.o bn_s_mp_invmod_fast.o bn_s_mp_invmod_slow.o bn_s_mp_karatsuba_mul.o \
-bn_s_mp_karatsuba_sqr.o bn_s_mp_montgomery_reduce_fast.o bn_s_mp_mul_digs.o bn_s_mp_mul_digs_fast.o \
-bn_s_mp_mul_high_digs.o bn_s_mp_mul_high_digs_fast.o bn_s_mp_prime_is_divisible.o \
-bn_s_mp_rand_jenkins.o bn_s_mp_rand_platform.o bn_s_mp_reverse.o bn_s_mp_sqr.o bn_s_mp_sqr_fast.o \
-bn_s_mp_sub.o bn_s_mp_toom_mul.o bn_s_mp_toom_sqr.o
+bn_mp_to_unsigned_bin_n.o bn_mp_todecimal_fast.o bn_mp_toradix.o bn_mp_toradix_n.o \
+bn_mp_unsigned_bin_size.o bn_mp_xor.o bn_mp_zero.o bn_prime_tab.o bn_s_mp_add.o bn_s_mp_balance_mul.o \
+bn_s_mp_exptmod.o bn_s_mp_exptmod_fast.o bn_s_mp_get_bit.o bn_s_mp_invmod_fast.o bn_s_mp_invmod_slow.o \
+bn_s_mp_karatsuba_mul.o bn_s_mp_karatsuba_sqr.o bn_s_mp_montgomery_reduce_fast.o bn_s_mp_mul_digs.o \
+bn_s_mp_mul_digs_fast.o bn_s_mp_mul_high_digs.o bn_s_mp_mul_high_digs_fast.o \
+bn_s_mp_prime_is_divisible.o bn_s_mp_rand_jenkins.o bn_s_mp_rand_platform.o bn_s_mp_reverse.o \
+bn_s_mp_sqr.o bn_s_mp_sqr_fast.o bn_s_mp_sub.o bn_s_mp_toom_mul.o bn_s_mp_toom_sqr.o
 
 #END_INS
 

diff --git a/makefile.unix b/makefile.unix
@@ -53,13 +53,13 @@ bn_mp_set.o bn_mp_set_double.o bn_mp_set_i32.o bn_mp_set_i64.o bn_mp_set_l.o bn_
 bn_mp_set_u32.o bn_mp_set_u64.o bn_mp_set_ul.o bn_mp_set_ull.o bn_mp_shrink.o bn_mp_signed_bin_size.o \
 bn_mp_signed_rsh.o bn_mp_sqr.o bn_mp_sqrmod.o bn_mp_sqrt.o bn_mp_sqrtmod_prime.o bn_mp_sub.o bn_mp_sub_d.o \
 bn_mp_submod.o bn_mp_to_signed_bin.o bn_mp_to_signed_bin_n.o bn_mp_to_unsigned_bin.o \
-bn_mp_to_unsigned_bin_n.o bn_mp_toradix.o bn_mp_toradix_n.o bn_mp_unsigned_bin_size.o bn_mp_xor.o \
-bn_mp_zero.o bn_prime_tab.o bn_s_mp_add.o bn_s_mp_balance_mul.o bn_s_mp_exptmod.o bn_s_mp_exptmod_fast.o \
-bn_s_mp_get_bit.o bn_s_mp_invmod_fast.o bn_s_mp_invmod_slow.o bn_s_mp_karatsuba_mul.o \
-bn_s_mp_karatsuba_sqr.o bn_s_mp_montgomery_reduce_fast.o bn_s_mp_mul_digs.o bn_s_mp_mul_digs_fast.o \
-bn_s_mp_mul_high_digs.o bn_s_mp_mul_high_digs_fast.o bn_s_mp_prime_is_divisible.o \
-bn_s_mp_rand_jenkins.o bn_s_mp_rand_platform.o bn_s_mp_reverse.o bn_s_mp_sqr.o bn_s_mp_sqr_fast.o \
-bn_s_mp_sub.o bn_s_mp_toom_mul.o bn_s_mp_toom_sqr.o
+bn_mp_to_unsigned_bin_n.o bn_mp_todecimal_fast.o bn_mp_toradix.o bn_mp_toradix_n.o \
+bn_mp_unsigned_bin_size.o bn_mp_xor.o bn_mp_zero.o bn_prime_tab.o bn_s_mp_add.o bn_s_mp_balance_mul.o \
+bn_s_mp_exptmod.o bn_s_mp_exptmod_fast.o bn_s_mp_get_bit.o bn_s_mp_invmod_fast.o bn_s_mp_invmod_slow.o \
+bn_s_mp_karatsuba_mul.o bn_s_mp_karatsuba_sqr.o bn_s_mp_montgomery_reduce_fast.o bn_s_mp_mul_digs.o \
+bn_s_mp_mul_digs_fast.o bn_s_mp_mul_high_digs.o bn_s_mp_mul_high_digs_fast.o \
+bn_s_mp_prime_is_divisible.o bn_s_mp_rand_jenkins.o bn_s_mp_rand_platform.o bn_s_mp_reverse.o \
+bn_s_mp_sqr.o bn_s_mp_sqr_fast.o bn_s_mp_sub.o bn_s_mp_toom_mul.o bn_s_mp_toom_sqr.o
 
 HEADERS_PUB=tommath.h
 HEADERS=tommath_private.h tommath_class.h tommath_superclass.h $(HEADERS_PUB)

diff --git a/tommath_class.h b/tommath_class.h
@@ -141,6 +141,7 @@
 #   define BN_MP_TO_SIGNED_BIN_N_C
 #   define BN_MP_TO_UNSIGNED_BIN_C
 #   define BN_MP_TO_UNSIGNED_BIN_N_C
+#   define BN_MP_TODECIMAL_FAST_C
 #   define BN_MP_TORADIX_C
 #   define BN_MP_TORADIX_N_C
 #   define BN_MP_UNSIGNED_BIN_SIZE_C
@@ -1070,6 +1071,28 @@
 #   define BN_MP_UNSIGNED_BIN_SIZE_C
 #endif
 
+#if defined(BN_MP_TODECIMAL_FAST_C)
+#   define BN_MP_2EXPT_C
+#   define BN_MP_ADD_C
+#   define BN_MP_ADD_D_C
+#   define BN_MP_BARRETT_TODECIMAL_C
+#   define BN_MP_BARRETT_TODECIMAL_REC_C
+#   define BN_MP_CMP_C
+#   define BN_MP_DIV_2D_C
+#   define BN_MP_GET_INT_C
+#   define BN_MP_GET_LONG_C
+#   define BN_MP_INIT_C
+#   define BN_MP_INIT_COPY_C
+#   define BN_MP_INIT_MULTI_C
+#   define BN_MP_INIT_SET_C
+#   define BN_MP_MUL_2_C
+#   define BN_MP_MUL_C
+#   define BN_MP_NEG_C
+#   define BN_MP_SQR_C
+#   define BN_MP_SUB_C
+#   define BN_MP_SUB_D_C
+#endif
+
 #if defined(BN_MP_TORADIX_C)
 #   define BN_MP_CLEAR_C
 #   define BN_MP_DIV_D_C
-Original file line number
+Diff line change
@@ Expand Up / @@ -832,6 +832,10 @@ @@
     			RelativePath="bn_mp_to_unsigned_bin_n.c"
     			>
     		</File>
+    		<File
+    			RelativePath="bn_mp_todecimal_fast.c"
+    			>
+    		</File>
     		<File
     			RelativePath="bn_mp_toradix.c"
     			>
@@ Expand Down @@