C 代码 包含 R8COL 的实用程序例程，即 双精度实 MxN 数组，被视为 N 列向量， 每个长度为M.rar

02 March 2022 04:19:46 PM R8COL_TEST C version Test the R8COL library. I4_LOG_10_TEST I4_LOG_10: whole part of log base 10, X, I4_LOG_10 0 0 1 0 2 0 3 0 9 0 10 1 11 1 99 1 101 2 -1 0 -2 0 -3 0 -9 0 I4_MAX_TEST I4_MAX returns the maximum of two I4's. A B C=I4_MAX(A,B) -57 92 92 66 12 66 -17 -87 -17 -49 -78 -49 -92 27 27 -88 -10 -10 -20 51 51 60 -100 60 80 -30 80 -81 -98 -81 I4_MIN_TEST I4_MIN returns the minimum of two I4's. A B C=I4_MIN(A,B) -57 92 -57 66 12 12 -17 -87 -87 -49 -78 -78 -92 27 -92 -88 -10 -88 -20 51 -20 60 -100 -100 80 -30 -30 -81 -98 -98 I4_POWER_TEST I4_POWER computes I^J I J I4_POWER(I,J) 0 1 0 1 2 1 2 3 8 3 3 27 10 3 1000 -1 4 1 -2 5 -32 I4_UNIFORM_AB_TEST I4_UNIFORM_AB computes pseudorandom values in an interval [A,B]. The lower endpoint A = -100 The upper endpoint B = 200 The initial seed is 123456789 1 -35 2 187 3 149 4 69 5 25 6 -81 7 -23 8 -67 9 -87 10 90 11 -82 12 35 13 20 14 127 15 139 16 -100 17 170 18 5 19 -72 20 -96 I4VEC_INDICATOR1_NEW_TEST I4VEC_INDICATOR1_NEW returns an indicator vector. The indicator1 vector: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 I4VEC_PRINT_TEST I4VEC_PRINT prints an I4VEC Here is the I4VEC: 0: 91 1: 92 2: 93 3: 94 I4VEC_TRANSPOSE_PRINT_TEST I4VEC_TRANSPOSE_PRINT prints an I4VEC with 5 entries to a row, and an optional title. Output from I4VEC_PRINT: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 10: 11 11: 12 Now call I4VEC_TRANSPOSE_PRINT with a short title: My array: 1 2 3 4 5 6 7 8 9 10 11 12 Now call I4VEC_TRANSPOSE_PRINT with no title: 1 2 3 4 5 6 7 8 9 10 11 12 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 perm0_check returns ierror = 1 Permutation 2: 4 1 3 0 2 perm0_check returns ierror = 0 Permutation 3: 0 2 1 3 2 PERM0_CHECK - Warning! Permutation is missing value 4 perm0_check returns ierror = 1 R8_UNIFORM_AB_TEST R8_UNIFORM_AB produces a random real in a given range. Using range 10.000000 <= A <= 25.000000. I A 0 10.001996 1 13.542136 2 22.679357 3 21.955784 4 10.855646 5 20.834593 6 21.996377 7 18.110503 8 18.220583 9 23.346688 R8COL_FIND_TEST R8COL_FIND finds a column in a table matching a given set of data. R8COL_FIND returns COL = 3 R8COL_INDICATOR_NEW_TEST R8COL_INDICATOR_NEW returns an indicator matrix; Indicator matrix: Col: 0 1 2 3 Row 0: 11 12 13 14 1: 21 22 23 24 2: 31 32 33 34 3: 41 42 43 44 4: 51 52 53 54 R8COL_INSERT_TEST R8COL_INSERT inserts new columns into a sorted R8COL. The unsorted matrix: Col: 0 1 2 3 Row 0: 2 4 1 3 1: 6 8 5 7 2: 10 12 9 11 The sorted matrix: Col: 0 1 2 3 Row 0: 1 2 3 4 1: 5 6 7 8 2: 9 10 11 12 New column: 0: 3 1: 7 2: 11 The data was already in column 3 New column: 0: 3 1: 4 2: 18 The updated matrix: Col: 0 1 2 3 Row 0: 1 2 3 3 1: 5 6 4 7 2: 9 10 18 11 R8COL_MAX_TEST R8COL_MAX computes maximums of an R8COL; The array: Col: 0 1 2 3 Row 0: 1 2 3 4 1: 5 6 7 8 2: 9 10 11 12 Column maximums: 0: 9 1: 10 2: 11 3: 12 R8COL_MEAN_TEST R8COL_MEAN computes means of an R8COL; The array: Col: 0 1 2 3 Row 0: 1 2 3 4 1: 5 6 7 8 2: 9 10 11 12 The column means: 0: 5 1: 6 2: 7 3: 8 R8COL_MIN_TEST R8COL_MIN computes minimums of an R8COL; The array: Col: 0 1 2 3 Row 0: 1 2 3 4 1: 5 6 7 8 2: 9 10 11 12 Column minimums: 0: 1 1: 2 2: 3 3: 4 R8COL_PERMUTE_TEST R8COL_PERMUTE permutes an R8COL in place. A (unpermuted): Col: 0 1 2 3 4 Row 0: 11 12 13 14 15 1: 21 22 23 24 25 2: 31 32 33 34 35 The (column) permutation vector: 0: 1 1: 3 2: 4 3: 0 4: 2 A (permuted): Col: