diff options
Diffstat (limited to 'ui_sub.c')
-rw-r--r-- | ui_sub.c | 35 |
1 files changed, 35 insertions, 0 deletions
diff --git a/ui_sub.c b/ui_sub.c new file mode 100644 index 0000000..a4effc4 --- /dev/null +++ b/ui_sub.c @@ -0,0 +1,35 @@ +/* mpc_ui_sub -- Subtract a complex number from an unsigned long int. + +Copyright (C) 2005 Marc Helbling, Paul Zimmermann. + +This file is part of the MPC Library. + +The MPC Library is free software; you can redistribute it and/or modify +it under the terms of the GNU Lesser General Public License as published by +the Free Software Foundation; either version 2.1 of the License, or (at your +option) any later version. + +The MPC Library is distributed in the hope that it will be useful, but +WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY +or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the MPC Library; see the file COPYING.LIB. If not, write to +the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, +MA 02111-1307, USA. */ + +#include "gmp.h" +#include "mpfr.h" +#include "mpc.h" + +int +mpc_ui_sub (mpc_ptr a, unsigned long int b, mpc_srcptr c, mpc_rnd_t rnd) +{ + int inex_re, inex_im; + + inex_re = mpfr_ui_sub (MPC_RE (a), b, MPC_RE (c), MPC_RND_RE (rnd)); + inex_im = mpfr_neg (MPC_IM (a), MPC_IM (c), MPC_RND_IM (rnd)); + + return MPC_INEX (inex_re, inex_im); +} |