gcc/testsuite/g77.f-torture/execute/f90-intrinsic-mathematical.f

   1 c  f90-intrinsic-mathematical.f
   2 c
   3 c Test Fortran 90 intrinsic mathematical functions - Section 13.10.3 and
   4 c 13.13
   5 c     David Billinghurst <David.Billinghurst@riotinto.com>
   6 c
   7 c Notes:
   8 c  * g77 does not fully comply with F90.  Noncompliances noted in comments.
   9 c  * Section 13.12: Specific names for intrinsic functions tested in
  10 c intrinsic77.f
  11
  12       logical fail
  13       common /flags/ fail
  14       fail = .false.
  15
  16 c     ACOS - Section 13.13.3
  17       call c_r(ACOS(0.54030231),1.0,'ACOS(real)')
  18       call c_d(ACOS(0.54030231d0),1.d0,'ACOS(double)')
  19
  20 c     ASIN - Section 13.13.12
  21       call c_r(ASIN(0.84147098),1.0,'ASIN(real)')
  22       call c_d(ASIN(0.84147098d0),1.d0,'ASIN(double)')
  23
  24 c     ATAN - Section 13.13.14
  25       call c_r(ATAN(1.5574077),1.0,'ATAN(real)')
  26       call c_d(ATAN(1.5574077d0),1.d0,'ATAN(double)')
  27
  28 c     ATAN2 - Section 13.13.15
  29       call c_r(ATAN2(1.5574077,1.),1.0,'ATAN2(real)')
  30       call c_d(ATAN2(1.5574077d0,1.d0),1.d0,'ATAN2(double)')
  31
  32 c     COS - Section 13.13.22
  33       call c_r(COS(1.0),0.54030231,'COS(real)')
  34       call c_d(COS(1.d0),0.54030231d0,'COS(double)')
  35       call c_c(COS((1.,0.)),(0.54030231,0.),'COS(complex)')
  36       call c_z(COS((1.d0,0.d0)),(0.54030231d0,0.d0),
  37      $     'COS(double complex)')
  38
  39 c     COSH - Section 13.13.23
  40       call c_r(COSH(1.0),1.5430806,'COSH(real)')
  41       call c_d(COSH(1.d0),1.5430806d0,'COSH(double)')
  42
  43 c     EXP - Section 13.13.34
  44       call c_r(EXP(1.0),2.7182818,'EXP(real)')
  45       call c_d(EXP(1.d0),2.7182818d0,'EXP(double)')
  46       call c_c(EXP((1.,0.)),(2.7182818,0.),'EXP(complex)')
  47       call c_z(EXP((1.d0,0.d0)),(2.7182818d0,0.d0),
  48      $     'EXP(double complex)')
  49
  50 c     LOG - Section 13.13.59
  51       call c_r(LOG(10.0),2.3025851,'LOG(real)')
  52       call c_d(LOG(10.d0),2.3025851d0,'LOG(double)')
  53       call c_c(LOG((10.,0.)),(2.3025851,0.),'LOG(complex)')
  54       call c_z(LOG((10.d0,0.)),(2.3025851d0,0.d0),
  55      $     'LOG(double complex)')
  56
  57 c     LOG10 - Section 13.13.60
  58       call c_r(LOG10(10.0),1.0,'LOG10(real)')
  59       call c_d(LOG10(10.d0),1.d0,'LOG10(double)')
  60
  61 c     SIN - Section 13.13.97
  62       call c_r(SIN(1.0),0.84147098,'SIN(real)')
  63       call c_d(SIN(1.d0),0.84147098d0,'SIN(double)')
  64       call c_c(SIN((1.,0.)),(0.84147098,0.),'SIN(complex)')
  65       call c_z(SIN((1.d0,0.d0)),(0.84147098d0,0.d0),
  66      $     'SIN(double complex)')
  67
  68 c     SINH - Section 13.13.98
  69       call c_r(SINH(1.0),1.175201,'SINH(real)')
  70       call c_d(SINH(1.d0),1.175201d0,'SINH(double)')
  71
  72 c     SQRT - Section 13.13.102
  73       call c_r(SQRT(4.0),2.0,'SQRT(real)')
  74       call c_d(SQRT(4.d0),2.d0,'SQRT(double)')
  75       call c_c(SQRT((4.,0.)),(2.,0.),'SQRT(complex)')
  76       call c_z(SQRT((4.d0,0.)),(2.d0,0.),
  77      $     'SQRT(double complex)')
  78
  79 c     TAN - Section 13.13.105
  80       call c_r(TAN(1.0),1.5574077,'TAN(real)')
  81       call c_d(TAN(1.d0),1.5574077d0,'TAN(double)')
  82
  83 c     TANH - Section 13.13.106
  84       call c_r(TANH(1.0),0.76159416,'TANH(real)')
  85       call c_d(TANH(1.d0),0.76159416d0,'TANH(double)')
  86
  87       if ( fail ) call abort()
  88       end
  89
  90       subroutine failure(label)
  91 c     Report failure and set flag
  92       character*(*) label
  93       logical fail
  94       common /flags/ fail
  95       write(6,'(a,a,a)') 'Test ',label,' FAILED'
  96       fail = .true.
  97       end
  98
  99       subroutine c_r(a,b,label)
 100 c     Check if REAL a equals b, and fail otherwise
 101       real a, b
 102       character*(*) label
 103       if ( abs(a-b) .gt. 1.0e-5 ) then
 104          call failure(label)
 105          write(6,*) 'Got ',a,' expected ', b
 106       end if
 107       end
 108
 109       subroutine c_d(a,b,label)
 110 c     Check if DOUBLE PRECISION a equals b, and fail otherwise
 111       double precision a, b
 112       character*(*) label
 113       if ( abs(a-b) .gt. 1.0d-5 ) then
 114          call failure(label)
 115          write(6,*) 'Got ',a,' expected ', b
 116       end if
 117       end
 118
 119       subroutine c_c(a,b,label)
 120 c     Check if COMPLEX a equals b, and fail otherwise
 121       complex a, b
 122       character*(*) label
 123       if ( abs(a-b) .gt. 1.0e-5 ) then
 124          call failure(label)
 125          write(6,*) 'Got ',a,' expected ', b
 126       end if
 127       end
 128
 129       subroutine c_z(a,b,label)
 130 c     Check if COMPLEX a equals b, and fail otherwise
 131       double complex a, b
 132       character*(*) label
 133       if ( abs(a-b) .gt. 1.0d-5 ) then
 134          call failure(label)
 135          write(6,*) 'Got ',a,' expected ', b
 136       end if
 137       end