Skip to content

AFNI/NIfTI Server

Sections
Personal tools
You are here: Home » AFNI » Documentation
 

Doxygen Source Code Documentation

Main Page   Alphabetical List   Data Structures   File List   Data Fields   Globals   Search  

cdf_65.c

Go to the documentation of this file.
00001 #include "cdflib.h"
00002 double gsumln(double *a,double *b)
00003 /*
00004 -----------------------------------------------------------------------
00005           EVALUATION OF THE FUNCTION LN(GAMMA(A + B))
00006           FOR 1 .LE. A .LE. 2  AND  1 .LE. B .LE. 2
00007 -----------------------------------------------------------------------
00008 */
00009 {
00010 static double gsumln,x,T1,T2;
00011 /*
00012      ..
00013      .. Executable Statements ..
00014 */
00015     x = *a+*b-2.e0;
00016     if(x > 0.25e0) goto S10;
00017     T1 = 1.0e0+x;
00018     gsumln = gamln1(&T1);
00019     return gsumln;
00020 S10:
00021     if(x > 1.25e0) goto S20;
00022     gsumln = gamln1(&x)+alnrel(&x);
00023     return gsumln;
00024 S20:
00025     T2 = x-1.0e0;
00026     gsumln = gamln1(&T2)+log(x*(1.0e0+x));
00027     return gsumln;
00028 } /* END */
 
Powered by Plone

This site conforms to the following standards: