42 extern void bessel_jn_r16 (gfc_array_r16 * const restrict ret, int n1,
47 bessel_jn_r16 (gfc_array_r16 * const restrict ret, int n1, int n2, GFC_REAL_16 x)  in bessel_jn_r16()  argument
58       size_t size = n2 < n1 ? 0 : n2-n1+1;   in bessel_jn_r16()
64   if (unlikely (n2 < n1))  in bessel_jn_r16()
68       && GFC_DESCRIPTOR_EXTENT(ret,0) != (n2-n1+1))  in bessel_jn_r16()
70 		  "(%ld vs. %ld)", (long int) n2-n1,  in bessel_jn_r16()
78       for (i = 1; i <= n2-n1; i++)  in bessel_jn_r16()
84   ret->base_addr[(n2-n1)*stride] = last1;  in bessel_jn_r16()
86   if (n1 == n2)  in bessel_jn_r16()
90   ret->base_addr[(n2-n1-1)*stride] = last2;  in bessel_jn_r16()
92   if (n1 + 1 == n2)  in bessel_jn_r16()
97   for (i = n2-n1-2; i >= 0; i--)  in bessel_jn_r16()
99       ret->base_addr[i*stride] = x2rev * (i+1+n1) * last2 - last1;  in bessel_jn_r16()
109 				     int n1, int n2, GFC_REAL_16 x);
113 bessel_yn_r16 (gfc_array_r16 * const restrict ret, int n1, int n2,  in bessel_yn_r16()  argument
125       size_t size = n2 < n1 ? 0 : n2-n1+1;   in bessel_yn_r16()
131   if (unlikely (n2 < n1))  in bessel_yn_r16()
135       && GFC_DESCRIPTOR_EXTENT(ret,0) != (n2-n1+1))  in bessel_yn_r16()
137 		  "(%ld vs. %ld)", (long int) n2-n1,  in bessel_yn_r16()
144       for (i = 0; i <= n2-n1; i++)  in bessel_yn_r16()
153   last1 = MATHFUNC(yn) (n1, x);  in bessel_yn_r16()
156   if (n1 == n2)  in bessel_yn_r16()
159   last2 = MATHFUNC(yn) (n1 + 1, x);  in bessel_yn_r16()
162   if (n1 + 1 == n2)  in bessel_yn_r16()
167   for (i = 2; i <= n2 - n1; i++)  in bessel_yn_r16()
177 	  ret->base_addr[i*stride] = x2rev * (i-1+n1) * last2 - last1;  in bessel_yn_r16()