26 December 2023 01:53:36 PM gsl_test(): C version Test gsl(), the GNU Scientific Library. gsl_eigen_nonsymm_test(): gsl_eigen_nonysmm() computes the eigenvalues and eigenvectors of a nonsymmetric matrix. Eigenvalue(0) = -6.41391 + 0i Eigenvector(0) = 0.0998822 + 0i 0.111251 + 0i -0.292501 + 0i -0.944505 + 0i Eigenvalue(1) = 5.54555 + 3.08545i Eigenvector(1) = 0.0430757 + 0.00968662i -0.0709124 + 0.138917i 0.516595 + -0.0160059i 0.839574 + 0.0413888i Eigenvalue(2) = 5.54555 + -3.08545i Eigenvector(2) = 0.0430757 + -0.00968662i -0.0709124 + -0.138917i 0.516595 + 0.0160059i 0.839574 + -0.0413888i Eigenvalue(3) = 2.3228 + 0i Eigenvector(3) = -0.144933 + 0i 0.356601 + 0i 0.919369 + 0i 0.0811836 + 0i gsl_fft_complex_test(): gsl_fft_complex() computes the fast Fourier transform of a complex vector of data. Input data (a symmetric pulse): 0 1.000000e+00 0.000000e+00 1 1.000000e+00 0.000000e+00 2 1.000000e+00 0.000000e+00 3 1.000000e+00 0.000000e+00 4 1.000000e+00 0.000000e+00 5 1.000000e+00 0.000000e+00 6 1.000000e+00 0.000000e+00 7 1.000000e+00 0.000000e+00 8 1.000000e+00 0.000000e+00 9 1.000000e+00 0.000000e+00 10 1.000000e+00 0.000000e+00 11 0.000000e+00 0.000000e+00 12 0.000000e+00 0.000000e+00 13 0.000000e+00 0.000000e+00 14 0.000000e+00 0.000000e+00 15 0.000000e+00 0.000000e+00 16 0.000000e+00 0.000000e+00 17 0.000000e+00 0.000000e+00 18 0.000000e+00 0.000000e+00 19 0.000000e+00 0.000000e+00 20 0.000000e+00 0.000000e+00 21 0.000000e+00 0.000000e+00 22 0.000000e+00 0.000000e+00 23 0.000000e+00 0.000000e+00 24 0.000000e+00 0.000000e+00 25 0.000000e+00 0.000000e+00 26 0.000000e+00 0.000000e+00 27 0.000000e+00 0.000000e+00 28 0.000000e+00 0.000000e+00 29 0.000000e+00 0.000000e+00 30 0.000000e+00 0.000000e+00 31 0.000000e+00 0.000000e+00 32 0.000000e+00 0.000000e+00 33 0.000000e+00 0.000000e+00 34 0.000000e+00 0.000000e+00 35 0.000000e+00 0.000000e+00 36 0.000000e+00 0.000000e+00 37 0.000000e+00 0.000000e+00 38 0.000000e+00 0.000000e+00 39 0.000000e+00 0.000000e+00 40 0.000000e+00 0.000000e+00 41 0.000000e+00 0.000000e+00 42 0.000000e+00 0.000000e+00 43 0.000000e+00 0.000000e+00 44 0.000000e+00 0.000000e+00 45 0.000000e+00 0.000000e+00 46 0.000000e+00 0.000000e+00 47 0.000000e+00 0.000000e+00 48 0.000000e+00 0.000000e+00 49 0.000000e+00 0.000000e+00 50 0.000000e+00 0.000000e+00 51 0.000000e+00 0.000000e+00 52 0.000000e+00 0.000000e+00 53 0.000000e+00 0.000000e+00 54 0.000000e+00 0.000000e+00 55 0.000000e+00 0.000000e+00 56 0.000000e+00 0.000000e+00 57 0.000000e+00 0.000000e+00 58 0.000000e+00 0.000000e+00 59 0.000000e+00 0.000000e+00 60 0.000000e+00 0.000000e+00 61 0.000000e+00 0.000000e+00 62 0.000000e+00 0.000000e+00 63 0.000000e+00 0.000000e+00 64 0.000000e+00 0.000000e+00 65 0.000000e+00 0.000000e+00 66 0.000000e+00 0.000000e+00 67 0.000000e+00 0.000000e+00 68 0.000000e+00 0.000000e+00 69 0.000000e+00 0.000000e+00 70 0.000000e+00 0.000000e+00 71 0.000000e+00 0.000000e+00 72 0.000000e+00 0.000000e+00 73 0.000000e+00 0.000000e+00 74 0.000000e+00 0.000000e+00 75 0.000000e+00 0.000000e+00 76 0.000000e+00 0.000000e+00 77 0.000000e+00 0.000000e+00 78 0.000000e+00 0.000000e+00 79 0.000000e+00 0.000000e+00 80 0.000000e+00 0.000000e+00 81 0.000000e+00 0.000000e+00 82 0.000000e+00 0.000000e+00 83 0.000000e+00 0.000000e+00 84 0.000000e+00 0.000000e+00 85 0.000000e+00 0.000000e+00 86 0.000000e+00 0.000000e+00 87 0.000000e+00 0.000000e+00 88 0.000000e+00 0.000000e+00 89 0.000000e+00 0.000000e+00 90 0.000000e+00 0.000000e+00 91 0.000000e+00 0.000000e+00 92 0.000000e+00 0.000000e+00 93 0.000000e+00 0.000000e+00 94 0.000000e+00 0.000000e+00 95 0.000000e+00 0.000000e+00 96 0.000000e+00 0.000000e+00 97 0.000000e+00 0.000000e+00 98 0.000000e+00 0.000000e+00 99 0.000000e+00 0.000000e+00 100 0.000000e+00 0.000000e+00 101 0.000000e+00 0.000000e+00 102 0.000000e+00 0.000000e+00 103 0.000000e+00 0.000000e+00 104 0.000000e+00 0.000000e+00 105 0.000000e+00 0.000000e+00 106 0.000000e+00 0.000000e+00 107 0.000000e+00 0.000000e+00 108 0.000000e+00 0.000000e+00 109 0.000000e+00 0.000000e+00 110 0.000000e+00 0.000000e+00 111 0.000000e+00 0.000000e+00 112 0.000000e+00 0.000000e+00 113 0.000000e+00 0.000000e+00 114 0.000000e+00 0.000000e+00 115 0.000000e+00 0.000000e+00 116 0.000000e+00 0.000000e+00 117 0.000000e+00 0.000000e+00 118 1.000000e+00 0.000000e+00 119 1.000000e+00 0.000000e+00 120 1.000000e+00 0.000000e+00 121 1.000000e+00 0.000000e+00 122 1.000000e+00 0.000000e+00 123 1.000000e+00 0.000000e+00 124 1.000000e+00 0.000000e+00 125 1.000000e+00 0.000000e+00 126 1.000000e+00 0.000000e+00 127 1.000000e+00 0.000000e+00 Output data: 0 1.856155e+00 0.000000e+00 1 1.775235e+00 1.177569e-16 2 1.545075e+00 0.000000e+00 3 1.201145e+00 2.453269e-17 4 7.952850e-01 -4.906539e-17 5 3.863025e-01 -6.869155e-17 6 2.955767e-02 -4.906539e-18 7 -2.324518e-01 -9.813078e-18 8 -3.767087e-01 2.943923e-17 9 -4.023202e-01 -5.151866e-17 10 -3.288421e-01 -1.226635e-17 11 -1.908223e-01 -3.679904e-17 12 -2.984509e-02 -5.397193e-17 13 1.141877e-01 -1.054906e-16 14 2.107341e-01 -1.496494e-16 15 2.437478e-01 -2.526868e-16 16 2.133883e-01 9.813078e-18 17 1.341899e-01 -3.618572e-17 18 3.033355e-02 -4.906539e-18 19 -7.075117e-02 -5.887847e-17 20 -1.449419e-01 2.453269e-18 21 -1.766837e-01 -1.226635e-17 22 -1.618772e-01 -9.813078e-18 23 -1.079124e-01 1.471962e-17 24 -3.103786e-02 3.925231e-17 25 4.815112e-02 2.698596e-17 26 1.099404e-01 2.759928e-17 27 1.401760e-01 -2.943923e-17 28 1.333280e-01 7.359808e-18 29 9.330968e-02 -1.079439e-16 30 3.198024e-02 -1.557826e-16 31 -3.418842e-02 -1.711155e-16 32 -8.838835e-02 0.000000e+00 33 -1.176209e-01 1.103971e-17 34 -1.157154e-01 5.642520e-17 35 -8.454075e-02 -1.717289e-17 36 -3.319203e-02 3.679904e-17 37 2.456534e-02 -1.226635e-17 38 7.390118e-02 -5.519856e-18 39 1.026526e-01 -1.226635e-18 40 1.042612e-01 9.813078e-18 41 7.921898e-02 -1.717289e-17 42 3.471627e-02 7.359808e-18 43 -1.736727e-02 -2.207943e-17 44 -6.358048e-02 -5.397193e-17 45 -9.232361e-02 -6.133174e-17 46 -9.671760e-02 -7.114481e-17 47 -7.621953e-02 -6.133174e-17 48 -3.661165e-02 0.000000e+00 49 1.159674e-02 -3.066587e-18 50 5.592184e-02 -4.906539e-18 51 8.510068e-02 -3.679904e-18 52 9.192081e-02 -4.906539e-18 53 7.498108e-02 -2.453269e-18 54 3.895842e-02 1.962616e-17 55 -6.664175e-03 -4.906539e-18 56 -5.006798e-02 0.000000e+00 57 -8.013033e-02 -2.207943e-17 58 -8.924785e-02 1.471962e-17 59 -7.524106e-02 9.813078e-18 60 -4.186758e-02 3.925231e-17 61 2.175052e-03 5.887847e-17 62 4.549549e-02 6.869155e-17 63 7.692872e-02 7.359808e-17 64 8.838835e-02 0.000000e+00 65 7.692872e-02 -1.962616e-17 66 4.549549e-02 -1.962616e-17 67 2.175052e-03 4.906539e-18 68 -4.186758e-02 9.813078e-18 69 -7.524106e-02 9.813078e-18 70 -8.924785e-02 1.471962e-17 71 -8.013033e-02 0.000000e+00 72 -5.006798e-02 -9.813078e-18 73 -6.664175e-03 7.359808e-18 74 3.895842e-02 7.359808e-18 75 7.498108e-02 2.453269e-18 76 9.192081e-02 4.906539e-18 77 8.510068e-02 -1.226635e-17 78 5.592184e-02 -2.207943e-17 79 1.159674e-02 -5.642520e-17 80 -3.661165e-02 -9.813078e-18 81 -7.621953e-02 -3.618572e-17 82 -9.671760e-02 -2.453269e-17 83 -9.232361e-02 -2.943923e-17 84 -6.358048e-02 1.717289e-17 85 -1.736727e-02 -2.453269e-18 86 3.471627e-02 9.813078e-18 87 7.921898e-02 -4.906539e-18 88 1.042612e-01 0.000000e+00 89 1.026526e-01 -2.453269e-18 90 7.390118e-02 1.839952e-18 91 2.456534e-02 9.813078e-18 92 -3.319203e-02 -4.661212e-17 93 -8.454075e-02 9.813078e-18 94 -1.157154e-01 3.802568e-17 95 -1.176209e-01 7.298477e-17 96 -8.838835e-02 0.000000e+00 97 -3.418842e-02 -6.991818e-17 98 3.198024e-02 -7.605135e-17 99 9.330968e-02 7.359808e-18 100 1.333280e-01 -3.679904e-17 101 1.401760e-01 3.189250e-17 102 1.099404e-01 1.533293e-17 103 4.815112e-02 -4.783875e-17 104 -3.103786e-02 -2.943923e-17 105 -1.079124e-01 -2.698596e-17 106 -1.618772e-01 1.717289e-17 107 -1.766837e-01 -2.453269e-18 108 -1.449419e-01 6.378501e-17 109 -7.075117e-02 1.005840e-16 110 3.033355e-02 1.251167e-16 111 1.341899e-01 2.134344e-16 112 2.133883e-01 0.000000e+00 113 2.437478e-01 3.618572e-17 114 2.107341e-01 -4.906539e-18 115 1.141877e-01 3.311914e-17 116 -2.984509e-02 2.453269e-17 117 -1.908223e-01 -2.207943e-17 118 -3.288421e-01 3.925231e-17 119 -4.023202e-01 -2.453269e-17 120 -3.767087e-01 -3.925231e-17 121 -2.324518e-01 8.586443e-17 122 2.955767e-02 1.471962e-17 123 3.863025e-01 1.079439e-16 124 7.952850e-01 3.925231e-17 125 1.201145e+00 1.962616e-16 126 1.545075e+00 8.831770e-17 127 1.775235e+00 2.600466e-16 gsl_integration_qng_test(): gsl_integration_qng() uses the Gauss-Kronrod rule to estimate the integral of f(x) over [a,b]. Integral estimate = 1.1547 Error estimate = 0.0122407 Number of function evaluations was 87 Exact integral = 1.1547 Exact error = 8.61352e-09 gsl_linalg_LU_test(): gsl_linalg_LU_decomp() computes LU factors of a matrix A. gsl_linalg_LU_solve() solves linear systems A*x=b. Computed solution vector X: -4.05205 -12.6056 1.66091 8.69377 gsl_midpoint_test(): Test the GSL implicit midpoint method ODE solver. parameters: mu = 10 t0 = 0 y0 = (1,0) tstop = 100 n = 100 1.00000e+00 -1.45695e+00 -1.15464e+01 2.00000e+00 -1.95608e+00 6.90648e-02 3.00000e+00 -1.88481e+00 7.36425e-02 4.00000e+00 -1.80842e+00 7.93679e-02 5.00000e+00 -1.72550e+00 8.68268e-02 6.00000e+00 -1.63383e+00 9.71293e-02 7.00000e+00 -1.52952e+00 1.12714e-01 8.00000e+00 -1.40454e+00 1.40375e-01 9.00000e+00 -1.23605e+00 2.10147e-01 1.00000e+01 -8.53914e-01 9.01981e-01 1.10000e+01 1.99274e+00 -6.69456e-02 1.20000e+01 1.92381e+00 -7.10573e-02 1.30000e+01 1.85032e+00 -7.61105e-02 1.40000e+01 1.77113e+00 -8.25399e-02 1.50000e+01 1.68453e+00 -9.11203e-02 1.60000e+01 1.58767e+00 -1.03412e-01 1.70000e+01 1.47528e+00 -1.23180e-01 1.80000e+01 1.33506e+00 -1.62838e-01 1.90000e+01 1.12080e+00 -3.05564e-01 2.00000e+01 -1.59836e+00 -9.82324e+00 2.10000e+01 -1.96149e+00 6.87434e-02 2.20000e+01 -1.89057e+00 7.32478e-02 2.30000e+01 -1.81462e+00 7.88664e-02 2.40000e+01 -1.73228e+00 8.61589e-02 2.50000e+01 -1.64141e+00 9.61766e-02 2.60000e+01 -1.53829e+00 1.11196e-01 2.70000e+01 -1.41542e+00 1.37411e-01 2.80000e+01 -1.25215e+00 2.00816e-01 2.90000e+01 -9.15983e-01 6.98110e-01 3.00000e+01 1.99798e+00 -6.66424e-02 3.10000e+01 1.92936e+00 -7.07054e-02 3.20000e+01 1.85626e+00 -7.56719e-02 3.30000e+01 1.77758e+00 -8.19715e-02 3.40000e+01 1.69164e+00 -9.03419e-02 3.50000e+01 1.59573e+00 -1.02252e-01 3.60000e+01 1.48486e+00 -1.21185e-01 3.70000e+01 1.34764e+00 -1.58254e-01 3.80000e+01 1.14378e+00 -2.81637e-01 3.90000e+01 -5.76491e-01 -1.26640e+01 4.00000e+01 -1.96686e+00 6.84268e-02 4.10000e+01 -1.89629e+00 7.28598e-02 4.20000e+01 -1.82078e+00 7.83749e-02 4.30000e+01 -1.73901e+00 8.55073e-02 4.40000e+01 -1.64891e+00 9.52531e-02 4.50000e+01 -1.54695e+00 1.09741e-01 4.60000e+01 -1.42608e+00 1.34631e-01 4.70000e+01 -1.26756e+00 1.92562e-01 4.80000e+01 -9.65239e-01 5.67652e-01 4.90000e+01 2.00319e+00 -6.62233e-02 5.00000e+01 1.93489e+00 -7.03588e-02 5.10000e+01 1.86218e+00 -7.52413e-02 5.20000e+01 1.78398e+00 -8.14155e-02 5.30000e+01 1.69869e+00 -8.95843e-02 5.40000e+01 1.60370e+00 -1.01132e-01 5.50000e+01 1.49428e+00 -1.19286e-01 5.60000e+01 1.35988e+00 -1.54033e-01 5.70000e+01 1.16506e+00 -2.61943e-01 5.80000e+01 1.43008e-01 -6.13092e+00 5.90000e+01 -1.97221e+00 6.81150e-02 6.00000e+01 -1.90199e+00 7.24784e-02 6.10000e+01 -1.82691e+00 7.78931e-02 6.20000e+01 -1.74569e+00 8.48712e-02 6.30000e+01 -1.65634e+00 9.43572e-02 6.40000e+01 -1.55550e+00 1.08344e-01 6.50000e+01 -1.43653e+00 1.32017e-01 6.60000e+01 -1.28236e+00 1.85204e-01 6.70000e+01 -1.00604e+00 4.78776e-01 6.80000e+01 2.00833e+00 -6.45142e-02 6.90000e+01 1.94039e+00 -7.00178e-02 7.00000e+01 1.86806e+00 -7.48185e-02 7.10000e+01 1.79034e+00 -8.08715e-02 7.20000e+01 1.70568e+00 -8.88466e-02 7.30000e+01 1.61159e+00 -1.00050e-01 7.40000e+01 1.50356e+00 -1.17478e-01 7.50000e+01 1.37179e+00 -1.50133e-01 7.60000e+01 1.18492e+00 -2.45459e-01 7.70000e+01 4.84678e-01 -3.05785e+00 7.80000e+01 -1.97754e+00 6.78075e-02 7.90000e+01 -1.90765e+00 7.21033e-02 8.00000e+01 -1.83299e+00 7.74208e-02 8.10000e+01 -1.75231e+00 8.42500e-02 8.20000e+01 -1.66370e+00 9.34878e-02 8.30000e+01 -1.56394e+00 1.07002e-01 8.40000e+01 -1.44678e+00 1.29552e-01 8.50000e+01 -1.29661e+00 1.78598e-01 8.60000e+01 -1.04095e+00 4.15153e-01 8.70000e+01 2.01299e+00 -4.86475e-02 8.80000e+01 1.94587e+00 -6.96820e-02 8.90000e+01 1.87391e+00 -7.44032e-02 9.00000e+01 1.79666e+00 -8.03389e-02 9.10000e+01 1.71262e+00 -8.81282e-02 9.20000e+01 1.61939e+00 -9.90043e-02 9.30000e+01 1.51270e+00 -1.15753e-01 9.40000e+01 1.38342e+00 -1.46515e-01 9.50000e+01 1.20360e+00 -2.31461e-01 9.60000e+01 6.68080e-01 -1.79244e+00 9.70000e+01 -1.98284e+00 6.75046e-02 9.80000e+01 -1.91329e+00 7.17344e-02 9.90000e+01 -1.83904e+00 7.69575e-02 1.00000e+02 -1.75889e+00 8.36432e-02 Solution data stored in 'gsl_midpoint_data.txt'. gsl_multiroot_fsolver_test(): Demonstrate the ability to find a root of a set of nonlinear equations. In this case, we have two functions in two unknowns, and the only root is X = (1,1). iter = 0 x = -10.000 -5.000 f(x) = 1.100e+01 -1.050e+03 iter = 1 x = -10.000 -5.000 f(x) = 1.100e+01 -1.050e+03 iter = 2 x = -3.976 24.827 f(x) = 4.976e+00 9.020e+01 iter = 3 x = -3.976 24.827 f(x) = 4.976e+00 9.020e+01 iter = 4 x = -3.976 24.827 f(x) = 4.976e+00 9.020e+01 iter = 5 x = -1.274 -5.680 f(x) = 2.274e+00 -7.302e+01 iter = 6 x = -1.274 -5.680 f(x) = 2.274e+00 -7.302e+01 iter = 7 x = 0.249 0.298 f(x) = 7.511e-01 2.359e+00 iter = 8 x = 0.249 0.298 f(x) = 7.511e-01 2.359e+00 iter = 9 x = 1.000 0.878 f(x) = -2.653e-10 -1.218e+00 iter = 10 x = 1.000 0.989 f(x) = -2.353e-11 -1.080e-01 iter = 11 x = 1.000 1.000 f(x) = 0.000e+00 0.000e+00 status = success gsl_poly_eval_test(): gsl_poly_eval() evaluates a real polynomial. p(x) = (x-1)*(x-3)*(x-6) = x^3 - 10x^2 + 27x - 18. x p(x) 0 -18 0.5 -6.88 1 0 1.5 3.38 2 4 2.5 2.62 3 0 3.5 -3.12 4 -6 4.5 -7.88 5 -8 5.5 -5.62 6 0 6.5 9.62 7 24 7.5 43.9 8 70 8.5 103 9 144 9.5 193 10 252 gsl_qrng_niederreiter_2_test: gsl_qrng_alloc sets aside space for a sequence; gsl_qrng_niederreiter_2 requests the Niederreiter_2 sequence; gsl_qrng_get gets the next entry of the requested sequence; Determine the first 25 points of the Niederreiter2 quasi-random sequence in 2 dimensions. I X(I) 0 0 0 1 0.5 0.5 2 0.75 0.25 3 0.25 0.75 4 0.375 0.375 5 0.875 0.875 6 0.625 0.125 7 0.125 0.625 8 0.1875 0.3125 9 0.6875 0.8125 10 0.9375 0.0625 11 0.4375 0.5625 12 0.3125 0.1875 13 0.8125 0.6875 14 0.5625 0.4375 15 0.0625 0.9375 16 0.09375 0.46875 17 0.59375 0.96875 18 0.84375 0.21875 19 0.34375 0.71875 20 0.46875 0.09375 21 0.96875 0.59375 22 0.71875 0.34375 23 0.21875 0.84375 24 0.15625 0.15625 gsl_qrng_sobol_test(): gsl_qrng_alloc() sets aside space for a sequence; gsl_qrng_sobol() requests the Sobol sequence; gsl_qrng_get() gets the next entry of the requested sequence; Determine the first 25 points of the Sobol quasi-random sequence in 2 dimensions. I X(I) 0 0.5 0.5 1 0.75 0.25 2 0.25 0.75 3 0.375 0.375 4 0.875 0.875 5 0.625 0.125 6 0.125 0.625 7 0.1875 0.3125 8 0.6875 0.8125 9 0.9375 0.0625 10 0.4375 0.5625 11 0.3125 0.1875 12 0.8125 0.6875 13 0.5625 0.4375 14 0.0625 0.9375 15 0.09375 0.46875 16 0.59375 0.96875 17 0.84375 0.21875 18 0.34375 0.71875 19 0.46875 0.09375 20 0.96875 0.59375 21 0.71875 0.34375 22 0.21875 0.84375 23 0.15625 0.15625 24 0.65625 0.65625 gsl_sf_bessel_J0_test(): gsl_sf_bessel_J0() evaluates the J0 Bessel function. X Exact Computed -5 -0.177597 -0.177597 -4 -0.39715 -0.39715 -3 -0.260052 -0.260052 -2 0.223891 0.223891 -1 0.765198 0.765198 0 1 1 1 0.765198 0.765198 2 0.223891 0.223891 3 -0.260052 -0.260052 4 -0.39715 -0.39715 5 -0.177597 -0.177597 6 0.150645 0.150645 7 0.300079 0.300079 8 0.171651 0.171651 9 -0.0903336 -0.0903336 10 -0.245936 -0.245936 11 -0.17119 -0.17119 12 0.0476893 0.0476893 13 0.206926 0.206926 14 0.171073 0.171073 15 -0.0142245 -0.0142245 gsl_sf_bessel_J1_test(): gsl_sf_bessel_J1() evaluates the J1 Bessel function. X Exact Computed -5 0.327579 0.327579 -4 0.0660433 0.0660433 -3 -0.339059 -0.339059 -2 -0.576725 -0.576725 -1 -0.440051 -0.440051 0 0 0 1 0.440051 0.440051 2 0.576725 0.576725 3 0.339059 0.339059 4 -0.0660433 -0.0660433 5 -0.327579 -0.327579 6 -0.276684 -0.276684 7 -0.00468282 -0.00468282 8 0.234636 0.234636 9 0.245312 0.245312 10 0.0434727 0.0434727 11 -0.176785 -0.176785 12 -0.223447 -0.223447 13 -0.0703181 -0.0703181 14 0.133375 0.133375 15 0.205104 0.205104 gsl_sf_coupling_3j_test(): gsl_sf_coupling_3j() returns values of the Wigner 3J coefficient. J1 J2 J3 M1 M2 M3 THREE_J 1 4.5 3.5 1 -3.5 2.5 0.2788866755113585 0.2788866755113585 2 4.5 3.5 1 -3.5 2.5 -0.09534625892455922 -0.09534625892455932 3 4.5 3.5 1 -3.5 2.5 -0.06741998624632421 -0.06741998624632435 4 4.5 3.5 1 -3.5 2.5 0.1533110351679666 0.1533110351679665 5 4.5 3.5 1 -3.5 2.5 -0.156446554693686 -0.1564465546936862 6 4.5 3.5 1 -3.5 2.5 0.1099450412156551 0.1099450412156551 7 4.5 3.5 1 -3.5 2.5 -0.05536235693131719 -0.05536235693131723 8 4.5 3.5 1 -3.5 2.5 0.01799835451137786 0.01799835451137783 gsl_sf_coupling_6j_test(): gsl_sf_coupling_6j() returns values of the Wigner 6J coefficient. J1 J2 J3 J4 J5 J6 SIX_J 1 8 7 6.5 7.5 7.5 0.034909051383733 0.03490905138373299 2 8 7 6.5 7.5 7.5 -0.03743025039659792 -0.03743025039659792 3 8 7 6.5 7.5 7.5 0.0189086639095956 0.0189086639095956 4 8 7 6.5 7.5 7.5 0.007342448254928643 0.007342448254928649 5 8 7 6.5 7.5 7.5 -0.02358935185081794 -0.02358935185081793 6 8 7 6.5 7.5 7.5 0.01913476955215437 0.01913476955215436 7 8 7 6.5 7.5 7.5 0.001288017397724172 0.001288017397724189 8 8 7 6.5 7.5 7.5 -0.01930018366290527 -0.01930018366290528 9 8 7 6.5 7.5 7.5 0.01677305949382889 0.0167730594938289 10 8 7 6.5 7.5 7.5 0.005501147274850949 0.005501147274850945 11 8 7 6.5 7.5 7.5 -0.02135439790896831 -0.02135439790896832 12 8 7 6.5 7.5 7.5 0.003460364451435387 0.003460364451435383 13 8 7 6.5 7.5 7.5 0.02520950054795585 0.02520950054795584 14 8 7 6.5 7.5 7.5 0.01483990561221713 0.01483990561221713 15 8 7 6.5 7.5 7.5 0.002708577680633186 0.002708577680633186 gsl_sf_coupling_9j_test(): gsl_sf_coupling_9j() returns values of the Wigner 9J coefficient. SOMETHING IS SERIOUSLY WRONG HERE. ALL RESULTS ARE 0. JVB, 15 February 2022 J1 J2 J3 J4 J5 J6 J7 J8 J9 NINE_J 2 16 14 13 15 15 12 20 12 1 8 7 6.5 7.5 7.5 6 10 6 0.0004270039294528318 0 3 16 14 13 15 15 12 20 12 1.5 8 7 6.5 7.5 7.5 6 10 6 -0.001228915451058514 0 4 16 14 13 15 15 12 20 12 2 8 7 6.5 7.5 7.5 6 10 6 -0.0001944260688400887 0 2 6 4 8 3 6 7 4 4 1 3 2 4 1.5 3 3.5 2 2 0.003338419923885592 0 3 6 4 8 3 6 7 4 4 1.5 3 2 4 1.5 3 3.5 2 2 -0.0007958936865080434 0 4 6 4 8 3 6 7 4 4 2 3 2 4 1.5 3 3.5 2 2 -0.004338208690251972 0 1 1 2 4 2 3 3 1 3 0.5 0.5 1 2 1 1.5 1.5 0.5 1.5 0.05379143536399187 0 2 1 2 4 2 3 3 1 3 1 0.5 1 2 1 1.5 1.5 0.5 1.5 0.006211299937499411 0 3 1 2 4 2 3 3 1 3 1.5 0.5 1 2 1 1.5 1.5 0.5 1.5 0.03042903097250921 0 gsl_sf_dawson_test(): gsl_sf_dawson() evaluates Dawson's integral. X Exact Computed 0 0 0 0.1 0.099336 0.099336 0.2 0.194751 0.194751 0.3 0.282632 0.282632 0.4 0.359943 0.359943 0.5 0.424436 0.424436 0.6 0.474763 0.474763 0.7 0.510504 0.510504 0.8 0.532102 0.532102 0.9 0.540724 0.540724 1 0.53808 0.53808 1.1 0.526207 0.526207 1.2 0.507273 0.507273 1.3 0.483398 0.483398 1.4 0.456507 0.456507 1.5 0.428249 0.428249 1.6 0.39994 0.39994 1.7 0.372559 0.372559 1.8 0.346773 0.346773 1.9 0.322974 0.322974 2 0.30134 0.30134 gsl_sf_hyperg_2F1_test(): Test gsl_sf_hyperg_2F1(), which evaluates the hypergeometric function 2F1 for real parameters and arguments. A B C X Hyper_2F1(A,B,C;X) -2.500000 3.300000 6.700000 0.250000 0.7235612934899779 0.7235612934899777 -0.500000 1.100000 6.700000 0.250000 0.9791110934527796 0.9791110934527796 0.500000 1.100000 6.700000 0.250000 1.021657814008856 1.021657814008857 2.500000 3.300000 6.700000 0.250000 1.405156320011213 1.405156320011212 -2.500000 3.300000 6.700000 0.550000 0.4696143163982161 0.4696143163982158 -0.500000 1.100000 6.700000 0.550000 0.9529619497744632 0.9529619497744686 0.500000 1.100000 6.700000 0.550000 1.051281421394799 1.051281421394799 2.500000 3.300000 6.700000 0.550000 2.399906290477786 2.399906290477785 -2.500000 3.300000 6.700000 0.850000 0.2910609592841472 0.2910609592841466 -0.500000 1.100000 6.700000 0.850000 0.9253696791037318 0.9253696791037372 0.500000 1.100000 6.700000 0.850000 1.0865504094807 1.086550409480699 2.500000 3.300000 6.700000 0.850000 5.738156552618904 5.738156552618895 3.300000 6.700000 -5.500000 0.250000 15090.66974870461 15090.6697487046 1.100000 6.700000 -0.500000 0.250000 -104.3117006736435 -104.3117006736435 1.100000 6.700000 0.500000 0.250000 21.17505070776881 21.17505070776881 3.300000 6.700000 4.500000 0.250000 4.194691581903192 4.194691581903191 3.300000 6.700000 -5.500000 0.550000 10170777974.04881 10170777974.04814 1.100000 6.700000 -0.500000 0.550000 -24708.63532248916 -24708.63532248926 1.100000 6.700000 0.500000 0.550000 1372.230454838499 1372.230454838499 3.300000 6.700000 4.500000 0.550000 58.09272870639465 58.09272870639469 3.300000 6.700000 -5.500000 0.850000 5.868208761512417e+18 5.868208761512327e+18 1.100000 6.700000 -0.500000 0.850000 -446350101.47296 -446350101.4729599 1.100000 6.700000 0.500000 0.850000 5383505.756129573 5383505.756129567 3.300000 6.700000 4.500000 0.850000 20396.91377601966 20396.91377601964 gsl_test(): Normal end of execution. 26 December 2023 01:53:36 PM