X-Git-Url: https://git.distorted.org.uk/u/mdw/catacomb/blobdiff_plain/ddcb6aa7914ebf7ccb2f82baba4451a954abf3d4..52cdaca9b7b8e15871a3de8b5b571c18dcc580f8:/mpx-ksqr.c diff --git a/mpx-ksqr.c b/mpx-ksqr.c index 4aa84dd..25dbb70 100644 --- a/mpx-ksqr.c +++ b/mpx-ksqr.c @@ -1,6 +1,6 @@ /* -*-c-*- * - * $Id: mpx-ksqr.c,v 1.5 2000/10/08 12:11:01 mdw Exp $ + * $Id: mpx-ksqr.c,v 1.6 2000/10/08 15:48:35 mdw Exp $ * * Karatsuba-based squaring algorithm * @@ -30,6 +30,9 @@ /*----- Revision history --------------------------------------------------* * * $Log: mpx-ksqr.c,v $ + * Revision 1.6 2000/10/08 15:48:35 mdw + * Rename Karatsuba constants now that we have @gfx_kmul@ too. + * * Revision 1.5 2000/10/08 12:11:01 mdw * Use @mpx_ueq@ instead of @MPX_UCMP@. * @@ -55,13 +58,13 @@ #include #include "mpx.h" -#include "mpx-kmac.h" +#include "karatsuba.h" /*----- Tweakables --------------------------------------------------------*/ #ifdef TEST_RIG -# undef KARATSUBA_CUTOFF -# define KARATSUBA_CUTOFF 2 +# undef MPK_THRESH +# define MPK_THRESH 2 #endif /*----- Main code ---------------------------------------------------------*/ @@ -82,7 +85,7 @@ * * The destination must be twice as large as the argument. The * scratch space must be twice as large as the argument, plus - * the magic number @KARATSUBA_SLOP@. + * the magic number @MPK_SLOP@. */ void mpx_ksqr(mpw *dv, mpw *dvl, @@ -102,7 +105,7 @@ void mpx_ksqr(mpw *dv, mpw *dvl, MPX_SHRINK(av, avl); - if (avl - av <= KARATSUBA_CUTOFF) { + if (avl - av <= MPK_THRESH) { mpx_usqr(dv, dvl, av, avl); return; } @@ -136,12 +139,12 @@ void mpx_ksqr(mpw *dv, mpw *dvl, mpw *rdv = tdv + m; UADD2(sv, svm, av, avm, avm, avl); - if (m > KARATSUBA_CUTOFF) + if (m > MPK_THRESH) mpx_ksqr(tdv, rdv + m + 4, sv, svm + 1, ssv, svl); else mpx_usqr(tdv, rdv + m + 4, sv, svm + 1); - if (m > KARATSUBA_CUTOFF) + if (m > MPK_THRESH) mpx_ksqr(sv, ssv, avm, avl, ssv, svl); else mpx_usqr(sv, ssv, avm, avl); @@ -149,7 +152,7 @@ void mpx_ksqr(mpw *dv, mpw *dvl, UADD(rdv, sv, svm + 1); USUB(tdv, sv, svn); - if (m > KARATSUBA_CUTOFF) + if (m > MPK_THRESH) mpx_ksqr(sv, ssv, av, avm, ssv, svl); else mpx_usqr(sv, ssv, av, avm);