1 |
/* specfunc/gamma.c |
/* cdf/gamma.c |
2 |
* |
* |
3 |
* Copyright (C) 1996, 1997, 1998, 1999, 2000 Gerard Jungman |
* Copyright (C) 1996, 1997, 1998, 1999, 2000 Gerard Jungman |
4 |
* |
* |
14 |
* |
* |
15 |
* You should have received a copy of the GNU General Public License |
* You should have received a copy of the GNU General Public License |
16 |
* along with this program; if not, write to the Free Software |
* along with this program; if not, write to the Free Software |
17 |
* Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. |
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307, USA. |
18 |
*/ |
*/ |
19 |
|
|
20 |
/* Author: G. Jungman */ |
/* Author: G. Jungman */ |
1051 |
result->err += 2.0 * GSL_DBL_EPSILON * result->val; |
result->err += 2.0 * GSL_DBL_EPSILON * result->val; |
1052 |
return GSL_SUCCESS; |
return GSL_SUCCESS; |
1053 |
} |
} |
1054 |
else if(x < GSL_SF_GAMMA_XMAX) { |
else if(x < GSL_CDF_GAMMA_XMAX) { |
1055 |
/* We do not want to exponentiate the logarithm |
/* We do not want to exponentiate the logarithm |
1056 |
* if x is large because of the inevitable |
* if x is large because of the inevitable |
1057 |
* inflation of the error. So we carefully |
* inflation of the error. So we carefully |