summaryrefslogtreecommitdiff
path: root/gcc/testsuite/gfortran.dg/nearest_3.f90
diff options
context:
space:
mode:
Diffstat (limited to 'gcc/testsuite/gfortran.dg/nearest_3.f90')
-rw-r--r--gcc/testsuite/gfortran.dg/nearest_3.f90339
1 files changed, 339 insertions, 0 deletions
diff --git a/gcc/testsuite/gfortran.dg/nearest_3.f90 b/gcc/testsuite/gfortran.dg/nearest_3.f90
new file mode 100644
index 000000000..7d6831670
--- /dev/null
+++ b/gcc/testsuite/gfortran.dg/nearest_3.f90
@@ -0,0 +1,339 @@
+! { dg-do run }
+! { dg-add-options ieee }
+! { dg-skip-if "NaN not supported" { spu-*-* } { "*" } { "" } }
+!
+! PR fortran/34209
+!
+! Test run-time implementation of NEAREST
+!
+program test
+ implicit none
+ real(4), volatile :: r4
+ real(8), volatile :: r8
+
+! Single precision with single-precision sign
+
+ r4 = 0.0_4
+ ! 0+ > 0
+ if (nearest(r4, 1.0) &
+ <= r4) &
+ call abort()
+ ! 0++ > 0+
+ if (nearest(nearest(r4, 1.0), 1.0) &
+ <= nearest(r4, 1.0)) &
+ call abort()
+ ! 0+++ > 0++
+ if (nearest(nearest(nearest(r4, 1.0), 1.0), 1.0) &
+ <= nearest(nearest(r4, 1.0), 1.0)) &
+ call abort()
+ ! 0+- = 0
+ if (nearest(nearest(r4, 1.0), -1.0) &
+ /= r4) &
+ call abort()
+ ! 0++- = 0+
+ if (nearest(nearest(nearest(r4, 1.0), 1.0), -1.0) &
+ /= nearest(r4, 1.0)) &
+ call abort()
+ ! 0++-- = 0
+ if (nearest(nearest(nearest(nearest(r4, 1.0), 1.0), -1.0), -1.0) &
+ /= r4) &
+ call abort()
+
+ ! 0- < 0
+ if (nearest(r4, -1.0) &
+ >= r4) &
+ call abort()
+ ! 0-- < 0+
+ if (nearest(nearest(r4, -1.0), -1.0) &
+ >= nearest(r4, -1.0)) &
+ call abort()
+ ! 0--- < 0--
+ if (nearest(nearest(nearest(r4, -1.0), -1.0), -1.0) &
+ >= nearest(nearest(r4, -1.0), -1.0)) &
+ call abort()
+ ! 0-+ = 0
+ if (nearest(nearest(r4, -1.0), 1.0) &
+ /= r4) &
+ call abort()
+ ! 0--+ = 0-
+ if (nearest(nearest(nearest(r4, -1.0), -1.0), 1.0) &
+ /= nearest(r4, -1.0)) &
+ call abort()
+ ! 0--++ = 0
+ if (nearest(nearest(nearest(nearest(r4, -1.0), -1.0), 1.0), 1.0) &
+ /= r4) &
+ call abort()
+
+ r4 = 42.0_4
+ ! 42++ > 42+
+ if (nearest(nearest(r4, 1.0), 1.0) &
+ <= nearest(r4, 1.0)) &
+ call abort()
+ ! 42-- < 42-
+ if (nearest(nearest(r4, -1.0), -1.0) &
+ >= nearest(r4, -1.0)) &
+ call abort()
+ ! 42-+ = 42
+ if (nearest(nearest(r4, -1.0), 1.0) &
+ /= r4) &
+ call abort()
+ ! 42+- = 42
+ if (nearest(nearest(r4, 1.0), -1.0) &
+ /= r4) &
+ call abort()
+
+ r4 = 0.0
+ ! INF+ = INF
+ if (nearest(1.0/r4, 1.0) /= 1.0/r4) call abort()
+ ! -INF- = -INF
+ if (nearest(-1.0/r4, -1.0) /= -1.0/r4) call abort()
+ ! NAN- = NAN
+ if (.not.isnan(nearest(0.0/r4, 1.0))) call abort()
+ ! NAN+ = NAN
+ if (.not.isnan(nearest(0.0/r4, -1.0))) call abort()
+
+! Double precision with single-precision sign
+
+ r8 = 0.0_8
+ ! 0+ > 0
+ if (nearest(r8, 1.0) &
+ <= r8) &
+ call abort()
+ ! 0++ > 0+
+ if (nearest(nearest(r8, 1.0), 1.0) &
+ <= nearest(r8, 1.0)) &
+ call abort()
+ ! 0+++ > 0++
+ if (nearest(nearest(nearest(r8, 1.0), 1.0), 1.0) &
+ <= nearest(nearest(r8, 1.0), 1.0)) &
+ call abort()
+ ! 0+- = 0
+ if (nearest(nearest(r8, 1.0), -1.0) &
+ /= r8) &
+ call abort()
+ ! 0++- = 0+
+ if (nearest(nearest(nearest(r8, 1.0), 1.0), -1.0) &
+ /= nearest(r8, 1.0)) &
+ call abort()
+ ! 0++-- = 0
+ if (nearest(nearest(nearest(nearest(r8, 1.0), 1.0), -1.0), -1.0) &
+ /= r8) &
+ call abort()
+
+ ! 0- < 0
+ if (nearest(r8, -1.0) &
+ >= r8) &
+ call abort()
+ ! 0-- < 0+
+ if (nearest(nearest(r8, -1.0), -1.0) &
+ >= nearest(r8, -1.0)) &
+ call abort()
+ ! 0--- < 0--
+ if (nearest(nearest(nearest(r8, -1.0), -1.0), -1.0) &
+ >= nearest(nearest(r8, -1.0), -1.0)) &
+ call abort()
+ ! 0-+ = 0
+ if (nearest(nearest(r8, -1.0), 1.0) &
+ /= r8) &
+ call abort()
+ ! 0--+ = 0-
+ if (nearest(nearest(nearest(r8, -1.0), -1.0), 1.0) &
+ /= nearest(r8, -1.0)) &
+ call abort()
+ ! 0--++ = 0
+ if (nearest(nearest(nearest(nearest(r8, -1.0), -1.0), 1.0), 1.0) &
+ /= r8) &
+ call abort()
+
+ r8 = 42.0_8
+ ! 42++ > 42+
+ if (nearest(nearest(r8, 1.0), 1.0) &
+ <= nearest(r8, 1.0)) &
+ call abort()
+ ! 42-- < 42-
+ if (nearest(nearest(r8, -1.0), -1.0) &
+ >= nearest(r8, -1.0)) &
+ call abort()
+ ! 42-+ = 42
+ if (nearest(nearest(r8, -1.0), 1.0) &
+ /= r8) &
+ call abort()
+ ! 42+- = 42
+ if (nearest(nearest(r8, 1.0), -1.0) &
+ /= r8) &
+ call abort()
+
+ r4 = 0.0
+ ! INF+ = INF
+ if (nearest(1.0/r4, 1.0) /= 1.0/r4) call abort()
+ ! -INF- = -INF
+ if (nearest(-1.0/r4, -1.0) /= -1.0/r4) call abort()
+ ! NAN- = NAN
+ if (.not.isnan(nearest(0.0/r4, 1.0))) call abort()
+ ! NAN+ = NAN
+ if (.not.isnan(nearest(0.0/r4, -1.0))) call abort()
+
+
+! Single precision with double-precision sign
+
+ r4 = 0.0_4
+ ! 0+ > 0
+ if (nearest(r4, 1.0d0) &
+ <= r4) &
+ call abort()
+ ! 0++ > 0+
+ if (nearest(nearest(r4, 1.0d0), 1.0d0) &
+ <= nearest(r4, 1.0d0)) &
+ call abort()
+ ! 0+++ > 0++
+ if (nearest(nearest(nearest(r4, 1.0d0), 1.0d0), 1.0d0) &
+ <= nearest(nearest(r4, 1.0d0), 1.0d0)) &
+ call abort()
+ ! 0+- = 0
+ if (nearest(nearest(r4, 1.0d0), -1.0d0) &
+ /= r4) &
+ call abort()
+ ! 0++- = 0+
+ if (nearest(nearest(nearest(r4, 1.0d0), 1.0d0), -1.0d0) &
+ /= nearest(r4, 1.0d0)) &
+ call abort()
+ ! 0++-- = 0
+ if (nearest(nearest(nearest(nearest(r4, 1.0d0), 1.0d0), -1.0d0), -1.0d0) &
+ /= r4) &
+ call abort()
+
+ ! 0- < 0
+ if (nearest(r4, -1.0d0) &
+ >= r4) &
+ call abort()
+ ! 0-- < 0+
+ if (nearest(nearest(r4, -1.0d0), -1.0d0) &
+ >= nearest(r4, -1.0d0)) &
+ call abort()
+ ! 0--- < 0--
+ if (nearest(nearest(nearest(r4, -1.0d0), -1.0d0), -1.0d0) &
+ >= nearest(nearest(r4, -1.0d0), -1.0d0)) &
+ call abort()
+ ! 0-+ = 0
+ if (nearest(nearest(r4, -1.0d0), 1.0d0) &
+ /= r4) &
+ call abort()
+ ! 0--+ = 0-
+ if (nearest(nearest(nearest(r4, -1.0d0), -1.0d0), 1.0d0) &
+ /= nearest(r4, -1.0d0)) &
+ call abort()
+ ! 0--++ = 0
+ if (nearest(nearest(nearest(nearest(r4, -1.0d0), -1.0d0), 1.0d0), 1.0d0) &
+ /= r4) &
+ call abort()
+
+ r4 = 42.0_4
+ ! 42++ > 42+
+ if (nearest(nearest(r4, 1.0d0), 1.0d0) &
+ <= nearest(r4, 1.0d0)) &
+ call abort()
+ ! 42-- < 42-
+ if (nearest(nearest(r4, -1.0d0), -1.0d0) &
+ >= nearest(r4, -1.0d0)) &
+ call abort()
+ ! 42-+ = 42
+ if (nearest(nearest(r4, -1.0d0), 1.0d0) &
+ /= r4) &
+ call abort()
+ ! 42+- = 42
+ if (nearest(nearest(r4, 1.0d0), -1.0d0) &
+ /= r4) &
+ call abort()
+
+ r4 = 0.0
+ ! INF+ = INF
+ if (nearest(1.0d0/r4, 1.0d0) /= 1.0d0/r4) call abort()
+ ! -INF- = -INF
+ if (nearest(-1.0d0/r4, -1.0d0) /= -1.0d0/r4) call abort()
+ ! NAN- = NAN
+ if (.not.isnan(nearest(0.0/r4, 1.0d0))) call abort()
+ ! NAN+ = NAN
+ if (.not.isnan(nearest(0.0/r4, -1.0d0))) call abort()
+
+! Double precision with double-precision sign
+
+ r8 = 0.0_8
+ ! 0+ > 0
+ if (nearest(r8, 1.0d0) &
+ <= r8) &
+ call abort()
+ ! 0++ > 0+
+ if (nearest(nearest(r8, 1.0d0), 1.0d0) &
+ <= nearest(r8, 1.0d0)) &
+ call abort()
+ ! 0+++ > 0++
+ if (nearest(nearest(nearest(r8, 1.0d0), 1.0d0), 1.0d0) &
+ <= nearest(nearest(r8, 1.0d0), 1.0d0)) &
+ call abort()
+ ! 0+- = 0
+ if (nearest(nearest(r8, 1.0d0), -1.0d0) &
+ /= r8) &
+ call abort()
+ ! 0++- = 0+
+ if (nearest(nearest(nearest(r8, 1.0d0), 1.0d0), -1.0d0) &
+ /= nearest(r8, 1.0d0)) &
+ call abort()
+ ! 0++-- = 0
+ if (nearest(nearest(nearest(nearest(r8, 1.0d0), 1.0d0), -1.0d0), -1.0d0) &
+ /= r8) &
+ call abort()
+
+ ! 0- < 0
+ if (nearest(r8, -1.0d0) &
+ >= r8) &
+ call abort()
+ ! 0-- < 0+
+ if (nearest(nearest(r8, -1.0d0), -1.0d0) &
+ >= nearest(r8, -1.0d0)) &
+ call abort()
+ ! 0--- < 0--
+ if (nearest(nearest(nearest(r8, -1.0d0), -1.0d0), -1.0d0) &
+ >= nearest(nearest(r8, -1.0d0), -1.0d0)) &
+ call abort()
+ ! 0-+ = 0
+ if (nearest(nearest(r8, -1.0d0), 1.0d0) &
+ /= r8) &
+ call abort()
+ ! 0--+ = 0-
+ if (nearest(nearest(nearest(r8, -1.0d0), -1.0d0), 1.0d0) &
+ /= nearest(r8, -1.0d0)) &
+ call abort()
+ ! 0--++ = 0
+ if (nearest(nearest(nearest(nearest(r8, -1.0d0), -1.0d0), 1.0d0), 1.0d0) &
+ /= r8) &
+ call abort()
+
+ r8 = 42.0_8
+ ! 42++ > 42+
+ if (nearest(nearest(r8, 1.0d0), 1.0d0) &
+ <= nearest(r8, 1.0d0)) &
+ call abort()
+ ! 42-- < 42-
+ if (nearest(nearest(r8, -1.0d0), -1.0d0) &
+ >= nearest(r8, -1.0d0)) &
+ call abort()
+ ! 42-+ = 42
+ if (nearest(nearest(r8, -1.0d0), 1.0d0) &
+ /= r8) &
+ call abort()
+ ! 42+- = 42
+ if (nearest(nearest(r8, 1.0d0), -1.0d0) &
+ /= r8) &
+ call abort()
+
+ r4 = 0.0
+ ! INF+ = INF
+ if (nearest(1.0d0/r4, 1.0d0) /= 1.0d0/r4) call abort()
+ ! -INF- = -INF
+ if (nearest(-1.0d0/r4, -1.0d0) /= -1.0d0/r4) call abort()
+ ! NAN- = NAN
+ if (.not.isnan(nearest(0.0/r4, 1.0d0))) call abort()
+ ! NAN+ = NAN
+ if (.not.isnan(nearest(0.0/r4, -1.0d0))) call abort()
+
+end program test
generated by cgit v1.2.3 (git 2.25.1) at 2025-04-26 02:00:25 +0000