gsl 1.12.109 → 1.14.5

Sign up to get free protection for your applications and to get access to all the features.
Files changed (511) hide show
  1. data/AUTHORS +6 -0
  2. data/COPYING +339 -0
  3. data/ChangeLog +556 -0
  4. data/{README.rdoc → README} +3 -0
  5. data/Rakefile +54 -10
  6. data/THANKS +17 -0
  7. data/VERSION +1 -2
  8. data/examples/alf/alf.gp +15 -0
  9. data/examples/alf/alf.rb +32 -0
  10. data/examples/blas/blas.rb +13 -0
  11. data/examples/blas/dnrm2.rb +16 -0
  12. data/examples/blas/level1.rb +81 -0
  13. data/examples/blas/level2.rb +11 -0
  14. data/examples/blas/level3.rb +12 -0
  15. data/examples/bspline.rb +57 -0
  16. data/examples/cdf.rb +16 -0
  17. data/examples/cheb.rb +21 -0
  18. data/examples/combination.rb +23 -0
  19. data/examples/complex/RC-lpf.rb +47 -0
  20. data/examples/complex/add.rb +36 -0
  21. data/examples/complex/coerce.rb +14 -0
  22. data/examples/complex/complex.rb +25 -0
  23. data/examples/complex/fpmi.rb +70 -0
  24. data/examples/complex/functions.rb +77 -0
  25. data/examples/complex/michelson.rb +36 -0
  26. data/examples/complex/mul.rb +28 -0
  27. data/examples/complex/oscillator.rb +17 -0
  28. data/examples/complex/set.rb +37 -0
  29. data/examples/const/physconst.rb +151 -0
  30. data/examples/const/travel.rb +45 -0
  31. data/examples/deriv/demo.rb +13 -0
  32. data/examples/deriv/deriv.rb +36 -0
  33. data/examples/deriv/diff.rb +35 -0
  34. data/examples/dht.rb +42 -0
  35. data/examples/dirac.rb +56 -0
  36. data/examples/eigen/eigen.rb +34 -0
  37. data/examples/eigen/herm.rb +22 -0
  38. data/examples/eigen/narray.rb +9 -0
  39. data/examples/eigen/nonsymm.rb +37 -0
  40. data/examples/eigen/nonsymmv.rb +43 -0
  41. data/examples/eigen/qhoscillator.gp +35 -0
  42. data/examples/eigen/qhoscillator.rb +90 -0
  43. data/examples/eigen/vander.rb +41 -0
  44. data/examples/fft/fft.rb +17 -0
  45. data/examples/fft/fft2.rb +17 -0
  46. data/examples/fft/forward.rb +25 -0
  47. data/examples/fft/forward2.rb +26 -0
  48. data/examples/fft/radix2.rb +18 -0
  49. data/examples/fft/real-halfcomplex.rb +33 -0
  50. data/examples/fft/real-halfcomplex2.rb +30 -0
  51. data/examples/fft/realradix2.rb +19 -0
  52. data/examples/fft/sunspot.dat +256 -0
  53. data/examples/fft/sunspot.rb +16 -0
  54. data/examples/fit/expdata.dat +20 -0
  55. data/examples/fit/expfit.rb +31 -0
  56. data/examples/fit/gaussfit.rb +29 -0
  57. data/examples/fit/gaussian_2peaks.rb +34 -0
  58. data/examples/fit/hillfit.rb +40 -0
  59. data/examples/fit/lognormal.rb +26 -0
  60. data/examples/fit/lorentzfit.rb +22 -0
  61. data/examples/fit/multifit.rb +72 -0
  62. data/examples/fit/ndlinear.rb +133 -0
  63. data/examples/fit/nonlinearfit.rb +89 -0
  64. data/examples/fit/plot.gp +36 -0
  65. data/examples/fit/polyfit.rb +9 -0
  66. data/examples/fit/powerfit.rb +21 -0
  67. data/examples/fit/sigmoidfit.rb +40 -0
  68. data/examples/fit/sinfit.rb +22 -0
  69. data/examples/fit/wlinear.rb +46 -0
  70. data/examples/fresnel.rb +11 -0
  71. data/examples/function/function.rb +36 -0
  72. data/examples/function/log.rb +7 -0
  73. data/examples/function/min.rb +33 -0
  74. data/examples/function/sin.rb +10 -0
  75. data/examples/function/synchrotron.rb +18 -0
  76. data/examples/gallery/butterfly.rb +7 -0
  77. data/examples/gallery/cayley.rb +12 -0
  78. data/examples/gallery/cornu.rb +23 -0
  79. data/examples/gallery/eight.rb +11 -0
  80. data/examples/gallery/koch.rb +40 -0
  81. data/examples/gallery/lemniscate.rb +11 -0
  82. data/examples/gallery/polar.rb +11 -0
  83. data/examples/gallery/rgplot/cossin.rb +35 -0
  84. data/examples/gallery/rgplot/rgplot.replaced +0 -0
  85. data/examples/gallery/rgplot/roesller.rb +55 -0
  86. data/examples/gallery/roesller.rb +39 -0
  87. data/examples/gallery/scarabaeus.rb +14 -0
  88. data/examples/histogram/cauchy.rb +27 -0
  89. data/examples/histogram/cauchy.sh +2 -0
  90. data/examples/histogram/exponential.rb +19 -0
  91. data/examples/histogram/gauss.rb +16 -0
  92. data/examples/histogram/gsl-histogram.rb +40 -0
  93. data/examples/histogram/histo2d.rb +31 -0
  94. data/examples/histogram/histo3d.rb +34 -0
  95. data/examples/histogram/histogram-pdf.rb +27 -0
  96. data/examples/histogram/histogram.rb +26 -0
  97. data/examples/histogram/integral.rb +28 -0
  98. data/examples/histogram/poisson.rb +27 -0
  99. data/examples/histogram/power.rb +25 -0
  100. data/examples/histogram/rebin.rb +17 -0
  101. data/examples/histogram/smp.dat +5 -0
  102. data/examples/histogram/xexp.rb +21 -0
  103. data/examples/integration/ahmed.rb +21 -0
  104. data/examples/integration/cosmology.rb +75 -0
  105. data/examples/integration/friedmann.gp +16 -0
  106. data/examples/integration/friedmann.rb +35 -0
  107. data/examples/integration/gamma-zeta.rb +35 -0
  108. data/examples/integration/integration.rb +22 -0
  109. data/examples/integration/qag.rb +8 -0
  110. data/examples/integration/qag2.rb +14 -0
  111. data/examples/integration/qag3.rb +8 -0
  112. data/examples/integration/qagi.rb +28 -0
  113. data/examples/integration/qagi2.rb +49 -0
  114. data/examples/integration/qagiu.rb +29 -0
  115. data/examples/integration/qagp.rb +20 -0
  116. data/examples/integration/qags.rb +14 -0
  117. data/examples/integration/qawc.rb +18 -0
  118. data/examples/integration/qawf.rb +41 -0
  119. data/examples/integration/qawo.rb +29 -0
  120. data/examples/integration/qaws.rb +30 -0
  121. data/examples/integration/qng.rb +17 -0
  122. data/examples/interp/demo.gp +20 -0
  123. data/examples/interp/demo.rb +45 -0
  124. data/examples/interp/interp.rb +37 -0
  125. data/examples/interp/points +10 -0
  126. data/examples/interp/spline.rb +20 -0
  127. data/examples/jacobi/deriv.rb +40 -0
  128. data/examples/jacobi/integrate.rb +34 -0
  129. data/examples/jacobi/interp.rb +43 -0
  130. data/examples/jacobi/jacobi.rb +11 -0
  131. data/examples/linalg/HH.rb +15 -0
  132. data/examples/linalg/HH_narray.rb +13 -0
  133. data/examples/linalg/LQ_solve.rb +73 -0
  134. data/examples/linalg/LU.rb +84 -0
  135. data/examples/linalg/LU2.rb +31 -0
  136. data/examples/linalg/LU_narray.rb +24 -0
  137. data/examples/linalg/PTLQ.rb +47 -0
  138. data/examples/linalg/QR.rb +18 -0
  139. data/examples/linalg/QRPT.rb +47 -0
  140. data/examples/linalg/QR_solve.rb +78 -0
  141. data/examples/linalg/QR_solve_narray.rb +13 -0
  142. data/examples/linalg/SV.rb +16 -0
  143. data/examples/linalg/SV_narray.rb +12 -0
  144. data/examples/linalg/SV_solve.rb +49 -0
  145. data/examples/linalg/chol.rb +29 -0
  146. data/examples/linalg/chol_narray.rb +15 -0
  147. data/examples/linalg/complex.rb +57 -0
  148. data/examples/linalg/invert_narray.rb +10 -0
  149. data/examples/math/const.rb +67 -0
  150. data/examples/math/elementary.rb +35 -0
  151. data/examples/math/functions.rb +41 -0
  152. data/examples/math/inf_nan.rb +34 -0
  153. data/examples/math/minmax.rb +22 -0
  154. data/examples/math/power.rb +18 -0
  155. data/examples/math/test.rb +31 -0
  156. data/examples/matrix/a.dat +0 -0
  157. data/examples/matrix/add.rb +45 -0
  158. data/examples/matrix/b.dat +4 -0
  159. data/examples/matrix/cat.rb +31 -0
  160. data/examples/matrix/colvectors.rb +24 -0
  161. data/examples/matrix/complex.rb +41 -0
  162. data/examples/matrix/det.rb +29 -0
  163. data/examples/matrix/diagonal.rb +23 -0
  164. data/examples/matrix/get_all.rb +159 -0
  165. data/examples/matrix/hilbert.rb +31 -0
  166. data/examples/matrix/iterator.rb +19 -0
  167. data/examples/matrix/matrix.rb +57 -0
  168. data/examples/matrix/minmax.rb +53 -0
  169. data/examples/matrix/mul.rb +39 -0
  170. data/examples/matrix/rand.rb +20 -0
  171. data/examples/matrix/read.rb +29 -0
  172. data/examples/matrix/rowcol.rb +47 -0
  173. data/examples/matrix/set.rb +41 -0
  174. data/examples/matrix/set_all.rb +100 -0
  175. data/examples/matrix/view.rb +32 -0
  176. data/examples/matrix/view_all.rb +148 -0
  177. data/examples/matrix/write.rb +23 -0
  178. data/examples/min.rb +29 -0
  179. data/examples/monte/miser.rb +47 -0
  180. data/examples/monte/monte.rb +47 -0
  181. data/examples/monte/plain.rb +47 -0
  182. data/examples/monte/vegas.rb +46 -0
  183. data/examples/multimin/bundle.rb +66 -0
  184. data/examples/multimin/cqp.rb +109 -0
  185. data/examples/multimin/fdfminimizer.rb +40 -0
  186. data/examples/multimin/fminimizer.rb +41 -0
  187. data/examples/multiroot/demo.rb +36 -0
  188. data/examples/multiroot/fdfsolver.rb +50 -0
  189. data/examples/multiroot/fsolver.rb +33 -0
  190. data/examples/multiroot/fsolver2.rb +32 -0
  191. data/examples/multiroot/fsolver3.rb +26 -0
  192. data/examples/narray/histogram.rb +14 -0
  193. data/examples/narray/mandel.rb +27 -0
  194. data/examples/narray/narray.rb +28 -0
  195. data/examples/narray/narray2.rb +44 -0
  196. data/examples/narray/sf.rb +26 -0
  197. data/examples/ntuple/create.rb +17 -0
  198. data/examples/ntuple/project.rb +31 -0
  199. data/examples/odeiv/binarysystem.gp +23 -0
  200. data/examples/odeiv/binarysystem.rb +104 -0
  201. data/examples/odeiv/demo.gp +24 -0
  202. data/examples/odeiv/demo.rb +69 -0
  203. data/examples/odeiv/demo2.gp +26 -0
  204. data/examples/odeiv/duffing.rb +45 -0
  205. data/examples/odeiv/frei1.rb +109 -0
  206. data/examples/odeiv/frei2.rb +76 -0
  207. data/examples/odeiv/legendre.rb +52 -0
  208. data/examples/odeiv/odeiv.rb +32 -0
  209. data/examples/odeiv/odeiv2.rb +45 -0
  210. data/examples/odeiv/oscillator.rb +42 -0
  211. data/examples/odeiv/sedov.rb +97 -0
  212. data/examples/odeiv/whitedwarf.gp +40 -0
  213. data/examples/odeiv/whitedwarf.rb +158 -0
  214. data/examples/ool/conmin.rb +100 -0
  215. data/examples/ool/gencan.rb +99 -0
  216. data/examples/ool/pgrad.rb +100 -0
  217. data/examples/ool/spg.rb +100 -0
  218. data/examples/pdf/bernoulli.rb +5 -0
  219. data/examples/pdf/beta.rb +7 -0
  220. data/examples/pdf/binomiral.rb +10 -0
  221. data/examples/pdf/cauchy.rb +6 -0
  222. data/examples/pdf/chisq.rb +8 -0
  223. data/examples/pdf/exponential.rb +7 -0
  224. data/examples/pdf/exppow.rb +6 -0
  225. data/examples/pdf/fdist.rb +7 -0
  226. data/examples/pdf/flat.rb +7 -0
  227. data/examples/pdf/gamma.rb +8 -0
  228. data/examples/pdf/gauss-tail.rb +5 -0
  229. data/examples/pdf/gauss.rb +6 -0
  230. data/examples/pdf/geometric.rb +5 -0
  231. data/examples/pdf/gumbel.rb +6 -0
  232. data/examples/pdf/hypergeometric.rb +11 -0
  233. data/examples/pdf/landau.rb +5 -0
  234. data/examples/pdf/laplace.rb +7 -0
  235. data/examples/pdf/logarithmic.rb +5 -0
  236. data/examples/pdf/logistic.rb +6 -0
  237. data/examples/pdf/lognormal.rb +6 -0
  238. data/examples/pdf/neg-binomiral.rb +10 -0
  239. data/examples/pdf/pareto.rb +7 -0
  240. data/examples/pdf/pascal.rb +10 -0
  241. data/examples/pdf/poisson.rb +5 -0
  242. data/examples/pdf/rayleigh-tail.rb +6 -0
  243. data/examples/pdf/rayleigh.rb +6 -0
  244. data/examples/pdf/tdist.rb +6 -0
  245. data/examples/pdf/weibull.rb +8 -0
  246. data/examples/permutation/ex1.rb +22 -0
  247. data/examples/permutation/permutation.rb +16 -0
  248. data/examples/poly/bell.rb +6 -0
  249. data/examples/poly/bessel.rb +6 -0
  250. data/examples/poly/cheb.rb +6 -0
  251. data/examples/poly/cheb_II.rb +6 -0
  252. data/examples/poly/cubic.rb +9 -0
  253. data/examples/poly/demo.rb +20 -0
  254. data/examples/poly/eval.rb +28 -0
  255. data/examples/poly/eval_derivs.rb +14 -0
  256. data/examples/poly/fit.rb +21 -0
  257. data/examples/poly/hermite.rb +6 -0
  258. data/examples/poly/poly.rb +13 -0
  259. data/examples/poly/quadratic.rb +25 -0
  260. data/examples/random/diffusion.rb +34 -0
  261. data/examples/random/gaussian.rb +9 -0
  262. data/examples/random/generator.rb +27 -0
  263. data/examples/random/hdsobol.rb +21 -0
  264. data/examples/random/poisson.rb +9 -0
  265. data/examples/random/qrng.rb +19 -0
  266. data/examples/random/randomwalk.rb +37 -0
  267. data/examples/random/randomwalk2d.rb +19 -0
  268. data/examples/random/rayleigh.rb +36 -0
  269. data/examples/random/rng.rb +33 -0
  270. data/examples/random/rngextra.rb +14 -0
  271. data/examples/roots/bisection.rb +25 -0
  272. data/examples/roots/brent.rb +43 -0
  273. data/examples/roots/demo.rb +30 -0
  274. data/examples/roots/newton.rb +46 -0
  275. data/examples/roots/recombination.gp +12 -0
  276. data/examples/roots/recombination.rb +61 -0
  277. data/examples/roots/steffenson.rb +48 -0
  278. data/examples/sf/ShiChi.rb +6 -0
  279. data/examples/sf/SiCi.rb +6 -0
  280. data/examples/sf/airy_Ai.rb +8 -0
  281. data/examples/sf/airy_Bi.rb +8 -0
  282. data/examples/sf/bessel_IK.rb +12 -0
  283. data/examples/sf/bessel_JY.rb +13 -0
  284. data/examples/sf/beta_inc.rb +9 -0
  285. data/examples/sf/clausen.rb +6 -0
  286. data/examples/sf/dawson.rb +5 -0
  287. data/examples/sf/debye.rb +9 -0
  288. data/examples/sf/dilog.rb +6 -0
  289. data/examples/sf/ellint.rb +6 -0
  290. data/examples/sf/expint.rb +8 -0
  291. data/examples/sf/fermi.rb +10 -0
  292. data/examples/sf/gamma_inc_P.rb +9 -0
  293. data/examples/sf/gegenbauer.rb +8 -0
  294. data/examples/sf/hyperg.rb +7 -0
  295. data/examples/sf/laguerre.rb +19 -0
  296. data/examples/sf/lambertW.rb +5 -0
  297. data/examples/sf/legendre_P.rb +10 -0
  298. data/examples/sf/lngamma.rb +5 -0
  299. data/examples/sf/psi.rb +54 -0
  300. data/examples/sf/sphbessel.gp +27 -0
  301. data/examples/sf/sphbessel.rb +30 -0
  302. data/examples/sf/synchrotron.rb +5 -0
  303. data/examples/sf/transport.rb +10 -0
  304. data/examples/sf/zetam1.rb +5 -0
  305. data/examples/siman.rb +44 -0
  306. data/examples/sort/heapsort.rb +23 -0
  307. data/examples/sort/heapsort_vector_complex.rb +21 -0
  308. data/examples/sort/sort.rb +23 -0
  309. data/examples/sort/sort2.rb +16 -0
  310. data/examples/stats/mean.rb +17 -0
  311. data/examples/stats/statistics.rb +18 -0
  312. data/examples/stats/test.rb +9 -0
  313. data/examples/sum.rb +34 -0
  314. data/examples/tamu_anova.rb +18 -0
  315. data/examples/vector/a.dat +0 -0
  316. data/examples/vector/add.rb +56 -0
  317. data/examples/vector/b.dat +4 -0
  318. data/examples/vector/c.dat +3 -0
  319. data/examples/vector/collect.rb +26 -0
  320. data/examples/vector/compare.rb +28 -0
  321. data/examples/vector/complex.rb +51 -0
  322. data/examples/vector/complex_get_all.rb +85 -0
  323. data/examples/vector/complex_set_all.rb +131 -0
  324. data/examples/vector/complex_view_all.rb +77 -0
  325. data/examples/vector/connect.rb +22 -0
  326. data/examples/vector/decimate.rb +38 -0
  327. data/examples/vector/diff.rb +31 -0
  328. data/examples/vector/filescan.rb +17 -0
  329. data/examples/vector/floor.rb +23 -0
  330. data/examples/vector/get_all.rb +82 -0
  331. data/examples/vector/gnuplot.rb +38 -0
  332. data/examples/vector/graph.rb +28 -0
  333. data/examples/vector/histogram.rb +22 -0
  334. data/examples/vector/linspace.rb +24 -0
  335. data/examples/vector/log.rb +17 -0
  336. data/examples/vector/logic.rb +33 -0
  337. data/examples/vector/logspace.rb +25 -0
  338. data/examples/vector/minmax.rb +47 -0
  339. data/examples/vector/mul.rb +49 -0
  340. data/examples/vector/narray.rb +46 -0
  341. data/examples/vector/read.rb +29 -0
  342. data/examples/vector/set.rb +35 -0
  343. data/examples/vector/set_all.rb +121 -0
  344. data/examples/vector/smpv.dat +15 -0
  345. data/examples/vector/test.rb +43 -0
  346. data/examples/vector/test_gslblock.rb +58 -0
  347. data/examples/vector/vector.rb +110 -0
  348. data/examples/vector/view.rb +35 -0
  349. data/examples/vector/view_all.rb +73 -0
  350. data/examples/vector/where.rb +29 -0
  351. data/examples/vector/write.rb +24 -0
  352. data/examples/vector/zip.rb +34 -0
  353. data/examples/wavelet/ecg.dat +256 -0
  354. data/examples/wavelet/wavelet1.rb +50 -0
  355. data/ext/extconf.rb +9 -0
  356. data/ext/gsl.c +10 -1
  357. data/ext/histogram.c +6 -2
  358. data/ext/integration.c +39 -0
  359. data/ext/matrix_complex.c +1 -1
  360. data/ext/multiset.c +214 -0
  361. data/ext/nmf.c +4 -0
  362. data/ext/nmf_wrap.c +3 -0
  363. data/ext/vector_complex.c +1 -1
  364. data/ext/vector_double.c +3 -3
  365. data/ext/vector_source.c +6 -6
  366. data/include/rb_gsl.h +7 -0
  367. data/include/rb_gsl_common.h +6 -0
  368. data/rdoc/alf.rdoc +77 -0
  369. data/rdoc/blas.rdoc +269 -0
  370. data/rdoc/bspline.rdoc +42 -0
  371. data/rdoc/changes.rdoc +164 -0
  372. data/rdoc/cheb.rdoc +99 -0
  373. data/rdoc/cholesky_complex.rdoc +46 -0
  374. data/rdoc/combi.rdoc +125 -0
  375. data/rdoc/complex.rdoc +210 -0
  376. data/rdoc/const.rdoc +546 -0
  377. data/rdoc/dht.rdoc +122 -0
  378. data/rdoc/diff.rdoc +133 -0
  379. data/rdoc/ehandling.rdoc +50 -0
  380. data/rdoc/eigen.rdoc +401 -0
  381. data/rdoc/fft.rdoc +535 -0
  382. data/rdoc/fit.rdoc +284 -0
  383. data/rdoc/function.rdoc +94 -0
  384. data/rdoc/graph.rdoc +137 -0
  385. data/rdoc/hist.rdoc +409 -0
  386. data/rdoc/hist2d.rdoc +279 -0
  387. data/rdoc/hist3d.rdoc +112 -0
  388. data/rdoc/index.rdoc +62 -0
  389. data/rdoc/integration.rdoc +398 -0
  390. data/rdoc/interp.rdoc +231 -0
  391. data/rdoc/intro.rdoc +27 -0
  392. data/rdoc/linalg.rdoc +681 -0
  393. data/rdoc/linalg_complex.rdoc +88 -0
  394. data/rdoc/math.rdoc +276 -0
  395. data/rdoc/matrix.rdoc +1093 -0
  396. data/rdoc/min.rdoc +189 -0
  397. data/rdoc/monte.rdoc +234 -0
  398. data/rdoc/multimin.rdoc +312 -0
  399. data/rdoc/multiroot.rdoc +293 -0
  400. data/rdoc/narray.rdoc +173 -0
  401. data/rdoc/ndlinear.rdoc +247 -0
  402. data/rdoc/nonlinearfit.rdoc +348 -0
  403. data/rdoc/ntuple.rdoc +88 -0
  404. data/rdoc/odeiv.rdoc +378 -0
  405. data/rdoc/perm.rdoc +221 -0
  406. data/rdoc/poly.rdoc +335 -0
  407. data/rdoc/qrng.rdoc +90 -0
  408. data/rdoc/randist.rdoc +233 -0
  409. data/rdoc/ref.rdoc +93 -0
  410. data/rdoc/rng.rdoc +203 -0
  411. data/rdoc/rngextra.rdoc +11 -0
  412. data/rdoc/roots.rdoc +305 -0
  413. data/rdoc/screenshot.rdoc +40 -0
  414. data/rdoc/sf.rdoc +1622 -0
  415. data/rdoc/siman.rdoc +89 -0
  416. data/rdoc/sort.rdoc +94 -0
  417. data/rdoc/start.rdoc +16 -0
  418. data/rdoc/stats.rdoc +219 -0
  419. data/rdoc/sum.rdoc +65 -0
  420. data/rdoc/tensor.rdoc +251 -0
  421. data/rdoc/tut.rdoc +5 -0
  422. data/rdoc/use.rdoc +177 -0
  423. data/rdoc/vector.rdoc +1243 -0
  424. data/rdoc/vector_complex.rdoc +347 -0
  425. data/rdoc/wavelet.rdoc +218 -0
  426. data/setup.rb +1585 -0
  427. data/tests/blas/amax.rb +14 -0
  428. data/tests/blas/asum.rb +16 -0
  429. data/tests/blas/axpy.rb +25 -0
  430. data/tests/blas/copy.rb +23 -0
  431. data/tests/blas/dot.rb +23 -0
  432. data/tests/bspline.rb +53 -0
  433. data/tests/cdf.rb +1388 -0
  434. data/tests/cheb.rb +112 -0
  435. data/tests/combination.rb +123 -0
  436. data/tests/complex.rb +17 -0
  437. data/tests/const.rb +24 -0
  438. data/tests/deriv.rb +85 -0
  439. data/tests/dht/dht1.rb +17 -0
  440. data/tests/dht/dht2.rb +23 -0
  441. data/tests/dht/dht3.rb +23 -0
  442. data/tests/dht/dht4.rb +23 -0
  443. data/tests/diff.rb +78 -0
  444. data/tests/eigen/eigen.rb +220 -0
  445. data/tests/eigen/gen.rb +105 -0
  446. data/tests/eigen/genherm.rb +66 -0
  447. data/tests/eigen/gensymm.rb +68 -0
  448. data/tests/eigen/nonsymm.rb +53 -0
  449. data/tests/eigen/nonsymmv.rb +53 -0
  450. data/tests/eigen/symm-herm.rb +74 -0
  451. data/tests/err.rb +58 -0
  452. data/tests/fit.rb +124 -0
  453. data/tests/gsl_test.rb +118 -0
  454. data/tests/gsl_test2.rb +107 -0
  455. data/tests/histo.rb +12 -0
  456. data/tests/integration/integration1.rb +72 -0
  457. data/tests/integration/integration2.rb +71 -0
  458. data/tests/integration/integration3.rb +71 -0
  459. data/tests/integration/integration4.rb +71 -0
  460. data/tests/interp.rb +45 -0
  461. data/tests/linalg/HH.rb +64 -0
  462. data/tests/linalg/LU.rb +47 -0
  463. data/tests/linalg/QR.rb +77 -0
  464. data/tests/linalg/SV.rb +24 -0
  465. data/tests/linalg/TDN.rb +116 -0
  466. data/tests/linalg/TDS.rb +122 -0
  467. data/tests/linalg/bidiag.rb +73 -0
  468. data/tests/linalg/cholesky.rb +20 -0
  469. data/tests/linalg/linalg.rb +158 -0
  470. data/tests/matrix/matrix_nmf_test.rb +39 -0
  471. data/tests/matrix/matrix_test.rb +48 -0
  472. data/tests/min.rb +99 -0
  473. data/tests/monte/miser.rb +31 -0
  474. data/tests/monte/vegas.rb +45 -0
  475. data/tests/multifit/test_2dgauss.rb +112 -0
  476. data/tests/multifit/test_brown.rb +90 -0
  477. data/tests/multifit/test_enso.rb +246 -0
  478. data/tests/multifit/test_filip.rb +155 -0
  479. data/tests/multifit/test_gauss.rb +97 -0
  480. data/tests/multifit/test_longley.rb +110 -0
  481. data/tests/multifit/test_multifit.rb +52 -0
  482. data/tests/multimin.rb +139 -0
  483. data/tests/multiroot.rb +131 -0
  484. data/tests/multiset.rb +52 -0
  485. data/tests/odeiv.rb +353 -0
  486. data/tests/poly/poly.rb +242 -0
  487. data/tests/poly/special.rb +65 -0
  488. data/tests/qrng.rb +131 -0
  489. data/tests/quartic.rb +29 -0
  490. data/tests/randist.rb +134 -0
  491. data/tests/rng.rb +305 -0
  492. data/tests/roots.rb +76 -0
  493. data/tests/run-test.sh +17 -0
  494. data/tests/sf/gsl_test_sf.rb +249 -0
  495. data/tests/sf/test_airy.rb +83 -0
  496. data/tests/sf/test_bessel.rb +306 -0
  497. data/tests/sf/test_coulomb.rb +17 -0
  498. data/tests/sf/test_dilog.rb +25 -0
  499. data/tests/sf/test_gamma.rb +209 -0
  500. data/tests/sf/test_hyperg.rb +356 -0
  501. data/tests/sf/test_legendre.rb +227 -0
  502. data/tests/sf/test_mathieu.rb +59 -0
  503. data/tests/sf/test_sf.rb +839 -0
  504. data/tests/stats.rb +174 -0
  505. data/tests/sum.rb +98 -0
  506. data/tests/sys.rb +323 -0
  507. data/tests/tensor.rb +419 -0
  508. data/tests/vector/vector_complex_test.rb +101 -0
  509. data/tests/vector/vector_test.rb +141 -0
  510. data/tests/wavelet.rb +142 -0
  511. metadata +596 -15
