1! { dg-do compile } 2! 3! Check whether BESXY functions take scalars and 4! arrays as arguments (PR31760). 5! 6PROGRAM test_erf 7 REAL :: r = 0.0, ra(2) = (/ 0.0, 1.0 /) 8 9 r = BESJ0(r) 10 r = BESJ1(r) 11 r = BESJN(0, r) 12 13 r = BESY0(r) 14 r = BESY1(r) 15 r = BESYN(0, r) 16 17 ra = BESJ0(ra) 18 ra = BESJ1(ra) 19 ra = BESJN(0, ra) 20 21 ra = BESY0(ra) 22 ra = BESY1(ra) 23 ra = BESYN(0, ra) 24 25 r = BESSEL_J0(r) 26 r = BESSEL_J1(r) 27 r = BESSEL_JN(0, r) 28 29 r = BESSEL_Y0(r) 30 r = BESSEL_Y1(r) 31 r = BESSEL_YN(0, r) 32 33 ra = BESSEL_J0(ra) 34 ra = BESSEL_J1(ra) 35 ra = BESSEL_JN(0, ra) 36 37 ra = BESSEL_Y0(ra) 38 ra = BESSEL_Y1(ra) 39 ra = BESSEL_YN(0, ra) 40 41END PROGRAM 42