From 8e70c271ec0d77772ffa04b7985d1599f064e963 Mon Sep 17 00:00:00 2001 From: "Kaveh R. Ghazi" Date: Fri, 2 Oct 2009 03:52:05 +0000 Subject: re PR fortran/33197 (Fortran 2008: math functions) PR fortran/33197 * gfortran.h (HAVE_mpc_arc): Define. * simplify.c (gfc_simplify_acos): Handle complex acos. (gfc_simplify_acosh): Likewise for acosh. (gfc_simplify_asin): Likewise for asin. (gfc_simplify_asinh): Likewise for asinh. (gfc_simplify_atan): Likewise for atan. (gfc_simplify_atanh): Likewise for atanh. From-SVN: r152394 --- gcc/fortran/simplify.c | 41 ++++++++++++++++++++++++++++++++++++++--- 1 file changed, 38 insertions(+), 3 deletions(-) (limited to 'gcc/fortran/simplify.c') diff --git a/gcc/fortran/simplify.c b/gcc/fortran/simplify.c index a22dd00..9856b2e 100644 --- a/gcc/fortran/simplify.c +++ b/gcc/fortran/simplify.c @@ -745,16 +745,21 @@ gfc_simplify_acos (gfc_expr *x) &x->where); return &gfc_bad_expr; } + result = gfc_constant_result (x->ts.type, x->ts.kind, &x->where); + mpfr_acos (result->value.real, x->value.real, GFC_RND_MODE); break; case BT_COMPLEX: +#ifdef HAVE_mpc_arc + result = gfc_constant_result (x->ts.type, x->ts.kind, &x->where); + mpc_acos (result->value.complex, x->value.complex, GFC_MPC_RND_MODE); + break; +#else return NULL; +#endif default: gfc_internal_error ("in gfc_simplify_acos(): Bad type"); } - result = gfc_constant_result (x->ts.type, x->ts.kind, &x->where); - - mpfr_acos (result->value.real, x->value.real, GFC_RND_MODE); return range_check (result, "ACOS"); } @@ -781,7 +786,13 @@ gfc_simplify_acosh (gfc_expr *x) mpfr_acosh (result->value.real, x->value.real, GFC_RND_MODE); break; case BT_COMPLEX: +#ifdef HAVE_mpc_arc + result = gfc_constant_result (x->ts.type, x->ts.kind, &x->where); + mpc_acosh (result->value.complex, x->value.complex, GFC_MPC_RND_MODE); + break; +#else return NULL; +#endif default: gfc_internal_error ("in gfc_simplify_acosh(): Bad type"); } @@ -1043,7 +1054,13 @@ gfc_simplify_asin (gfc_expr *x) mpfr_asin (result->value.real, x->value.real, GFC_RND_MODE); break; case BT_COMPLEX: +#ifdef HAVE_mpc_arc + result = gfc_constant_result (x->ts.type, x->ts.kind, &x->where); + mpc_asin (result->value.complex, x->value.complex, GFC_MPC_RND_MODE); + break; +#else return NULL; +#endif default: gfc_internal_error ("in gfc_simplify_asin(): Bad type"); } @@ -1067,7 +1084,13 @@ gfc_simplify_asinh (gfc_expr *x) mpfr_asinh (result->value.real, x->value.real, GFC_RND_MODE); break; case BT_COMPLEX: +#ifdef HAVE_mpc_arc + result = gfc_constant_result (x->ts.type, x->ts.kind, &x->where); + mpc_asinh (result->value.complex, x->value.complex, GFC_MPC_RND_MODE); + break; +#else return NULL; +#endif default: gfc_internal_error ("in gfc_simplify_asinh(): Bad type"); } @@ -1091,7 +1114,13 @@ gfc_simplify_atan (gfc_expr *x) mpfr_atan (result->value.real, x->value.real, GFC_RND_MODE); break; case BT_COMPLEX: +#ifdef HAVE_mpc_arc + result = gfc_constant_result (x->ts.type, x->ts.kind, &x->where); + mpc_atan (result->value.complex, x->value.complex, GFC_MPC_RND_MODE); + break; +#else return NULL; +#endif default: gfc_internal_error ("in gfc_simplify_atan(): Bad type"); } @@ -1123,7 +1152,13 @@ gfc_simplify_atanh (gfc_expr *x) mpfr_atanh (result->value.real, x->value.real, GFC_RND_MODE); break; case BT_COMPLEX: +#ifdef HAVE_mpc_arc + result = gfc_constant_result (x->ts.type, x->ts.kind, &x->where); + mpc_atanh (result->value.complex, x->value.complex, GFC_MPC_RND_MODE); + break; +#else return NULL; +#endif default: gfc_internal_error ("in gfc_simplify_atanh(): Bad type"); } -- cgit v1.1