Annotation of OpenXM_contrib/gmp/demos/expr/exprfa.c, Revision 1.1
1.1 ! ohara 1: /* mpf expression evaluation */
! 2:
! 3: /*
! 4: Copyright 2000, 2001 Free Software Foundation, Inc.
! 5:
! 6: This file is part of the GNU MP Library.
! 7:
! 8: The GNU MP Library is free software; you can redistribute it and/or modify
! 9: it under the terms of the GNU Lesser General Public License as published by
! 10: the Free Software Foundation; either version 2.1 of the License, or (at your
! 11: option) any later version.
! 12:
! 13: The GNU MP Library is distributed in the hope that it will be useful, but
! 14: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
! 15: or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
! 16: License for more details.
! 17:
! 18: You should have received a copy of the GNU Lesser General Public License
! 19: along with the GNU MP Library; see the file COPYING.LIB. If not, write to
! 20: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
! 21: MA 02111-1307, USA.
! 22: */
! 23:
! 24: /* Future: Bitwise "&", "|" and "&" could be done, if desired. Not sure
! 25: those functions would be much value though. */
! 26:
! 27:
! 28: #include <ctype.h>
! 29: #include <stdio.h>
! 30:
! 31: #include "gmp.h"
! 32: #include "expr-impl.h"
! 33:
! 34:
! 35: /* Change this to "#define TRACE(x) x" to get some traces. */
! 36: #define TRACE(x)
! 37:
! 38:
! 39: static size_t
! 40: e_mpf_number (mpf_ptr res, __gmp_const char *e, size_t elen, int base)
! 41: {
! 42: char *edup;
! 43: size_t i, ret, extra=0;
! 44: int mant_base, exp_base;
! 45:
! 46: TRACE (printf ("mpf_number base=%d \"%.*s\"\n", base, (int) elen, e));
! 47:
! 48: /* mpf_set_str doesn't currently accept 0x for hex in base==0, so do it
! 49: here instead. FIXME: Would prefer to let mpf_set_str handle this. */
! 50: if (base == 0 && elen >= 2 && e[0] == '0' && (e[1] == 'x' || e[1] == 'X'))
! 51: {
! 52: base = 16;
! 53: extra = 2;
! 54: e += extra;
! 55: elen -= extra;
! 56: }
! 57:
! 58: if (base == 0)
! 59: mant_base = 10;
! 60: else if (base < 0)
! 61: mant_base = -base;
! 62: else
! 63: mant_base = base;
! 64:
! 65: /* exponent in decimal if base is negative */
! 66: if (base < 0)
! 67: exp_base = 10;
! 68: else if (base == 0)
! 69: exp_base = 10;
! 70: else
! 71: exp_base = base;
! 72:
! 73: #define IS_EXPONENT(c) \
! 74: (c == '@' || (base <= 10 && base >= -10 && (e[i] == 'e' || e[i] == 'E')))
! 75:
! 76: i = 0;
! 77: for (;;)
! 78: {
! 79: if (i >= elen)
! 80: goto parsed;
! 81: if (e[i] == '.')
! 82: break;
! 83: if (IS_EXPONENT (e[i]))
! 84: goto exponent;
! 85: if (! isasciidigit_in_base (e[i], mant_base))
! 86: goto parsed;
! 87: i++;
! 88: }
! 89:
! 90: /* fraction */
! 91: i++;
! 92: for (;;)
! 93: {
! 94: if (i >= elen)
! 95: goto parsed;
! 96: if (IS_EXPONENT (e[i]))
! 97: goto exponent;
! 98: if (! isasciidigit_in_base (e[i], mant_base))
! 99: goto parsed;
! 100: i++;
! 101: }
! 102:
! 103: exponent:
! 104: i++;
! 105: if (i >= elen)
! 106: goto parsed;
! 107: if (e[i] == '-')
! 108: i++;
! 109: for (;;)
! 110: {
! 111: if (i >= elen)
! 112: goto parsed;
! 113: if (! isasciidigit_in_base (e[i], exp_base))
! 114: break;
! 115: i++;
! 116: }
! 117:
! 118: parsed:
! 119: TRACE (printf (" parsed i=%u \"%.*s\"\n", i, (int) i, e));
! 120:
! 121: edup = (*__gmp_allocate_func) (i+1);
! 122: memcpy (edup, e, i);
! 123: edup[i] = '\0';
! 124:
! 125: if (mpf_set_str (res, edup, base) == 0)
! 126: ret = i + extra;
! 127: else
! 128: ret = 0;
! 129:
! 130: (*__gmp_free_func) (edup, i+1);
! 131: return ret;
! 132: }
! 133:
! 134: static int
! 135: e_mpf_ulong_p (mpf_srcptr f)
! 136: {
! 137: return mpf_integer_p (f) && mpf_fits_ulong_p (f);
! 138: }
! 139:
! 140: /* Don't want to change the precision of w, can only do an actual swap when
! 141: w and x have the same precision. */
! 142: static void
! 143: e_mpf_set_or_swap (mpf_ptr w, mpf_ptr x)
! 144: {
! 145: if (mpf_get_prec (w) == mpf_get_prec (x))
! 146: mpf_swap (w, x);
! 147: else
! 148: mpf_set (w, x);
! 149: }
! 150:
! 151:
! 152: int
! 153: mpf_expr_a (__gmp_const struct mpexpr_operator_t *table,
! 154: mpf_ptr res, int base, unsigned long prec,
! 155: __gmp_const char *e, size_t elen,
! 156: mpf_srcptr var[26])
! 157: {
! 158: struct mpexpr_parse_t p;
! 159:
! 160: p.table = table;
! 161: p.res = (mpX_ptr) res;
! 162: p.base = base;
! 163: p.prec = prec;
! 164: p.e = e;
! 165: p.elen = elen;
! 166: p.var = (mpX_srcptr *) var;
! 167:
! 168: p.mpX_clear = (mpexpr_fun_one_t) mpf_clear;
! 169: p.mpX_ulong_p = (mpexpr_fun_i_unary_t) e_mpf_ulong_p;
! 170: p.mpX_get_ui = (mpexpr_fun_get_ui_t) mpf_get_ui;
! 171: p.mpX_init = (mpexpr_fun_unary_ui_t) mpf_init;
! 172: p.mpX_number = (mpexpr_fun_number_t) e_mpf_number;
! 173: p.mpX_set = (mpexpr_fun_unary_t) mpf_set;
! 174: p.mpX_set_or_swap = (mpexpr_fun_unary_t) e_mpf_set_or_swap;
! 175: p.mpX_set_si = (mpexpr_fun_set_si_t) mpf_set_si;
! 176: p.mpX_swap = (mpexpr_fun_swap_t) mpf_swap;
! 177:
! 178: return mpexpr_evaluate (&p);
! 179: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to exprfa.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / demos / expr