Lines Matching refs:n1
38 extern void bessel_jn_r8 (gfc_array_r8 * const restrict ret, int n1,
43 bessel_jn_r8 (gfc_array_r8 * const restrict ret, int n1, int n2, GFC_REAL_8 x) in bessel_jn_r8() argument
54 size_t size = n2 < n1 ? 0 : n2-n1+1; in bessel_jn_r8()
60 if (unlikely (n2 < n1)) in bessel_jn_r8()
64 && GFC_DESCRIPTOR_EXTENT(ret,0) != (n2-n1+1)) in bessel_jn_r8()
66 "(%ld vs. %ld)", (long int) n2-n1, in bessel_jn_r8()
74 for (i = 1; i <= n2-n1; i++) in bessel_jn_r8()
80 ret->base_addr[(n2-n1)*stride] = last1; in bessel_jn_r8()
82 if (n1 == n2) in bessel_jn_r8()
86 ret->base_addr[(n2-n1-1)*stride] = last2; in bessel_jn_r8()
88 if (n1 + 1 == n2) in bessel_jn_r8()
93 for (i = n2-n1-2; i >= 0; i--) in bessel_jn_r8()
95 ret->base_addr[i*stride] = x2rev * (i+1+n1) * last2 - last1; in bessel_jn_r8()
105 int n1, int n2, GFC_REAL_8 x);
109 bessel_yn_r8 (gfc_array_r8 * const restrict ret, int n1, int n2, in bessel_yn_r8() argument
121 size_t size = n2 < n1 ? 0 : n2-n1+1; in bessel_yn_r8()
127 if (unlikely (n2 < n1)) in bessel_yn_r8()
131 && GFC_DESCRIPTOR_EXTENT(ret,0) != (n2-n1+1)) in bessel_yn_r8()
133 "(%ld vs. %ld)", (long int) n2-n1, in bessel_yn_r8()
140 for (i = 0; i <= n2-n1; i++) in bessel_yn_r8()
149 last1 = MATHFUNC(yn) (n1, x); in bessel_yn_r8()
152 if (n1 == n2) in bessel_yn_r8()
155 last2 = MATHFUNC(yn) (n1 + 1, x); in bessel_yn_r8()
158 if (n1 + 1 == n2) in bessel_yn_r8()
163 for (i = 2; i <= n2 - n1; i++) in bessel_yn_r8()
173 ret->base_addr[i*stride] = x2rev * (i-1+n1) * last2 - last1; in bessel_yn_r8()