Blame mpz/pow_ui.c

Packit 5c3484
/* mpz_pow_ui -- mpz raised to ulong.
Packit 5c3484
Packit 5c3484
Copyright 2001, 2008 Free Software Foundation, Inc.
Packit 5c3484
Packit 5c3484
This file is part of the GNU MP Library.
Packit 5c3484
Packit 5c3484
The GNU MP Library is free software; you can redistribute it and/or modify
Packit 5c3484
it under the terms of either:
Packit 5c3484
Packit 5c3484
  * the GNU Lesser General Public License as published by the Free
Packit 5c3484
    Software Foundation; either version 3 of the License, or (at your
Packit 5c3484
    option) any later version.
Packit 5c3484
Packit 5c3484
or
Packit 5c3484
Packit 5c3484
  * the GNU General Public License as published by the Free Software
Packit 5c3484
    Foundation; either version 2 of the License, or (at your option) any
Packit 5c3484
    later version.
Packit 5c3484
Packit 5c3484
or both in parallel, as here.
Packit 5c3484
Packit 5c3484
The GNU MP Library is distributed in the hope that it will be useful, but
Packit 5c3484
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
Packit 5c3484
or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
Packit 5c3484
for more details.
Packit 5c3484
Packit 5c3484
You should have received copies of the GNU General Public License and the
Packit 5c3484
GNU Lesser General Public License along with the GNU MP Library.  If not,
Packit 5c3484
see https://www.gnu.org/licenses/.  */
Packit 5c3484
Packit 5c3484
#include "gmp.h"
Packit 5c3484
#include "gmp-impl.h"
Packit 5c3484
Packit 5c3484
void
Packit 5c3484
mpz_pow_ui (mpz_ptr r, mpz_srcptr b, unsigned long int e)
Packit 5c3484
{
Packit 5c3484
  /* We test some small exponents here, mainly to avoid the overhead of
Packit 5c3484
     mpz_n_pow_ui for small bases and exponents.  */
Packit 5c3484
  switch (e)
Packit 5c3484
    {
Packit 5c3484
    case 0:
Packit 5c3484
      mpz_set_ui (r, 1);
Packit 5c3484
      break;
Packit 5c3484
    case 1:
Packit 5c3484
      mpz_set (r, b);
Packit 5c3484
      break;
Packit 5c3484
    case 2:
Packit 5c3484
      mpz_mul (r, b, b);
Packit 5c3484
      break;
Packit 5c3484
    default:
Packit 5c3484
      mpz_n_pow_ui (r, PTR(b), (mp_size_t) SIZ(b), e);
Packit 5c3484
    }
Packit 5c3484
}