|
Packit |
5c3484 |
/* mpf expression evaluation
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
Copyright 2000-2002, 2004 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 |
|
|
Packit |
5c3484 |
/* Future: Bitwise "&", "|" and "&" could be done, if desired. Not sure
|
|
Packit |
5c3484 |
those functions would be much value though. */
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#include <ctype.h>
|
|
Packit |
5c3484 |
#include <stdio.h>
|
|
Packit |
5c3484 |
#include <string.h>
|
|
Packit |
5c3484 |
|
|
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 size_t
|
|
Packit |
5c3484 |
e_mpf_number (mpf_ptr res, const char *e, size_t elen, int base)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
char *edup;
|
|
Packit |
5c3484 |
size_t i, ret, extra=0;
|
|
Packit |
5c3484 |
int mant_base, exp_base;
|
|
Packit |
5c3484 |
void *(*allocate_func) (size_t);
|
|
Packit |
5c3484 |
void (*free_func) (void *, size_t);
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
TRACE (printf ("mpf_number base=%d \"%.*s\"\n", base, (int) elen, e));
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* mpf_set_str doesn't currently accept 0x for hex in base==0, so do it
|
|
Packit |
5c3484 |
here instead. FIXME: Would prefer to let mpf_set_str handle this. */
|
|
Packit |
5c3484 |
if (base == 0 && elen >= 2 && e[0] == '0' && (e[1] == 'x' || e[1] == 'X'))
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
base = 16;
|
|
Packit |
5c3484 |
extra = 2;
|
|
Packit |
5c3484 |
e += extra;
|
|
Packit |
5c3484 |
elen -= extra;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (base == 0)
|
|
Packit |
5c3484 |
mant_base = 10;
|
|
Packit |
5c3484 |
else if (base < 0)
|
|
Packit |
5c3484 |
mant_base = -base;
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
mant_base = base;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* exponent in decimal if base is negative */
|
|
Packit |
5c3484 |
if (base < 0)
|
|
Packit |
5c3484 |
exp_base = 10;
|
|
Packit |
5c3484 |
else if (base == 0)
|
|
Packit |
5c3484 |
exp_base = 10;
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
exp_base = base;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
#define IS_EXPONENT(c) \
|
|
Packit |
5c3484 |
(c == '@' || (base <= 10 && base >= -10 && (e[i] == 'e' || e[i] == 'E')))
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
i = 0;
|
|
Packit |
5c3484 |
for (;;)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
if (i >= elen)
|
|
Packit |
5c3484 |
goto parsed;
|
|
Packit |
5c3484 |
if (e[i] == '.')
|
|
Packit |
5c3484 |
break;
|
|
Packit |
5c3484 |
if (IS_EXPONENT (e[i]))
|
|
Packit |
5c3484 |
goto exponent;
|
|
Packit |
5c3484 |
if (! isasciidigit_in_base (e[i], mant_base))
|
|
Packit |
5c3484 |
goto parsed;
|
|
Packit |
5c3484 |
i++;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* fraction */
|
|
Packit |
5c3484 |
i++;
|
|
Packit |
5c3484 |
for (;;)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
if (i >= elen)
|
|
Packit |
5c3484 |
goto parsed;
|
|
Packit |
5c3484 |
if (IS_EXPONENT (e[i]))
|
|
Packit |
5c3484 |
goto exponent;
|
|
Packit |
5c3484 |
if (! isasciidigit_in_base (e[i], mant_base))
|
|
Packit |
5c3484 |
goto parsed;
|
|
Packit |
5c3484 |
i++;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
exponent:
|
|
Packit |
5c3484 |
i++;
|
|
Packit |
5c3484 |
if (i >= elen)
|
|
Packit |
5c3484 |
goto parsed;
|
|
Packit |
5c3484 |
if (e[i] == '-')
|
|
Packit |
5c3484 |
i++;
|
|
Packit |
5c3484 |
for (;;)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
if (i >= elen)
|
|
Packit |
5c3484 |
goto parsed;
|
|
Packit |
5c3484 |
if (! isasciidigit_in_base (e[i], exp_base))
|
|
Packit |
5c3484 |
break;
|
|
Packit |
5c3484 |
i++;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
parsed:
|
|
Packit |
5c3484 |
TRACE (printf (" parsed i=%u \"%.*s\"\n", i, (int) i, e));
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
mp_get_memory_functions (&allocate_func, NULL, &free_func);
|
|
Packit |
5c3484 |
edup = (*allocate_func) (i+1);
|
|
Packit |
5c3484 |
memcpy (edup, e, i);
|
|
Packit |
5c3484 |
edup[i] = '\0';
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
if (mpf_set_str (res, edup, base) == 0)
|
|
Packit |
5c3484 |
ret = i + extra;
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
ret = 0;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
(*free_func) (edup, i+1);
|
|
Packit |
5c3484 |
return ret;
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
static int
|
|
Packit |
5c3484 |
e_mpf_ulong_p (mpf_srcptr f)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
return mpf_integer_p (f) && mpf_fits_ulong_p (f);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
/* Don't want to change the precision of w, can only do an actual swap when
|
|
Packit |
5c3484 |
w and x have the same precision. */
|
|
Packit |
5c3484 |
static void
|
|
Packit |
5c3484 |
e_mpf_set_or_swap (mpf_ptr w, mpf_ptr x)
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
if (mpf_get_prec (w) == mpf_get_prec (x))
|
|
Packit |
5c3484 |
mpf_swap (w, x);
|
|
Packit |
5c3484 |
else
|
|
Packit |
5c3484 |
mpf_set (w, x);
|
|
Packit |
5c3484 |
}
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
int
|
|
Packit |
5c3484 |
mpf_expr_a (const struct mpexpr_operator_t *table,
|
|
Packit |
5c3484 |
mpf_ptr res, int base, unsigned long prec,
|
|
Packit |
5c3484 |
const char *e, size_t elen,
|
|
Packit |
5c3484 |
mpf_srcptr var[26])
|
|
Packit |
5c3484 |
{
|
|
Packit |
5c3484 |
struct mpexpr_parse_t p;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
p.table = table;
|
|
Packit |
5c3484 |
p.res = (mpX_ptr) res;
|
|
Packit |
5c3484 |
p.base = base;
|
|
Packit |
5c3484 |
p.prec = prec;
|
|
Packit |
5c3484 |
p.e = e;
|
|
Packit |
5c3484 |
p.elen = elen;
|
|
Packit |
5c3484 |
p.var = (mpX_srcptr *) var;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
p.mpX_clear = (mpexpr_fun_one_t) mpf_clear;
|
|
Packit |
5c3484 |
p.mpX_ulong_p = (mpexpr_fun_i_unary_t) e_mpf_ulong_p;
|
|
Packit |
5c3484 |
p.mpX_get_ui = (mpexpr_fun_get_ui_t) mpf_get_ui;
|
|
Packit |
5c3484 |
p.mpX_init = (mpexpr_fun_unary_ui_t) mpf_init2;
|
|
Packit |
5c3484 |
p.mpX_number = (mpexpr_fun_number_t) e_mpf_number;
|
|
Packit |
5c3484 |
p.mpX_set = (mpexpr_fun_unary_t) mpf_set;
|
|
Packit |
5c3484 |
p.mpX_set_or_swap = (mpexpr_fun_unary_t) e_mpf_set_or_swap;
|
|
Packit |
5c3484 |
p.mpX_set_si = (mpexpr_fun_set_si_t) mpf_set_si;
|
|
Packit |
5c3484 |
p.mpX_swap = (mpexpr_fun_swap_t) mpf_swap;
|
|
Packit |
5c3484 |
|
|
Packit |
5c3484 |
return mpexpr_evaluate (&p);
|
|
Packit |
5c3484 |
}
|