File: nearest_6.f90

package info (click to toggle)
gcc-arm-none-eabi 15%3A14.2.rel1-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 1,099,328 kB
  • sloc: cpp: 3,627,108; ansic: 2,571,498; ada: 834,230; f90: 235,082; makefile: 79,231; asm: 74,984; xml: 51,692; exp: 39,736; sh: 33,298; objc: 15,629; python: 15,069; fortran: 14,429; pascal: 7,003; awk: 5,070; perl: 3,106; ml: 285; lisp: 253; lex: 204; haskell: 135
file content (26 lines) | stat: -rw-r--r-- 820 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 
! { dg-do run }
! PR fortran/109186 - Verify that NEAREST produces same results at
!                     compile-time and run-time for corner cases
! Reported by John Harper

program p
  implicit none
  integer,  parameter :: sp = selected_real_kind (6)
  integer,  parameter :: dp = selected_real_kind (13)
  real(sp), parameter :: x1 = huge (1._sp), t1 = tiny (1._sp)
  real(dp), parameter :: x2 = huge (1._dp), t2 = tiny (1._dp)
  real(sp), volatile  :: y1, z1
  real(dp), volatile  :: y2, z2
  y1 = x1
  z1 = nearest (y1, -1._sp)
  if (nearest (x1,  -1._sp) /= z1) stop 1
  y2 = x2
  z2 = nearest (y2, -1._dp)
  if (nearest (x2,  -1._dp) /= z2) stop 2
  y1 = t1
  z1 = nearest (y1,  1._sp)
  if (nearest (t1,   1._sp) /= z1) stop 3
  y2 = t2
  z2 = nearest (y2,  1._dp)
  if (nearest (t2,   1._dp) /= z2) stop 4
end