Home | History | Annotate | Line # | Download | only in expr 
      1 /* mpf expression evaluation
      2 
      3 Copyright 2000-2002 Free Software Foundation, Inc.
      4 
      5 This file is part of the GNU MP Library.
      6 
      7 The GNU MP Library is free software; you can redistribute it and/or modify
      8 it under the terms of either:
      9 
     10   * the GNU Lesser General Public License as published by the Free
     11     Software Foundation; either version 3 of the License, or (at your
     12     option) any later version.
     13 
     14 or
     15 
     16   * the GNU General Public License as published by the Free Software
     17     Foundation; either version 2 of the License, or (at your option) any
     18     later version.
     19 
     20 or both in parallel, as here.
     21 
     22 The GNU MP Library is distributed in the hope that it will be useful, but
     23 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
     24 or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
     25 for more details.
     26 
     27 You should have received copies of the GNU General Public License and the
     28 GNU Lesser General Public License along with the GNU MP Library.  If not,
     29 see https://www.gnu.org/licenses/.  */
     30 
     31 #include <stdio.h>
     32 #include <string.h>
     33 #include "gmp.h"
     34 #include "expr-impl.h"
     35 
     36 
     37 /* Change this to "#define TRACE(x) x" to get some traces. */
     38 #define TRACE(x)
     39 
     40 
     41 static int
     42 e_mpf_sgn (mpf_srcptr x)
     43 {
     44   return mpf_sgn (x);
     45 }
     46 
     47 
     48 static const struct mpexpr_operator_t  _mpf_expr_standard_table[] = {
     49 
     50   { "**",  (mpexpr_fun_t) mpf_pow_ui,
     51     MPEXPR_TYPE_BINARY_UI | MPEXPR_TYPE_RIGHTASSOC,                   220 },
     52 
     53   { "!",   (mpexpr_fun_t) e_mpf_sgn,
     54     MPEXPR_TYPE_LOGICAL_NOT | MPEXPR_TYPE_PREFIX,                     210 },
     55   { "-",   (mpexpr_fun_t) mpf_neg,
     56     MPEXPR_TYPE_UNARY | MPEXPR_TYPE_PREFIX,                           210 },
     57 
     58   { "*",   (mpexpr_fun_t) mpf_mul,           MPEXPR_TYPE_BINARY,      200 },
     59   { "/",   (mpexpr_fun_t) mpf_div,           MPEXPR_TYPE_BINARY,      200 },
     60 
     61   { "+",   (mpexpr_fun_t) mpf_add,           MPEXPR_TYPE_BINARY,      190 },
     62   { "-",   (mpexpr_fun_t) mpf_sub,           MPEXPR_TYPE_BINARY,      190 },
     63 
     64   { "<<",  (mpexpr_fun_t) mpf_mul_2exp,      MPEXPR_TYPE_BINARY_UI,   180 },
     65   { ">>",  (mpexpr_fun_t) mpf_div_2exp,      MPEXPR_TYPE_BINARY_UI,   180 },
     66 
     67   { "<=",  (mpexpr_fun_t) mpf_cmp,           MPEXPR_TYPE_CMP_LE,      170 },
     68   { "<",   (mpexpr_fun_t) mpf_cmp,           MPEXPR_TYPE_CMP_LT,      170 },
     69   { ">=",  (mpexpr_fun_t) mpf_cmp,           MPEXPR_TYPE_CMP_GE,      170 },
     70   { ">",   (mpexpr_fun_t) mpf_cmp,           MPEXPR_TYPE_CMP_GT,      170 },
     71 
     72   { "==",  (mpexpr_fun_t) mpf_cmp,           MPEXPR_TYPE_CMP_EQ,      160 },
     73   { "!=",  (mpexpr_fun_t) mpf_cmp,           MPEXPR_TYPE_CMP_NE,      160 },
     74 
     75   { "&&",  (mpexpr_fun_t) e_mpf_sgn,         MPEXPR_TYPE_LOGICAL_AND, 120 },
     76   { "||",  (mpexpr_fun_t) e_mpf_sgn,         MPEXPR_TYPE_LOGICAL_OR,  110 },
     77 
     78   { ":",   NULL,                             MPEXPR_TYPE_COLON,       101 },
     79   { "?",   (mpexpr_fun_t) e_mpf_sgn,         MPEXPR_TYPE_QUESTION,    100 },
     80 
     81   { ")",   NULL,                             MPEXPR_TYPE_CLOSEPAREN,    4 },
     82   { "(",   NULL,                             MPEXPR_TYPE_OPENPAREN,     3 },
     83   { ",",   NULL,                             MPEXPR_TYPE_ARGSEP,        2 },
     84   { "$",   NULL,                             MPEXPR_TYPE_VARIABLE,      1 },
     85 
     86   { "abs",      (mpexpr_fun_t) mpf_abs,          MPEXPR_TYPE_UNARY        },
     87   { "ceil",     (mpexpr_fun_t) mpf_ceil,         MPEXPR_TYPE_UNARY        },
     88   { "cmp",      (mpexpr_fun_t) mpf_cmp,          MPEXPR_TYPE_I_BINARY     },
     89   { "eq",       (mpexpr_fun_t) mpf_eq,           MPEXPR_TYPE_I_TERNARY_UI },
     90   { "floor",    (mpexpr_fun_t) mpf_floor,        MPEXPR_TYPE_UNARY        },
     91   { "integer_p",(mpexpr_fun_t) mpf_integer_p,    MPEXPR_TYPE_I_UNARY      },
     92   { "max",   (mpexpr_fun_t) mpf_cmp, MPEXPR_TYPE_MAX | MPEXPR_TYPE_PAIRWISE },
     93   { "min",   (mpexpr_fun_t) mpf_cmp, MPEXPR_TYPE_MIN | MPEXPR_TYPE_PAIRWISE },
     94   { "reldiff",  (mpexpr_fun_t) mpf_reldiff,      MPEXPR_TYPE_BINARY       },
     95   { "sgn",      (mpexpr_fun_t) e_mpf_sgn,        MPEXPR_TYPE_I_UNARY      },
     96   { "sqrt",     (mpexpr_fun_t) mpf_sqrt,         MPEXPR_TYPE_UNARY        },
     97   { "trunc",    (mpexpr_fun_t) mpf_trunc,        MPEXPR_TYPE_UNARY        },
     98 
     99   { NULL }
    100 };
    101 
    102 const struct mpexpr_operator_t * const mpf_expr_standard_table
    103 = _mpf_expr_standard_table;
    104 
    105 
    106 int
    107 mpf_expr (mpf_ptr res, int base, const char *e, ...)
    108 {
    109   mpf_srcptr  var[MPEXPR_VARIABLES];
    110   va_list     ap;
    111   int         ret;
    112   va_start (ap, e);
    113 
    114   TRACE (printf ("mpf_expr(): base %d, %s\n", base, e));
    115   ret = mpexpr_va_to_var ((void **) var, ap);
    116   va_end (ap);
    117 
    118   if (ret != MPEXPR_RESULT_OK)
    119     return ret;
    120 
    121   return mpf_expr_a (mpf_expr_standard_table, res, base,
    122 		     mpf_get_prec (res), e, strlen(e), var);
    123 }
    124