/[gsl]/gsl/cdf/gamma_func.c
ViewVC logotype

Diff of /gsl/cdf/gamma_func.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1.2.1 by jstover, Tue Dec 31 16:16:15 2002 UTC revision 1.1.2.2 by jstover, Thu Jan 2 23:19:18 2003 UTC
# Line 1  Line 1 
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   *   *
# Line 14  Line 14 
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 */
# Line 1051  gamma_xgthalf(const double x, gsl_cdf_re Line 1051  gamma_xgthalf(const double x, gsl_cdf_re
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

Legend:
Removed from v.1.1.2.1  
changed lines
  Added in v.1.1.2.2

savannah-hackers-public@gnu.org
ViewVC Help
Powered by ViewVC 1.1.26