Blame demos/expr/exprq.c

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
}