From 659a5a81045af899310dca5f70907da287c56d19 Mon Sep 17 00:00:00 2001 From: Mattias Andrée Date: Wed, 27 Apr 2016 14:35:13 +0200 Subject: Optimise zpowu and zmodpowu MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: Mattias Andrée --- src/zmodpowu.c | 6 ++++-- src/zpowu.c | 6 ++++-- 2 files changed, 8 insertions(+), 4 deletions(-) diff --git a/src/zmodpowu.c b/src/zmodpowu.c index c9c8af5..fd5e925 100644 --- a/src/zmodpowu.c +++ b/src/zmodpowu.c @@ -27,9 +27,11 @@ zmodpowu(z_t a, z_t b, unsigned long long int c, z_t d) zset(td, d); zsetu(a, 1); - for (; c; c >>= 1) { + if (c & 1) + zmodmul(a, a, tb, td); + while (c >>= 1) { + zmodsqr(tb, tb, td); if (c & 1) zmodmul(a, a, tb, td); - zmodsqr(tb, tb, td); } } diff --git a/src/zpowu.c b/src/zpowu.c index 9618010..3d58aa4 100644 --- a/src/zpowu.c +++ b/src/zpowu.c @@ -23,10 +23,12 @@ zpowu(z_t a, z_t b, unsigned long long int c) zabs(tb, b); zsetu(a, 1); - for (; c; c >>= 1) { + if (c & 1) + zmul_impl(a, a, tb); + while (c >>= 1) { + zsqr_impl(tb, tb); if (c & 1) zmul_impl(a, a, tb); - zsqr_impl(tb, tb); } if (neg) -- cgit v1.2.3-70-g09d2