Viewing file:  thyperbolic.c (7.14 KB) -rw-r--r--
Select action/file-type:
 (+) |  (+) |  (+) | Code (+) | Session (+) |  (+) | SDB (+) |  (+) |  (+) |  (+) |  (+) |  (+) |
/* Test file for hyperbolic function : mpfr_cosh, mpfr_sinh, mpfr_tanh, mpfr_acosh, mpfr_asinh, mpfr_atanh.
Copyright 2001-2004, 2006-2020 Free Software Foundation, Inc.
Contributed by the AriC and Caramba projects, INRIA.
This file is part of the GNU MPFR Library.
The GNU MPFR Library is free software; you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as published by
the Free Software Foundation; either version 3 of the License, or (at your
option) any later version.
The GNU MPFR Library is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
License for more details.
You should have received a copy of the GNU Lesser General Public License
along with the GNU MPFR Library; see the file COPYING.LESSER. If not, see
https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc.,
51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA. */
#include "mpfr-test.h"
static int
check_NAN (void)
{
mpfr_t t, ch,sh,th,ach,ash,ath;
int tester;
int fail = 0;
mpfr_init2(t,200);
mpfr_init2(ch,200);
mpfr_init2(sh,200);
mpfr_init2(th,200);
mpfr_init2(ach,200);
mpfr_init2(ash,200);
mpfr_init2(ath,200);
MPFR_SET_NAN(t);
/******cosh********/
tester=mpfr_cosh(ch,t,MPFR_RNDD);
if (!MPFR_IS_NAN(ch) || tester!=0)
{
printf("cosh NAN \n");
fail = 1;
goto clean_up;
}
/******sinh********/
tester=mpfr_sinh(sh,t,MPFR_RNDD);
if (!MPFR_IS_NAN(sh) || tester!=0)
{
printf("sinh NAN \n");
fail = 1;
goto clean_up;
}
/******tanh********/
tester=mpfr_tanh(th,t,MPFR_RNDD);
if (!MPFR_IS_NAN(th) || tester!=0)
{
printf("tanh NAN \n");
fail = 1;
goto clean_up;
}
/******acosh********/
tester=mpfr_acosh(ach,t,MPFR_RNDD);
if (!MPFR_IS_NAN(ach) || tester!=0)
{
printf("acosh NAN \n");
fail = 1;
goto clean_up;
}
/******asinh********/
tester=mpfr_asinh(ash,t,MPFR_RNDD);
if (!MPFR_IS_NAN(ash) || tester!=0)
{
printf("asinh NAN \n");
fail = 1;
goto clean_up;
}
/******atanh********/
tester=mpfr_atanh(ath,t,MPFR_RNDD);
if (!MPFR_IS_NAN(ath) || tester!=0)
{
printf("atanh NAN \n");
fail = 1;
goto clean_up;
}
clean_up:
mpfr_clear(t);
mpfr_clear(ch);
mpfr_clear(sh);
mpfr_clear(th);
mpfr_clear(ach);
mpfr_clear(ash);
mpfr_clear(ath);
return fail;
}
static int
check_zero (void)
{
mpfr_t t, ch,sh,th,ach,ash,ath;
int tester;
int fail = 0;
mpfr_init2(t,200);
mpfr_init2(ch,200);
mpfr_init2(sh,200);
mpfr_init2(th,200);
mpfr_init2(ach,200);
mpfr_init2(ash,200);
mpfr_init2(ath,200);
mpfr_set_ui(t,0,MPFR_RNDD);
/******cosh********/
tester = mpfr_cosh (ch, t, MPFR_RNDD);
if (mpfr_cmp_ui(ch, 1) || tester)
{
printf("cosh(0) \n");
fail = 1;
goto clean_up;
}
/******sinh********/
tester = mpfr_sinh (sh, t, MPFR_RNDD);
if (!MPFR_IS_ZERO(sh) || tester)
{
printf("sinh(0) \n");
fail = 1;
goto clean_up;
}
/******tanh********/
tester = mpfr_tanh (th, t, MPFR_RNDD);
if (!MPFR_IS_ZERO(th) || tester)
{
printf("tanh(0) \n");
fail = 1;
goto clean_up;
}
/******acosh********/
tester=mpfr_acosh(ach,t,MPFR_RNDD);
if (!MPFR_IS_NAN(ach) || tester)
{
printf("acosh(0) \n");
fail = 1;
goto clean_up;
}
/******asinh********/
tester=mpfr_asinh(ash,t,MPFR_RNDD);
if (!MPFR_IS_ZERO(ash) || tester)
{
printf("asinh(0) \n");
fail = 1;
goto clean_up;
}
/******atanh********/
tester=mpfr_atanh(ath,t,MPFR_RNDD);
if (!MPFR_IS_ZERO(ath) || tester)
{
printf("atanh(0) \n");
fail = 1;
goto clean_up;
}
clean_up:
mpfr_clear(t);
mpfr_clear(ch);
mpfr_clear(sh);
mpfr_clear(th);
mpfr_clear(ach);
mpfr_clear(ash);
mpfr_clear(ath);
return fail;
}
static int
check_INF (void)
{
mpfr_t t, ch, sh, th, ach, ash, ath;
int tester;
int fail = 0;
mpfr_init2 (t, 200);
mpfr_init2 (ch, 200);
mpfr_init2 (sh, 200);
mpfr_init2 (th, 200);
mpfr_init2 (ach, 200);
mpfr_init2 (ash, 200);
mpfr_init2 (ath, 200);
MPFR_SET_INF(t);
if(MPFR_IS_NEG (t))
MPFR_CHANGE_SIGN(t);
/******cosh********/
tester = mpfr_cosh(ch,t,MPFR_RNDD);
if (!MPFR_IS_INF(ch) || MPFR_IS_NEG (ch) || tester!=0)
{
printf("cosh(INF) \n");
fail = 1;
goto clean_up;
}
/******sinh********/
tester=mpfr_sinh(sh,t,MPFR_RNDD);
if (!MPFR_IS_INF(sh) || MPFR_IS_NEG (sh) || tester!=0)
{
printf("sinh(INF) \n");
fail = 1;
goto clean_up;
}
/******tanh********/
tester=mpfr_tanh(th,t,MPFR_RNDD);
if (mpfr_cmp_ui(th,1) != 0 || tester!=0)
{
printf("tanh(INF) \n");
fail = 1;
goto clean_up;
}
/******acosh********/
tester=mpfr_acosh(ach,t,MPFR_RNDD);
if (!MPFR_IS_INF(ach) || MPFR_IS_NEG (ach) || tester!=0)
{
printf("acosh(INF) \n");
fail = 1;
goto clean_up;
}
/******asinh********/
tester=mpfr_asinh(ash,t,MPFR_RNDD);
if (!MPFR_IS_INF(ash) || MPFR_IS_NEG (ash) || tester!=0)
{
printf("asinh(INF) \n");
fail = 1;
goto clean_up;
}
/******atanh********/
tester = mpfr_atanh (ath, t, MPFR_RNDD);
if (!MPFR_IS_NAN(ath) || tester != 0)
{
printf("atanh(INF) \n");
fail = 1;
goto clean_up;
}
MPFR_CHANGE_SIGN(t);
/******cosh********/
tester=mpfr_cosh(ch,t,MPFR_RNDD);
if (!MPFR_IS_INF(ch) || MPFR_IS_NEG (ch) || tester!=0)
{
printf("cosh(-INF) \n");
fail = 1;
goto clean_up;
}
/******sinh********/
tester=mpfr_sinh(sh,t,MPFR_RNDD);
if (!MPFR_IS_INF(sh) || MPFR_IS_POS (sh) || tester!=0)
{
printf("sinh(-INF) \n");
fail = 1;
goto clean_up;
}
/******tanh********/
tester=mpfr_tanh(th,t,MPFR_RNDD);
if (!mpfr_cmp_ui(th,-1) || tester!=0)
{
printf("tanh(-INF) \n");
fail = 1;
goto clean_up;
}
/******acosh********/
tester=mpfr_acosh(ach,t,MPFR_RNDD);
if (!MPFR_IS_NAN(ach) || tester!=0)
{
printf("acosh(-INF) \n");
fail = 1;
goto clean_up;
}
/******asinh********/
tester=mpfr_asinh(ash,t,MPFR_RNDD);
if (!MPFR_IS_INF(ash) || MPFR_IS_POS (ash) || tester!=0)
{
printf("asinh(-INF) \n");
fail = 1;
goto clean_up;
}
/******atanh********/
tester = mpfr_atanh (ath, t, MPFR_RNDD);
if (!MPFR_IS_NAN(ath) || tester != 0)
{
printf("atanh(-INF) \n");
fail = 1;
goto clean_up;
}
clean_up:
mpfr_clear(t);
mpfr_clear(ch);
mpfr_clear(sh);
mpfr_clear(th);
mpfr_clear(ach);
mpfr_clear(ash);
mpfr_clear(ath);
return fail;
}
int
main(void)
{
tests_start_mpfr ();
if (check_zero ())
{
printf ("Error in evaluation at 0\n");
exit (1);
}
if (check_INF ())
{
printf ("Error in evaluation of INF\n");
exit (1);
}
if (check_NAN ())
{
printf ("Error in evaluation of NAN\n");
exit (1);
}
tests_end_mpfr ();
return 0;
}
|
