erfc.f   [plain text]

c { dg-do run }
c============================================== test.f
      real x, y
      real(kind=8) x1, y1
      x=0.
      y = erfc(x)
      if (y .ne. 1.) call abort

      x=1.1
      y = erfc(x)
      if (abs(y - .1197949) .ge. 1.e-6) call abort

c modified from x=10, y .gt. 1.5e-44 to avoid lack of -mieee on Alphas.
      x=8
      y = erfc(x)
      if (y .gt. 1.2e-28) call abort

      x1=0.
      y1 = erfc(x1)
      if (y1 .ne. 1.) call abort

      x1=1.1d0
      y1 = erfc(x1)
      if (abs(y1 - .1197949d0) .ge. 1.d-6) call abort

      x1=10
      y1 = erfc(x1)
      if (y1 .gt. 1.5d-44) call abort
      end
c=================================================
!output:
!  0.  1.875
!  1.10000002  1.48958981
!  10.  5.00220949E-06
!
!The values should be:
!erfc(0)=1
!erfc(1.1)= 0.1197949
!erfc(10)<1.543115467311259E-044