1 files changed, 41 insertions, 0 deletions
diff --git a/gcc/testsuite/gfortran.dg/besxy.f90 b/gcc/testsuite/gfortran.dg/besxy.f90
new file mode 100644
index 000000000..5cd5c8a96
--- /dev/null
+++ b/gcc/testsuite/gfortran.dg/besxy.f90
@@ -0,0 +1,41 @@
+! { dg-do compile }
+!
+! Check whether BESXY functions take scalars and
+! arrays as arguments (PR31760).
+!
+PROGRAM test_erf
+  REAL :: r = 0.0, ra(2) = (/ 0.0, 1.0 /)
+
+  r  = BESJ0(r)
+  r  = BESJ1(r)
+  r  = BESJN(0, r)
+
+  r  = BESY0(r)
+  r  = BESY1(r)
+  r  = BESYN(0, r)
+
+  ra = BESJ0(ra)
+  ra = BESJ1(ra)
+  ra = BESJN(0, ra)
+
+  ra = BESY0(ra)
+  ra = BESY1(ra)
+  ra = BESYN(0, ra)
+
+  r  = BESSEL_J0(r)
+  r  = BESSEL_J1(r)
+  r  = BESSEL_JN(0, r)
+
+  r  = BESSEL_Y0(r)
+  r  = BESSEL_Y1(r)
+  r  = BESSEL_YN(0, r)
+
+  ra = BESSEL_J0(ra)
+  ra = BESSEL_J1(ra)
+  ra = BESSEL_JN(0, ra)
+
+  ra = BESSEL_Y0(ra)
+  ra = BESSEL_Y1(ra)
+  ra = BESSEL_YN(0, ra)
+
+END PROGRAM