hiro/arb
1
0
Fork 0
mirror of https://github.com/vale981/arb synced 2025-03-04 17:01:40 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

add simple wrappers for erfc, erfi, expint

Browse source
This commit is contained in:
Fredrik Johansson 2015-01-21 17:28:16 +01:00
parent de03803f57
commit 3eb1cf06cd
4 changed files with 137 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
acb_hypgeom.h
View file

@ -77,6 +77,12 @@ void acb_hypgeom_gamma_upper_singular(acb_t res, long s, const acb_t z, int modi
void acb_hypgeom_gamma_upper(acb_t res, const acb_t s, const acb_t z, int modified, long prec);
void acb_hypgeom_expint(acb_t res, const acb_t s, const acb_t z, long prec);
void acb_hypgeom_erfc(acb_t res, const acb_t z, long prec);
void acb_hypgeom_erfi(acb_t res, const acb_t z, long prec);
#ifdef __cplusplus
}
#endif

57
acb_hypgeom/erfc.c Normal file
View file

@ -0,0 +1,57 @@
/*=============================================================================
This file is part of ARB.
ARB is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
ARB 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 General Public License for more details.
You should have received a copy of the GNU General Public License
along with ARB; if not, write to the Free Software
Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
=============================================================================*/
/******************************************************************************
Copyright (C) 2015 Fredrik Johansson
******************************************************************************/
#include "acb_hypgeom.h"
void
acb_hypgeom_erfc(acb_t res, const acb_t z, long prec)
{
if (arb_is_positive(acb_realref(z)))
{
acb_t t, u;
acb_init(t);
acb_init(u);
acb_set_si(t, 1);
acb_mul_2exp_si(t, t, -1);
acb_mul(u, z, z, prec);
acb_hypgeom_gamma_upper(res, t, u, 0, prec);
arb_const_sqrt_pi(acb_realref(t), prec);
acb_div_arb(res, res, acb_realref(t), prec);
acb_clear(t);
acb_clear(u);
}
else
{
acb_hypgeom_erf(res, z, prec);
acb_sub_ui(res, res, 1, prec);
acb_neg(res, res);
}
}

36
acb_hypgeom/erfi.c Normal file
View file

@ -0,0 +1,36 @@
/*=============================================================================
This file is part of ARB.
ARB is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
ARB 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 General Public License for more details.
You should have received a copy of the GNU General Public License
along with ARB; if not, write to the Free Software
Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
=============================================================================*/
/******************************************************************************
Copyright (C) 2015 Fredrik Johansson
******************************************************************************/
#include "acb_hypgeom.h"
void
acb_hypgeom_erfi(acb_t res, const acb_t z, long prec)
{
acb_mul_onei(res, z);
acb_hypgeom_erf(res, res, prec);
acb_mul_onei(res, res);
acb_neg(res, res);
}

38
acb_hypgeom/expint.c Normal file
View file

@ -0,0 +1,38 @@
/*=============================================================================
This file is part of ARB.
ARB is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
ARB 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 General Public License for more details.
You should have received a copy of the GNU General Public License
along with ARB; if not, write to the Free Software
Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
=============================================================================*/
/******************************************************************************
Copyright (C) 2015 Fredrik Johansson
******************************************************************************/
#include "acb_hypgeom.h"
void
acb_hypgeom_expint(acb_t res, const acb_t s, const acb_t z, long prec)
{
acb_t t;
acb_init(t);
acb_sub_ui(t, s, 1, prec);
acb_neg(t, t);
acb_hypgeom_gamma_upper(res, t, z, 1, prec);
acb_clear(t);
}