Blame mpq/inv.c

Packit 5c3484
/* mpq_inv(dest,src) -- invert a rational number, i.e. set DEST to SRC
Packit 5c3484
   with the numerator and denominator swapped.
Packit 5c3484
Packit 5c3484
Copyright 1991, 1994, 1995, 2000, 2001, 2012 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
mpq_inv (mpq_ptr dest, mpq_srcptr src)
Packit 5c3484
{
Packit 5c3484
  mp_size_t num_size = SIZ(NUM(src));
Packit 5c3484
  mp_size_t den_size = SIZ(DEN(src));
Packit 5c3484
Packit 5c3484
  if (num_size < 0)
Packit 5c3484
    {
Packit 5c3484
      num_size = -num_size;
Packit 5c3484
      den_size = -den_size;
Packit 5c3484
    }
Packit 5c3484
  else if (UNLIKELY (num_size == 0))
Packit 5c3484
    DIVIDE_BY_ZERO;
Packit 5c3484
Packit 5c3484
  SIZ(DEN(dest)) = num_size;
Packit 5c3484
  SIZ(NUM(dest)) = den_size;
Packit 5c3484
Packit 5c3484
  /* If dest == src we may just swap the numerator and denominator;
Packit 5c3484
     we ensured that the new denominator is positive.  */
Packit 5c3484
Packit 5c3484
  if (dest == src)
Packit 5c3484
    {
Packit 5c3484
      MP_PTR_SWAP (PTR(NUM(dest)), PTR(DEN(dest)));
Packit 5c3484
      MP_SIZE_T_SWAP (ALLOC(NUM(dest)), ALLOC(DEN(dest)));
Packit 5c3484
    }
Packit 5c3484
  else
Packit 5c3484
    {
Packit 5c3484
      mp_ptr dp;
Packit 5c3484
Packit 5c3484
      den_size = ABS (den_size);
Packit 5c3484
      dp = MPZ_NEWALLOC (NUM(dest), den_size);
Packit 5c3484
      MPN_COPY (dp, PTR(DEN(src)), den_size);
Packit 5c3484
Packit 5c3484
      dp = MPZ_NEWALLOC (DEN(dest), num_size);
Packit 5c3484
      MPN_COPY (dp, PTR(NUM(src)), num_size);
Packit 5c3484
    }
Packit 5c3484
}