diff options
| author | Tobias Schlüter <tobias.schlueter@physik.uni-muenchen.de> | 2004-10-07 01:16:43 +0200 |
|---|---|---|
| committer | Tobias Schlüter <tobi@gcc.gnu.org> | 2004-10-07 01:16:43 +0200 |
| commit | c45abcaf8fe7c588468d6c2c9fb0b2fec24e3cf7 (patch) | |
| tree | d32fe02580ebb8b9e04104f199396953a5678b02 /gcc | |
| parent | 9e749a8b426c8907c288e006388e85c015231788 (diff) | |
| download | gcc-c45abcaf8fe7c588468d6c2c9fb0b2fec24e3cf7.zip gcc-c45abcaf8fe7c588468d6c2c9fb0b2fec24e3cf7.tar.gz gcc-c45abcaf8fe7c588468d6c2c9fb0b2fec24e3cf7.tar.bz2 | |
* simplify.c (gfc_simplify_abs): Use mpfr_hypot for CABS.
From-SVN: r88647
Diffstat (limited to 'gcc')
| -rw-r--r-- | gcc/fortran/ChangeLog | 2 | ||||
| -rw-r--r-- | gcc/fortran/simplify.c | 13 |
2 files changed, 4 insertions, 11 deletions
diff --git a/gcc/fortran/ChangeLog b/gcc/fortran/ChangeLog index d52914f..276e586 100644 --- a/gcc/fortran/ChangeLog +++ b/gcc/fortran/ChangeLog @@ -4,6 +4,8 @@ * simplify.c (twos_complement): New function. (gfc_simplify_ishft, gfc_simplify_ishftc): Revise. + * simplify.c (gfc_simplify_abs): Use mpfr_hypot for CABS. + 2004-10-06 Paul Brook <paul@codesourcery.com> * trans-stmt.c (gfc_trans_simple_do): New function. diff --git a/gcc/fortran/simplify.c b/gcc/fortran/simplify.c index 0ac9d61..253f686 100644 --- a/gcc/fortran/simplify.c +++ b/gcc/fortran/simplify.c @@ -176,7 +176,6 @@ gfc_expr * gfc_simplify_abs (gfc_expr * e) { gfc_expr *result; - mpfr_t a, b; if (e->expr_type != EXPR_CONSTANT) return NULL; @@ -203,17 +202,9 @@ gfc_simplify_abs (gfc_expr * e) result = gfc_constant_result (BT_REAL, e->ts.kind, &e->where); gfc_set_model_kind (e->ts.kind); - mpfr_init (a); - mpfr_init (b); - /* FIXME: Possible numerical problems. */ - mpfr_mul (a, e->value.complex.r, e->value.complex.r, GFC_RND_MODE); - mpfr_mul (b, e->value.complex.i, e->value.complex.i, GFC_RND_MODE); - mpfr_add (a, a, b, GFC_RND_MODE); - mpfr_sqrt (result->value.real, a, GFC_RND_MODE); - - mpfr_clear (a); - mpfr_clear (b); + mpfr_hypot (result->value.real, e->value.complex.r, + e->value.complex.i, GFC_RND_MODE); result = range_check (result, "CABS"); break; |