Re: [Bug-gsl] gsl_sf_hyperg_1F1

[Weibin Li]

I compile with
gcc bug_hyper.c -lgsl -lm -o bugh

In the output, the error seems larger than the main value, e.g., the
first tens of output
*****************************************************************************************************
38.8687 -6.651329595131e+18     3.776536780306e+19
38.8687 -6.440866203123e+18     3.776536780987e+19
38.8687 -5.845910332028e+18     3.776536781668e+19
38.8687 -7.409159267905e+18     3.776536782349e+19
38.8687 -6.87604313712e+18      3.776536783029e+19
38.8687 -7.337048036355e+18     3.77653678371e+19
38.8687 -7.790012654018e+18     3.776536784391e+19
38.8687 -6.082386334844e+18     3.776536785071e+19
38.8687 -5.935589739723e+18     3.776536785752e+19
38.8687 -7.11026720186e+18      3.776536786433e+19
38.8687 -7.645859192076e+18     3.776536787113e+19
38.8687 -8.182429272082e+18     3.776536787794e+19
38.8687 -4.881283642078e+18     3.776536788475e+19
38.8687 -4.776660818405e+18     3.776536789156e+19
38.8687 -5.10451703914e+18      3.776536789836e+19
38.8687 -6.999912261672e+18     3.776536790517e+19
38.8687 -9.998028827588e+18     3.776536791198e+19
38.8687 -7.923642832751e+18     3.776536791878e+19
38.8687 -6.13227395409e+18      3.776536792559e+19
38.8687 -7.342701081958e+18     3.77653679324e+19
38.8687 -8.774127765439e+18     3.77653679392e+19
38.8687 -7.097891724944e+18     3.776536794601e+19
38.8687 -6.575969198714e+18     3.776536795282e+19
38.8687 -7.13604838894e+18      3.776536795962e+19
38.8687 -8.460246152981e+18     3.776536796643e+19
38.8687 -5.136538170138e+18     3.776536797324e+19
38.8687 -8.009477006139e+18     3.776536798005e+19
38.8687 -4.825528765175e+18     3.776536798685e+19
38.8687 -6.373481350919e+18     3.776536799366e+19
38.8687 -6.434091393798e+18     3.776536800047e+19
38.8687 -8.213791928209e+18     3.776536800727e+19
38.8687 -8.598940254748e+18     3.776536801408e+19
38.8687 -5.095280478498e+18     3.776536802089e+19
38.8687 -8.493031884986e+18     3.776536802769e+19
38.8687 -8.312410653034e+18     3.77653680345e+19
38.8687 -7.229368958853e+18     3.776536804131e+19
38.8687 -5.767473774847e+18     3.776536804812e+19
38.8687 -3.905511944514e+18     3.776536805492e+19
38.8687 -5.799374618919e+18     3.776536806173e+19
38.8687 -5.860381873821e+18     3.776536806854e+19
38.8687 -7.296747576845e+18     3.776536807534e+19
38.8687 -7.247582347993e+18     3.776536808215e+19
38.8687 -7.422721893279e+18     3.776536808896e+19
38.8687 -8.105539444036e+18     3.776536809576e+19
38.8687 -6.076738251137e+18     3.776536810257e+19
38.8687 -6.864860945364e+18     3.776536810938e+19
38.8687 -7.453084201213e+18     3.776536811619e+19
38.8687 -7.438670308571e+18     3.776536812299e+19
38.8687 -8.65478307854e+18      3.77653681298e+19
38.8687 -6.520998476832e+18     3.776536813661e+19
38.8687 -5.883689107301e+18     3.776536814341e+19
38.8687 -7.451158908856e+18     3.776536815022e+19
38.8687 -5.989263565292e+18     3.776536815703e+19
38.8687 -8.563064953954e+18     3.776536816383e+19
38.8687 -7.017148186973e+18     3.776536817064e+19
38.8687 -6.893653490311e+18     3.776536817745e+19
38.8687 -5.623495276948e+18     3.776536818426e+19
38.8687 -7.370414825849e+18     3.776536819106e+19
38.8687 -8.037443520789e+18     3.776536819787e+19
38.8687 -6.733180438873e+18     3.776536820468e+19
38.8687 -5.716431795564e+18     3.776536821148e+19
38.8687 -4.740281441476e+18     3.776536821829e+19
38.8687 -8.827657073298e+18     3.77653682251e+19
38.8687 -7.426134675734e+18     3.77653682319e+19
38.8687 -6.973240151763e+18     3.776536823871e+19
38.8687 -5.298255360332e+18     3.776536824552e+19
38.8687 -5.802182514645e+18     3.776536825232e+19
38.8687 -6.105419408952e+18     3.776536825913e+19
38.8687 -7.246035882174e+18     3.776536826594e+19
38.8687 -9.301167993379e+18     3.776536827275e+19
38.8687 -7.185591375806e+18     3.776536827955e+19
38.8687 -7.576075779048e+18     3.776536828636e+19
38.8687 -8.48791412194e+18      3.776536829317e+19
38.8687 -8.355980811712e+18     3.776536829997e+19
38.8687 -4.679452434633e+18     3.776536830678e+19
38.8687 -7.281209291425e+18     3.776536831359e+19
38.8687 -7.846481326287e+18     3.776536832039e+19
38.8687 -7.560879476736e+18     3.77653683272e+19
38.8687 -6.329833255143e+18     3.776536833401e+19
38.8687 -7.003356313895e+18     3.776536834082e+19
38.8687 -4.950067740148e+18     3.776536834762e+19
*****************************************************************************
> At Mon, 30 Nov 2009 15:26:11 +0100,
> Weibin Li wrote:
>   
>> I experienced bugs with gsl_sf_hyperg_1F1. The version  is gsl_1.12, but
>> the testing is also done with gsl_1.13, using a  mac with version 10.5.8
>> and hp-workstation with gnome_2.24.1.
>>     
>
> Thanks for the email.  Can you sent the exact output you get from the
> program as well, just so we can be sure we're seeing the same problem.
>
>