@@ -0,0 +1,227 @@
1
+ #!/usr/bin/env ruby
2
+ require("gsl")
3
+ require("../gsl_test2.rb")
4
+ require("./gsl_test_sf.rb")
5
+ include GSL::Test
6
+ include Math
7
+ include GSL::Test::Sf
8
+ s = 0
9
+ m = GSL::MODE_DEFAULT
10
+ TEST_SF(s, "GSL::Sf::legendre_P1_e", "(-0.5)", -0.5, TEST_TOL0, GSL::SUCCESS)
11
+ TEST_SF(s, "GSL::Sf::legendre_P1_e", "( 0.5)", 0.5, TEST_TOL0, GSL::SUCCESS)
12
+ TEST_SF(s, "GSL::Sf::legendre_P2_e", "(0.0)", -0.5 , TEST_TOL0, GSL::SUCCESS)
13
+ TEST_SF(s, "GSL::Sf::legendre_P2_e", "(0.5)", -0.125, TEST_TOL0, GSL::SUCCESS)
14
+ TEST_SF(s, "GSL::Sf::legendre_P2_e", "(1.0)", 1.0 , TEST_TOL0, GSL::SUCCESS)
15
+ TEST_SF(s, "GSL::Sf::legendre_P2_e", "(100.0)", 14999.5 , TEST_TOL0, GSL::SUCCESS)
16
+ TEST_SF(s, "GSL::Sf::legendre_P3_e", "( -0.5)", 0.4375, TEST_TOL0, GSL::SUCCESS)
17
+ TEST_SF(s, "GSL::Sf::legendre_P3_e", "( 0.5)", -0.4375, TEST_TOL0, GSL::SUCCESS)
18
+ TEST_SF(s, "GSL::Sf::legendre_P3_e", "( 1.0)", 1.0 , TEST_TOL0, GSL::SUCCESS)
19
+ TEST_SF(s, "GSL::Sf::legendre_P3_e", "(100.0)", 2.49985e+06, TEST_TOL0, GSL::SUCCESS)
20
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(1, -0.5)", -0.5, TEST_TOL0, GSL::SUCCESS)
21
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(1, 1.0e-8)", 1.0e-08, TEST_TOL0, GSL::SUCCESS)
22
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(1, 0.5)", 0.5, TEST_TOL0, GSL::SUCCESS)
23
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(1, 1.0)", 1.0, TEST_TOL0, GSL::SUCCESS)
24
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(10, -0.5)", -0.1882286071777345, TEST_TOL0, GSL::SUCCESS)
25
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(10, 1.0e-8)", -0.24609374999999864648, TEST_TOL0, GSL::SUCCESS)
26
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(10, 0.5)", -0.18822860717773437500, TEST_TOL0, GSL::SUCCESS)
27
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(10, 1.0)", 1.0, TEST_TOL0, GSL::SUCCESS)
28
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(99, -0.5)", 0.08300778172138770477, TEST_TOL0, GSL::SUCCESS)
29
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(99, 1.0e-8)", -7.958923738716563193e-08, TEST_TOL0, GSL::SUCCESS)
30
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(99, 0.5)", -0.08300778172138770477, TEST_TOL0, GSL::SUCCESS)
31
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(99, 0.999)", -0.3317727359254778874, TEST_TOL2, GSL::SUCCESS)
32
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(99, 1.0)", 1.0, TEST_TOL0, GSL::SUCCESS)
33
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(1000, -0.5)", -0.019168251091650277878, TEST_TOL2, GSL::SUCCESS)
34
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(1000, 1.0e-8)", 0.02522501817709828, TEST_TOL2, GSL::SUCCESS)
35
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(1000, 0.5)", -0.019168251091650277878, TEST_TOL2, GSL::SUCCESS)
36
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(1000, 1.0)", 1.0, TEST_TOL0, GSL::SUCCESS)
37
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(4000, -0.5)", -0.009585404456573080972, TEST_TOL2, GSL::SUCCESS)
38
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(4000, 0.5)", -0.009585404456573080972, TEST_TOL2, GSL::SUCCESS)
39
+ TEST_SF(s, "GSL::Sf::legendre_Pl_e", "(4000, 1.0)", 1.0, TEST_TOL0, GSL::SUCCESS)
40
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(10, 0, -0.5)", -0.18822860717773437500, TEST_TOL0, GSL::SUCCESS)
41
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(10, 0, 1.0e-08)", -0.24609374999999864648, TEST_TOL0, GSL::SUCCESS)
42
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(10, 0, 0.5)", -0.18822860717773437500, TEST_TOL0, GSL::SUCCESS)
43
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(10, 1, -0.5)", -2.0066877394361256516, TEST_TOL0, GSL::SUCCESS)
44
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(10, 1, 1.0e-08)", -2.7070312499999951725e-07, TEST_TOL0, GSL::SUCCESS)
45
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(10, 1, 0.5)", 2.0066877394361256516, TEST_TOL0, GSL::SUCCESS)
46
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(10, 5, -0.5)", -30086.169706116174977, TEST_TOL0, GSL::SUCCESS)
47
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(10, 5, 1.0e-08)", -0.0025337812499999964949, TEST_TOL0, GSL::SUCCESS)
48
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(10, 5, 0.5)", 30086.169706116174977, TEST_TOL0, GSL::SUCCESS)
49
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(10, 5, 0.999)", -0.5036411489013270406, TEST_TOL1, GSL::SUCCESS)
50
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(100, 5, -0.5)", -6.617107444248382171e+08, TEST_TOL0, GSL::SUCCESS)
51
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(100, 5, 1.0e-08)", 817.8987598063712851, TEST_TOL0, GSL::SUCCESS)
52
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(100, 5, 0.5)", 6.617107444248382171e+08, TEST_TOL0, GSL::SUCCESS)
53
+ TEST_SF(s, "GSL::Sf::legendre_Plm_e", "(100, 5, 0.999)", -1.9831610803806212189e+09, TEST_TOL2, GSL::SUCCESS)
54
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(10, 0, -0.5)", -0.24332702369300133776, TEST_TOL0, GSL::SUCCESS)
55
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(10, 0, 0.5)", -0.24332702369300133776, TEST_TOL0, GSL::SUCCESS)
56
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(10, 0, 0.999)", 1.2225754122797385990, TEST_TOL1, GSL::SUCCESS)
57
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(10, 5, -0.5)", -0.3725739049803293972, TEST_TOL0, GSL::SUCCESS)
58
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(10, 5, 1.0e-08)", -3.1377233589376792243e-08, TEST_TOL0, GSL::SUCCESS)
59
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(10, 5, 0.5)", 0.3725739049803293972, TEST_TOL0, GSL::SUCCESS)
60
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(10, 5, 0.999)", -6.236870674727370094e-06, TEST_TOL2, GSL::SUCCESS)
61
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(10, 10, -0.5)", 0.12876871185785724117, TEST_TOL1, GSL::SUCCESS)
62
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(10, 10, 0.5)", 0.12876871185785724117, TEST_TOL1, GSL::SUCCESS)
63
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(10, 10, 0.999)", 1.7320802307583118647e-14, TEST_TOL2, GSL::SUCCESS)
64
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(200, 1, -0.5)", 0.3302975570099492931, TEST_TOL1, GSL::SUCCESS)
65
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(200, 1, 0.5)", -0.3302975570099492931, TEST_TOL1, GSL::SUCCESS)
66
+ TEST_SF(s, "GSL::Sf::legendre_sphPlm_e", "(200, 1, 0.999)", -1.4069792055546256912, TEST_TOL2, GSL::SUCCESS)
67
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(0.0, -0.5)", 0.8573827581049917129, TEST_TOL0, GSL::SUCCESS)
68
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(0.0, 0.5)", 0.8573827581049917129, TEST_TOL0, GSL::SUCCESS)
69
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(0.0, 2.0)", 0.6062611623284649811, TEST_TOL0, GSL::SUCCESS)
70
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(0.0, 100.0)", 0.07979045091636735635, TEST_TOL0, GSL::SUCCESS)
71
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(10.0, -0.5)", 5.345484922591867188e+08, TEST_TOL1, GSL::SUCCESS)
72
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(10.0, 0.5)", 15137.910380385258370, TEST_TOL1, GSL::SUCCESS)
73
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(10.0, 2.0)", 0.4992680691891618544, TEST_TOL1, GSL::SUCCESS)
74
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(10.0, 100.0)", -0.07272008163718195685, TEST_TOL2, GSL::SUCCESS)
75
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(200.0, -1.0e-3)", 1.3347639529084185010e+136, TEST_TOL2, GSL::SUCCESS)
76
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(200.0, 1.0e-8)", 1.0928098010940058507e+136, TEST_TOL2, GSL::SUCCESS)
77
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(200.0, 0.5)", 3.895546021611205442e+90, TEST_TOL2, GSL::SUCCESS)
78
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(200.0, 10.0)", -0.04308567180833581268, TEST_TOL3, GSL::SUCCESS)
79
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(200.0, 100.0)", -0.04694669186576399194, TEST_TOL3, GSL::SUCCESS)
80
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(200.0, 1000.0)", 0.023698140704121273277, TEST_TOL3, GSL::SUCCESS)
81
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(200.0, 1.0e+8)", -0.00006790983312124277891, TEST_TOL3, GSL::SUCCESS)
82
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(1.0e+8, 1.1)", 1.1599311133054742944, TEST_SQRT_TOL0, GSL::SUCCESS)
83
+ TEST_SF(s, "GSL::Sf::conicalP_half_e", "(1.0e+8, 100.0)", 0.07971967557381557875, TEST_SQRT_TOL0, GSL::SUCCESS)
84
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(0.0, -0.5)", 1.7956982494514644808, TEST_TOL0, GSL::SUCCESS)
85
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(0.0, 0.5)", 0.8978491247257322404, TEST_TOL0, GSL::SUCCESS)
86
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(0.0, 2.0)", 0.7984204253272901551, TEST_TOL0, GSL::SUCCESS)
87
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(0.0, 100.0)", 0.4227531369388072584, TEST_TOL0, GSL::SUCCESS)
88
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(10.0, -0.5)", 5.345484922591867181e+07, TEST_TOL1, GSL::SUCCESS)
89
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(10.0, 0.5)", 1513.7910356104985334, TEST_TOL1, GSL::SUCCESS)
90
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(10.0, 2.0)", 0.03439243987215615642, TEST_TOL1, GSL::SUCCESS)
91
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(10.0, 100.0)", 0.003283756665952609624, TEST_TOL2, GSL::SUCCESS)
92
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(200.0, -0.5)", 1.7699538115312304280e+179, TEST_TOL2, GSL::SUCCESS)
93
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(200.0, 1.0e-8)", 5.464049005470029253e+133, TEST_TOL2, GSL::SUCCESS)
94
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(200.0, 0.5)", 1.9477730108056027211e+88, TEST_TOL2, GSL::SUCCESS)
95
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(200.0, 10.0)", 0.0012462575917716355362, TEST_TOL2, GSL::SUCCESS)
96
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(200.0, 100.0)", -0.0003225881344802625149, TEST_TOL2, GSL::SUCCESS)
97
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(200.0, 1000.0)", -0.00004330652890886567623, TEST_TOL3, GSL::SUCCESS)
98
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(200.0, 1.0e+8)", 2.0943091278037078483e-07, TEST_TOL3, GSL::SUCCESS)
99
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(1.0e+8, 1.1)", 2.092320445620989618e-09, 16.0*TEST_SQRT_TOL0, GSL::SUCCESS)
100
+ TEST_SF(s, "GSL::Sf::conicalP_mhalf_e", "(1.0e+8, 100.0)", -3.359967833599016923e-11, 256.0*TEST_SQRT_TOL0, GSL::SUCCESS)
101
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(0.0, -0.5)", 1.3728805006183501647, TEST_TOL0, GSL::SUCCESS)
102
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(0.0, 0.5)", 1.0731820071493643751, TEST_TOL0, GSL::SUCCESS)
103
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(0.0, 2.0)", 0.9012862993604472987, TEST_TOL0, GSL::SUCCESS)
104
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(0.0, 100.0)", 0.30091748588199264556, TEST_TOL0, GSL::SUCCESS)
105
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(10.0, -0.5)", 1.6795592815421804669e+08, TEST_TOL1, GSL::SUCCESS)
106
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(10.0, 0.5)", 4826.034132009618240, TEST_TOL1, GSL::SUCCESS)
107
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(10.0, 2.0)", 0.18798468917758716146, TEST_TOL2, GSL::SUCCESS)
108
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(10.0, 100.0)", -0.008622130749987962529, TEST_TOL2, GSL::SUCCESS)
109
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(200.0, -0.5)", 2.502194818646823e+180, TEST_TOL4, GSL::SUCCESS)
110
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(1000.0, 100.0)", 0.0017908817653497715844, TEST_TOL3, GSL::SUCCESS)
111
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(1000.0, 1000.0)", -0.0006566893804926284301, TEST_TOL3, GSL::SUCCESS)
112
+ TEST_SF(s, "GSL::Sf::conicalP_0_e", "(1000.0, 1.0e+8)", 2.3167213561756390068e-06, TEST_TOL4, GSL::SUCCESS)
113
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(0.0, -0.5)", 0.4939371126656998499, TEST_TOL1, GSL::SUCCESS)
114
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(0.0, 0.5)", 0.14933621085538265636, TEST_TOL1, GSL::SUCCESS)
115
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(0.0, 2.0)", -0.13666874968871549533, TEST_TOL1, GSL::SUCCESS)
116
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(0.0, 100.0)", -0.10544528203156629098, TEST_TOL2, GSL::SUCCESS)
117
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(10.0, -0.5)", 1.7253802958788312520e+09, TEST_TOL2, GSL::SUCCESS)
118
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(10.0, 0.5)", 46781.02294059967988, TEST_TOL1, GSL::SUCCESS)
119
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(10.0, 2.0)", 0.26613342643657444400, TEST_TOL2, GSL::SUCCESS)
120
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(10.0, 100.0)", -0.23281959695501029796, TEST_TOL2, GSL::SUCCESS)
121
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(200.0, -0.999)", 2.71635193199341135e+270, TEST_TOL2, GSL::SUCCESS)
122
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(200.0, -0.9)", 4.2952493176812905e+234, TEST_TOL2, GSL::SUCCESS)
123
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(200.0, -0.5)", 5.01159205956053439e+182, TEST_TOL3, GSL::SUCCESS)
124
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(200.0, 0.999)", 195733.0396081538, TEST_TOL2, GSL::SUCCESS)
125
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(200.0, 10.0)", -2.9272610662414349553, TEST_TOL2, GSL::SUCCESS)
126
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(1000.0, 100.0)", -1.7783258105862399857, TEST_TOL6, GSL::SUCCESS)
127
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(1000.0, 1000.0)", 0.4535161075156427179, TEST_TOL4, GSL::SUCCESS)
128
+ TEST_SF(s, "GSL::Sf::conicalP_1_e", "(1000.0, 1.0e+8)", 0.0009983414549874888478, TEST_SQRT_TOL0, GSL::SUCCESS)
129
+ TEST_SF(s, "GSL::Sf::conicalP_sph_reg_e", "(2, 1.0, -0.5)", 1.6406279287008789526, TEST_TOL0, GSL::SUCCESS)
130
+ TEST_SF(s, "GSL::Sf::conicalP_sph_reg_e", "(10, 1.0, -0.5)", 0.000029315266725049129448, TEST_TOL1, GSL::SUCCESS)
131
+ TEST_SF(s, "GSL::Sf::conicalP_sph_reg_e", "(20, 1.0, -0.5)", 7.335769429462034431e-15, TEST_TOL1, GSL::SUCCESS)
132
+ TEST_SF(s, "GSL::Sf::conicalP_sph_reg_e", "(30, 1.0, -0.5)", 1.3235612394267378871e-26, TEST_TOL2, GSL::SUCCESS)
133
+ TEST_SF(s, "GSL::Sf::conicalP_sph_reg_e", "(10, 1.0, 0.5)", 2.7016087199857873954e-10, TEST_TOL1, GSL::SUCCESS)
134
+ TEST_SF(s, "GSL::Sf::conicalP_sph_reg_e", "(20, 1.0, 0.5)", 1.1782569701435933399e-24, TEST_TOL1, GSL::SUCCESS)
135
+ TEST_SF(s, "GSL::Sf::conicalP_sph_reg_e", "(30, 1.0, 0.5)", 3.636240588303797919e-41, TEST_TOL1, GSL::SUCCESS)
136
+ TEST_SF(s, "GSL::Sf::conicalP_sph_reg_e", "(10, 1.0, 2.0)", 2.4934929626284934483e-10, TEST_TOL1, GSL::SUCCESS)
137
+ TEST_SF(s, "GSL::Sf::conicalP_sph_reg_e", "(20, 1.0, 2.0)", 1.1284762488012616191e-24, TEST_TOL2, GSL::SUCCESS)
138
+ TEST_SF(s, "GSL::Sf::conicalP_sph_reg_e", "(30, 100.0, 100.0)", -1.6757772087159526048e-64, TEST_TOL6, GSL::SUCCESS)
139
+ TEST_SF(s, "GSL::Sf::conicalP_cyl_reg_e", "(2, 1.0, -0.5)", 2.2048510472375258708, TEST_TOL0, GSL::SUCCESS)
140
+ TEST_SF(s, "GSL::Sf::conicalP_cyl_reg_e", "(10, 1.0, -0.5)", 0.00007335034531618655690, TEST_TOL1, GSL::SUCCESS)
141
+ TEST_SF(s, "GSL::Sf::conicalP_cyl_reg_e", "(20, 1.0, -0.5)", 2.5419860619212164696e-14, TEST_TOL1, GSL::SUCCESS)
142
+ TEST_SF(s, "GSL::Sf::conicalP_cyl_reg_e", "(30, 1.0, -0.5)", 5.579714972260536827e-26, TEST_TOL2, GSL::SUCCESS)
143
+ TEST_SF(s, "GSL::Sf::conicalP_cyl_reg_e", "(10, 1.0, 0.5)", 1.1674078819646475282e-09, TEST_TOL0, GSL::SUCCESS)
144
+ TEST_SF(s, "GSL::Sf::conicalP_cyl_reg_e", "(20, 1.0, 0.5)", 7.066408031229072207e-24, TEST_TOL1, GSL::SUCCESS)
145
+ TEST_SF(s, "GSL::Sf::conicalP_cyl_reg_e", "(30, 1.0, 0.5)", 2.6541973286862588488e-40, TEST_TOL1, GSL::SUCCESS)
146
+ TEST_SF(s, "GSL::Sf::conicalP_cyl_reg_e", "(10, 1.0, 2.0)", 1.0736109751890863051e-09, TEST_TOL2, GSL::SUCCESS)
147
+ TEST_SF(s, "GSL::Sf::conicalP_cyl_reg_e", "(20, 1.0, 2.0)", 6.760965304863386741e-24, TEST_TOL2, GSL::SUCCESS)
148
+ TEST_SF(s, "GSL::Sf::conicalP_cyl_reg_e", "(30, 100.0, 100.0)", -4.268753482520651007e-63, TEST_TOL4, GSL::SUCCESS)
149
+ TEST_SF(s, "GSL::Sf::legendre_H3d_0_e", "(1.0e-06, 1.0e-06)", 0.9999999999998333333 , TEST_TOL0, GSL::SUCCESS)
150
+ TEST_SF(s, "GSL::Sf::legendre_H3d_0_e", "(1.0, 0.0)", 1.0 , TEST_TOL0, GSL::SUCCESS)
151
+ TEST_SF(s, "GSL::Sf::legendre_H3d_0_e", "(1.0, 1.0)", 0.7160229153604338713 , TEST_TOL0, GSL::SUCCESS)
152
+ TEST_SF(s, "GSL::Sf::legendre_H3d_0_e", "(1.0, 100.0)", -3.767437313149604566e-44 , TEST_TOL2, GSL::SUCCESS)
153
+ TEST_SF(s, "GSL::Sf::legendre_H3d_0_e", "(1.0, 500.0)", -6.665351935878582205e-218, TEST_TOL2, GSL::SUCCESS)
154
+ TEST_SF(s, "GSL::Sf::legendre_H3d_0_e", "(100.0, 1.0)", -0.004308757035378200029 , TEST_TOL0, GSL::SUCCESS)
155
+ TEST_SF(s, "GSL::Sf::legendre_H3d_0_e", "(100.0, 10.0)", 7.508054627912986427e-07 , TEST_TOL0, GSL::SUCCESS)
156
+ TEST_SF(s, "GSL::Sf::legendre_H3d_0_e", "(1000.0, 1.0)", 0.0007036067909088818319 , TEST_TOL0, GSL::SUCCESS)
157
+ TEST_SF(s, "GSL::Sf::legendre_H3d_0_e", "(1.0e+08, 1.0)", 7.927485371429105968e-09 , TEST_TOL3, GSL::SUCCESS)
158
+ TEST_SF(s, "GSL::Sf::legendre_H3d_0_e", "(1.0e+08, 100.0)", -3.627118904186918957e-52 , 32.0*TEST_SQRT_TOL0, GSL::SUCCESS)
159
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(1.0e-06, 1.0e-06)", 3.333333333334222222e-07, TEST_TOL0, GSL::SUCCESS)
160
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(1.0, 1.0e-10)", 4.714045207910316829e-11, TEST_TOL0, GSL::SUCCESS)
161
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(1.0, 1.0)", 0.3397013994799344639, TEST_TOL0, GSL::SUCCESS)
162
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(1.0, 100.0)", -7.200624449531811272e-44, TEST_TOL2, GSL::SUCCESS)
163
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(1.0, 500.0)", 4.192260336821728677e-218, TEST_TOL2, GSL::SUCCESS)
164
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(100.0, 0.01)", 0.30117664944267412324 , TEST_TOL1, GSL::SUCCESS)
165
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(100.0, 1.0)", -0.007393833425336299309 , TEST_TOL0, GSL::SUCCESS)
166
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(100.0, 10.0)", -5.031062029821254982e-07 , TEST_TOL0, GSL::SUCCESS)
167
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(1000.0, 0.001)", 0.30116875865090396421 , TEST_TOL0, GSL::SUCCESS)
168
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(1000.0, 1.0)", -0.0004776144516074971885 , TEST_TOL0, GSL::SUCCESS)
169
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(1.0e+08, 1.0e-08)", 0.30116867893975679722 , TEST_TOL1, GSL::SUCCESS)
170
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(1.0e+08, 1.0)", 3.0921097047369081582e-09, TEST_TOL4, GSL::SUCCESS)
171
+ TEST_SF(s, "GSL::Sf::legendre_H3d_1_e", "(1.0e+08, 100.0)", -6.496142701296286936e-52 , 32.0*TEST_SQRT_TOL0, GSL::SUCCESS)
172
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1.0e-06, 1.0e-06)", 1.1544011544013627977e-32, TEST_TOL2, GSL::SUCCESS)
173
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1.0, 1.0e-10)", 2.0224912016958766992e-52, TEST_TOL2, GSL::SUCCESS)
174
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1.0, 1.0)", 0.011498635037491577728, TEST_TOL1, GSL::SUCCESS)
175
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1.0, 5.0)", 0.0020696945662545205776, TEST_TOL4, GSL::SUCCESS)
176
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1.0, 7.0)", -0.0017555303787488993676, TEST_TOL4, GSL::SUCCESS)
177
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1.0, 10.0)", 0.00008999979724504887101, TEST_TOL2, GSL::SUCCESS)
178
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1.0, 100.0)", -4.185397793298567945e-44, TEST_TOL2, GSL::SUCCESS)
179
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1.0, 500.0)", 1.4235113901091961263e-217, TEST_TOL3, GSL::SUCCESS)
180
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 100.0, 0.001)", 9.642762597222417946e-10, TEST_TOL2, GSL::SUCCESS)
181
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 100.0, 0.002)", 3.0821201254308036109e-08, TEST_TOL2, GSL::SUCCESS)
182
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 100.0, 0.01)", 0.00009281069019005840532, TEST_TOL1, GSL::SUCCESS)
183
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 100.0, 1.0)", -0.008043100696178624653, TEST_TOL2, GSL::SUCCESS)
184
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 100.0, 10.0)", -3.927678432813974207e-07, TEST_TOL3, GSL::SUCCESS)
185
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1000.0, 0.001)", 0.00009256365284253254503, TEST_TOL1, GSL::SUCCESS)
186
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1000.0, 0.01)", -0.05553733815473079983, TEST_TOL0, GSL::SUCCESS)
187
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1.0e+08, 1.0e-08)", 0.00009256115861125841299, TEST_TOL2, GSL::SUCCESS)
188
+ TEST_SF(s, "GSL::Sf::legendre_H3d_e", "(5, 1.0e+08, 100.0)", -6.496143209092860765e-52 , 128.0*TEST_SQRT_TOL0, GSL::SUCCESS)
189
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "(-0.9999847412109375)", -5.8917472200477175158028143531855, TEST_TOL4, GSL::SUCCESS)
190
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "(-0.5)", -0.5493061443340548457, TEST_TOL0, GSL::SUCCESS)
191
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "(-1e-10)", -1.000000000000000000e-10, TEST_TOL0, GSL::SUCCESS)
192
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "(0.0)", 0.0, TEST_TOL0, GSL::SUCCESS)
193
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "(1e-10)", 1.000000000000000000e-10, TEST_TOL0, GSL::SUCCESS)
194
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "(0.9999847412109375)", 5.8917472200477175158028143531855, TEST_TOL4, GSL::SUCCESS)
195
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "( 1.0000152587890625)", 5.8917548494422489138325509750429, TEST_TOL4, GSL::SUCCESS)
196
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "( 1.5)", 0.8047189562170501873, TEST_TOL0, GSL::SUCCESS)
197
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "( 9.99)", 0.1004364599660005447, TEST_TOL0, GSL::SUCCESS)
198
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "( 10.0)", 0.1003353477310755806, TEST_TOL0, GSL::SUCCESS)
199
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "( 10.01)", 0.1002344395571710243, TEST_TOL0, GSL::SUCCESS)
200
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "( 100)", 0.010000333353334762015, TEST_TOL0, GSL::SUCCESS)
201
+ TEST_SF(s, "GSL::Sf::legendre_Q0_e", "( 1e10)", 1.000000000000000000e-10, TEST_TOL0, GSL::SUCCESS)
202
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(-0.9999847412109375)", 4.8916573191196772369, TEST_TOL4, GSL::SUCCESS)
203
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(-0.5)", -0.7253469278329725772, TEST_TOL0, GSL::SUCCESS)
204
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(-0.01)", -0.9998999966664666524, TEST_TOL0, GSL::SUCCESS)
205
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(-1e-10)", -0.999999999999999999, TEST_TOL0, GSL::SUCCESS)
206
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(0.0)", -1.0, TEST_TOL0, GSL::SUCCESS)
207
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(1e-10)", -0.999999999999999999, TEST_TOL0, GSL::SUCCESS)
208
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(0.0001)", -0.9999999899999999667, TEST_TOL0, GSL::SUCCESS)
209
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(0.01)", -0.9998999966664666524, TEST_TOL0, GSL::SUCCESS)
210
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(0.5)", -0.7253469278329725772, TEST_TOL0, GSL::SUCCESS)
211
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(0.9999847412109375)", 4.8916573191196772369, TEST_TOL4, GSL::SUCCESS)
212
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "(1.0000152587890625)", 4.8918447504867045145, TEST_TOL4, GSL::SUCCESS)
213
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "( 1.5)", 0.20707843432557528095, TEST_TOL0, GSL::SUCCESS)
214
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "( 9.99)", 3.360235060345441639e-3, TEST_TOL0, GSL::SUCCESS)
215
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "( 10.0)", 3.353477310755806357e-3, TEST_TOL0, GSL::SUCCESS)
216
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "( 10.01)", 3.346739967281953346e-3, TEST_TOL0, GSL::SUCCESS)
217
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "( 100.0)", 3.333533347620158821e-5, TEST_TOL0, GSL::SUCCESS)
218
+ TEST_SF(s, "GSL::Sf::legendre_Q1_e", "( 1e10)", 3.333333333333333333e-21, TEST_TOL0, GSL::SUCCESS)
219
+ TEST_SF(s, "GSL::Sf::legendre_Ql_e", "(10, -0.5)", -0.29165813966586752393, TEST_TOL0, GSL::SUCCESS)
220
+ TEST_SF(s, "GSL::Sf::legendre_Ql_e", "(10, 0.5)", 0.29165813966586752393, TEST_TOL0, GSL::SUCCESS)
221
+ TEST_SF(s, "GSL::Sf::legendre_Ql_e", "(10, 1.5)", 0.000014714232718207477406, TEST_TOL0, GSL::SUCCESS)
222
+ TEST_SF(s, "GSL::Sf::legendre_Ql_e", "(100, -0.5)", -0.09492507395207282096, TEST_TOL1, GSL::SUCCESS)
223
+ TEST_SF(s, "GSL::Sf::legendre_Ql_e", "(100, 0.5)", 0.09492507395207282096, TEST_TOL1, GSL::SUCCESS)
224
+ TEST_SF(s, "GSL::Sf::legendre_Ql_e", "(100, 1.5)", 1.1628163435044121988e-43, TEST_TOL2, GSL::SUCCESS)
225
+ TEST_SF(s, "GSL::Sf::legendre_Ql_e", "(1000, -0.5)", -0.030105074974005303500, TEST_TOL1, GSL::SUCCESS)
226
+ TEST_SF(s, "GSL::Sf::legendre_Ql_e", "(1000, 0.5)", 0.030105074974005303500, TEST_TOL1, GSL::SUCCESS)
227
+ TEST_SF(s, "GSL::Sf::legendre_Ql_e", "(1000, 1.1)", 1.0757258447825356443e-194, TEST_TOL3, GSL::SUCCESS)
@@ -0,0 +1,59 @@
1
+ #!/usr/bin/env ruby
2
+ require("gsl")
3
+ require("../gsl_test2.rb")
4
+ require("./gsl_test_sf.rb")
5
+ include GSL::Test
6
+ include Math
7
+ include GSL::Test::Sf
8
+ s = 0
9
+
10
+ TEST_SNGL = 1.0e-06
11
+
12
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 0.0, 0.0)", 0.7071067811865475, TEST_SNGL, GSL::SUCCESS)
13
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 0.0, GSL::M_PI_2)",
14
+ 0.7071067811865475, TEST_SNGL, GSL::SUCCESS);
15
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 5.0, 0.0)",
16
+ 0.04480018165188902, TEST_SNGL, GSL::SUCCESS);
17
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 5.0, GSL::M_PI_2)",
18
+ 1.334848674698019, TEST_SNGL, GSL::SUCCESS);
19
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 10.0, 0.0)",
20
+ 0.007626517570935782, TEST_SNGL, GSL::SUCCESS);
21
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 10.0, GSL::M_PI_2)",
22
+ 1.468660470712856, TEST_SNGL, GSL::SUCCESS);
23
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 15.0, 0.0)",
24
+ 0.001932508315204592, TEST_SNGL, GSL::SUCCESS);
25
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 15.0, GSL::M_PI_2)",
26
+ 1.550108146686649, TEST_SNGL, GSL::SUCCESS);
27
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 20.0, 0.0)",
28
+ 0.0006037438292242197, TEST_SNGL, GSL::SUCCESS);
29
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 20.0, GSL::M_PI_2)",
30
+ 1.609890857395926, TEST_SNGL, GSL::SUCCESS);
31
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 25.0, 0.0)",
32
+ 0.0002158630184146612, TEST_SNGL, GSL::SUCCESS);
33
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(0, 25.0, GSL::M_PI_2)",
34
+ 1.657510298323475, TEST_SNGL, GSL::SUCCESS);
35
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(1, 0.0, 0.0)",
36
+ 1.00000000, TEST_SNGL, GSL::SUCCESS);
37
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(1, 5.0, 0.0)",
38
+ 0.2565428793223637, TEST_SNGL, GSL::SUCCESS);
39
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(1, 10.0, 0.0)",
40
+ 0.05359874774717657, TEST_SNGL, GSL::SUCCESS);
41
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(1, 15.0, 0.0)",
42
+ 0.01504006645382623, TEST_SNGL, GSL::SUCCESS);
43
+ TEST_SF(s, "GSL::Sf::mathieu_se_e", "(5, 0.0, GSL::M_PI_2)",
44
+ 1.0000000, TEST_SNGL, GSL::SUCCESS);
45
+ TEST_SF(s, "GSL::Sf::mathieu_se_e", "(5, 5.0, GSL::M_PI_2)",
46
+ 0.9060779302023551, TEST_SNGL, GSL::SUCCESS);
47
+ TEST_SF(s, "GSL::Sf::mathieu_se_e", "(5, 10.0, GSL::M_PI_2)",
48
+ 0.8460384335355106, TEST_SNGL, GSL::SUCCESS);
49
+ TEST_SF(s, "GSL::Sf::mathieu_se_e", "(5, 15.0, GSL::M_PI_2)",
50
+ 0.837949340012484, TEST_SNGL, GSL::SUCCESS);
51
+ TEST_SF(s, "GSL::Sf::mathieu_se_e", "(5, 20.0, GSL::M_PI_2)",
52
+ 0.8635431218533667, TEST_SNGL, GSL::SUCCESS);
53
+ TEST_SF(s, "GSL::Sf::mathieu_se_e", "(5, 25.0, GSL::M_PI_2)",
54
+ 0.8992683245108413, TEST_SNGL, GSL::SUCCESS);
55
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(10, 0.0, 0.0)",
56
+ 1.00000000, TEST_SNGL, GSL::SUCCESS);
57
+ TEST_SF(s, "GSL::Sf::mathieu_ce_e", "(10, 0.0, GSL::M_PI_2)",
58
+ -1.00000000, TEST_SNGL, GSL::SUCCESS);
59
+
@@ -0,0 +1,839 @@
1
+ #!/usr/bin/env ruby
2
+ require("gsl")
3
+ require("../gsl_test2.rb")
4
+ require("./gsl_test_sf.rb")
5
+ include GSL::Test
6
+ include Math
7
+ include GSL::Test::Sf
8
+ s = 0
9
+ m = GSL::MODE_DEFAULT
10
+ mode = GSL::MODE_DEFAULT
11
+ GSL::set_error_handler_off
12
+ TEST_SF(s, "GSL::Sf::clausen_e", "(M_PI/20.0)", 0.4478882448133546, TEST_TOL0, GSL::SUCCESS)
13
+ TEST_SF(s, "GSL::Sf::clausen_e", "(M_PI/6.0)", 0.8643791310538927, TEST_TOL0, GSL::SUCCESS)
14
+ TEST_SF(s, "GSL::Sf::clausen_e", "(M_PI/3.0)", 1.0149416064096535, TEST_TOL0, GSL::SUCCESS)
15
+ TEST_SF(s, "GSL::Sf::clausen_e", "( 2.0*M_PI + M_PI/3.0)", 1.0149416064096535, TEST_TOL0, GSL::SUCCESS)
16
+ TEST_SF(s, "GSL::Sf::clausen_e", "(100.0*M_PI + M_PI/3.0)", 1.0149416064096535, TEST_TOL0, GSL::SUCCESS)
17
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(0, 1, 1, 0, 1, -1)", sqrt(1.0/2.0), TEST_TOL0, GSL::SUCCESS)
18
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(1, 1, 2, 1, -1, 0)", sqrt(1.0/6.0), TEST_TOL0, GSL::SUCCESS)
19
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(2, 4, 6, 0, 2, -2)", sqrt(8.0/105.0), TEST_TOL0, GSL::SUCCESS)
20
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(4, 4, 8, 0, 0, 0)", sqrt(2.0/35.0), TEST_TOL0, GSL::SUCCESS)
21
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(4, 4, 8, 2, -2, 0)", 2.0/3.0*sqrt(2.0/35.0), TEST_TOL2, GSL::SUCCESS)
22
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(4, 4, 8, 4, -4, 0)", 1.0/(3.0*sqrt(70.0)), TEST_TOL2, GSL::SUCCESS)
23
+ # TEST_SF(s, "GSL::Sf::coupling_3j_e", "(-1, 1, 2, 1, -1, 0)", GSL::NAN, GSL::NAN, GSL::EDOM)
24
+ # TEST_SF(s, "GSL::Sf::coupling_3j_e", "(1, -1, 2, 1, -1, 0)", GSL::NAN, GSL::NAN, GSL::EDOM)
25
+ # TEST_SF(s, "GSL::Sf::coupling_3j_e", "(1, 1, -2, 1, -1, 0)", GSL::NAN, GSL::NAN, GSL::EDOM)
26
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(1, 1, 2, 2, -1, 0)", 0, 0, GSL::SUCCESS)
27
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(1, 1, 2, 1, -2, 0)", 0, 0, GSL::SUCCESS)
28
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(1, 1, 2, 1, -1, 3)", 0, 0, GSL::SUCCESS)
29
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(1, 1, 3, 1, -1, 0)", 0, 0, GSL::SUCCESS)
30
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(1, 4, 2, 1, -1, 0)", 0, 0, GSL::SUCCESS)
31
+ TEST_SF(s, "GSL::Sf::coupling_3j_e", "(4, 1, 2, 1, -1, 0)", 0, 0, GSL::SUCCESS)
32
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, 2, 4, 2, 2, 2)", 1.0/6.0, TEST_TOL0, GSL::SUCCESS)
33
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(4, 4, 2, 4, 4, 4)", -1.0/10.0, TEST_TOL0, GSL::SUCCESS)
34
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(4, 4, 2, 4, 4, 2)", 1.0/6.0, TEST_TOL0, GSL::SUCCESS)
35
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(4, 4, 2, 2, 2, 2)", -0.5/sqrt(5.0), TEST_TOL0, GSL::SUCCESS)
36
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(4, 4, 4, 2, 2, 2)", sqrt(7.0/3.0)/10.0, TEST_TOL0, GSL::SUCCESS)
37
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(6, 6, 6, 4, 4, 4)", -sqrt(3.0/5.0)/14.0, TEST_TOL0, GSL::SUCCESS)
38
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(6, 6, 6, 4, 4, 2)", -sqrt(3.0/5.0)/7.0, TEST_TOL0, GSL::SUCCESS)
39
+ # TEST_SF(s, "GSL::Sf::coupling_6j_e", "(-2, 2, 4, 2, 2, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
40
+ # TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, -2, 4, 2, 2, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
41
+ # TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, 2, -4, 2, 2, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
42
+ # TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, 2, 4, -2, 2, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
43
+ # TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, 2, 4, 2, -2, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
44
+ # TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, 2, 4, 2, 2, -2)", GSL::NAN, GSL::NAN, GSL::EDOM)
45
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, 2, 4, 2, 2, 7)", 0, 0, GSL::SUCCESS)
46
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, 2, 4, 2, 7, 2)", 0, 0, GSL::SUCCESS)
47
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, 2, 4, 7, 2, 2)", 0, 0, GSL::SUCCESS)
48
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, 2, 7, 2, 2, 2)", 0, 0, GSL::SUCCESS)
49
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(2, 7, 4, 2, 2, 2)", 0, 0, GSL::SUCCESS)
50
+ TEST_SF(s, "GSL::Sf::coupling_6j_e", "(7, 2, 4, 2, 2, 2)", 0, 0, GSL::SUCCESS)
51
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, 3, 2, 1, 1, 2)", -sqrt(1.0/6.0)/10.0, TEST_TOL2, GSL::SUCCESS)
52
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(8, 4, 10, 7, 3, 8, 1, 1, 2)", sqrt(7.0/3.0)/60.0, TEST_TOL2, GSL::SUCCESS)
53
+ # TEST_SF(s, "GSL::Sf::coupling_9j_e", "(-4, 2, 4, 3, 3, 2, 1, 1, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
54
+ # TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, -2, 4, 3, 3, 2, 1, 1, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
55
+ # TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, -4, 3, 3, 2, 1, 1, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
56
+ # TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, -3, 3, 2, 1, 1, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
57
+ # TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, -3, 2, 1, 1, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
58
+ # TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, 3, -2, 1, 1, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
59
+ # TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, 3, 2, -1, 1, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
60
+ # TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, 3, 2, 1, -1, 2)", GSL::NAN, GSL::NAN, GSL::EDOM)
61
+ # TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, 3, 2, 1, 1, -2)", GSL::NAN, GSL::NAN, GSL::EDOM)
62
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(10, 2, 4, 3, 3, 2, 1, 1, 2)", 0, 0, GSL::SUCCESS)
63
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 10, 4, 3, 3, 2, 1, 1, 2)", 0, 0, GSL::SUCCESS)
64
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 10, 3, 3, 2, 1, 1, 2)", 0, 0, GSL::SUCCESS)
65
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 10, 3, 2, 1, 1, 2)", 0, 0, GSL::SUCCESS)
66
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, 10, 2, 1, 1, 2)", 0, 0, GSL::SUCCESS)
67
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, 3, 10, 1, 1, 2)", 0, 0, GSL::SUCCESS)
68
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, 3, 2, 10, 1, 2)", 0, 0, GSL::SUCCESS)
69
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, 3, 2, 1, 10, 2)", 0, 0, GSL::SUCCESS)
70
+ TEST_SF(s, "GSL::Sf::coupling_9j_e", "(4, 2, 4, 3, 3, 2, 1, 1, 10)", 0, 0, GSL::SUCCESS)
71
+ TEST_SF(s, "GSL::Sf::dawson_e", "(1.0e-15)", 1.0e-15, TEST_TOL0, GSL::SUCCESS)
72
+ TEST_SF(s, "GSL::Sf::dawson_e", "(0.5)", 0.4244363835020222959, TEST_TOL0, GSL::SUCCESS)
73
+ TEST_SF(s, "GSL::Sf::dawson_e", "(2.0)", 0.30134038892379196603, TEST_TOL0, GSL::SUCCESS)
74
+ TEST_SF(s, "GSL::Sf::dawson_e", "(1000.0)", 0.0005000002500003750009, TEST_TOL0, GSL::SUCCESS)
75
+ TEST_SF(s, "GSL::Sf::debye_1_e", "(0.1)", 0.975277750004723276, TEST_TOL0, GSL::SUCCESS)
76
+ TEST_SF(s, "GSL::Sf::debye_1_e", "(1.0)", 0.777504634112248239, TEST_TOL0, GSL::SUCCESS)
77
+ TEST_SF(s, "GSL::Sf::debye_1_e", "(10.0)", 0.164443465679946027, TEST_TOL0, GSL::SUCCESS)
78
+ TEST_SF(s, "GSL::Sf::debye_2_e", "(0.1)", 0.967083287045302664, TEST_TOL0, GSL::SUCCESS)
79
+ TEST_SF(s, "GSL::Sf::debye_2_e", "(1.0)", 0.70787847562782924, TEST_TOL0, GSL::SUCCESS)
80
+ TEST_SF(s, "GSL::Sf::debye_2_e", "(10.0)", 0.0479714980201218708, TEST_TOL0, GSL::SUCCESS)
81
+ TEST_SF(s, "GSL::Sf::debye_3_e", "(0.1)", 0.962999940487211048, TEST_TOL0, GSL::SUCCESS)
82
+ TEST_SF(s, "GSL::Sf::debye_3_e", "(1.0)", 0.674415564077814667, TEST_TOL0, GSL::SUCCESS)
83
+ TEST_SF(s, "GSL::Sf::debye_3_e", "(10.0)", 0.0192957656903454886, TEST_TOL0, GSL::SUCCESS)
84
+ TEST_SF(s, "GSL::Sf::debye_4_e", "(0.1)", 0.960555486124335944, TEST_TOL0, GSL::SUCCESS)
85
+ TEST_SF(s, "GSL::Sf::debye_4_e", "(1.0)", 0.654874068886737049, TEST_TOL0, GSL::SUCCESS)
86
+ TEST_SF(s, "GSL::Sf::debye_4_e", "(10.0)", 0.00967367556027115896, TEST_TOL0, GSL::SUCCESS)
87
+ x = 0.2*DBL_MAX
88
+ TEST_SF(s, "GSL::Sf::multiply_e", "(-3.0,2.0)", -6.0, TEST_TOL0, GSL::SUCCESS)
89
+ TEST_SF(s, "GSL::Sf::multiply_e", "(#{x}, 1.0/#{x})", 1.0, TEST_TOL0, GSL::SUCCESS)
90
+ TEST_SF(s, "GSL::Sf::multiply_e", "(#{x}, 0.2)", 0.04*DBL_MAX, TEST_TOL1, GSL::SUCCESS)
91
+ TEST_SF(s, "GSL::Sf::multiply_e", "(#{x}, 4.0)", 0.8*DBL_MAX, TEST_TOL1, GSL::SUCCESS)
92
+ TEST_SF(s, "GSL::Sf::ellint_Kcomp_e", "( 0.99, #{mode})", 3.3566005233611923760, TEST_TOL0, GSL::SUCCESS)
93
+ TEST_SF(s, "GSL::Sf::ellint_Kcomp_e", "( 0.50, #{mode})", 1.6857503548125960429, TEST_TOL0, GSL::SUCCESS)
94
+ TEST_SF(s, "GSL::Sf::ellint_Kcomp_e", "(0.010, #{mode})", 1.5708355989121522360, TEST_TOL0, GSL::SUCCESS)
95
+ TEST_SF(s, "GSL::Sf::ellint_Ecomp_e", "(0.99, #{mode})", 1.0284758090288040010, TEST_TOL0, GSL::SUCCESS)
96
+ TEST_SF(s, "GSL::Sf::ellint_Ecomp_e", "(0.50, #{mode})", 1.4674622093394271555, TEST_TOL0, GSL::SUCCESS)
97
+ TEST_SF(s, "GSL::Sf::ellint_Ecomp_e", "(0.01, #{mode})", 1.5707570561503852873, TEST_TOL0, GSL::SUCCESS)
98
+ TEST_SF(s, "GSL::Sf::ellint_F_e", "(M_PI/3.0, 0.99, #{mode})", 1.3065333392738766762, TEST_TOL0, GSL::SUCCESS)
99
+ TEST_SF(s, "GSL::Sf::ellint_F_e", "(M_PI/3.0, 0.50, #{mode})", 1.0895506700518854093, TEST_TOL0, GSL::SUCCESS)
100
+ TEST_SF(s, "GSL::Sf::ellint_F_e", "(M_PI/3.0, 0.01, #{mode})", 1.0472129063770918952, TEST_TOL0, GSL::SUCCESS)
101
+ TEST_SF(s, "GSL::Sf::ellint_E_e", "(M_PI/3.0, 0.99, #{mode})", 0.8704819220377943536, TEST_TOL0, GSL::SUCCESS)
102
+ TEST_SF(s, "GSL::Sf::ellint_E_e", "(M_PI/3.0, 0.50, #{mode})", 1.0075555551444720293, TEST_TOL0, GSL::SUCCESS)
103
+ TEST_SF(s, "GSL::Sf::ellint_E_e", "(M_PI/3.0, 0.01, #{mode})", 1.0471821963889481104, TEST_TOL0, GSL::SUCCESS)
104
+ TEST_SF(s, "GSL::Sf::ellint_P_e", "(M_PI/3.0, 0.99, 0.5, #{mode})", 1.1288726598764099882, TEST_TOL0, GSL::SUCCESS)
105
+ TEST_SF(s, "GSL::Sf::ellint_P_e", "(M_PI/3.0, 0.50, 0.5, #{mode})", 0.9570574331323584890, TEST_TOL0, GSL::SUCCESS)
106
+ TEST_SF(s, "GSL::Sf::ellint_P_e", "(M_PI/3.0, 0.01, 0.5, #{mode})", 0.9228868127118118465, TEST_TOL0, GSL::SUCCESS)
107
+ TEST_SF(s, "GSL::Sf::ellint_RF_e", "(5.0e-11, 1.0e-10, 1.0, #{mode})", 12.36441982979439, TEST_TOL0, GSL::SUCCESS)
108
+ TEST_SF(s, "GSL::Sf::ellint_RF_e", "(1.0, 2.0, 3.0, #{mode})", 0.7269459354689082, TEST_TOL0, GSL::SUCCESS)
109
+ TEST_SF(s, "GSL::Sf::ellint_RD_e", "(5.0e-11, 1.0e-10, 1.0, #{mode})", 34.0932594919337362, TEST_TOL0, GSL::SUCCESS)
110
+ TEST_SF(s, "GSL::Sf::ellint_RD_e", "(1.0, 2.0, 3.0, #{mode})", 0.2904602810289906, TEST_TOL0, GSL::SUCCESS)
111
+ TEST_SF(s, "GSL::Sf::ellint_RC_e", "(1.0, 2.0, #{mode})", 0.7853981633974482, TEST_TOL0, GSL::SUCCESS)
112
+ TEST_SF(s, "GSL::Sf::ellint_RJ_e", "(2.0, 3.0, 4.0, 5.0, #{mode})", 0.1429757966715675, TEST_TOL0, GSL::SUCCESS)
113
+ TEST_SF(s, "GSL::Sf::erfc_e", "(-10.0)", 2.0, TEST_TOL0, GSL::SUCCESS)
114
+ TEST_SF(s, "GSL::Sf::erfc_e", "(-5.0000002)", 1.9999999999984625433, TEST_TOL0, GSL::SUCCESS)
115
+ TEST_SF(s, "GSL::Sf::erfc_e", "(-5.0)", 1.9999999999984625402, TEST_TOL0, GSL::SUCCESS)
116
+ TEST_SF(s, "GSL::Sf::erfc_e", "(-1.0)", 1.8427007929497148693, TEST_TOL0, GSL::SUCCESS)
117
+ TEST_SF(s, "GSL::Sf::erfc_e", "(-0.5)", 1.5204998778130465377, TEST_TOL0, GSL::SUCCESS)
118
+ TEST_SF(s, "GSL::Sf::erfc_e", "(1.0)", 0.15729920705028513066, TEST_TOL0, GSL::SUCCESS)
119
+ TEST_SF(s, "GSL::Sf::erfc_e", "(3.0)", 0.000022090496998585441373, TEST_TOL1, GSL::SUCCESS)
120
+ TEST_SF(s, "GSL::Sf::erfc_e", "(7.0)", 4.183825607779414399e-23, TEST_TOL2, GSL::SUCCESS)
121
+ TEST_SF(s, "GSL::Sf::erfc_e", "(10.0)", 2.0884875837625447570e-45, TEST_TOL2, GSL::SUCCESS)
122
+ TEST_SF(s, "GSL::Sf::log_erfc_e", "(-1.0)", log(1.842700792949714869), TEST_TOL0, GSL::SUCCESS)
123
+ TEST_SF(s, "GSL::Sf::log_erfc_e", "(-0.1)", 0.106576400586522485015, TEST_TOL0, GSL::SUCCESS)
124
+ TEST_SF(s, "GSL::Sf::log_erfc_e", "(-1e-10)", 1.1283791670318505967e-10, TEST_TOL0, GSL::SUCCESS)
125
+ TEST_SF(s, "GSL::Sf::log_erfc_e", "(0.0)", log(1.0), TEST_TOL0, GSL::SUCCESS)
126
+ TEST_SF(s, "GSL::Sf::log_erfc_e", "(1e-10)", -1.128379167159174551e-10, TEST_TOL0, GSL::SUCCESS)
127
+ TEST_SF(s, "GSL::Sf::log_erfc_e", "(0.001)", -0.0011290158896213548027, TEST_TOL0, GSL::SUCCESS)
128
+ TEST_SF(s, "GSL::Sf::log_erfc_e", "(0.1)", -0.119304973737395598329, TEST_TOL0, GSL::SUCCESS)
129
+ TEST_SF(s, "GSL::Sf::log_erfc_e", "(1.0)", log(0.15729920705028513066), TEST_TOL0, GSL::SUCCESS)
130
+ TEST_SF(s, "GSL::Sf::log_erfc_e", "(10.0)", log(2.0884875837625447570e-45), TEST_TOL0, GSL::SUCCESS)
131
+ TEST_SF(s, "GSL::Sf::erf_e", "(-10.0)", -1.0000000000000000000, TEST_TOL0, GSL::SUCCESS)
132
+ TEST_SF(s, "GSL::Sf::erf_e", "(0.5)", 0.5204998778130465377, TEST_TOL0, GSL::SUCCESS)
133
+ TEST_SF(s, "GSL::Sf::erf_e", "(1.0)", 0.8427007929497148693, TEST_TOL0, GSL::SUCCESS)
134
+ TEST_SF(s, "GSL::Sf::erf_e", "(10.0)", 1.0000000000000000000, TEST_TOL0, GSL::SUCCESS)
135
+ TEST_SF(s, "GSL::Sf::erf_Z_e", "(1.0)", 0.24197072451914334980, TEST_TOL0, GSL::SUCCESS)
136
+ TEST_SF(s, "GSL::Sf::erf_Q_e", "(10.0)", 7.619853024160526066e-24, TEST_TOL2, GSL::SUCCESS)
137
+ TEST_SF(s, "GSL::Sf::hazard_e", "(-20.0)", 5.5209483621597631896e-88, TEST_TOL2, GSL::SUCCESS)
138
+ TEST_SF(s, "GSL::Sf::hazard_e", "(-10.0)", 7.6945986267064193463e-23, TEST_TOL2, GSL::SUCCESS)
139
+ TEST_SF(s, "GSL::Sf::hazard_e", "(-1.0)", 0.28759997093917836123, TEST_TOL0, GSL::SUCCESS)
140
+ TEST_SF(s, "GSL::Sf::hazard_e", "( 0.0)", 0.79788456080286535588, TEST_TOL0, GSL::SUCCESS)
141
+ TEST_SF(s, "GSL::Sf::hazard_e", "( 1.0)", 1.5251352761609812091, TEST_TOL0, GSL::SUCCESS)
142
+ TEST_SF(s, "GSL::Sf::hazard_e", "(10.0)", 10.098093233962511963, TEST_TOL2, GSL::SUCCESS)
143
+ TEST_SF(s, "GSL::Sf::hazard_e", "(20.0)", 20.049753068527850542, TEST_TOL2, GSL::SUCCESS)
144
+ TEST_SF(s, "GSL::Sf::hazard_e", "(30.0)", 30.033259667433677037, TEST_TOL2, GSL::SUCCESS)
145
+ TEST_SF(s, "GSL::Sf::hazard_e", "(50.0)", 50.019984031905639809, TEST_TOL0, GSL::SUCCESS)
146
+ TEST_SF(s, "GSL::Sf::hazard_e", "(80.0)", 80.012496096798234468, TEST_TOL0, GSL::SUCCESS)
147
+ TEST_SF(s, "GSL::Sf::hazard_e", "(150.0)", 150.00666607420571802, TEST_TOL0, GSL::SUCCESS)
148
+ TEST_SF(s, "GSL::Sf::hazard_e", "(300.0)", 300.00333325926337415, TEST_TOL0, GSL::SUCCESS)
149
+ TEST_SF(s, "GSL::Sf::hazard_e", "(900.0)", 900.00111110836764382, TEST_TOL0, GSL::SUCCESS)
150
+ TEST_SF(s, "GSL::Sf::hazard_e", "(1001.0)", 1001.0009989990049990, TEST_TOL0, GSL::SUCCESS)
151
+ TEST_SF(s, "GSL::Sf::hazard_e", "(2000.0)", 2000.0004999997500003, TEST_TOL0, GSL::SUCCESS)
152
+ TEST_SF(s, "GSL::Sf::exp_e", "(-10.0)", exp(-10.0), TEST_TOL0, GSL::SUCCESS)
153
+ TEST_SF(s, "GSL::Sf::exp_e", "( 10.0)", exp( 10.0), TEST_TOL0, GSL::SUCCESS)
154
+ TEST_SF(s, "GSL::Sf::exp_err_e", "(-10.0, TEST_TOL1)", exp(-10.0), TEST_TOL1, GSL::SUCCESS)
155
+ TEST_SF(s, "GSL::Sf::exp_err_e", "( 10.0, TEST_TOL1)", exp( 10.0), TEST_TOL1, GSL::SUCCESS)
156
+ x = 0
157
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "(-10.0, 1.0e-06)", 1.0e-06*exp(-10.0), TEST_TOL0, GSL::SUCCESS)
158
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "(-10.0, 2.0)", 2.0*exp(-10.0), TEST_TOL0, GSL::SUCCESS)
159
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "(-10.0, -2.0)", -2.0*exp(-10.0), TEST_TOL0, GSL::SUCCESS)
160
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "( 10.0, 1.0e-06)", 1.0e-06*exp( 10.0), TEST_TOL0, GSL::SUCCESS)
161
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "( 10.0, -2.0)", -2.0*exp( 10.0), TEST_TOL0, GSL::SUCCESS)
162
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "(#{x}, 1.00001)", 1.00001*exp(x), TEST_TOL3, GSL::SUCCESS)
163
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "(#{x}, 1.000001)", 1.000001*exp(x), TEST_TOL3, GSL::SUCCESS)
164
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "(#{x}, 1.000000001)", 1.000000001*exp(x), TEST_TOL3, GSL::SUCCESS)
165
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "(#{x}, 100.0)", 100.0*exp(x), TEST_TOL3, GSL::SUCCESS)
166
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "(#{x}, 1.0e+20)", 1.0e+20*exp(x), TEST_TOL3, GSL::SUCCESS)
167
+ TEST_SF(s, "GSL::Sf::exp_mult_e", "(#{x}, exp(-#{x})*exp(M_LN2))", 2.0, TEST_TOL4, GSL::SUCCESS )
168
+ TEST_SF(s, "GSL::Sf::exp_mult_err_e", "(-10.0, TEST_SQRT_TOL0, 2.0, TEST_SQRT_TOL0)", 2.0*exp(-10.0), TEST_SQRT_TOL0, GSL::SUCCESS)
169
+ TEST_SF(s, "GSL::Sf::exp_mult_err_e", "(#{x}, TEST_SQRT_TOL0*#{x}, exp(-#{x})*exp(M_LN2), TEST_SQRT_TOL0*#{x})", 2.0, TEST_SQRT_TOL0, GSL::SUCCESS )
170
+ TEST_SF(s, "GSL::Sf::expm1_e", "(-10.0)", exp(-10.0)-1.0, TEST_TOL0, GSL::SUCCESS)
171
+ TEST_SF(s, "GSL::Sf::expm1_e", "(-0.001)", -0.00099950016662500845, TEST_TOL0, GSL::SUCCESS)
172
+ TEST_SF(s, "GSL::Sf::expm1_e", "(-1.0e-8)", -1.0e-08 + 0.5e-16, TEST_TOL0, GSL::SUCCESS)
173
+ TEST_SF(s, "GSL::Sf::expm1_e", "( 1.0e-8)", 1.0e-08 + 0.5e-16, TEST_TOL0, GSL::SUCCESS)
174
+ TEST_SF(s, "GSL::Sf::expm1_e", "( 0.001)", 0.0010005001667083417, TEST_TOL0, GSL::SUCCESS)
175
+ TEST_SF(s, "GSL::Sf::expm1_e", "( 10.0)", exp(10.0)-1.0, TEST_TOL0, GSL::SUCCESS)
176
+ TEST_SF(s, "GSL::Sf::exprel_e", "(-10.0)", 0.0999954600070237515, TEST_TOL0, GSL::SUCCESS)
177
+ TEST_SF(s, "GSL::Sf::exprel_e", "(-0.001)", 0.9995001666250084, TEST_TOL0, GSL::SUCCESS)
178
+ TEST_SF(s, "GSL::Sf::exprel_e", "(-1.0e-8)", 1.0 - 0.5e-08, TEST_TOL0, GSL::SUCCESS)
179
+ TEST_SF(s, "GSL::Sf::exprel_e", "( 1.0e-8)", 1.0 + 0.5e-08, TEST_TOL0, GSL::SUCCESS)
180
+ TEST_SF(s, "GSL::Sf::exprel_e", "( 0.001)", 1.0005001667083417, TEST_TOL0, GSL::SUCCESS)
181
+ TEST_SF(s, "GSL::Sf::exprel_e", "( 10.0)", 2202.5465794806716517, TEST_TOL0, GSL::SUCCESS)
182
+ TEST_SF(s, "GSL::Sf::exprel_2_e", "(-10.0)", 0.18000090799859524970, TEST_TOL0, GSL::SUCCESS)
183
+ TEST_SF(s, "GSL::Sf::exprel_2_e", "(-0.001)", 0.9996667499833361107, TEST_TOL0, GSL::SUCCESS)
184
+ TEST_SF(s, "GSL::Sf::exprel_2_e", "(-1.0e-8)", 0.9999999966666666750, TEST_TOL0, GSL::SUCCESS)
185
+ TEST_SF(s, "GSL::Sf::exprel_2_e", "( 1.0e-8)", 1.0000000033333333417, TEST_TOL0, GSL::SUCCESS)
186
+ TEST_SF(s, "GSL::Sf::exprel_2_e", "( 0.001)", 1.0003334166833361115, TEST_TOL0, GSL::SUCCESS)
187
+ TEST_SF(s, "GSL::Sf::exprel_2_e", "( 10.0)", 440.3093158961343303, TEST_TOL0, GSL::SUCCESS)
188
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, -1000.0)", 0.00299400600000000000, TEST_TOL0, GSL::SUCCESS)
189
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, -100.0)", 0.02940600000000000000, TEST_TOL0, GSL::SUCCESS)
190
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, -10.0)", 0.24599972760042142509, TEST_TOL0, GSL::SUCCESS)
191
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, -3.0)", 0.5444917625849191238, TEST_TOL0, GSL::SUCCESS)
192
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, -0.001)", 0.9997500499916678570, TEST_TOL0, GSL::SUCCESS)
193
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, -1.0e-8)", 0.9999999975000000050, TEST_TOL0, GSL::SUCCESS)
194
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, 1.0e-8)", 1.0000000025000000050, TEST_TOL0, GSL::SUCCESS)
195
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, 0.001)", 1.0002500500083345240, TEST_TOL0, GSL::SUCCESS)
196
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, 3.0)", 2.5745637607083706091, TEST_TOL0, GSL::SUCCESS)
197
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, 3.1)", 2.6772417068460206247, TEST_TOL0, GSL::SUCCESS)
198
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, 10.0)", 131.79279476884029910, TEST_TOL1, GSL::SUCCESS)
199
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(3, 100.0)", 1.6128702850896812690e+38, TEST_TOL2, GSL::SUCCESS)
200
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, -1000.0)", 0.04766231609253975959, TEST_TOL0, GSL::SUCCESS)
201
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, -100.0)", 0.3348247572345889317, TEST_TOL0, GSL::SUCCESS)
202
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, -10.0)", 0.8356287051853286482, TEST_TOL0, GSL::SUCCESS)
203
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, -3.0)", 0.9443881609152163615, TEST_TOL0, GSL::SUCCESS)
204
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, -1.0)", 0.980762245565660617, TEST_TOL0, GSL::SUCCESS)
205
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, -1.0e-8)", 1.0 -1.0e-8/51.0, TEST_TOL0, GSL::SUCCESS)
206
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, 1.0e-8)", 1.0 +1.0e-8/51.0, TEST_TOL0, GSL::SUCCESS)
207
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, 1.0)", 1.01999216583666790, TEST_TOL0, GSL::SUCCESS)
208
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, 3.0)", 1.0624205757460368307, TEST_TOL0, GSL::SUCCESS)
209
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, 48.0)", 7.499573876877194416, TEST_TOL0, GSL::SUCCESS)
210
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, 50.1)", 9.311803306230992272, TEST_TOL4, GSL::SUCCESS)
211
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, 100.0)", 8.175664432485807634e+07, TEST_TOL4, GSL::SUCCESS)
212
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(50, 500.0)", 4.806352370663185330e+146, TEST_TOL3, GSL::SUCCESS)
213
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, -1000.0)", 0.3334815803127619256, TEST_TOL0, GSL::SUCCESS)
214
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, -100.0)", 0.8335646217536183909, TEST_TOL0, GSL::SUCCESS)
215
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, -10.0)", 0.9804297803131823066, TEST_TOL0, GSL::SUCCESS)
216
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, -3.0)", 0.9940475488850672997, TEST_TOL0, GSL::SUCCESS)
217
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, -1.0)", 0.9980079602383488808, TEST_TOL0, GSL::SUCCESS)
218
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, -1.0e-8)", 1.0 -1.0e-8/501.0, TEST_TOL0, GSL::SUCCESS)
219
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, 1.0e-8)", 1.0 +1.0e-8/501.0, TEST_TOL0, GSL::SUCCESS)
220
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, 1.0)", 1.0019999920160634252, TEST_TOL0, GSL::SUCCESS)
221
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, 3.0)", 1.0060240236632444934, TEST_TOL0, GSL::SUCCESS)
222
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, 48.0)", 1.1059355517981272174, TEST_TOL0, GSL::SUCCESS)
223
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, 100.0)", 1.2492221464878287204, TEST_TOL1, GSL::SUCCESS)
224
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, 500.0)", 28.363019877927630858, TEST_TOL2, GSL::SUCCESS)
225
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, 1000.0)", 2.4037563160335300322e+68, TEST_TOL4, GSL::SUCCESS)
226
+ TEST_SF(s, "GSL::Sf::exprel_n_e", "(500, 1600.0)", 7.899293535320607403e+226, TEST_TOL4, GSL::SUCCESS)
227
+ TEST_SF(s, "GSL::Sf::expint_E1_e", "(-1.0)", -1.8951178163559367555, TEST_TOL0, GSL::SUCCESS)
228
+ TEST_SF(s, "GSL::Sf::expint_E1_e", "(1.0e-10)", 22.448635265138923980, TEST_TOL0, GSL::SUCCESS)
229
+ TEST_SF(s, "GSL::Sf::expint_E1_e", "(1.0e-05)", 10.935719800043695615, TEST_TOL0, GSL::SUCCESS)
230
+ TEST_SF(s, "GSL::Sf::expint_E1_e", "(0.1)", 1.82292395841939066610, TEST_TOL0, GSL::SUCCESS)
231
+ TEST_SF(s, "GSL::Sf::expint_E1_e", "(1.0)", 0.21938393439552027368, TEST_TOL0, GSL::SUCCESS)
232
+ TEST_SF(s, "GSL::Sf::expint_E1_e", "(10.0)", 4.156968929685324277e-06, TEST_TOL1, GSL::SUCCESS)
233
+ TEST_SF(s, "GSL::Sf::expint_E1_e", "(50.0)", 3.783264029550459019e-24, TEST_TOL2, GSL::SUCCESS)
234
+ TEST_SF(s, "GSL::Sf::expint_E1_e", "(300.0)", 1.710384276804510115e-133, TEST_TOL2, GSL::SUCCESS)
235
+ TEST_SF(s, "GSL::Sf::expint_E2_e", "(-1.0)", 0.8231640121031084799, TEST_TOL1, GSL::SUCCESS)
236
+ TEST_SF(s, "GSL::Sf::expint_E2_e", "(1.0/4294967296.0)", 0.9999999947372139168, TEST_TOL0, GSL::SUCCESS)
237
+ TEST_SF(s, "GSL::Sf::expint_E2_e", "(1.0/65536.0)", 0.9998243233207178845, TEST_TOL0, GSL::SUCCESS)
238
+ TEST_SF(s, "GSL::Sf::expint_E2_e", "(0.1)", 0.7225450221940205066, TEST_TOL0, GSL::SUCCESS)
239
+ TEST_SF(s, "GSL::Sf::expint_E2_e", "(1.0)", 0.14849550677592204792, TEST_TOL0, GSL::SUCCESS)
240
+ TEST_SF(s, "GSL::Sf::expint_E2_e", "(10.0)", 3.830240465631608762e-06, TEST_TOL1, GSL::SUCCESS)
241
+ TEST_SF(s, "GSL::Sf::expint_E2_e", "(50.0)", 3.711783318868827367e-24, TEST_TOL2, GSL::SUCCESS)
242
+ TEST_SF(s, "GSL::Sf::expint_E2_e", "(300.0)", 1.7047391998483433998e-133, TEST_TOL2, GSL::SUCCESS)
243
+
244
+ if GSL_VERSION >= "1.9.90"
245
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(1, -1.0)", -1.8951178163559367555, TEST_TOL2, GSL::SUCCESS)
246
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(1, 1.0e-10)", 22.448635265138923980, TEST_TOL2, GSL::SUCCESS)
247
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(1, 1.0e-5)", 10.93571980004369561, TEST_TOL2, GSL::SUCCESS)
248
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(1,0.1)", 1.82292395841939066610, TEST_TOL0, GSL::SUCCESS);
249
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(1,1.0)", 0.21938393439552027368, TEST_TOL0, GSL::SUCCESS);
250
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(1,10.0)", 4.156968929685324277e-06, TEST_TOL1, GSL::SUCCESS);
251
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(1,50.0)", 3.783264029550459019e-24, TEST_TOL2, GSL::SUCCESS);
252
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(1,300.0)", 1.710384276804510115e-133, TEST_TOL2, GSL::SUCCESS);
253
+
254
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(2,-1.0)", 0.8231640121031084799, TEST_TOL1, GSL::SUCCESS);
255
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(2,0.0)", 1.0, TEST_TOL0, GSL::SUCCESS);
256
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(2,1.0/4294967296.0)", 0.9999999947372139168, TEST_TOL0, GSL::SUCCESS);
257
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(2,1.0/65536.0)", 0.9998243233207178845, TEST_TOL0, GSL::SUCCESS);
258
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(2,0.1)", 0.7225450221940205066, TEST_TOL0, GSL::SUCCESS);
259
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(2,1.0)", 0.14849550677592204792, TEST_TOL0, GSL::SUCCESS);
260
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(2,10.0)", 3.830240465631608762e-06, TEST_TOL1, GSL::SUCCESS);
261
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(2,50.0)", 3.711783318868827367e-24, TEST_TOL2, GSL::SUCCESS);
262
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(2,300.0)", 1.7047391998483433998e-133, TEST_TOL2, GSL::SUCCESS);
263
+
264
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(3,0.0)", 0.5, TEST_TOL0, GSL::SUCCESS);
265
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(3,1.0/4294967296.0)", 0.499999999767169356972, TEST_TOL1, GSL::SUCCESS);
266
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(3,1.0/65536.0)", 0.4999847426094515610, TEST_TOL0, GSL::SUCCESS);
267
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(3,0.1)", 0.4162914579082787612543, TEST_TOL0, GSL::SUCCESS);
268
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(3,1.0)", 0.10969196719776013683858, TEST_TOL1, GSL::SUCCESS);
269
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(3,10.0)",0.000003548762553084381959981, TEST_TOL1, GSL::SUCCESS);
270
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(3,50.0)", 3.6429094264752049812e-24, TEST_TOL2, GSL::SUCCESS);
271
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(3,300.0)",1.699131143349179084e-133, TEST_TOL2, GSL::SUCCESS);
272
+
273
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(10,0.0)", 0.111111111111111111, TEST_TOL0, GSL::SUCCESS);
274
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(10,1.0/4294967296.0)", 0.111111111082007280658, TEST_TOL2, GSL::SUCCESS);
275
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(10,1.0/65536.0)", 0.11110920377910896018606, TEST_TOL1, GSL::SUCCESS);
276
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(10,0.1)", 0.099298432000896813567905, TEST_TOL1, GSL::SUCCESS);
277
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(10,1.0)", 0.036393994031416401634164534, TEST_TOL1, GSL::SUCCESS);
278
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(10,10.0)", 0.00000232530265702821081778968, TEST_TOL1, GSL::SUCCESS);
279
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(10,50.0)", 3.223296586749110919572e-24, TEST_TOL2, GSL::SUCCESS);
280
+ TEST_SF(s, "GSL::Sf::expint_En_e", "(10,300.0)", 1.6608815083360041367294736e-133, TEST_TOL2, GSL::SUCCESS);
281
+ end
282
+
283
+ TEST_SF(s, "GSL::Sf::expint_Ei_e", "(-1.0)", -0.21938393439552027368, TEST_TOL0, GSL::SUCCESS)
284
+ TEST_SF(s, "GSL::Sf::expint_Ei_e", "(1.0/4294967296.0)", -21.603494112783886397, TEST_TOL0, GSL::SUCCESS)
285
+ TEST_SF(s, "GSL::Sf::expint_Ei_e", "(1.0)", 1.8951178163559367555, TEST_TOL0, GSL::SUCCESS)
286
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(-10000.0)", -0.00010001000200060024012, TEST_TOL0, GSL::SUCCESS)
287
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(-1000.0)", -0.0010010020060241207251, TEST_TOL0, GSL::SUCCESS)
288
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(-10.0)", -0.11314702047341077803, TEST_TOL0, GSL::SUCCESS)
289
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(-1.0)", -0.69717488323506606877, TEST_TOL0, GSL::SUCCESS)
290
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(1.0e-10)", 22.448635267383787506, TEST_TOL0, GSL::SUCCESS)
291
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(1.0e-05)", 10.935829157788483865, TEST_TOL0, GSL::SUCCESS)
292
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(0.1)", 2.0146425447084516791, TEST_TOL0, GSL::SUCCESS)
293
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(1.0)", 0.59634736232319407434, TEST_TOL0, GSL::SUCCESS)
294
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(10.0)", 0.091563333939788081876, TEST_TOL0, GSL::SUCCESS)
295
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(50.0)", 0.019615109930114870365, TEST_TOL0, GSL::SUCCESS)
296
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(300.0)", 0.0033222955652707070644, TEST_TOL0, GSL::SUCCESS)
297
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(1000.0)", 0.00099900199402388071500, TEST_TOL0, GSL::SUCCESS)
298
+ TEST_SF(s, "GSL::Sf::expint_E1_scaled_e", "(10000.0)", 0.000099990001999400239880, TEST_TOL0, GSL::SUCCESS)
299
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(-10000.0)", -0.00010002000600240120072, TEST_TOL3, GSL::SUCCESS)
300
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(-1000.0)", -0.0010020060241207250807, TEST_TOL3, GSL::SUCCESS)
301
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(-10.0)", -0.13147020473410778034, TEST_TOL1, GSL::SUCCESS)
302
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(-1.0)", 0.30282511676493393123, TEST_TOL1, GSL::SUCCESS)
303
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(1.0/4294967296.0)", 0.99999999497004455927, TEST_TOL0, GSL::SUCCESS)
304
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(1.0/65536.0)", 0.99983957954556245453, TEST_TOL0, GSL::SUCCESS)
305
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(0.1)", 0.79853574552915483209, TEST_TOL0, GSL::SUCCESS)
306
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(1.0)", 0.40365263767680592566, TEST_TOL0, GSL::SUCCESS)
307
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(10.0)", 0.084366660602119181239, TEST_TOL1, GSL::SUCCESS)
308
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(50.0)", 0.019244503494256481735, TEST_TOL2, GSL::SUCCESS)
309
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(300.0)", 0.0033113304187878806691, TEST_TOL0, GSL::SUCCESS)
310
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(1000.0)", 0.00099800597611928500004, TEST_TOL0, GSL::SUCCESS)
311
+ TEST_SF(s, "GSL::Sf::expint_E2_scaled_e", "(10000.0)", 0.000099980005997601199281, TEST_TOL0, GSL::SUCCESS)
312
+ TEST_SF(s, "GSL::Sf::expint_Ei_scaled_e", "(-1000.0)", -0.00099900199402388071500, TEST_TOL0, GSL::SUCCESS)
313
+ TEST_SF(s, "GSL::Sf::expint_Ei_scaled_e", "(-1.0)", -0.59634736232319407434, TEST_TOL0, GSL::SUCCESS)
314
+ TEST_SF(s, "GSL::Sf::expint_Ei_scaled_e", "(1.0/4294967296.0)", -21.603494107753930958, TEST_TOL0, GSL::SUCCESS)
315
+ TEST_SF(s, "GSL::Sf::expint_Ei_scaled_e", "(1.0)", 0.69717488323506606877, TEST_TOL0, GSL::SUCCESS)
316
+ TEST_SF(s, "GSL::Sf::expint_Ei_scaled_e", "(1000.0)", 0.0010010020060241207251, TEST_TOL0, GSL::SUCCESS)
317
+ TEST_SF(s, "GSL::Sf::Shi_e", "(-1.0)", -1.0572508753757285146, TEST_TOL0, GSL::SUCCESS)
318
+ TEST_SF(s, "GSL::Sf::Shi_e", "(1.0/4294967296.0)", 2.3283064365386962891e-10, TEST_TOL0, GSL::SUCCESS)
319
+ TEST_SF(s, "GSL::Sf::Shi_e", "(1.0/65536.0)", 0.00001525878906269737298, TEST_TOL0, GSL::SUCCESS)
320
+ TEST_SF(s, "GSL::Sf::Shi_e", "(0.1)", 0.1000555722250569955, TEST_TOL0, GSL::SUCCESS)
321
+ TEST_SF(s, "GSL::Sf::Shi_e", "(1.0)", 1.0572508753757285146, TEST_TOL0, GSL::SUCCESS)
322
+ TEST_SF(s, "GSL::Sf::Shi_e", "(10.0)", 1246.1144901994233444, TEST_TOL1, GSL::SUCCESS)
323
+ TEST_SF(s, "GSL::Sf::Shi_e", "(50.0)", 5.292818448565845482e+19, TEST_TOL2, GSL::SUCCESS)
324
+ TEST_SF(s, "GSL::Sf::Shi_e", "(300.0)", 3.248241254044332895e+127, TEST_TOL2, GSL::SUCCESS)
325
+ TEST_SF(s, "GSL::Sf::Chi_e", "(-1.0)", 0.8378669409802082409, TEST_TOL0, GSL::SUCCESS)
326
+ TEST_SF(s, "GSL::Sf::Chi_e", "(1.0/4294967296.0)", -21.603494113016717041, TEST_TOL0, GSL::SUCCESS)
327
+ TEST_SF(s, "GSL::Sf::Chi_e", "(1.0/65536.0)", -10.513139223999384429, TEST_TOL0, GSL::SUCCESS)
328
+ TEST_SF(s, "GSL::Sf::Chi_e", "(1.0/8.0)", -1.4983170827635760646, TEST_TOL0, GSL::SUCCESS)
329
+ TEST_SF(s, "GSL::Sf::Chi_e", "(1.0)", 0.8378669409802082409, TEST_TOL0, GSL::SUCCESS)
330
+ TEST_SF(s, "GSL::Sf::Chi_e", "(10.0)", 1246.1144860424544147, TEST_TOL1, GSL::SUCCESS)
331
+ TEST_SF(s, "GSL::Sf::Chi_e", "(50.0)", 5.292818448565845482e+19, TEST_TOL2, GSL::SUCCESS)
332
+ TEST_SF(s, "GSL::Sf::Chi_e", "(300.0)", 3.248241254044332895e+127, TEST_TOL2, GSL::SUCCESS)
333
+ TEST_SF(s, "GSL::Sf::expint_3_e", "(1.0e-10)", 1.0e-10, TEST_TOL0, GSL::SUCCESS)
334
+ TEST_SF(s, "GSL::Sf::expint_3_e", "(1.0e-05)", 9.9999999999999975e-06, TEST_TOL0, GSL::SUCCESS)
335
+ TEST_SF(s, "GSL::Sf::expint_3_e", "(0.1)", 0.09997500714119079665122, TEST_TOL0, GSL::SUCCESS)
336
+ TEST_SF(s, "GSL::Sf::expint_3_e", "(0.5)", 0.48491714311363971332427, TEST_TOL0, GSL::SUCCESS)
337
+ TEST_SF(s, "GSL::Sf::expint_3_e", "(1.0)", 0.80751118213967145285833, TEST_TOL0, GSL::SUCCESS)
338
+ TEST_SF(s, "GSL::Sf::expint_3_e", "(2.0)", 0.89295351429387631138208, TEST_TOL0, GSL::SUCCESS)
339
+ TEST_SF(s, "GSL::Sf::expint_3_e", "(5.0)", 0.89297951156924921121856, TEST_TOL0, GSL::SUCCESS)
340
+ TEST_SF(s, "GSL::Sf::expint_3_e", "(10.0)", 0.89297951156924921121856, TEST_TOL0, GSL::SUCCESS)
341
+ TEST_SF(s, "GSL::Sf::expint_3_e", "(100.0)", 0.89297951156924921121856, TEST_TOL0, GSL::SUCCESS)
342
+ TEST_SF(s, "GSL::Sf::Si_e", "(-1.0)", -0.9460830703671830149, TEST_TOL0, GSL::SUCCESS)
343
+ TEST_SF(s, "GSL::Sf::Si_e", "(1.0e-10)", 1.0e-10, TEST_TOL0, GSL::SUCCESS)
344
+ TEST_SF(s, "GSL::Sf::Si_e", "(1.0e-05)", 9.999999999944444444e-06, TEST_TOL0, GSL::SUCCESS)
345
+ TEST_SF(s, "GSL::Sf::Si_e", "(0.1)", 0.09994446110827695016, TEST_TOL0, GSL::SUCCESS)
346
+ TEST_SF(s, "GSL::Sf::Si_e", "(1.0)", 0.9460830703671830149, TEST_TOL0, GSL::SUCCESS)
347
+ TEST_SF(s, "GSL::Sf::Si_e", "(10.0)", 1.6583475942188740493, TEST_TOL0, GSL::SUCCESS)
348
+ TEST_SF(s, "GSL::Sf::Si_e", "(50.0)", 1.5516170724859358947, TEST_TOL0, GSL::SUCCESS)
349
+ TEST_SF(s, "GSL::Sf::Si_e", "(300.0)", 1.5708810882137495193, TEST_TOL0, GSL::SUCCESS)
350
+ TEST_SF(s, "GSL::Sf::Si_e", "(1.0e+20)", 1.5707963267948966192, TEST_TOL0, GSL::SUCCESS)
351
+ TEST_SF(s, "GSL::Sf::Ci_e", "(1.0/4294967296.0)", -21.603494113016717041, TEST_TOL0, GSL::SUCCESS)
352
+ TEST_SF(s, "GSL::Sf::Ci_e", "(1.0/65536.0)", -10.513139224115799751, TEST_TOL0, GSL::SUCCESS)
353
+ TEST_SF(s, "GSL::Sf::Ci_e", "(1.0/8.0)", -1.5061295845296396649, TEST_TOL0, GSL::SUCCESS)
354
+ TEST_SF(s, "GSL::Sf::Ci_e", "(1.0)", 0.3374039229009681347, TEST_TOL0, GSL::SUCCESS)
355
+ TEST_SF(s, "GSL::Sf::Ci_e", "(10.0)", -0.04545643300445537263, TEST_TOL0, GSL::SUCCESS)
356
+ TEST_SF(s, "GSL::Sf::Ci_e", "(50.0)", -0.005628386324116305440, TEST_TOL0, GSL::SUCCESS)
357
+ TEST_SF(s, "GSL::Sf::Ci_e", "(300.0)", -0.003332199918592111780, TEST_TOL0, GSL::SUCCESS)
358
+ TEST_SF(s, "GSL::Sf::Ci_e", "(65536.0)", 0.000010560248837656279453, TEST_TOL0, GSL::SUCCESS)
359
+ TEST_SF(s, "GSL::Sf::Ci_e", "(4294967296.0)", -1.0756463261957757485e-10, TEST_SQRT_TOL0, GSL::SUCCESS)
360
+ TEST_SF(s, "GSL::Sf::Ci_e", "(1099511627776.0)", -3.689865584710764214e-13, 1024.0*TEST_SQRT_TOL0, GSL::SUCCESS)
361
+ TEST_SF(s, "GSL::Sf::atanint_e", "(1.0e-10)", 1.0e-10, TEST_TOL0, GSL::SUCCESS)
362
+ TEST_SF(s, "GSL::Sf::atanint_e", "(1.0e-05)", 9.99999999988888888889e-06, TEST_TOL0, GSL::SUCCESS)
363
+ TEST_SF(s, "GSL::Sf::atanint_e", "(0.1)", 0.09988928686033618404, TEST_TOL0, GSL::SUCCESS)
364
+ TEST_SF(s, "GSL::Sf::atanint_e", "(1.0)", 0.91596559417721901505, TEST_TOL0, GSL::SUCCESS)
365
+ TEST_SF(s, "GSL::Sf::atanint_e", "(2.0)", 1.57601540344632342236, TEST_TOL0, GSL::SUCCESS)
366
+ TEST_SF(s, "GSL::Sf::atanint_e", "(10.0)", 3.71678149306806859029, TEST_TOL0, GSL::SUCCESS)
367
+ TEST_SF(s, "GSL::Sf::atanint_e", "(50.0)", 6.16499047850274874222, TEST_TOL0, GSL::SUCCESS)
368
+ TEST_SF(s, "GSL::Sf::atanint_e", "(300.0)", 8.96281388924518959990, TEST_TOL0, GSL::SUCCESS)
369
+ TEST_SF(s, "GSL::Sf::atanint_e", "(1.0e+5)", 18.084471031038661920, TEST_TOL0, GSL::SUCCESS)
370
+ TEST_SF(s, "GSL::Sf::fermi_dirac_m1_e", "(-10.0)", 0.00004539786870243439450, TEST_TOL0, GSL::SUCCESS)
371
+ TEST_SF(s, "GSL::Sf::fermi_dirac_m1_e", "( -1.0)", 0.26894142136999512075, TEST_TOL0, GSL::SUCCESS)
372
+ TEST_SF(s, "GSL::Sf::fermi_dirac_m1_e", "( 1.0)", 0.7310585786300048793, TEST_TOL0, GSL::SUCCESS)
373
+ # 2008/Oct/17 YT
374
+ # This test fails in Darwin 9.5.0, gcc4.0.1
375
+ # expected: 0.9999546021312975
376
+ # obtained: 0.9999546021312976 2.217778082901727e-19 2.21788e-19
377
+ # fracdiff: 5.551367143362476e-17
378
+ TEST_SF(s, "GSL::Sf::fermi_dirac_m1_e", "( 10.0)", 0.9999546021312975656, TEST_TOL0, GSL::SUCCESS)
379
+ TEST_SF(s, "GSL::Sf::fermi_dirac_0_e", "(-10.0)", 0.00004539889921686464677, TEST_TOL0, GSL::SUCCESS)
380
+ TEST_SF(s, "GSL::Sf::fermi_dirac_0_e", "( -1.0)", 0.31326168751822283405, TEST_TOL0, GSL::SUCCESS)
381
+ TEST_SF(s, "GSL::Sf::fermi_dirac_0_e", "( 1.0)", 1.3132616875182228340, TEST_TOL0, GSL::SUCCESS)
382
+ TEST_SF(s, "GSL::Sf::fermi_dirac_0_e", "( 10.0)", 10.000045398899216865, TEST_TOL0, GSL::SUCCESS)
383
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "(-10.0)", 0.00004539941448447633524, TEST_TOL0, GSL::SUCCESS)
384
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( -2.0)", 0.13101248471442377127, TEST_TOL0, GSL::SUCCESS)
385
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( -1.0)", 0.3386479964034521798, TEST_TOL0, GSL::SUCCESS)
386
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( -0.4)", 0.5825520806897909028, TEST_TOL0, GSL::SUCCESS)
387
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( 0.4)", 1.1423819861584355337, TEST_TOL0, GSL::SUCCESS)
388
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( 1.0)", 1.8062860704447742567, TEST_TOL0, GSL::SUCCESS)
389
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( 1.5)", 2.5581520872227806402, TEST_TOL0, GSL::SUCCESS)
390
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( 2.5)", 4.689474797599761667, TEST_TOL0, GSL::SUCCESS)
391
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( 10.0)", 51.64488866743374196, TEST_TOL0, GSL::SUCCESS)
392
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( 12.0)", 73.64492792264531092, TEST_TOL0, GSL::SUCCESS)
393
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( 20.0)", 201.64493406478707282, TEST_TOL0, GSL::SUCCESS)
394
+ TEST_SF(s, "GSL::Sf::fermi_dirac_1_e", "( 50.0)", 1251.6449340668482264, TEST_TOL0, GSL::SUCCESS)
395
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "(-10.0)", 0.00004539967212174776662, TEST_TOL0, GSL::SUCCESS)
396
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( -2.0)", 0.13313272938565030508, TEST_TOL0, GSL::SUCCESS)
397
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( -1.0)", 0.3525648792978077590, TEST_TOL0, GSL::SUCCESS)
398
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( -0.4)", 0.6229402647001272120, TEST_TOL0, GSL::SUCCESS)
399
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( 0.4)", 1.2915805581060844533, TEST_TOL0, GSL::SUCCESS)
400
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( 1.0)", 2.1641656128127008622, TEST_TOL0, GSL::SUCCESS)
401
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( 1.5)", 3.247184513920792475, TEST_TOL0, GSL::SUCCESS)
402
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( 2.5)", 6.797764392735056317, TEST_TOL0, GSL::SUCCESS)
403
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( 10.0)", 183.11605273482105278, TEST_TOL0, GSL::SUCCESS)
404
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( 12.0)", 307.73921494638635166, TEST_TOL0, GSL::SUCCESS)
405
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( 20.0)", 1366.2320146723590157, TEST_TOL0, GSL::SUCCESS)
406
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "( 50.0)", 20915.580036675744655, TEST_TOL0, GSL::SUCCESS)
407
+ TEST_SF(s, "GSL::Sf::fermi_dirac_2_e", "(200.0)", 1.3336623201467029786e+06, TEST_TOL0, GSL::SUCCESS)
408
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "(-10.0)", 0.00004539847236080549532, TEST_TOL0, GSL::SUCCESS)
409
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( -2.0)", 0.12366562180120994266, TEST_TOL0, GSL::SUCCESS)
410
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( -1.0)", 0.29402761761145122022, TEST_TOL0, GSL::SUCCESS)
411
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( -0.4)", 0.4631755336886027800, TEST_TOL0, GSL::SUCCESS)
412
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( 0.4)", 0.7654084737661656915, TEST_TOL0, GSL::SUCCESS)
413
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( 1.0)", 1.0270571254743506890, TEST_TOL0, GSL::SUCCESS)
414
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( 1.5)", 1.2493233478527122008, TEST_TOL0, GSL::SUCCESS)
415
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( 2.5)", 1.6663128834358313625, TEST_TOL0, GSL::SUCCESS)
416
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( 10.0)", 3.552779239536617160, TEST_TOL0, GSL::SUCCESS)
417
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( 12.0)", 3.897268231925439359, TEST_TOL0, GSL::SUCCESS)
418
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( 20.0)", 5.041018507535328603, TEST_TOL0, GSL::SUCCESS)
419
+ TEST_SF(s, "GSL::Sf::fermi_dirac_mhalf_e", "( 50.0)", 7.977530858581869960, TEST_TOL1, GSL::SUCCESS)
420
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "(-10.0)", 0.00004539920105264132755, TEST_TOL1, GSL::SUCCESS)
421
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( -2.0)", 0.12929851332007559106, TEST_TOL0, GSL::SUCCESS)
422
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( -1.0)", 0.3277951592607115477, TEST_TOL0, GSL::SUCCESS)
423
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( -0.4)", 0.5522452153690688947, TEST_TOL0, GSL::SUCCESS)
424
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( 0.4)", 1.0386797503389389277, TEST_TOL0, GSL::SUCCESS)
425
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( 1.0)", 1.5756407761513002308, TEST_TOL0, GSL::SUCCESS)
426
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( 1.5)", 2.1448608775831140360, TEST_TOL0, GSL::SUCCESS)
427
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( 2.5)", 3.606975377950373251, TEST_TOL0, GSL::SUCCESS)
428
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( 10.0)", 24.084656964637653615, TEST_TOL0, GSL::SUCCESS)
429
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( 12.0)", 31.540203287044242593, TEST_TOL0, GSL::SUCCESS)
430
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( 20.0)", 67.49151222165892049, TEST_TOL0, GSL::SUCCESS)
431
+ TEST_SF(s, "GSL::Sf::fermi_dirac_half_e", "( 50.0)", 266.09281252136259343, TEST_TOL1, GSL::SUCCESS)
432
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "(-10.0)", 0.00004539956540456176333, TEST_TOL0, GSL::SUCCESS)
433
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( -2.0)", 0.13224678225177236685, TEST_TOL0, GSL::SUCCESS)
434
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( -1.0)", 0.3466747947990574170, TEST_TOL0, GSL::SUCCESS)
435
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( -0.4)", 0.6056120213305040910, TEST_TOL0, GSL::SUCCESS)
436
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( 0.4)", 1.2258236403963668282, TEST_TOL0, GSL::SUCCESS)
437
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( 1.0)", 2.0022581487784644573, TEST_TOL0, GSL::SUCCESS)
438
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( 1.5)", 2.9277494127932173068, TEST_TOL0, GSL::SUCCESS)
439
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( 2.5)", 5.768879312210516582, TEST_TOL0, GSL::SUCCESS)
440
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( 10.0)", 101.00510084332600020, TEST_TOL2, GSL::SUCCESS)
441
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( 12.0)", 156.51518642795728036, TEST_TOL1, GSL::SUCCESS)
442
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( 20.0)", 546.5630100657601959, TEST_TOL1, GSL::SUCCESS)
443
+ TEST_SF(s, "GSL::Sf::fermi_dirac_3half_e", "( 50.0)", 5332.353566687145552, TEST_TOL1, GSL::SUCCESS)
444
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(3, -2.0)", 0.1342199155038680215, TEST_TOL0, GSL::SUCCESS)
445
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(3, 0.0)", 0.9470328294972459176, TEST_TOL0, GSL::SUCCESS)
446
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(3, 0.1)", 1.0414170610956165759, TEST_TOL0, GSL::SUCCESS)
447
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(3, 1.0)", 2.3982260822489407070, TEST_TOL0, GSL::SUCCESS)
448
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(3, 3.0)", 12.621635313399690724, TEST_TOL1, GSL::SUCCESS)
449
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(3, 100.0)", 4.174893231066566793e+06, TEST_TOL1, GSL::SUCCESS)
450
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(3, 500.0)", 2.604372285319088354e+09, TEST_TOL1, GSL::SUCCESS)
451
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(5, -2.0)", 0.13505242246823676478, TEST_TOL0, GSL::SUCCESS)
452
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(5, 0.0)", 0.9855510912974351041, TEST_TOL0, GSL::SUCCESS)
453
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(5, 0.1)", 1.0876519750101492782, TEST_TOL0, GSL::SUCCESS)
454
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(5, 1.0)", 2.6222337848692390539, TEST_TOL0, GSL::SUCCESS)
455
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(5, 3.0)", 17.008801618012113022, TEST_TOL1, GSL::SUCCESS)
456
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(5, 100.0)", 1.3957522531334869874e+09, TEST_TOL1, GSL::SUCCESS)
457
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(5, 500.0)", 2.1705672808114817955e+13, TEST_TOL2, GSL::SUCCESS)
458
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(7, -2.0)", 0.1352641105671255851, TEST_TOL0, GSL::SUCCESS)
459
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(7, 0.0)", 0.9962330018526478992, TEST_TOL0, GSL::SUCCESS)
460
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(7, 0.1)", 1.1005861815180315485, TEST_TOL0, GSL::SUCCESS)
461
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(7, 1.0)", 2.6918878172003129203, TEST_TOL0, GSL::SUCCESS)
462
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(7, 3.0)", 19.033338976999367642, TEST_TOL2, GSL::SUCCESS)
463
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(7, 10.0)", 5654.530932873610014, TEST_TOL1, GSL::SUCCESS)
464
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(7, 50.0)", 1.005005069985066278e+09, TEST_TOL2, GSL::SUCCESS)
465
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(7, 500.0)", 9.691690268341569514e+16, TEST_TOL3, GSL::SUCCESS)
466
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(9, -2.0)", 0.1353174385330242691, TEST_TOL0, GSL::SUCCESS)
467
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(9, 0.0)", 0.9990395075982715656, TEST_TOL0, GSL::SUCCESS)
468
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(9, 0.1)", 1.1039997234712941212, TEST_TOL0, GSL::SUCCESS)
469
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(9, 1.0)", 2.7113648898129249947, TEST_TOL0, GSL::SUCCESS)
470
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(9, 3.0)", 19.768544008138602223, TEST_TOL2, GSL::SUCCESS)
471
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(9, 10.0)", 10388.990167312912478, TEST_TOL2, GSL::SUCCESS)
472
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(9, 50.0)", 2.85466960802601649e+10, TEST_TOL1, GSL::SUCCESS)
473
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(9, 500.0)", 2.69273849842695876e+20, 2*TEST_TOL1, GSL::SUCCESS)
474
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(10, -2.0)", 0.13532635396712288092, TEST_TOL0, GSL::SUCCESS)
475
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(10, 0.0)", 0.9995171434980607541, TEST_TOL0, GSL::SUCCESS)
476
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(10, 0.1)", 1.1045818238852612296, TEST_TOL0, GSL::SUCCESS)
477
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(10, 1.0)", 2.7147765350346120647, TEST_TOL0, GSL::SUCCESS)
478
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(10, 3.0)", 19.917151938411675171, TEST_TOL1, GSL::SUCCESS)
479
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(10, 10.0)", 12790.918595516495955, TEST_TOL2, GSL::SUCCESS)
480
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(10, 50.0)", 1.3147703201869657654e+11, TEST_TOL2, GSL::SUCCESS)
481
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(10, 500.0)", 1.2241331244469204398e+22, TEST_TOL2, GSL::SUCCESS)
482
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(11, -2.0)", 0.1353308162894847149, TEST_TOL0, GSL::SUCCESS)
483
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(11, 0.0)", 0.9997576851438581909, TEST_TOL0, GSL::SUCCESS)
484
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(11, 0.1)", 1.1048751811565850418, TEST_TOL0, GSL::SUCCESS)
485
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(11, 1.0)", 2.7165128749007313436, TEST_TOL0, GSL::SUCCESS)
486
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(11, 3.0)", 19.997483022044603065, TEST_TOL2, GSL::SUCCESS)
487
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(11, 10.0)", 14987.996005901818036, TEST_TOL2, GSL::SUCCESS)
488
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(11, 50.0)", 5.558322924078990628e+11, TEST_TOL2, GSL::SUCCESS)
489
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(11, 500.0)", 5.101293089606198280e+23, TEST_TOL2, GSL::SUCCESS)
490
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(20, -2.0)", 0.13533527450327238373, TEST_TOL0, GSL::SUCCESS)
491
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(20, 0.0)", 0.9999995232582155428, TEST_TOL0, GSL::SUCCESS)
492
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(20, 0.1)", 1.1051703357941368203, TEST_TOL0, GSL::SUCCESS)
493
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(20, 1.0)", 2.7182783069905721654, TEST_TOL0, GSL::SUCCESS)
494
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(20, 3.0)", 20.085345296028242734, TEST_TOL2, GSL::SUCCESS)
495
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(20, 10.0)", 21898.072920149606475, TEST_TOL2, GSL::SUCCESS)
496
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(20, 50.0)", 1.236873256595717618e+16, TEST_TOL2, GSL::SUCCESS)
497
+ TEST_SF(s, "GSL::Sf::fermi_dirac_int_e", "(20, 500.0)", 9.358938204369557277e+36, TEST_TOL2, GSL::SUCCESS)
498
+ TEST_SF(s, "GSL::Sf::gegenpoly_1_e", "(-0.2, 1.0)", -0.4, TEST_TOL0, GSL::SUCCESS)
499
+ TEST_SF(s, "GSL::Sf::gegenpoly_1_e", "( 0.0, 1.0)", 2.0, TEST_TOL0, GSL::SUCCESS)
500
+ TEST_SF(s, "GSL::Sf::gegenpoly_1_e", "( 1.0, 1.0)", 2.0, TEST_TOL0, GSL::SUCCESS)
501
+ TEST_SF(s, "GSL::Sf::gegenpoly_1_e", "( 1.0, 0.5)", 1.0, TEST_TOL0, GSL::SUCCESS)
502
+ TEST_SF(s, "GSL::Sf::gegenpoly_1_e", "( 5.0, 1.0)", 10.0, TEST_TOL0, GSL::SUCCESS)
503
+ TEST_SF(s, "GSL::Sf::gegenpoly_1_e", "( 100.0, 0.5)", 100.0, TEST_TOL0, GSL::SUCCESS)
504
+ TEST_SF(s, "GSL::Sf::gegenpoly_2_e", "(-0.2, 0.5)", 0.12, TEST_TOL0, GSL::SUCCESS)
505
+ TEST_SF(s, "GSL::Sf::gegenpoly_2_e", "( 0.0, 1.0)", 1.00, TEST_TOL0, GSL::SUCCESS)
506
+ TEST_SF(s, "GSL::Sf::gegenpoly_2_e", "( 1.0, 1.0)", 3.00, TEST_TOL0, GSL::SUCCESS)
507
+ TEST_SF(s, "GSL::Sf::gegenpoly_2_e", "( 1.0, 0.1)", -0.96, TEST_TOL0, GSL::SUCCESS)
508
+ TEST_SF(s, "GSL::Sf::gegenpoly_2_e", "( 5.0, 1.0)", 55.0, TEST_TOL0, GSL::SUCCESS)
509
+ TEST_SF(s, "GSL::Sf::gegenpoly_2_e", "( 100.0, 0.5)", 4950.0, TEST_TOL0, GSL::SUCCESS)
510
+ TEST_SF(s, "GSL::Sf::gegenpoly_3_e", "(-0.2, 0.5)", 0.112, TEST_TOL0, GSL::SUCCESS)
511
+ TEST_SF(s, "GSL::Sf::gegenpoly_3_e", "( 0.0, 1.0)", -2.0/3.0, TEST_TOL0, GSL::SUCCESS)
512
+ TEST_SF(s, "GSL::Sf::gegenpoly_3_e", "( 1.0, 1.0)", 4.000, TEST_TOL0, GSL::SUCCESS)
513
+ TEST_SF(s, "GSL::Sf::gegenpoly_3_e", "( 1.0, 0.1)", -0.392, TEST_TOL0, GSL::SUCCESS)
514
+ TEST_SF(s, "GSL::Sf::gegenpoly_3_e", "( 5.0, 1.0)", 220.000, TEST_TOL0, GSL::SUCCESS)
515
+ TEST_SF(s, "GSL::Sf::gegenpoly_3_e", "( 100.0, 0.5)", 161600.000, TEST_TOL0, GSL::SUCCESS)
516
+ TEST_SF(s, "GSL::Sf::gegenpoly_n_e", "(1, 1.0, 1.0)", 2.000 , TEST_TOL0, GSL::SUCCESS)
517
+ TEST_SF(s, "GSL::Sf::gegenpoly_n_e", "(10, 1.0, 1.0)", 11.000 , TEST_TOL0, GSL::SUCCESS)
518
+ TEST_SF(s, "GSL::Sf::gegenpoly_n_e", "(10, 1.0, 0.1)", -0.4542309376 , TEST_TOL0, GSL::SUCCESS)
519
+ TEST_SF(s, "GSL::Sf::gegenpoly_n_e", "(10, 5.0, 1.0)", 9.23780e+4 , TEST_TOL0, GSL::SUCCESS)
520
+ TEST_SF(s, "GSL::Sf::gegenpoly_n_e", "(10, 100.0, 0.5)", 1.5729338392690000e+13, TEST_TOL0, GSL::SUCCESS)
521
+ TEST_SF(s, "GSL::Sf::gegenpoly_n_e", "(1000, 100.0, 1.0)", 3.3353666135627322e+232, TEST_TOL1, GSL::SUCCESS)
522
+ TEST_SF(s, "GSL::Sf::gegenpoly_n_e", "(100, 2000.0, 1.0)", 5.8753432034937579e+202, TEST_TOL0, GSL::SUCCESS)
523
+ TEST_SF(s, "GSL::Sf::gegenpoly_n_e", "(103, 207.0, 2.0)", 1.4210272202235983e+145, TEST_TOL0, GSL::SUCCESS)
524
+ TEST_SF(s, "GSL::Sf::gegenpoly_n_e", "(103, -0.4, 0.3)", -1.64527498094522e-04, TEST_TOL1, GSL::SUCCESS)
525
+ TEST_SF(s, "GSL::Sf::laguerre_1_e", "(0.5, -1.0)", 2.5, TEST_TOL0, GSL::SUCCESS)
526
+ TEST_SF(s, "GSL::Sf::laguerre_1_e", "(0.5, 1.0)", 0.5, TEST_TOL0, GSL::SUCCESS)
527
+ TEST_SF(s, "GSL::Sf::laguerre_1_e", "(1.0, 1.0)", 1.0, TEST_TOL0, GSL::SUCCESS)
528
+ TEST_SF(s, "GSL::Sf::laguerre_2_e", "( 0.5, -1.0)", 4.875, TEST_TOL0, GSL::SUCCESS)
529
+ TEST_SF(s, "GSL::Sf::laguerre_2_e", "( 0.5, 1.0)", -0.125, TEST_TOL0, GSL::SUCCESS)
530
+ TEST_SF(s, "GSL::Sf::laguerre_2_e", "( 1.0, 1.0)", 0.5, TEST_TOL0, GSL::SUCCESS)
531
+ TEST_SF(s, "GSL::Sf::laguerre_2_e", "(-1.0, 1.0)", -0.5, TEST_TOL0, GSL::SUCCESS)
532
+ TEST_SF(s, "GSL::Sf::laguerre_2_e", "(-2.0, 1.0)", 0.5, TEST_TOL0, GSL::SUCCESS)
533
+ TEST_SF(s, "GSL::Sf::laguerre_2_e", "(-3.0, 1.0)", 2.5, TEST_TOL0, GSL::SUCCESS)
534
+ TEST_SF(s, "GSL::Sf::laguerre_3_e", "(0.5, -1.0)", 8.479166666666666667, TEST_TOL0, GSL::SUCCESS)
535
+ TEST_SF(s, "GSL::Sf::laguerre_3_e", "(0.5, 1.0)", -0.6041666666666666667, TEST_TOL0, GSL::SUCCESS)
536
+ TEST_SF(s, "GSL::Sf::laguerre_3_e", "(1.0, 1.0)", -0.16666666666666666667, TEST_TOL1, GSL::SUCCESS)
537
+ TEST_SF(s, "GSL::Sf::laguerre_3_e", "( 2.0, 1.0)", 2.3333333333333333333, TEST_TOL0, GSL::SUCCESS)
538
+ TEST_SF(s, "GSL::Sf::laguerre_3_e", "(-2.0, 1.0)", 1.0/3.0, TEST_TOL0, GSL::SUCCESS)
539
+ TEST_SF(s, "GSL::Sf::laguerre_3_e", "(-3.0, 1.0)", -1.0/6.0, TEST_TOL0, GSL::SUCCESS)
540
+ TEST_SF(s, "GSL::Sf::laguerre_3_e", "(-4.0, 1.0)", -8.0/3.0, TEST_TOL0, GSL::SUCCESS)
541
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(1, 0.5, 1.0)", 0.5, TEST_TOL0, GSL::SUCCESS)
542
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(2, 1.0, 1.0)", 0.5, TEST_TOL1, GSL::SUCCESS)
543
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(3, 2.0, 1.0)", 2.3333333333333333333, TEST_TOL1, GSL::SUCCESS)
544
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(4, 2.0, 0.5)", 6.752604166666666667, TEST_TOL1, GSL::SUCCESS)
545
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(90, 2.0, 0.5)", -48.79047157201507897, TEST_TOL1, GSL::SUCCESS)
546
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(90, 2.0, -100.0)", 2.5295879275042410902e+63, TEST_TOL2, GSL::SUCCESS)
547
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(90, 2.0, 100.0)", -2.0929042259546928670e+20, TEST_TOL1, GSL::SUCCESS)
548
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, 2.0, -0.5)", 2.2521795545919391405e+07, TEST_TOL2, GSL::SUCCESS)
549
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, 2.0, 0.5)", -28.764832945909097418, TEST_TOL2, GSL::SUCCESS)
550
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(1000, 2.0, -0.5)", 2.4399915170947549589e+21, TEST_TOL3, GSL::SUCCESS)
551
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100000, 2.0, 1.0)", 5107.73491348319, TEST_TOL4, GSL::SUCCESS)
552
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(1, -2.0, 1.0)", -2.0, TEST_TOL0, GSL::SUCCESS)
553
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(2, -2.0, 1.0)", 0.5, TEST_TOL0, GSL::SUCCESS)
554
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(3, -2.0, 1.0)", 1.0/3.0, TEST_TOL0, GSL::SUCCESS)
555
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(10, -2.0, 1.0)", -0.04654954805996472663, TEST_TOL2, GSL::SUCCESS)
556
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(10, -5.0, 1.0)", -0.0031385030864197530864, TEST_TOL2, GSL::SUCCESS)
557
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(10, -9.0, 1.0)", -2.480158730158730159e-06, TEST_TOL5, GSL::SUCCESS)
558
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(10, -11.0, 1.0)", 2.7182818011463844797, TEST_TOL2, GSL::SUCCESS)
559
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(10, -11.0, -1.0)", 0.3678794642857142857, TEST_TOL2, GSL::SUCCESS)
560
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -2.0, 1.0)", -0.0027339992019526273866, TEST_SQRT_TOL0, GSL::SUCCESS)
561
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -2.0, -1.0)", 229923.09193402028290, TEST_TOL5, GSL::SUCCESS)
562
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -10.0, 1.0)", 3.25966665871244092e-11, TEST_TOL6, GSL::SUCCESS)
563
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -10.0, -1.0)", 0.00016484365618205810025, TEST_TOL6, GSL::SUCCESS)
564
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -20.0, 1.0)", 5.09567630343671251e-21, TEST_TOL3, GSL::SUCCESS)
565
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -30.0, 1.0)", 3.46063150272466192e-34, TEST_TOL1, GSL::SUCCESS)
566
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -50.0, 1.0)", 1.20981872933162889e-65, TEST_TOL1, GSL::SUCCESS)
567
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -50.0, -1.0)", 8.60763477742332922e-65, TEST_TOL1, GSL::SUCCESS)
568
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -50.5, 1.0)", 4.84021010426688393e-31, TEST_TOL1, GSL::SUCCESS)
569
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -50.5, -1.0)", 8.49861345212160618e-33, TEST_TOL1, GSL::SUCCESS)
570
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -101.0, 1.0)", 2.7182818284590452354, TEST_TOL1, GSL::SUCCESS)
571
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -101.0, -1.0)", 0.3678794411714423216, TEST_TOL1, GSL::SUCCESS)
572
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -102.0, 1.0)", 271.8281828459045235, TEST_TOL1, GSL::SUCCESS)
573
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -102.0, -1.0)", 37.52370299948711680, TEST_TOL1, GSL::SUCCESS)
574
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -110.0, 1.0)", 1.0666955248998831554e+13, TEST_TOL1, GSL::SUCCESS)
575
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -110.0, -1.0)", 1.7028306108058225871e+12, TEST_TOL1, GSL::SUCCESS)
576
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -200.0, 1.0)", 7.47851889721356628e+58, TEST_TOL1, GSL::SUCCESS)
577
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -200.0, -1.0)", 2.73740299754732273e+58, TEST_TOL1, GSL::SUCCESS)
578
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -50.0, 10.0)", 4.504712811317745591e-21, TEST_SQRT_TOL0, GSL::SUCCESS)
579
+ TEST_SF(s, "GSL::Sf::laguerre_n_e", "(100, -50.0, -10.0)", 1.475165520610679937e-11, TEST_TOL1, GSL::SUCCESS)
580
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(0.0)", 0.0, TEST_TOL0, GSL::SUCCESS)
581
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(1.0)", 0.567143290409783872999969, TEST_TOL0, GSL::SUCCESS)
582
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(2.0)", 0.852605502013725491346472, TEST_TOL0, GSL::SUCCESS)
583
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(20.0)", 2.205003278024059970493066, TEST_TOL0, GSL::SUCCESS)
584
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(1000.0)", 5.24960285240159622712606, TEST_TOL0, GSL::SUCCESS)
585
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(1.0e+6)", 11.38335808614005262200016, TEST_TOL0, GSL::SUCCESS)
586
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(1.0e+12)", 24.43500440493491313826305, TEST_TOL0, GSL::SUCCESS)
587
+ # This fails in Darwin9.5.0, gcc4.0.1
588
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(1.0e+308)", 702.641362034106812081125, TEST_TOL0, GSL::SUCCESS)
589
+ # TEST_SF(s, "GSL::Sf::lambert_W0_e", "(-1.0/M_E - GSL::DBL_EPSILON)", -1.0, TEST_TOL0, GSL::EDOM)
590
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(-1.0/M_E + 1.0/(1024.0*1024.0*1024.0))", -0.999928845560308370714970, TEST_TOL0, GSL::SUCCESS)
591
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(-1.0/M_E + 1.0/(1024.0*1024.0))", -0.997724730359774141620354, TEST_TOL0, GSL::SUCCESS)
592
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(-1.0/M_E + 1.0/512.0)", -0.900335676696088773044678, TEST_TOL0, GSL::SUCCESS)
593
+ TEST_SF(s, "GSL::Sf::lambert_W0_e", "(-1.0/M_E + 0.25)", -0.1349044682661213545487599, TEST_TOL0, GSL::SUCCESS)
594
+ TEST_SF(s, "GSL::Sf::lambert_Wm1_e", "(0.0)", 0.0, TEST_TOL0, GSL::SUCCESS)
595
+ TEST_SF(s, "GSL::Sf::lambert_Wm1_e", "(1.0)", 0.567143290409783872999969, TEST_TOL0, GSL::SUCCESS)
596
+ TEST_SF(s, "GSL::Sf::lambert_Wm1_e", "(2.0)", 0.852605502013725491346472, TEST_TOL0, GSL::SUCCESS)
597
+ TEST_SF(s, "GSL::Sf::lambert_Wm1_e", "(20.0)", 2.205003278024059970493066, TEST_TOL0, GSL::SUCCESS)
598
+ # TEST_SF(s, "GSL::Sf::lambert_Wm1_e", "(-1.0/M_E - GSL::DBL_EPSILON)", -1.0, TEST_TOL0, GSL::EDOM)
599
+ TEST_SF(s, "GSL::Sf::lambert_Wm1_e", "(-1.0/M_E + 1.0/(1024.0*1024.0*1024.0))", -1.000071157815154608049055, TEST_TOL1, GSL::SUCCESS)
600
+ TEST_SF(s, "GSL::Sf::lambert_Wm1_e", "(-1.0/M_E + 1.0/(1024.0*1024.0))", -1.002278726118593023934693, TEST_TOL1, GSL::SUCCESS)
601
+ TEST_SF(s, "GSL::Sf::lambert_Wm1_e", "(-1.0/M_E + 1.0/512.0)", -1.106761200865743124599130, TEST_TOL1, GSL::SUCCESS)
602
+ TEST_SF(s, "GSL::Sf::lambert_Wm1_e", "(-1.0/M_E + 1.0/64.0)", -1.324240940341812125489772, TEST_TOL1, GSL::SUCCESS)
603
+ TEST_SF(s, "GSL::Sf::lambert_Wm1_e", "(-1.0/M_E + 0.25)", -3.345798131120112, TEST_TOL1, GSL::SUCCESS)
604
+ TEST_SF(s, "GSL::Sf::log_e", "(0.1)", -2.3025850929940456840, TEST_TOL0, GSL::SUCCESS)
605
+ TEST_SF(s, "GSL::Sf::log_e", "(1.1)", 0.09531017980432486004, TEST_TOL1, GSL::SUCCESS)
606
+ TEST_SF(s, "GSL::Sf::log_e", "(1000.0)", 6.907755278982137052, TEST_TOL0, GSL::SUCCESS)
607
+ TEST_SF(s, "GSL::Sf::log_abs_e", "(-0.1)", -2.3025850929940456840, TEST_TOL0, GSL::SUCCESS)
608
+ TEST_SF(s, "GSL::Sf::log_abs_e", "(-1.1)", 0.09531017980432486004, TEST_TOL1, GSL::SUCCESS)
609
+ TEST_SF(s, "GSL::Sf::log_abs_e", "(-1000.0)", 6.907755278982137052, TEST_TOL0, GSL::SUCCESS)
610
+ TEST_SF(s, "GSL::Sf::log_abs_e", "(0.1)", -2.3025850929940456840, TEST_TOL0, GSL::SUCCESS)
611
+ TEST_SF(s, "GSL::Sf::log_abs_e", "(1.1)", 0.09531017980432486004, TEST_TOL1, GSL::SUCCESS)
612
+ TEST_SF(s, "GSL::Sf::log_abs_e", "(1000.0)", 6.907755278982137052, TEST_TOL0, GSL::SUCCESS)
613
+ TEST_SF(s, "GSL::Sf::log_1plusx_e", "(1.0e-10)", 9.999999999500000000e-11, TEST_TOL0, GSL::SUCCESS)
614
+ TEST_SF(s, "GSL::Sf::log_1plusx_e", "(1.0e-8)", 9.999999950000000333e-09, TEST_TOL0, GSL::SUCCESS)
615
+ TEST_SF(s, "GSL::Sf::log_1plusx_e", "(1.0e-4)", 0.00009999500033330833533, TEST_TOL0, GSL::SUCCESS)
616
+ TEST_SF(s, "GSL::Sf::log_1plusx_e", "(0.1)", 0.09531017980432486004, TEST_TOL0, GSL::SUCCESS)
617
+ TEST_SF(s, "GSL::Sf::log_1plusx_e", "(0.49)", 0.3987761199573677730, TEST_TOL0, GSL::SUCCESS)
618
+ TEST_SF(s, "GSL::Sf::log_1plusx_e", "(-0.49)", -0.6733445532637655964, TEST_TOL0, GSL::SUCCESS)
619
+ TEST_SF(s, "GSL::Sf::log_1plusx_e", "(1.0)", M_LN2, TEST_TOL0, GSL::SUCCESS)
620
+ TEST_SF(s, "GSL::Sf::log_1plusx_e", "(-0.99)", -4.605170185988091368, TEST_TOL0, GSL::SUCCESS)
621
+ TEST_SF(s, "GSL::Sf::log_1plusx_mx_e", "(1.0e-10)", -4.999999999666666667e-21, TEST_TOL0, GSL::SUCCESS)
622
+ TEST_SF(s, "GSL::Sf::log_1plusx_mx_e", "(1.0e-8)", -4.999999966666666917e-17, TEST_TOL0, GSL::SUCCESS)
623
+ TEST_SF(s, "GSL::Sf::log_1plusx_mx_e", "(1.0e-4)", -4.999666691664666833e-09, TEST_TOL0, GSL::SUCCESS)
624
+ TEST_SF(s, "GSL::Sf::log_1plusx_mx_e", "(0.1)", -0.004689820195675139956, TEST_TOL0, GSL::SUCCESS)
625
+ TEST_SF(s, "GSL::Sf::log_1plusx_mx_e", "(0.49)", -0.09122388004263222704, TEST_TOL0, GSL::SUCCESS)
626
+ TEST_SF(s, "GSL::Sf::log_1plusx_mx_e", "(-0.49)", -0.18334455326376559639, TEST_TOL0, GSL::SUCCESS)
627
+ TEST_SF(s, "GSL::Sf::log_1plusx_mx_e", "(1.0)", M_LN2-1.0, TEST_TOL0, GSL::SUCCESS)
628
+ TEST_SF(s, "GSL::Sf::log_1plusx_mx_e", "(-0.99)", -3.615170185988091368, TEST_TOL0, GSL::SUCCESS)
629
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(2.0, 3)", 8.0, TEST_TOL0, GSL::SUCCESS)
630
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(-2.0, 3)", -8.0, TEST_TOL0, GSL::SUCCESS)
631
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(2.0, -3)", 1.0/8.0, TEST_TOL0, GSL::SUCCESS)
632
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(-2.0, -3)", -1.0/8.0, TEST_TOL0, GSL::SUCCESS)
633
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(10.0, 4)", 1.0e+4, TEST_TOL0, GSL::SUCCESS)
634
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(10.0, -4)", 1.0e-4, TEST_TOL0, GSL::SUCCESS)
635
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(-10.0, 4)", 1.0e+4, TEST_TOL0, GSL::SUCCESS)
636
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(-10.0, -4)", 1.0e-4, TEST_TOL0, GSL::SUCCESS)
637
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(10.0, 40)", 1.0e+40, TEST_TOL0, GSL::SUCCESS)
638
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(8.0, -40)", 7.523163845262640051e-37, TEST_TOL0, GSL::SUCCESS)
639
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(-10.0, 40)", 1.0e+40, TEST_TOL0, GSL::SUCCESS)
640
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(-8.0, -40)", 7.523163845262640051e-37, TEST_TOL0, GSL::SUCCESS)
641
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(10.0, 41)", 1.0e+41, TEST_TOL0, GSL::SUCCESS)
642
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(8.0, -41)", 9.403954806578300064e-38, TEST_TOL0, GSL::SUCCESS)
643
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(-10.0, 41)", -1.0e+41, TEST_TOL0, GSL::SUCCESS)
644
+ TEST_SF(s, "GSL::Sf::pow_int_e", "(-8.0, -41)", -9.403954806578300064e-38, TEST_TOL0, GSL::SUCCESS)
645
+ TEST_SF(s, "GSL::Sf::psi_int_e", "(1)", -0.57721566490153286060, TEST_TOL0, GSL::SUCCESS)
646
+ TEST_SF(s, "GSL::Sf::psi_int_e", "(2)", 0.42278433509846713939, TEST_TOL0, GSL::SUCCESS)
647
+ TEST_SF(s, "GSL::Sf::psi_int_e", "(3)", 0.92278433509846713939, TEST_TOL0, GSL::SUCCESS)
648
+ TEST_SF(s, "GSL::Sf::psi_int_e", "(4)", 1.2561176684318004727, TEST_TOL0, GSL::SUCCESS)
649
+ TEST_SF(s, "GSL::Sf::psi_int_e", "(5)", 1.5061176684318004727, TEST_TOL0, GSL::SUCCESS)
650
+ TEST_SF(s, "GSL::Sf::psi_int_e", "(100)", 4.600161852738087400, TEST_TOL0, GSL::SUCCESS)
651
+ TEST_SF(s, "GSL::Sf::psi_int_e", "(110)", 4.695928024251535633, TEST_TOL0, GSL::SUCCESS)
652
+ TEST_SF(s, "GSL::Sf::psi_int_e", "(5000)", 8.517093188082904107, TEST_TOL0, GSL::SUCCESS)
653
+ TEST_SF(s, "GSL::Sf::psi_e", "(5000.0)", 8.517093188082904107, TEST_TOL0, GSL::SUCCESS)
654
+ TEST_SF(s, "GSL::Sf::psi_e", "(5.0)", 1.5061176684318004727, TEST_TOL0, GSL::SUCCESS)
655
+ TEST_SF(s, "GSL::Sf::psi_e", "(-10.5)", 2.3982391295357816134, TEST_TOL0, GSL::SUCCESS)
656
+ TEST_SF(s, "GSL::Sf::psi_e", "(-100.5)", 4.615124601338064117, TEST_TOL2, GSL::SUCCESS)
657
+ TEST_SF(s, "GSL::Sf::psi_e", "(-1.0e+5-0.5)", 11.512935464924395337, 4.0*TEST_TOL4, GSL::SUCCESS)
658
+ TEST_SF(s, "GSL::Sf::psi_e", "(-262144.0-0.5)", 12.476653064769611581, 4.0*TEST_TOL4, GSL::SUCCESS)
659
+ TEST_SF(s, "GSL::Sf::psi_1piy_e", "(0.8)", -0.07088340212750589223, TEST_TOL1, GSL::SUCCESS)
660
+ TEST_SF(s, "GSL::Sf::psi_1piy_e", "(1.0)", 0.09465032062247697727, TEST_TOL0, GSL::SUCCESS)
661
+ TEST_SF(s, "GSL::Sf::psi_1piy_e", "(5.0)", 1.6127848446157465854, TEST_TOL2, GSL::SUCCESS)
662
+ TEST_SF(s, "GSL::Sf::psi_1piy_e", "(100.0)", 4.605178519404762003, TEST_TOL0, GSL::SUCCESS)
663
+ TEST_SF(s, "GSL::Sf::psi_1piy_e", "(2000.0)", 7.600902480375416216, TEST_TOL0, GSL::SUCCESS)
664
+ TEST_SF(s, "GSL::Sf::psi_1piy_e", "(-0.8)", -0.07088340212750589223, TEST_TOL1, GSL::SUCCESS)
665
+ TEST_SF(s, "GSL::Sf::psi_1piy_e", "(-1.0)", 0.09465032062247697727, TEST_TOL0, GSL::SUCCESS)
666
+ TEST_SF(s, "GSL::Sf::psi_1piy_e", "(-5.0)", 1.6127848446157465854, TEST_TOL2, GSL::SUCCESS)
667
+ TEST_SF(s, "GSL::Sf::psi_1piy_e", "(-100.0)", 4.605178519404762003, TEST_TOL0, GSL::SUCCESS)
668
+ TEST_SF(s, "GSL::Sf::psi_1piy_e", "(-2000.0)", 7.600902480375416216, TEST_TOL0, GSL::SUCCESS)
669
+ TEST_SF(s, "GSL::Sf::psi_1_int_e", "(1)", 1.6449340668482264364, TEST_TOL0, GSL::SUCCESS)
670
+ TEST_SF(s, "GSL::Sf::psi_1_int_e", "(2)", 0.64493406684822643647, TEST_TOL0, GSL::SUCCESS)
671
+ TEST_SF(s, "GSL::Sf::psi_1_int_e", "(3)", 0.39493406684822643647, TEST_TOL0, GSL::SUCCESS)
672
+ TEST_SF(s, "GSL::Sf::psi_1_int_e", "(4)", 0.28382295573711532536, TEST_TOL0, GSL::SUCCESS)
673
+ TEST_SF(s, "GSL::Sf::psi_1_int_e", "(1)", 1.6449340668482264365, TEST_TOL0, GSL::SUCCESS)
674
+ TEST_SF(s, "GSL::Sf::psi_1_int_e", "(5)", 0.22132295573711532536, TEST_TOL0, GSL::SUCCESS)
675
+ TEST_SF(s, "GSL::Sf::psi_1_int_e", "(100)", 0.010050166663333571395, TEST_TOL0, GSL::SUCCESS)
676
+ TEST_SF(s, "GSL::Sf::psi_1_int_e", "(110)", 0.009132356622022545705, TEST_TOL0, GSL::SUCCESS)
677
+ TEST_SF(s, "GSL::Sf::psi_1_int_e", "(500)", 0.0020020013333322666697, TEST_TOL0, GSL::SUCCESS)
678
+ if RB_GSL_VERSION >= "1.5"
679
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(1.0/32.0)", 1025.5728544782377089, TEST_TOL0, GSL::SUCCESS)
680
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(1.0)", 1.6449340668482264365, TEST_TOL0, GSL::SUCCESS)
681
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(5.0)", 0.22132295573711532536, TEST_TOL0, GSL::SUCCESS)
682
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(100.0)", 0.010050166663333571395, TEST_TOL0, GSL::SUCCESS)
683
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(110.0)", 0.009132356622022545705, TEST_TOL0, GSL::SUCCESS)
684
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(500.0)", 0.0020020013333322666697, TEST_TOL0, GSL::SUCCESS)
685
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(-1.0 - 1.0/128.0)", 16386.648472598746587, TEST_TOL0, GSL::SUCCESS)
686
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(-1.50)", 9.3792466449891237539, TEST_TOL0, GSL::SUCCESS)
687
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(-10.5)", 9.7787577398148123845, TEST_TOL0, GSL::SUCCESS)
688
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(-15.5)", 9.8071247184113896201, TEST_TOL0, GSL::SUCCESS)
689
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(-50.5)", 9.8499971860824842274, TEST_TOL0, GSL::SUCCESS)
690
+ TEST_SF(s, "GSL::Sf::psi_1_e", "(-1000.5)", 9.8686054001734414233, TEST_TOL0, GSL::SUCCESS)
691
+ end
692
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(1, 1)", 1.6449340668482264364, TEST_TOL0, GSL::SUCCESS)
693
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(1, 2)", 0.64493406684822643647, TEST_TOL0, GSL::SUCCESS)
694
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(1, 3)", 0.39493406684822643647, TEST_TOL0, GSL::SUCCESS)
695
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(1, 4)", 0.28382295573711532536, TEST_TOL0, GSL::SUCCESS)
696
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(1, 5)", 0.22132295573711532536, TEST_TOL0, GSL::SUCCESS)
697
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(1, 100)", 0.010050166663333571395, TEST_TOL0, GSL::SUCCESS)
698
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(1, 110)", 0.009132356622022545705, TEST_TOL0, GSL::SUCCESS)
699
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(1, 500)", 0.0020020013333322666697, TEST_TOL0, GSL::SUCCESS)
700
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(3, 5.0)", 0.021427828192755075022, TEST_TOL0, GSL::SUCCESS)
701
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(3, 500.0)", 1.6048063999872000683e-08, TEST_TOL0, GSL::SUCCESS)
702
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(10, 5.0)", -0.08675107579196581317, TEST_TOL1, GSL::SUCCESS)
703
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(10, 50.0)", -4.101091112731268288e-12, TEST_TOL0, GSL::SUCCESS)
704
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(0, -1.5)", 0.70315664064524318723, TEST_TOL0, GSL::SUCCESS)
705
+ TEST_SF(s, "GSL::Sf::psi_n_e", "(1, -1.5)", 9.3792466449891237539, TEST_TOL0, GSL::SUCCESS)
706
+ TEST_SF(s, "GSL::Sf::synchrotron_1_e", "(0.01)", 0.444972504114210632, TEST_TOL0, GSL::SUCCESS)
707
+ TEST_SF(s, "GSL::Sf::synchrotron_1_e", "(1.0)", 0.651422815355364504, TEST_TOL1, GSL::SUCCESS)
708
+ TEST_SF(s, "GSL::Sf::synchrotron_1_e", "(10.0)", 0.000192238264300868882, TEST_TOL1, GSL::SUCCESS)
709
+ TEST_SF(s, "GSL::Sf::synchrotron_1_e", "(100.0)", 4.69759366592220221e-43, TEST_TOL1, GSL::SUCCESS)
710
+ TEST_SF(s, "GSL::Sf::synchrotron_2_e", "(0.01)", 0.23098077342226277732, TEST_TOL2, GSL::SUCCESS)
711
+ TEST_SF(s, "GSL::Sf::synchrotron_2_e", "(1.0)", 0.4944750621042082670, TEST_TOL1, GSL::SUCCESS)
712
+ TEST_SF(s, "GSL::Sf::synchrotron_2_e", "(10.0)", 0.00018161187569530204281, TEST_TOL1, GSL::SUCCESS)
713
+ TEST_SF(s, "GSL::Sf::transport_2_e", "(1.0e-10)", 9.9999999999999999999e-11, TEST_TOL0, GSL::SUCCESS)
714
+ TEST_SF(s, "GSL::Sf::transport_2_e", "(1.0)", 0.97303256135517012845, TEST_TOL0, GSL::SUCCESS)
715
+ TEST_SF(s, "GSL::Sf::transport_2_e", "(3.0)", 2.41105004901695346199, TEST_TOL0, GSL::SUCCESS)
716
+ TEST_SF(s, "GSL::Sf::transport_2_e", "(10.0)", 3.28432911449795173575, TEST_TOL0, GSL::SUCCESS)
717
+ TEST_SF(s, "GSL::Sf::transport_2_e", "(100.0)", 3.28986813369645287294, TEST_TOL0, GSL::SUCCESS)
718
+ TEST_SF(s, "GSL::Sf::transport_2_e", "(1.0e+05)", 3.28986813369645287294, TEST_TOL0, GSL::SUCCESS)
719
+ TEST_SF(s, "GSL::Sf::transport_3_e", "(1.0e-10)", 4.999999999999999999997e-21, TEST_TOL0, GSL::SUCCESS)
720
+ TEST_SF(s, "GSL::Sf::transport_3_e", "(1.0)", 0.479841006572417499939, TEST_TOL0, GSL::SUCCESS)
721
+ TEST_SF(s, "GSL::Sf::transport_3_e", "(3.0)", 3.210604662942246772338, TEST_TOL0, GSL::SUCCESS)
722
+ TEST_SF(s, "GSL::Sf::transport_3_e", "(5.0)", 5.614386613842273228585, TEST_TOL0, GSL::SUCCESS)
723
+ TEST_SF(s, "GSL::Sf::transport_3_e", "(10.0)", 7.150322712008592975030, TEST_TOL0, GSL::SUCCESS)
724
+ TEST_SF(s, "GSL::Sf::transport_3_e", "(30.0)", 7.212341416160946511930, TEST_TOL0, GSL::SUCCESS)
725
+ TEST_SF(s, "GSL::Sf::transport_3_e", "(100.0)", 7.212341418957565712398, TEST_TOL0, GSL::SUCCESS)
726
+ TEST_SF(s, "GSL::Sf::transport_3_e", "(1.0e+05)", 7.212341418957565712398, TEST_TOL0, GSL::SUCCESS)
727
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(1.0e-10)", 3.33333333333333333333e-31, TEST_TOL0, GSL::SUCCESS)
728
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(1.0e-07)", 3.33333333333333166666e-22, TEST_TOL0, GSL::SUCCESS)
729
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(1.0e-04)", 3.33333333166666666726e-13, TEST_TOL0, GSL::SUCCESS)
730
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(0.1)", 0.000333166726172109903824, TEST_TOL0, GSL::SUCCESS)
731
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(1.0)", 0.31724404523442648241, TEST_TOL0, GSL::SUCCESS)
732
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(3.0)", 5.96482239737147652446, TEST_TOL0, GSL::SUCCESS)
733
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(5.0)", 15.3597843168821829816, TEST_TOL0, GSL::SUCCESS)
734
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(10.0)", 25.2736676770304417334, TEST_TOL0, GSL::SUCCESS)
735
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(30.0)", 25.9757575220840937469, TEST_TOL0, GSL::SUCCESS)
736
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(100.0)", 25.9757576090673165963, TEST_TOL1, GSL::SUCCESS)
737
+ TEST_SF(s, "GSL::Sf::transport_4_e", "(1.0e+05)", 25.9757576090673165963, TEST_TOL2, GSL::SUCCESS)
738
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(1.0e-10)", 2.49999999999999999999e-41, TEST_TOL0, GSL::SUCCESS)
739
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(1.0e-07)", 2.49999999999999861111e-29, TEST_TOL0, GSL::SUCCESS)
740
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(1.0e-04)", 2.49999999861111111163e-17, TEST_TOL0, GSL::SUCCESS)
741
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(0.1)", 0.000024986116317791487410, TEST_TOL0, GSL::SUCCESS)
742
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(1.0)", 0.236615879239094789259153, TEST_TOL0, GSL::SUCCESS)
743
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(3.0)", 12.77055769104415951115760, TEST_TOL0, GSL::SUCCESS)
744
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(5.0)", 50.26309221817518778543615, TEST_TOL0, GSL::SUCCESS)
745
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(10.0)", 116.3807454024207107698556, TEST_TOL0, GSL::SUCCESS)
746
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(30.0)", 124.4313279083858954839911, TEST_TOL0, GSL::SUCCESS)
747
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(100.0)", 124.4313306172043911597639, TEST_TOL0, GSL::SUCCESS)
748
+ TEST_SF(s, "GSL::Sf::transport_5_e", "(1.0e+05)", 124.43133061720439115976, TEST_TOL0, GSL::SUCCESS)
749
+ TEST_SF(s, "GSL::Sf::sin_e", "(-10.0)", 0.5440211108893698134, TEST_TOL0, GSL::SUCCESS)
750
+ TEST_SF(s, "GSL::Sf::sin_e", "(1.0)", 0.8414709848078965067, TEST_TOL0, GSL::SUCCESS)
751
+ TEST_SF(s, "GSL::Sf::sin_e", "(1000.0)", 0.8268795405320025603, TEST_TOL0, GSL::SUCCESS)
752
+ TEST_SF(s, "GSL::Sf::sin_e", "(1048576.75)", 0.8851545351115651914, TEST_TOL1, GSL::SUCCESS)
753
+ TEST_SF(s, "GSL::Sf::sin_e", "(62831853.75)", 0.6273955953485000827, TEST_TOL3, GSL::SUCCESS)
754
+ TEST_SF(s, "GSL::Sf::sin_e", "(1073741822.5)", -0.8284043541754465988, TEST_SQRT_TOL0, GSL::SUCCESS)
755
+ TEST_SF(s, "GSL::Sf::sin_e", "(1073741824.0)", -0.6173264150460421708, TEST_SQRT_TOL0, GSL::SUCCESS)
756
+ TEST_SF(s, "GSL::Sf::sin_e", "(1073741825.5)", 0.7410684679436226926, TEST_SQRT_TOL0, GSL::SUCCESS)
757
+ TEST_SF(s, "GSL::Sf::sin_e", "(1099511627776.0)", -0.4057050115328287198, 32.0*TEST_SQRT_TOL0, GSL::SUCCESS)
758
+ TEST_SF(s, "GSL::Sf::cos_e", "(-10.0)", -0.8390715290764524523, TEST_TOL0, GSL::SUCCESS)
759
+ TEST_SF(s, "GSL::Sf::cos_e", "(1.0)", 0.5403023058681397174, TEST_TOL0, GSL::SUCCESS)
760
+ TEST_SF(s, "GSL::Sf::cos_e", "(1000.0)", 0.5623790762907029911, TEST_TOL1, GSL::SUCCESS)
761
+ TEST_SF(s, "GSL::Sf::cos_e", "(1048576.75)", 0.4652971620066351799, TEST_TOL2, GSL::SUCCESS)
762
+ TEST_SF(s, "GSL::Sf::cos_e", "(62831853.75)", 0.7787006914966116436, TEST_TOL2, GSL::SUCCESS)
763
+ TEST_SF(s, "GSL::Sf::cos_e", "(1073741822.5)", -0.5601305436977716102, TEST_SQRT_TOL0, GSL::SUCCESS)
764
+ TEST_SF(s, "GSL::Sf::cos_e", "(1073741824.0)", 0.7867071229411881196, TEST_SQRT_TOL0, GSL::SUCCESS)
765
+ TEST_SF(s, "GSL::Sf::cos_e", "(1099511627776.0)", -0.9140040719915570023, 128.0*TEST_SQRT_TOL0, GSL::SUCCESS)
766
+ TEST_SF(s, "GSL::Sf::sinc_e", "(1.0/1024.0)", 0.9999984312693665404, TEST_TOL0, GSL::SUCCESS)
767
+ TEST_SF(s, "GSL::Sf::sinc_e", "(1.0/2.0)", 2.0/M_PI, TEST_TOL0, GSL::SUCCESS)
768
+ TEST_SF(s, "GSL::Sf::sinc_e", "(80.5)", 0.0039541600768172754, TEST_TOL0, GSL::SUCCESS)
769
+ TEST_SF(s, "GSL::Sf::sinc_e", "(100.5)", 0.0031672625490924445, TEST_TOL0, GSL::SUCCESS)
770
+ TEST_SF(s, "GSL::Sf::sinc_e", "(1.0e+06 + 0.5)", 3.18309727028927157e-07, TEST_TOL0, GSL::SUCCESS)
771
+ # TEST_SF(s, "GSL::Sf::sin_pi_x_e", "(1000.5)", 1.0, TEST_TOL0, GSL::SUCCESS)
772
+ # TEST_SF(s, "GSL::Sf::sin_pi_x_e", "(10000.0 + 1.0/65536.0)", 0.00004793689960306688455, TEST_TOL0, GSL::SUCCESS)
773
+ # TEST_SF(s, "GSL::Sf::sin_pi_x_e", "(1099511627776.0 + 1 + 0.125)", -0.3826834323650897717, TEST_TOL0, GSL::SUCCESS)
774
+ TEST_SF(s, "GSL::Sf::lnsinh_e", "(0.1)", -2.3009189815304652235, TEST_TOL0, GSL::SUCCESS)
775
+ TEST_SF(s, "GSL::Sf::lnsinh_e", "(1.0)", 0.16143936157119563361, TEST_TOL0, GSL::SUCCESS)
776
+ TEST_SF(s, "GSL::Sf::lnsinh_e", "(5.0)", 4.306807418479684201, TEST_TOL0, GSL::SUCCESS)
777
+ TEST_SF(s, "GSL::Sf::lnsinh_e", "(100.0)", 99.30685281944005469, TEST_TOL0, GSL::SUCCESS)
778
+ TEST_SF(s, "GSL::Sf::lncosh_e", "(0.125)", 0.007792239318898252791, TEST_TOL0, GSL::SUCCESS)
779
+ TEST_SF(s, "GSL::Sf::lncosh_e", "(1.0)", 0.4337808304830271870, TEST_TOL0, GSL::SUCCESS)
780
+ TEST_SF(s, "GSL::Sf::lncosh_e", "(5.0)", 4.306898218339271555, TEST_TOL0, GSL::SUCCESS)
781
+ TEST_SF(s, "GSL::Sf::lncosh_e", "(100.0)", 99.30685281944005469, TEST_TOL0, GSL::SUCCESS)
782
+ TEST_SF(s, "GSL::Sf::zeta_int_e", "(-61.0)", -3.30660898765775767257e+34, TEST_TOL0, GSL::SUCCESS)
783
+ TEST_SF(s, "GSL::Sf::zeta_int_e", "(-5.0)", -0.003968253968253968253968, TEST_TOL0, GSL::SUCCESS)
784
+ TEST_SF(s, "GSL::Sf::zeta_int_e", "( 5.0)", 1.0369277551433699263313655, TEST_TOL0, GSL::SUCCESS)
785
+ TEST_SF(s, "GSL::Sf::zeta_int_e", "(31.0)", 1.0000000004656629065033784, TEST_TOL0, GSL::SUCCESS)
786
+ TEST_SF(s, "GSL::Sf::zetam1_int_e", "(-61.0)", -3.30660898765775767257e+34, TEST_TOL0, GSL::SUCCESS)
787
+ TEST_SF(s, "GSL::Sf::zetam1_int_e", "(-5.0)", -1.003968253968253968253968, TEST_TOL0, GSL::SUCCESS)
788
+ TEST_SF(s, "GSL::Sf::zetam1_int_e", "( 5.0)", 0.0369277551433699263313655, TEST_TOL0, GSL::SUCCESS)
789
+ TEST_SF(s, "GSL::Sf::zetam1_int_e", "(31.0)", 0.0000000004656629065033784, TEST_TOL0, GSL::SUCCESS)
790
+ TEST_SF(s, "GSL::Sf::zeta_e", "(-151)", 8.195215221831378294e+143, TEST_TOL2, GSL::SUCCESS)
791
+ TEST_SF(s, "GSL::Sf::zeta_e", "(-51)", 9.68995788746359406565e+24, TEST_TOL1, GSL::SUCCESS)
792
+ TEST_SF(s, "GSL::Sf::zeta_e", "(-5)", -0.003968253968253968253968, TEST_TOL1, GSL::SUCCESS)
793
+ TEST_SF(s, "GSL::Sf::zeta_e", "(-0.5)", -0.207886224977354566017307, TEST_TOL1, GSL::SUCCESS)
794
+ TEST_SF(s, "GSL::Sf::zeta_e", "(-1e-10)", -0.49999999990810614668948, TEST_TOL1, GSL::SUCCESS)
795
+ TEST_SF(s, "GSL::Sf::zeta_e", "(0.0)", -0.5, TEST_TOL0, GSL::SUCCESS)
796
+ TEST_SF(s, "GSL::Sf::zeta_e", "(1e-10)", -0.50000000009189385333058, TEST_TOL0, GSL::SUCCESS)
797
+ TEST_SF(s, "GSL::Sf::zeta_e", "(0.5)", -1.460354508809586812889499, TEST_TOL0, GSL::SUCCESS)
798
+ TEST_SF(s, "GSL::Sf::zeta_e", "(1.0-1.0/1024.0)", -1023.4228554489429787, TEST_TOL0, GSL::SUCCESS)
799
+ TEST_SF(s, "GSL::Sf::zeta_e", "(1.0+1.0/1048576)", 1.0485765772157343441e+06, TEST_TOL0, GSL::SUCCESS)
800
+ TEST_SF(s, "GSL::Sf::zeta_e", "(5.0)", 1.036927755143369926331365, TEST_TOL0, GSL::SUCCESS)
801
+ TEST_SF(s, "GSL::Sf::zeta_e", "(25.5)", 1.000000021074106110269959, TEST_TOL0, GSL::SUCCESS)
802
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(0.5)", -2.460354508809586812889499, TEST_TOL0, GSL::SUCCESS)
803
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(2.0)", 0.64493406684822643647, TEST_TOL1, GSL::SUCCESS)
804
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(3.0)", 0.20205690315959428540, TEST_TOL1, GSL::SUCCESS)
805
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(5.0)", 0.0369277551433699263314, TEST_TOL1, GSL::SUCCESS)
806
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(9.5)", 0.0014125906121736622712, TEST_TOL1, GSL::SUCCESS)
807
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(10.5)", 0.000700842641736155219500, TEST_TOL1, GSL::SUCCESS)
808
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(12.5)", 0.000173751733643178193390, TEST_TOL1, GSL::SUCCESS)
809
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(13.5)", 0.000086686727462338155188, TEST_TOL1, GSL::SUCCESS)
810
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(15.5)", 0.000021619904246069108133, TEST_TOL1, GSL::SUCCESS)
811
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(16.5)", 0.000010803124900178547671, TEST_TOL0, GSL::SUCCESS)
812
+ TEST_SF(s, "GSL::Sf::zetam1_e", "(25.5)", 0.000000021074106110269959, TEST_TOL0, GSL::SUCCESS)
813
+ TEST_SF(s, "GSL::Sf::hzeta_e", "(2, 1.0)", 1.6449340668482264365, TEST_TOL0, GSL::SUCCESS)
814
+ TEST_SF(s, "GSL::Sf::hzeta_e", "(2, 10.0)", 0.1051663356816857461, TEST_TOL0, GSL::SUCCESS)
815
+ TEST_SF(s, "GSL::Sf::hzeta_e", "(5, 1.0)", 1.0369277551433699263, TEST_TOL0, GSL::SUCCESS)
816
+ TEST_SF(s, "GSL::Sf::hzeta_e", "(5, 10.0)", 0.000030413798676470276, TEST_TOL0, GSL::SUCCESS)
817
+ TEST_SF(s, "GSL::Sf::hzeta_e", "(9, 0.1)", 1.0000000004253980e+09, TEST_TOL0, GSL::SUCCESS)
818
+ TEST_SF(s, "GSL::Sf::hzeta_e", "(30, 0.5)", 1.0737418240000053e+09, TEST_TOL0, GSL::SUCCESS)
819
+ TEST_SF(s, "GSL::Sf::hzeta_e", "(30, 0.9)", 2.3589824880264765e+01, TEST_TOL1, GSL::SUCCESS)
820
+ TEST_SF(s, "GSL::Sf::hzeta_e", "(75, 0.25)", 1.4272476927059599e+45, TEST_TOL1, GSL::SUCCESS)
821
+ TEST_SF(s, "GSL::Sf::eta_int_e", "(-91)", -4.945598888750002040e+94, TEST_TOL0, GSL::SUCCESS)
822
+ TEST_SF(s, "GSL::Sf::eta_int_e", "(-51)", -4.363969073121683116e+40, TEST_TOL0, GSL::SUCCESS)
823
+ TEST_SF(s, "GSL::Sf::eta_int_e", "(-5)", 0.25, TEST_TOL0, GSL::SUCCESS)
824
+ TEST_SF(s, "GSL::Sf::eta_int_e", "(-1)", 0.25, TEST_TOL0, GSL::SUCCESS)
825
+ TEST_SF(s, "GSL::Sf::eta_int_e", "( 0)", 0.5, TEST_TOL0, GSL::SUCCESS)
826
+ TEST_SF(s, "GSL::Sf::eta_int_e", "( 5)", 0.9721197704469093059, TEST_TOL0, GSL::SUCCESS)
827
+ TEST_SF(s, "GSL::Sf::eta_int_e", "( 6)", 0.9855510912974351041, TEST_TOL0, GSL::SUCCESS)
828
+ TEST_SF(s, "GSL::Sf::eta_int_e", "( 20)", 0.9999990466115815221, TEST_TOL0, GSL::SUCCESS)
829
+ TEST_SF(s, "GSL::Sf::eta_int_e", "( 1000)", 1.0, TEST_TOL0, GSL::SUCCESS)
830
+ TEST_SF(s, "GSL::Sf::eta_e", "(-51.5)", -1.2524184036924703656e+41, TEST_TOL2, GSL::SUCCESS)
831
+ TEST_SF(s, "GSL::Sf::eta_e", "(-5)", 0.25, TEST_TOL0, GSL::SUCCESS)
832
+ TEST_SF(s, "GSL::Sf::eta_e", "(0.5)", 0.6048986434216303702, TEST_TOL0, GSL::SUCCESS)
833
+ TEST_SF(s, "GSL::Sf::eta_e", "(0.999)", 0.6929872789683383574, TEST_TOL0, GSL::SUCCESS)
834
+ TEST_SF(s, "GSL::Sf::eta_e", "(1.0)", 0.6931471805599453094, TEST_TOL0, GSL::SUCCESS)
835
+ TEST_SF(s, "GSL::Sf::eta_e", "(1.0+1.0e-10)", 0.6931471805759321998, TEST_TOL0, GSL::SUCCESS)
836
+ TEST_SF(s, "GSL::Sf::eta_e", "( 5)", 0.9721197704469093059, TEST_TOL0, GSL::SUCCESS)
837
+ TEST_SF(s, "GSL::Sf::eta_e", "( 5.2)", 0.9755278712546684682, TEST_TOL0, GSL::SUCCESS)
838
+ TEST_SF(s, "GSL::Sf::eta_e", "( 6)", 0.9855510912974351041, TEST_TOL0, GSL::SUCCESS)
839
+ TEST_SF(s, "GSL::Sf::eta_e", "( 20)", 0.9999990466115815221, TEST_TOL0, GSL::SUCCESS)