22 January 2017 12:05:13 PM ST_TO_CC_PRB C version Test the ST_TO_CC library. TEST01 Convert a sparse matrix from ST to CC format. ST: sparse triplet, I, J, A. CC: compressed column, I, CC, A. Sparse Triplet (ST) header information: Minimum row index I_MIN = 1 Maximum row index I_MAX = 5 Minimum col index J_MIN = 1 Maximum col index J_MAX = 5 Number of rows M = 5 Number of columns N = 5 Number of nonzeros NST = 12 The matrix in ST format: # I J A ---- ---- ---- -------------- 0 0 0 2 1 0 1 3 2 1 0 3 3 1 2 4 4 1 4 6 5 2 1 -1 6 2 2 -3 7 2 3 2 8 3 2 1 9 4 1 4 10 4 2 2 11 4 4 1 Number of CC values = 12 CC Matrix: # I J A ---- ---- ---- ---------------- 0 0 0 2 1 1 0 3 2 0 0 3 3 2 0 -1 4 4 0 4 5 1 1 4 6 2 1 -3 7 3 1 1 8 4 1 2 9 2 2 2 10 1 3 6 11 4 3 1 TEST02 Convert a sparse matrix from ST to CC format. ST: sparse triplet, I, J, A. CC: compressed column, I, CC, A. This matrix is read from the file 'west_st.txt' Sparse Triplet (ST) header information: Minimum row index I_MIN = 1 Maximum row index I_MAX = 67 Minimum col index J_MIN = 1 Maximum col index J_MAX = 67 Number of rows M = 67 Number of columns N = 67 Number of nonzeros NST = 299 The matrix in ST format: # I J A ---- ---- ---- -------------- 0 44 55 -1.863354 1 54 61 -1.863354 2 29 37 -1.567398 3 44 56 -1.490683 4 54 62 -1.490683 5 9 12 -1.265823 6 29 38 -1.253918 7 44 57 -1.118012 8 54 63 -1.118012 9 15 31 -1.05 10 16 32 -1.05 11 17 33 -1.05 12 18 34 -1.05 13 19 35 -1.05 14 24 31 -1.05 15 25 32 -1.05 16 26 33 -1.05 17 27 34 -1.05 18 28 35 -1.05 19 9 13 -1.012658 20 10 20 -1 21 11 21 -1 22 12 22 -1 23 13 23 -1 24 14 24 -1 25 30 49 -0.9722222 26 31 50 -0.9722222 27 32 51 -0.9722222 28 33 52 -0.9722222 29 34 53 -0.9722222 30 39 49 -0.9722222 31 40 50 -0.9722222 32 41 51 -0.9722222 33 42 52 -0.9722222 34 43 53 -0.9722222 35 35 25 -0.9583187 36 36 26 -0.9583187 37 37 27 -0.9583187 38 38 28 -0.9583187 39 49 55 -0.9444444 40 50 56 -0.9444444 41 51 57 -0.9444444 42 52 58 -0.9444444 43 53 59 -0.9444444 44 29 39 -0.9404389 45 20 1 -0.9159533 46 21 2 -0.9159533 47 22 3 -0.9159533 48 23 4 -0.9159533 49 0 7 -0.8341818 50 1 8 -0.8341818 51 2 9 -0.8341818 52 3 10 -0.8341818 53 45 43 -0.8242248 54 46 44 -0.8242248 55 47 45 -0.8242248 56 48 46 -0.8242248 57 4 1 -0.8 58 5 2 -0.8 59 6 3 -0.8 60 7 4 -0.8 61 8 5 -0.8 62 9 14 -0.7594937 63 44 58 -0.7453416 64 54 64 -0.7453416 65 29 40 -0.6269592 66 9 15 -0.5063291 67 44 59 -0.3726708 68 54 65 -0.3726708 69 0 17 -0.3361556 70 29 41 -0.3134796 71 1 17 -0.2939196 72 20 42 -0.2788416 73 4 0 -0.2788416 74 21 42 -0.2680186 75 5 0 -0.2680186 76 52 54 -0.2667757 77 53 54 -0.2630706 78 48 66 -0.2541193 79 9 16 -0.2531646 80 37 60 -0.2475675 81 47 66 -0.2421498 82 33 48 -0.2362845 83 42 48 -0.2362845 84 36 60 -0.2356469 85 22 42 -0.2323717 86 6 0 -0.2323717 87 32 48 -0.2303917 88 41 48 -0.2303917 89 17 30 -0.2286264 90 26 30 -0.2286264 91 16 30 -0.2232997 92 25 30 -0.2232997 93 2 17 -0.2214815 94 12 19 -0.2144206 95 11 19 -0.2140392 96 51 54 -0.2122056 97 38 60 -0.2074873 98 10 19 -0.2071759 99 15 30 -0.2070986 100 24 30 -0.2070986 101 35 60 -0.2069954 102 18 30 -0.2024528 103 27 30 -0.2024528 104 13 19 -0.1986768 105 31 48 -0.1947711 106 40 48 -0.1947711 107 46 66 -0.1918557 108 34 48 -0.18039 109 43 48 -0.18039 110 14 19 -0.1656874 111 30 48 -0.1581626 112 39 48 -0.1581626 113 23 42 -0.1575082 114 7 0 -0.1575082 115 50 54 -0.1514908 116 45 66 -0.1443354 117 19 30 -0.1385226 118 28 30 -0.1385226 119 3 17 -0.118986 120 49 54 -0.1064573 121 8 0 -0.06325978 122 8 6 0.01178291 123 28 0 0.03162989 124 7 6 0.047594391 125 39 54 0.053228639 126 40 54 0.075745419 127 27 0 0.078754112 128 24 36 0.08147449 129 6 6 0.088592619 130 30 36 0.090527207 131 28 36 0.092419088 132 25 36 0.097890154 133 19 19 0.099412464 134 34 36 0.1026879 135 41 54 0.1061028 136 31 36 0.1087668 137 26 36 0.1131608 138 27 36 0.1150555 139 26 0 0.1161859 140 5 6 0.1175679 141 18 19 0.1192061 142 15 19 0.1243055 143 32 36 0.1257342 144 33 36 0.1278394 145 16 19 0.1284235 146 17 19 0.1286524 147 43 54 0.1315353 148 42 54 0.1333878 149 25 0 0.1340093 150 4 6 0.1344622 151 24 0 0.1394208 152 30 43 0.25 153 31 44 0.25 154 32 45 0.25 155 33 46 0.25 156 34 47 0.25 157 10 12 0.3333333 158 11 13 0.3333333 159 12 14 0.3333333 160 13 15 0.3333333 161 14 16 0.3333333 162 24 1 0.4 163 25 2 0.4 164 26 3 0.4 165 27 4 0.4 166 28 5 0.4 167 4 12 0.4 168 4 7 0.4 169 5 13 0.4 170 5 8 0.4 171 6 14 0.4 172 6 9 0.4 173 7 10 0.4 174 7 15 0.4 175 8 11 0.4 176 8 16 0.4 177 49 61 0.4444444 178 50 62 0.4444444 179 51 63 0.4444444 180 52 64 0.4444444 181 53 65 0.4444444 182 15 25 0.45 183 16 26 0.45 184 17 27 0.45 185 18 28 0.45 186 19 29 0.45 187 39 55 0.4722222 188 40 56 0.4722222 189 41 57 0.4722222 190 42 58 0.4722222 191 43 59 0.4722222 192 39 25 0.5 193 40 26 0.5 194 41 27 0.5 195 42 28 0.5 196 43 29 0.5 197 49 43 0.5 198 50 44 0.5 199 51 45 0.5 200 52 46 0.5 201 53 47 0.5 202 3 15 0.5063291 203 15 20 0.6 204 16 21 0.6 205 17 22 0.6 206 18 23 0.6 207 19 24 0.6 208 23 40 0.6269592 209 24 37 0.65 210 25 38 0.65 211 26 39 0.65 212 27 40 0.65 213 28 41 0.65 214 14 18 0.6666667 215 30 37 0.7222222 216 31 38 0.7222222 217 32 39 0.7222222 218 33 40 0.7222222 219 34 41 0.7222222 220 38 58 0.7453416 221 48 64 0.7453416 222 2 14 0.7594937 223 22 39 0.9404389 224 59 31 0.5 225 59 32 0.5 226 59 33 0.5 227 59 34 0.5 228 59 35 0.5 229 29 42 1 230 44 60 1 231 54 66 1 232 55 18 1 233 56 10 1 234 56 11 1 235 56 7 1 236 56 8 1 237 56 9 1 238 57 12 1 239 57 13 1 240 57 14 1 241 57 15 1 242 57 16 1 243 58 20 1 244 58 21 1 245 58 22 1 246 58 23 1 247 58 24 1 248 59 31 0.5 249 59 32 0.5 250 59 33 0.5 251 59 34 0.5 252 59 35 0.5 253 60 1 1 254 60 2 1 255 60 3 1 256 60 4 1 257 60 5 1 258 61 37 1 259 61 38 1 260 61 39 1 261 61 40 1 262 61 41 1 263 62 49 1 264 62 50 1 265 62 51 1 266 62 52 1 267 62 53 1 268 63 25 1 269 63 26 1 270 63 27 1 271 63 28 1 272 63 29 1 273 64 55 1 274 64 56 1 275 64 57 1 276 64 58 1 277 64 59 1 278 65 43 1 279 65 44 1 280 65 45 1 281 65 46 1 282 65 47 1 283 66 61 1 284 66 62 1 285 66 63 1 286 66 64 1 287 66 65 1 288 9 17 1 289 1 13 1.012658 290 37 57 1.118012 291 47 63 1.118012 292 21 38 1.253918 293 0 12 1.265823 294 36 56 1.490683 295 46 62 1.490683 296 20 37 1.567398 297 35 55 1.863354 298 45 61 1.863354 Number of CC values = 294 CC Matrix: # I J A ---- ---- ---- ---------------- 0 4 0 -0.2788416 1 5 0 -0.2680186 2 6 0 -0.2323717 3 7 0 -0.1575082 4 8 0 -0.06325978 5 24 0 0.1394208 6 25 0 0.1340093 7 26 0 0.1161859 8 27 0 0.078754112 9 28 0 0.03162989 10 4 0 -0.8 11 20 0 -0.9159533 12 24 0 0.4 13 60 0 1 14 5 1 -0.8 15 21 1 -0.9159533 16 25 1 0.4 17 60 1 1 18 6 2 -0.8 19 22 2 -0.9159533 20 26 2 0.4 21 60 2 1 22 7 3 -0.8 23 23 3 -0.9159533 24 27 3 0.4 25 60 3 1 26 8 4 -0.8 27 28 4 0.4 28 60 4 1 29 4 5 0.1344622 30 5 5 0.1175679 31 6 5 0.088592619 32 7 5 0.047594391 33 8 5 0.01178291 34 0 6 -0.8341818 35 4 6 0.4 36 56 6 1 37 1 7 -0.8341818 38 5 7 0.4 39 56 7 1 40 2 8 -0.8341818 41 6 8 0.4 42 56 8 1 43 3 9 -0.8341818 44 7 9 0.4 45 56 9 1 46 8 10 0.4 47 56 10 1 48 0 11 1.265823 49 4 11 0.4 50 9 11 -1.265823 51 10 11 0.3333333 52 57 11 1 53 1 12 1.012658 54 5 12 0.4 55 9 12 -1.012658 56 11 12 0.3333333 57 57 12 1 58 2 13 0.7594937 59 6 13 0.4 60 9 13 -0.7594937 61 12 13 0.3333333 62 57 13 1 63 3 14 0.5063291 64 7 14 0.4 65 9 14 -0.5063291 66 13 14 0.3333333 67 57 14 1 68 8 15 0.4 69 9 15 -0.2531646 70 14 15 0.3333333 71 57 15 1 72 0 16 -0.3361556 73 1 16 -0.2939196 74 2 16 -0.2214815 75 3 16 -0.118986 76 9 16 1 77 14 17 0.6666667 78 55 17 1 79 10 18 -0.2071759 80 11 18 -0.2140392 81 12 18 -0.2144206 82 13 18 -0.1986768 83 14 18 -0.1656874 84 15 18 0.1243055 85 16 18 0.1284235 86 17 18 0.1286524 87 18 18 0.1192061 88 19 18 0.099412464 89 10 19 -1 90 15 19 0.6 91 58 19 1 92 11 20 -1 93 16 20 0.6 94 58 20 1 95 12 21 -1 96 17 21 0.6 97 58 21 1 98 13 22 -1 99 18 22 0.6 100 58 22 1 101 14 23 -1 102 19 23 0.6 103 58 23 1 104 15 24 0.45 105 35 24 -0.9583187 106 39 24 0.5 107 63 24 1 108 16 25 0.45 109 36 25 -0.9583187 110 40 25 0.5 111 63 25 1 112 17 26 0.45 113 37 26 -0.9583187 114 41 26 0.5 115 63 26 1 116 18 27 0.45 117 38 27 -0.9583187 118 42 27 0.5 119 63 27 1 120 19 28 0.45 121 43 28 0.5 122 63 28 1 123 15 29 -0.2070986 124 16 29 -0.2232997 125 17 29 -0.2286264 126 18 29 -0.2024528 127 19 29 -0.1385226 128 24 29 -0.2070986 129 25 29 -0.2232997 130 26 29 -0.2286264 131 27 29 -0.2024528 132 28 29 -0.1385226 133 15 30 -1.05 134 24 30 -1.05 135 59 30 1 136 16 31 -1.05 137 25 31 -1.05 138 59 31 1 139 17 32 -1.05 140 26 32 -1.05 141 59 32 1 142 18 33 -1.05 143 27 33 -1.05 144 59 33 1 145 19 34 -1.05 146 28 34 -1.05 147 59 34 1 148 24 35 0.08147449 149 25 35 0.097890154 150 26 35 0.1131608 151 27 35 0.1150555 152 28 35 0.092419088 153 30 35 0.090527207 154 31 35 0.1087668 155 32 35 0.1257342 156 33 35 0.1278394 157 34 35 0.1026879 158 20 36 1.567398 159 24 36 0.65 160 29 36 -1.567398 161 30 36 0.7222222 162 61 36 1 163 21 37 1.253918 164 25 37 0.65 165 29 37 -1.253918 166 31 37 0.7222222 167 61 37 1 168 22 38 0.9404389 169 26 38 0.65 170 29 38 -0.9404389 171 32 38 0.7222222 172 61 38 1 173 23 39 0.6269592 174 27 39 0.65 175 29 39 -0.6269592 176 33 39 0.7222222 177 61 39 1 178 28 40 0.65 179 29 40 -0.3134796 180 34 40 0.7222222 181 61 40 1 182 20 41 -0.2788416 183 21 41 -0.2680186 184 22 41 -0.2323717 185 23 41 -0.1575082 186 29 41 1 187 30 42 0.25 188 45 42 -0.8242248 189 49 42 0.5 190 65 42 1 191 31 43 0.25 192 46 43 -0.8242248 193 50 43 0.5 194 65 43 1 195 32 44 0.25 196 47 44 -0.8242248 197 51 44 0.5 198 65 44 1 199 33 45 0.25 200 48 45 -0.8242248 201 52 45 0.5 202 65 45 1 203 34 46 0.25 204 53 46 0.5 205 65 46 1 206 30 47 -0.1581626 207 31 47 -0.1947711 208 32 47 -0.2303917 209 33 47 -0.2362845 210 34 47 -0.18039 211 39 47 -0.1581626 212 40 47 -0.1947711 213 41 47 -0.2303917 214 42 47 -0.2362845 215 43 47 -0.18039 216 30 48 -0.9722222 217 39 48 -0.9722222 218 62 48 1 219 31 49 -0.9722222 220 40 49 -0.9722222 221 62 49 1 222 32 50 -0.9722222 223 41 50 -0.9722222 224 62 50 1 225 33 51 -0.9722222 226 42 51 -0.9722222 227 62 51 1 228 34 52 -0.9722222 229 43 52 -0.9722222 230 62 52 1 231 39 53 0.053228639 232 40 53 0.075745419 233 41 53 0.1061028 234 42 53 0.1333878 235 43 53 0.1315353 236 49 53 -0.1064573 237 50 53 -0.1514908 238 51 53 -0.2122056 239 52 53 -0.2667757 240 53 53 -0.2630706 241 35 54 1.863354 242 39 54 0.4722222 243 44 54 -1.863354 244 49 54 -0.9444444 245 64 54 1 246 36 55 1.490683 247 40 55 0.4722222 248 44 55 -1.490683 249 50 55 -0.9444444 250 64 55 1 251 37 56 1.118012 252 41 56 0.4722222 253 44 56 -1.118012 254 51 56 -0.9444444 255 64 56 1 256 38 57 0.7453416 257 42 57 0.4722222 258 44 57 -0.7453416 259 52 57 -0.9444444 260 64 57 1 261 43 58 0.4722222 262 44 58 -0.3726708 263 53 58 -0.9444444 264 64 58 1 265 35 59 -0.2069954 266 36 59 -0.2356469 267 37 59 -0.2475675 268 38 59 -0.2074873 269 44 59 1 270 45 60 1.863354 271 49 60 0.4444444 272 54 60 -1.863354 273 66 60 1 274 46 61 1.490683 275 50 61 0.4444444 276 54 61 -1.490683 277 66 61 1 278 47 62 1.118012 279 51 62 0.4444444 280 54 62 -1.118012 281 66 62 1 282 48 63 0.7453416 283 52 63 0.4444444 284 54 63 -0.7453416 285 66 63 1 286 53 64 0.4444444 287 54 64 -0.3726708 288 66 64 1 289 45 65 -0.1443354 290 46 65 -0.1918557 291 47 65 -0.2421498 292 48 65 -0.2541193 293 54 65 1 TEST03 Convert a sparse matrix from ST to CC format. ST: sparse triplet, I, J, A. CC: compressed column, I, CC, A. The ST matrix is read from the file 'west_st.txt' and the CC matrix is written to the files: 'west_icc.txt', 'west_ccc.txt', and 'west_acc.txt'. Sparse Triplet (ST) header information: Minimum row index I_MIN = 1 Maximum row index I_MAX = 67 Minimum col index J_MIN = 1 Maximum col index J_MAX = 67 Number of rows M = 67 Number of columns N = 67 Number of nonzeros NST = 299 Number of CC values = 294 TEST04 Convert a sparse matrix from ST to CC format. ST: sparse triplet, I, J, A. CC: compressed column, I, CC, A. The ST matrix is the Wathen finite element matrix. It has many repeated index pairs. To check, compare ACC*X - AST*X for a random X. Number of ST values = 576 Sparse Triplet (ST) header information: Minimum row index I_MIN = 0 Maximum row index I_MAX = 39 Minimum col index J_MIN = 0 Maximum col index J_MAX = 39 Number of rows M = 40 Number of columns N = 40 Number of nonzeros NST = 576 Number of CC values = 472 || ACC*X - AST*X|| = 1.43521e-12 ST_TO_CC_PRB Normal end of execution. 22 January 2017 12:05:13 PM