1/* Reference rational routines. 2 3Copyright 2001 Free Software Foundation, Inc. 4 5This file is part of the GNU MP Library. 6 7The GNU MP Library is free software; you can redistribute it and/or modify 8it under the terms of the GNU Lesser General Public License as published by 9the Free Software Foundation; either version 3 of the License, or (at your 10option) any later version. 11 12The GNU MP Library is distributed in the hope that it will be useful, but 13WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 14or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public 15License for more details. 16 17You should have received a copy of the GNU Lesser General Public License 18along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */ 19 20#include "gmp.h" 21#include "gmp-impl.h" 22#include "tests.h" 23 24 25void 26refmpq_add (mpq_ptr w, mpq_srcptr x, mpq_srcptr y) 27{ 28 mpz_mul (mpq_numref(w), mpq_numref(x), mpq_denref(y)); 29 mpz_addmul (mpq_numref(w), mpq_denref(x), mpq_numref(y)); 30 mpz_mul (mpq_denref(w), mpq_denref(x), mpq_denref(y)); 31 mpq_canonicalize (w); 32} 33 34void 35refmpq_sub (mpq_ptr w, mpq_srcptr x, mpq_srcptr y) 36{ 37 mpz_mul (mpq_numref(w), mpq_numref(x), mpq_denref(y)); 38 mpz_submul (mpq_numref(w), mpq_denref(x), mpq_numref(y)); 39 mpz_mul (mpq_denref(w), mpq_denref(x), mpq_denref(y)); 40 mpq_canonicalize (w); 41} 42