extended QGAMMA( extended arg ) { //============================ if( fabs( arg ) <= xmin ) { if( arg <= 0.0 ) { return( -xinf ); } else { return( xinf ); } } else if( fabs( arg ) >= big ) { return( xinf ); } return( __gamma( arg, &xinf ) ); }
single GAMMA( single arg ) { //========================== if( fabs( arg ) <= xmin ) { if( arg <= 0.0 ) { return( -xinf ); } else { return( xinf ); } } else if( fabs( arg ) >= big ) { return( xinf ); } return( __gamma( arg, &xinf ) ); }