← Back to branch summary

~ubuntu-branches/ubuntu/precise/eglibc/precise

« back to all changes in this revision

Viewing changes to sysdeps/ieee754/ldbl-128/e_lgammal_r.c

  • Committer: Package Import Robot
  • Author(s): Matthias Klose
  • Date: 2011-10-04 17:48:26 UTC
  • mfrom: (216.1.23 oneiric)
  • Revision ID: package-import@ubuntu.com-20111004174826-2cyb9ewn3ucymlsx
Tags: 2.13-20ubuntu5
https://launchpad.net/bugs/853688
libc6-dev: Don't break the current {gnat,gcj}-4.4-base versons. LP: #853688.

Show diffs side-by-side

added added

removed removed

Lines of Context:
sysdeps/ieee754/ldbl-128/e_lgammal_r.c
754
754
}
755
755
 
756
756
 
 
757
#ifdef __STDC__
757
758
long double
758
759
__ieee754_lgammal_r (long double x, int *signgamp)
 
760
#else
 
761
long double
 
762
__ieee754_lgammal_r (x, signgamp)
 
763
     long double x;
 
764
     int *signgamp;
 
765
#endif
759
766
{
760
767
  long double p, q, w, z, nx;
761
768
  int i, nn;
768
775
  if (x == 0.0L)
769
776
    {
770
777
      if (__signbitl (x))
771
 
        *signgamp = -1;
 
778
        *signgamp = -1;
772
779
    }
773
780
 
774
781
  if (x < 0.0L)
856
863
                {
857
864
                  z = x - 0.75L;
858
865
                  p = z * neval (z, RN1r75, NRN1r75)
859
 
                        / deval (z, RD1r75, NRD1r75);
 
866
                        / deval (z, RD1r75, NRD1r75);
860
867
                  p += lgam1r75b;
861
868
                  p += lgam1r75a;
862
869
                }
1033
1040
  q += neval (p, RASY, NRASY) / x;
1034
1041
  return (q);
1035
1042
}
1036
 
strong_alias (__ieee754_lgammal_r, __lgammal_r_finite)