|
Packit |
5c3484 |
/* mpq expression evaluation
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
Copyright 2000-2002 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 <stdio.h>
|
|
Packit |
5c3484 |
#include <string.h>
|
|
Packit |
5c3484 |
#include "gmp.h"
|
|
Packit |
5c3484 |
#include "expr-impl.h"
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* Change this to "#define TRACE(x) x" to get some traces. */
|
|
Packit |
5c3484 |
#define TRACE(x)
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
static void
|
|
Packit |
5c3484 |
e_mpq_pow_ui (mpq_ptr r, mpq_srcptr b, unsigned long e)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mpz_pow_ui (mpq_numref(r), mpq_numref(b), e);
|
|
Packit |
5c3484 |
mpz_pow_ui (mpq_denref(r), mpq_denref(b), e);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* Wrapped because mpq_sgn is a macro. */
|
|
Packit |
5c3484 |
static int
|
|
Packit |
5c3484 |
e_mpq_sgn (mpq_srcptr x)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
return mpq_sgn (x);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* Wrapped because mpq_equal only guarantees a non-zero return, whereas we
|
|
Packit |
5c3484 |
want 1 or 0 for == and !=. */
|
|
Packit |
5c3484 |
static int
|
|
Packit |
5c3484 |
e_mpq_equal (mpq_srcptr x, mpq_srcptr y)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
return mpq_equal (x, y) != 0;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
static int
|
|
Packit |
5c3484 |
e_mpq_notequal (mpq_srcptr x, mpq_srcptr y)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
return ! mpq_equal (x, y);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
static void
|
|
Packit |
5c3484 |
e_mpq_num (mpq_ptr w, mpq_srcptr x)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
if (w != x)
|
|
Packit |
5c3484 |
mpz_set (mpq_numref(w), mpq_numref(x));
|
|
Packit |
5c3484 |
mpz_set_ui (mpq_denref(w), 1L);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
static void
|
|
Packit |
5c3484 |
e_mpq_den (mpq_ptr w, mpq_srcptr x)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
if (w == x)
|
|
Packit |
5c3484 |
mpz_swap (mpq_numref(w), mpq_denref(w));
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
mpz_set (mpq_numref(w), mpq_denref(x));
|
|
Packit |
5c3484 |
mpz_set_ui (mpq_denref(w), 1L);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
static const struct mpexpr_operator_t _mpq_expr_standard_table[] = {
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ "**", (mpexpr_fun_t) e_mpq_pow_ui,
|
|
Packit |
5c3484 |
MPEXPR_TYPE_BINARY_UI | MPEXPR_TYPE_RIGHTASSOC, 220 },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ "!", (mpexpr_fun_t) e_mpq_sgn,
|
|
Packit |
5c3484 |
MPEXPR_TYPE_LOGICAL_NOT | MPEXPR_TYPE_PREFIX, 210 },
|
|
Packit |
5c3484 |
{ "-", (mpexpr_fun_t) mpq_neg,
|
|
Packit |
5c3484 |
MPEXPR_TYPE_UNARY | MPEXPR_TYPE_PREFIX, 210 },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ "*", (mpexpr_fun_t) mpq_mul, MPEXPR_TYPE_BINARY, 200 },
|
|
Packit |
5c3484 |
{ "/", (mpexpr_fun_t) mpq_div, MPEXPR_TYPE_BINARY, 200 },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ "+", (mpexpr_fun_t) mpq_add, MPEXPR_TYPE_BINARY, 190 },
|
|
Packit |
5c3484 |
{ "-", (mpexpr_fun_t) mpq_sub, MPEXPR_TYPE_BINARY, 190 },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ "<<", (mpexpr_fun_t) mpq_mul_2exp, MPEXPR_TYPE_BINARY_UI, 180 },
|
|
Packit |
5c3484 |
{ ">>", (mpexpr_fun_t) mpq_div_2exp, MPEXPR_TYPE_BINARY_UI, 180 },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ "<=", (mpexpr_fun_t) mpq_cmp, MPEXPR_TYPE_CMP_LE, 170 },
|
|
Packit |
5c3484 |
{ "<", (mpexpr_fun_t) mpq_cmp, MPEXPR_TYPE_CMP_LT, 170 },
|
|
Packit |
5c3484 |
{ ">=", (mpexpr_fun_t) mpq_cmp, MPEXPR_TYPE_CMP_GE, 170 },
|
|
Packit |
5c3484 |
{ ">", (mpexpr_fun_t) mpq_cmp, MPEXPR_TYPE_CMP_GT, 170 },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ "==", (mpexpr_fun_t) e_mpq_equal, MPEXPR_TYPE_I_BINARY, 160 },
|
|
Packit |
5c3484 |
{ "!=", (mpexpr_fun_t) e_mpq_notequal, MPEXPR_TYPE_I_BINARY, 160 },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ "&&", (mpexpr_fun_t) e_mpq_sgn, MPEXPR_TYPE_LOGICAL_AND, 120 },
|
|
Packit |
5c3484 |
{ "||", (mpexpr_fun_t) e_mpq_sgn, MPEXPR_TYPE_LOGICAL_OR, 110 },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ ":", NULL, MPEXPR_TYPE_COLON, 101 },
|
|
Packit |
5c3484 |
{ "?", (mpexpr_fun_t) e_mpq_sgn, MPEXPR_TYPE_QUESTION, 100 },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ ")", (mpexpr_fun_t) e_mpq_sgn, MPEXPR_TYPE_CLOSEPAREN, 4 },
|
|
Packit |
5c3484 |
{ "(", (mpexpr_fun_t) e_mpq_sgn, MPEXPR_TYPE_OPENPAREN, 3 },
|
|
Packit |
5c3484 |
{ ",", (mpexpr_fun_t) e_mpq_sgn, MPEXPR_TYPE_ARGSEP, 2 },
|
|
Packit |
5c3484 |
{ "$", NULL, MPEXPR_TYPE_VARIABLE, 1 },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ "abs", (mpexpr_fun_t) mpq_abs, MPEXPR_TYPE_UNARY },
|
|
Packit |
5c3484 |
{ "cmp", (mpexpr_fun_t) mpq_cmp, MPEXPR_TYPE_I_BINARY },
|
|
Packit |
5c3484 |
{ "den", (mpexpr_fun_t) e_mpq_den, MPEXPR_TYPE_UNARY },
|
|
Packit |
5c3484 |
{ "max", (mpexpr_fun_t) mpq_cmp, MPEXPR_TYPE_MAX | MPEXPR_TYPE_PAIRWISE },
|
|
Packit |
5c3484 |
{ "min", (mpexpr_fun_t) mpq_cmp, MPEXPR_TYPE_MIN | MPEXPR_TYPE_PAIRWISE },
|
|
Packit |
5c3484 |
{ "num", (mpexpr_fun_t) e_mpq_num, MPEXPR_TYPE_UNARY },
|
|
Packit |
5c3484 |
{ "sgn", (mpexpr_fun_t) e_mpq_sgn, MPEXPR_TYPE_I_UNARY },
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
{ NULL }
|
|
Packit |
5c3484 |
};
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
const struct mpexpr_operator_t * const mpq_expr_standard_table
|
|
Packit |
5c3484 |
= _mpq_expr_standard_table;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
int
|
|
Packit |
5c3484 |
mpq_expr (mpq_ptr res, int base, const char *e, ...)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
mpq_srcptr var[MPEXPR_VARIABLES];
|
|
Packit |
5c3484 |
va_list ap;
|
|
Packit |
5c3484 |
int ret;
|
|
Packit |
5c3484 |
va_start (ap, e);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
TRACE (printf ("mpq_expr(): base %d, %s\n", base, e));
|
|
Packit |
5c3484 |
ret = mpexpr_va_to_var ((void **) var, ap);
|
|
Packit |
5c3484 |
va_end (ap);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (ret != MPEXPR_RESULT_OK)
|
|
Packit |
5c3484 |
return ret;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
return mpq_expr_a (mpq_expr_standard_table, res, base, e, strlen(e), var);
|
|
Packit |
5c3484 |
}
|