C 代码包含许多实用程序例程(线性代数)， 使用双精度实数 （R8） 算法.rar

28 February 2022 11:25:07 AM r8lib_test(): C version GCC version 7.5.0 Test r8lib(). I4INT_TO_R8INT_TEST For data in an interval, I4INT_TO_R8INT converts an integer to a real; Integer interval: [1, 11] Real interval: [100, 200] R I(R) R(I(R)) 113.394 2 110 209.321 12 210 192.836 10 190 158.02 7 160 138.99 5 140 93.5954 0 90 118.485 3 120 99.2944 1 100 90.6978 0 90 167.416 8 170 PERM0_CHECK_TEST PERM0_CHECK checks a permutation of 0, ..., N-1. Permutation 1: 5 2 3 4 1 PERM0_CHECK - Warning! Permutation is missing value 0 Permutation 2: 4 1 3 0 2 Permutation 3: 0 2 1 3 2 PERM0_CHECK - Warning! Permutation is missing value 4 PERM0_UNIFORM_TEST PERM0_UNIFORM randomly selects a permutation of 0,...,N-1. 2 9 8 6 3 5 7 4 0 1 0 5 1 8 3 2 9 4 7 6 8 0 2 3 5 9 6 1 4 7 6 0 8 1 9 7 3 4 2 5 5 4 6 0 3 8 2 7 9 1 PERM1_CHECK_TEST PERM1_CHECK checks a permutation of 1, ..., N. Permutation 1: 5 2 3 4 1 Permutation 2: 4 1 3 0 2 PERM1_CHECK - Fatal error! Permutation is missing value 5 Permutation 3: 0 2 1 3 2 PERM1_CHECK - Fatal error! Permutation is missing value 4 PERM1_UNIFORM_TEST PERM1_UNIFORM randomly selects a permutation of 1,...,N. 3 10 9 7 4 6 8 5 1 2 1 6 2 9 4 3 10 5 8 7 9 1 3 4 6 10 7 2 5 8 7 1 9 2 10 8 4 5 3 6 6 5 7 1 4 9 3 8 10 2 R8_ABS_TEST R8_ABS returns the absolute value of an R8. X R8_ABS(X) -1.252654 1.252654 4.650541 4.650541 3.636074 3.636074 1.493564 1.493564 0.322457 0.322457 -2.471050 2.471050 -0.939378 0.939378 -2.120346 2.120346 -2.649368 2.649368 2.071726 2.071726 R8_ACOS_TEST R8_ACOS computes the arc-cosine of an angle. C R8_ACOS(C) ACOS(C) -1.16667 3.14159 -1 3.14159 3.14159 -0.833333 2.55591 2.55591 -0.666667 2.30052 2.30052 -0.5 2.0944 2.0944 -0.333333 1.91063 1.91063 -0.166667 1.73824 1.73824 0 1.5708 1.5708 0.166667 1.40335 1.40335 0.333333 1.23096 1.23096 0.5 1.0472 1.0472 0.666667 0.841069 0.841069 0.833333 0.585686 0.585686 1 0 0 1.16667 0 R8_ACOSH_TEST R8_ACOSH computes the arc-hyperbolic-cosine of an angle. X A=R8_ACOSH(X) COSH(A) 1 0 1 1.5 0.962424 1.5 2 1.31696 2 2.5 1.5668 2.5 3 1.76275 3 3.5 1.92485 3.5 4 2.06344 4 4.5 2.18464 4.5 5 2.29243 5 R8_ASINH_TEST R8_ASINH computes the inverse hyperbolic sine of a given value. X R8_ASINH(X) SINH(R8_ASINH(X)) 1.000000 0.881374 1.000000 1.200000 1.015973 1.200000 1.400000 1.137982 1.400000 1.600000 1.248983 1.600000 1.800000 1.350441 1.800000 2.000000 1.443635 2.000000 2.200000 1.529660 2.200000 2.400000 1.609438 2.400000 2.600000 1.683743 2.600000 2.800000 1.753229 2.800000 3.000000 1.818446 3.000000 R8_ATAN_TEST R8_ATAN computes the arc-tangent given Y and X; ATAN2 is the system version of this routine. X Y ATAN2(Y,X) R8_ATAN(Y,X) 1.000000 0.000000 0.000000 0.000000 1.000000 1.000000 0.785398 0.785398 0.000000 1.000000 1.570796 1.570796 -1.000000 1.000000 2.356194 2.356194 -1.000000 0.000000 3.141593 3.141593 -1.000000 -1.000000 -2.356194 3.926991 0.000000 -1.000000 -1.570796 4.712389 1.000000 -1.000000 -0.785398 5.497787 R8_ATANH_TEST R8_ATANH computes the inverse hyperbolic tangent of a given value. X R8_ATANH(X) TANH(R8_ATANH(X)) -0.2 -0.202733 -0.2 -0.1 -0.100335 -0.1 0 0 0 0.1 0.100335 0.1 0.2 0.202733 0.2 0.3 0.30952 0.3 0.4 0.423649 0.4 0.5 0.549306 0.5 0.6 0.693147 0.6 0.7 0.867301 0.7 0.8 1.09861 0.8 0.9 1.47222 0.9 R8_BIG_TEST R8_BIG returns a 'big' R8 value; R8_BIG = 1e+30 R8_CAS_TEST R8_CAS evaluates the casine of a number. X R8_CAS ( X ) 0.000000 1.000000 0.261799 1.224745 0.523599 1.366025 0.785398 1.414214 1.047198 1.366025 1.308997 1.224745 1.570796 1.000000 1.832596 0.707107 2.094395 0.366025 2.356194 0.000000 2.617994 -0.366025 2.879793 -0.707107 3.141593 -1.000000 R8_CEILING_TEST R8_CEILING rounds a value up. X R8_CEILING(X) -1.200000 -1.000000 -1.000000 -1.000000 -0.800000 0.000000 -0.600000 0.000000 -0.400000 0.000000 -0.200000 0.000000 0.000000 0.000000 0.200000 1.000000 0.400000 1.000000 0.600000 1.000000 0.800000 1.000000 1.000000 1.000000 1.200000 2.000000 R8_CHOOSE_TEST R8_CHOOSE evaluates C(N,K). N K CNK 0 0 1 1 0 1 1 1 1 2 0 1 2 1 2 2 2 1 3 0 1 3 1 3 3 2 3 3 3 1 4 0 1 4 1 4 4 2 6 4 3 4 4 4 1 5 0 1 5 1 5 5 2 10 5 3 10 5 4 5 5 5 1 R8_COSD_TEST R8_COSD computes the cosine of an angle given in degrees. ANGLE R8_COSD(ANGLE) 0.00 1 15.00 0.965926 30.00 0.866025 45.00 0.707107 60.00 0.5 75.00 0.258819 9