From 0d98678c1b2c3db9c6bd1860e9740c1928470c87 Mon Sep 17 00:00:00 2001 From: Mattias Andrée Date: Thu, 3 Mar 2016 10:45:50 +0100 Subject: Add new functions: 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/zpowu.c | 32 ++++++++++++++++++++++++++++++++ 1 file changed, 32 insertions(+) create mode 100644 src/zpowu.c (limited to 'src/zpowu.c') diff --git a/src/zpowu.c b/src/zpowu.c new file mode 100644 index 0000000..bc0c6dd --- /dev/null +++ b/src/zpowu.c @@ -0,0 +1,32 @@ +/* See LICENSE file for copyright and license details. */ +#include "internals" + +#include + +#define tb libzahl_tmp_pow_b + + +void +zpowu(z_t a, z_t b, unsigned long long int c) +{ + if (!c) { + if (zzero(b)) { + errno = EDOM; /* Indeterminate form: 0:th power of 0 */ + FAILURE_JUMP(); + } + zsetu(a, 1); + return; + } else if (zzero(b)) { + SET_SIGNUM(a, 0); + return; + } + + zset(tb, b); + zsetu(a, 1); + + for (; c; c >>= 1) { + if (c & 1) + zmul(a, a, tb); + zsqr(tb, tb); + } +} -- cgit v1.2.3-70-g09d2