15 September 2021 3:23:37.251 PM ST_TO_ccs_TEST FORTRAN90 version Test st_to_ccs. TEST01 Convert a sparse matrix from ST to CCS format. ST: sparse triplet, I, J, A. CCS: 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 ---- ---- ---- -------------- 1 1 1 2.0000000 2 1 2 3.0000000 3 2 1 3.0000000 4 2 3 4.0000000 5 2 5 6.0000000 6 3 2 -1.0000000 7 3 3 -3.0000000 8 3 4 2.0000000 9 4 3 1.0000000 10 5 2 4.0000000 11 5 3 2.0000000 12 5 5 1.0000000 Number of CCS values = 12 CCS Matrix: # I J A ---- ---- ---- ---------------- 1 1 1 2.0000000 2 2 1 3.0000000 3 1 2 3.0000000 4 3 2 -1.0000000 5 5 2 4.0000000 6 2 3 4.0000000 7 3 3 -3.0000000 8 4 3 1.0000000 9 5 3 2.0000000 10 3 4 2.0000000 11 2 5 6.0000000 12 5 5 1.0000000 TEST02 Convert a sparse matrix from ST to CCS format. ST: sparse triplet, I, J, A. CCS: 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 ---- ---- ---- -------------- 1 45 56 -1.8633540 2 55 62 -1.8633540 3 30 38 -1.5673980 4 45 57 -1.4906830 5 55 63 -1.4906830 6 10 13 -1.2658230 7 30 39 -1.2539180 8 45 58 -1.1180120 9 55 64 -1.1180120 10 16 32 -1.0500000 11 17 33 -1.0500000 12 18 34 -1.0500000 13 19 35 -1.0500000 14 20 36 -1.0500000 15 25 32 -1.0500000 16 26 33 -1.0500000 17 27 34 -1.0500000 18 28 35 -1.0500000 19 29 36 -1.0500000 20 10 14 -1.0126580 21 11 21 -1.0000000 22 12 22 -1.0000000 23 13 23 -1.0000000 24 14 24 -1.0000000 25 15 25 -1.0000000 26 31 50 -0.97222220 27 32 51 -0.97222220 28 33 52 -0.97222220 29 34 53 -0.97222220 30 35 54 -0.97222220 31 40 50 -0.97222220 32 41 51 -0.97222220 33 42 52 -0.97222220 34 43 53 -0.97222220 35 44 54 -0.97222220 36 36 26 -0.95831870 37 37 27 -0.95831870 38 38 28 -0.95831870 39 39 29 -0.95831870 40 50 56 -0.94444440 41 51 57 -0.94444440 42 52 58 -0.94444440 43 53 59 -0.94444440 44 54 60 -0.94444440 45 30 40 -0.94043890 46 21 2 -0.91595330 47 22 3 -0.91595330 48 23 4 -0.91595330 49 24 5 -0.91595330 50 1 8 -0.83418180 51 2 9 -0.83418180 52 3 10 -0.83418180 53 4 11 -0.83418180 54 46 44 -0.82422480 55 47 45 -0.82422480 56 48 46 -0.82422480 57 49 47 -0.82422480 58 5 2 -0.80000000 59 6 3 -0.80000000 60 7 4 -0.80000000 61 8 5 -0.80000000 62 9 6 -0.80000000 63 10 15 -0.75949370 64 45 59 -0.74534160 65 55 65 -0.74534160 66 30 41 -0.62695920 67 10 16 -0.50632910 68 45 60 -0.37267080 69 55 66 -0.37267080 70 1 18 -0.33615560 71 30 42 -0.31347960 72 2 18 -0.29391960 73 21 43 -0.27884160 74 5 1 -0.27884160 75 22 43 -0.26801860 76 6 1 -0.26801860 77 53 55 -0.26677570 78 54 55 -0.26307060 79 49 67 -0.25411930 80 10 17 -0.25316460 81 38 61 -0.24756750 82 48 67 -0.24214980 83 34 49 -0.23628450 84 43 49 -0.23628450 85 37 61 -0.23564690 86 23 43 -0.23237170 87 7 1 -0.23237170 88 33 49 -0.23039170 89 42 49 -0.23039170 90 18 31 -0.22862640 91 27 31 -0.22862640 92 17 31 -0.22329970 93 26 31 -0.22329970 94 3 18 -0.22148150 95 13 20 -0.21442060 96 12 20 -0.21403920 97 52 55 -0.21220560 98 39 61 -0.20748730 99 11 20 -0.20717590 100 16 31 -0.20709860 101 25 31 -0.20709860 102 36 61 -0.20699540 103 19 31 -0.20245280 104 28 31 -0.20245280 105 14 20 -0.19867680 106 32 49 -0.19477110 107 41 49 -0.19477110 108 47 67 -0.19185570 109 35 49 -0.18039000 110 44 49 -0.18039000 111 15 20 -0.16568740 112 31 49 -0.15816260 113 40 49 -0.15816260 114 24 43 -0.15750820 115 8 1 -0.15750820 116 51 55 -0.15149080 117 46 67 -0.14433540 118 20 31 -0.13852260 119 29 31 -0.13852260 120 4 18 -0.11898600 121 50 55 -0.10645730 122 9 1 -0.63259780E-01 123 9 7 0.11782910E-01 124 29 1 0.31629890E-01 125 8 7 0.47594391E-01 126 40 55 0.53228639E-01 127 41 55 0.75745419E-01 128 28 1 0.78754112E-01 129 25 37 0.81474490E-01 130 7 7 0.88592619E-01 131 31 37 0.90527207E-01 132 29 37 0.92419088E-01 133 26 37 0.97890154E-01 134 20 20 0.99412464E-01 135 35 37 0.10268790 136 42 55 0.10610280 137 32 37 0.10876680 138 27 37 0.11316080 139 28 37 0.11505550 140 27 1 0.11618590 141 6 7 0.11756790 142 19 20 0.11920610 143 16 20 0.12430550 144 33 37 0.12573420 145 34 37 0.12783940 146 17 20 0.12842350 147 18 20 0.12865240 148 44 55 0.13153530 149 43 55 0.13338780 150 26 1 0.13400930 151 5 7 0.13446220 152 25 1 0.13942080 153 31 44 0.25000000 154 32 45 0.25000000 155 33 46 0.25000000 156 34 47 0.25000000 157 35 48 0.25000000 158 11 13 0.33333330 159 12 14 0.33333330 160 13 15 0.33333330 161 14 16 0.33333330 162 15 17 0.33333330 163 25 2 0.40000000 164 26 3 0.40000000 165 27 4 0.40000000 166 28 5 0.40000000 167 29 6 0.40000000 168 5 13 0.40000000 169 5 8 0.40000000 170 6 14 0.40000000 171 6 9 0.40000000 172 7 15 0.40000000 173 7 10 0.40000000 174 8 11 0.40000000 175 8 16 0.40000000 176 9 12 0.40000000 177 9 17 0.40000000 178 50 62 0.44444440 179 51 63 0.44444440 180 52 64 0.44444440 181 53 65 0.44444440 182 54 66 0.44444440 183 16 26 0.45000000 184 17 27 0.45000000 185 18 28 0.45000000 186 19 29 0.45000000 187 20 30 0.45000000 188 40 56 0.47222220 189 41 57 0.47222220 190 42 58 0.47222220 191 43 59 0.47222220 192 44 60 0.47222220 193 40 26 0.50000000 194 41 27 0.50000000 195 42 28 0.50000000 196 43 29 0.50000000 197 44 30 0.50000000 198 50 44 0.50000000 199 51 45 0.50000000 200 52 46 0.50000000 201 53 47 0.50000000 202 54 48 0.50000000 203 4 16 0.50632910 204 16 21 0.60000000 205 17 22 0.60000000 206 18 23 0.60000000 207 19 24 0.60000000 208 20 25 0.60000000 209 24 41 0.62695920 210 25 38 0.65000000 211 26 39 0.65000000 212 27 40 0.65000000 213 28 41 0.65000000 214 29 42 0.65000000 215 15 19 0.66666670 216 31 38 0.72222220 217 32 39 0.72222220 218 33 40 0.72222220 219 34 41 0.72222220 220 35 42 0.72222220 221 39 59 0.74534160 222 49 65 0.74534160 223 3 15 0.75949370 224 23 40 0.94043890 225 60 32 0.50000000 226 60 33 0.50000000 227 60 34 0.50000000 228 60 35 0.50000000 229 60 36 0.50000000 230 30 43 1.0000000 231 45 61 1.0000000 232 55 67 1.0000000 233 56 19 1.0000000 234 57 11 1.0000000 235 57 12 1.0000000 236 57 8 1.0000000 237 57 9 1.0000000 238 57 10 1.0000000 239 58 13 1.0000000 240 58 14 1.0000000 241 58 15 1.0000000 242 58 16 1.0000000 243 58 17 1.0000000 244 59 21 1.0000000 245 59 22 1.0000000 246 59 23 1.0000000 247 59 24 1.0000000 248 59 25 1.0000000 249 60 32 0.50000000 250 60 33 0.50000000 251 60 34 0.50000000 252 60 35 0.50000000 253 60 36 0.50000000 254 61 2 1.0000000 255 61 3 1.0000000 256 61 4 1.0000000 257 61 5 1.0000000 258 61 6 1.0000000 259 62 38 1.0000000 260 62 39 1.0000000 261 62 40 1.0000000 262 62 41 1.0000000 263 62 42 1.0000000 264 63 50 1.0000000 265 63 51 1.0000000 266 63 52 1.0000000 267 63 53 1.0000000 268 63 54 1.0000000 269 64 26 1.0000000 270 64 27 1.0000000 271 64 28 1.0000000 272 64 29 1.0000000 273 64 30 1.0000000 274 65 56 1.0000000 275 65 57 1.0000000 276 65 58 1.0000000 277 65 59 1.0000000 278 65 60 1.0000000 279 66 44 1.0000000 280 66 45 1.0000000 281 66 46 1.0000000 282 66 47 1.0000000 283 66 48 1.0000000 284 67 62 1.0000000 285 67 63 1.0000000 286 67 64 1.0000000 287 67 65 1.0000000 288 67 66 1.0000000 289 10 18 1.0000000 290 2 14 1.0126580 291 38 58 1.1180120 292 48 64 1.1180120 293 22 39 1.2539180 294 1 13 1.2658230 295 37 57 1.4906830 296 47 63 1.4906830 297 21 38 1.5673980 298 36 56 1.8633540 299 46 62 1.8633540 Number of CCS values = 294 CCS Matrix: # I J A ---- ---- ---- ---------------- 1 5 1 -0.27884160 2 6 1 -0.26801860 3 7 1 -0.23237170 4 8 1 -0.15750820 5 9 1 -0.63259780E-01 6 25 1 0.13942080 7 26 1 0.13400930 8 27 1 0.11618590 9 28 1 0.78754112E-01 10 29 1 0.31629890E-01 11 5 2 -0.80000000 12 21 2 -0.91595330 13 25 2 0.40000000 14 61 2 1.0000000 15 6 3 -0.80000000 16 22 3 -0.91595330 17 26 3 0.40000000 18 61 3 1.0000000 19 7 4 -0.80000000 20 23 4 -0.91595330 21 27 4 0.40000000 22 61 4 1.0000000 23 8 5 -0.80000000 24 24 5 -0.91595330 25 28 5 0.40000000 26 61 5 1.0000000 27 9 6 -0.80000000 28 29 6 0.40000000 29 61 6 1.0000000 30 5 7 0.13446220 31 6 7 0.11756790 32 7 7 0.88592619E-01 33 8 7 0.47594391E-01 34 9 7 0.11782910E-01 35 1 8 -0.83418180 36 5 8 0.40000000 37 57 8 1.0000000 38 2 9 -0.83418180 39 6 9 0.40000000 40 57 9 1.0000000 41 3 10 -0.83418180 42 7 10 0.40000000 43 57 10 1.0000000 44 4 11 -0.83418180 45 8 11 0.40000000 46 57 11 1.0000000 47 9 12 0.40000000 48 57 12 1.0000000 49 1 13 1.2658230 50 5 13 0.40000000 51 10 13 -1.2658230 52 11 13 0.33333330 53 58 13 1.0000000 54 2 14 1.0126580 55 6 14 0.40000000 56 10 14 -1.0126580 57 12 14 0.33333330 58 58 14 1.0000000 59 3 15 0.75949370 60 7 15 0.40000000 61 10 15 -0.75949370 62 13 15 0.33333330 63 58 15 1.0000000 64 4 16 0.50632910 65 8 16 0.40000000 66 10 16 -0.50632910 67 14 16 0.33333330 68 58 16 1.0000000 69 9 17 0.40000000 70 10 17 -0.25316460 71 15 17 0.33333330 72 58 17 1.0000000 73 1 18 -0.33615560 74 2 18 -0.29391960 75 3 18 -0.22148150 76 4 18 -0.11898600 77 10 18 1.0000000 78 15 19 0.66666670 79 56 19 1.0000000 80 11 20 -0.20717590 81 12 20 -0.21403920 82 13 20 -0.21442060 83 14 20 -0.19867680 84 15 20 -0.16568740 85 16 20 0.12430550 86 17 20 0.12842350 87 18 20 0.12865240 88 19 20 0.11920610 89 20 20 0.99412464E-01 90 11 21 -1.0000000 91 16 21 0.60000000 92 59 21 1.0000000 93 12 22 -1.0000000 94 17 22 0.60000000 95 59 22 1.0000000 96 13 23 -1.0000000 97 18 23 0.60000000 98 59 23 1.0000000 99 14 24 -1.0000000 100 19 24 0.60000000 101 59 24 1.0000000 102 15 25 -1.0000000 103 20 25 0.60000000 104 59 25 1.0000000 105 16 26 0.45000000 106 36 26 -0.95831870 107 40 26 0.50000000 108 64 26 1.0000000 109 17 27 0.45000000 110 37 27 -0.95831870 111 41 27 0.50000000 112 64 27 1.0000000 113 18 28 0.45000000 114 38 28 -0.95831870 115 42 28 0.50000000 116 64 28 1.0000000 117 19 29 0.45000000 118 39 29 -0.95831870 119 43 29 0.50000000 120 64 29 1.0000000 121 20 30 0.45000000 122 44 30 0.50000000 123 64 30 1.0000000 124 16 31 -0.20709860 125 17 31 -0.22329970 126 18 31 -0.22862640 127 19 31 -0.20245280 128 20 31 -0.13852260 129 25 31 -0.20709860 130 26 31 -0.22329970 131 27 31 -0.22862640 132 28 31 -0.20245280 133 29 31 -0.13852260 134 16 32 -1.0500000 135 25 32 -1.0500000 136 60 32 1.0000000 137 17 33 -1.0500000 138 26 33 -1.0500000 139 60 33 1.0000000 140 18 34 -1.0500000 141 27 34 -1.0500000 142 60 34 1.0000000 143 19 35 -1.0500000 144 28 35 -1.0500000 145 60 35 1.0000000 146 20 36 -1.0500000 147 29 36 -1.0500000 148 60 36 1.0000000 149 25 37 0.81474490E-01 150 26 37 0.97890154E-01 151 27 37 0.11316080 152 28 37 0.11505550 153 29 37 0.92419088E-01 154 31 37 0.90527207E-01 155 32 37 0.10876680 156 33 37 0.12573420 157 34 37 0.12783940 158 35 37 0.10268790 159 21 38 1.5673980 160 25 38 0.65000000 161 30 38 -1.5673980 162 31 38 0.72222220 163 62 38 1.0000000 164 22 39 1.2539180 165 26 39 0.65000000 166 30 39 -1.2539180 167 32 39 0.72222220 168 62 39 1.0000000 169 23 40 0.94043890 170 27 40 0.65000000 171 30 40 -0.94043890 172 33 40 0.72222220 173 62 40 1.0000000 174 24 41 0.62695920 175 28 41 0.65000000 176 30 41 -0.62695920 177 34 41 0.72222220 178 62 41 1.0000000 179 29 42 0.65000000 180 30 42 -0.31347960 181 35 42 0.72222220 182 62 42 1.0000000 183 21 43 -0.27884160 184 22 43 -0.26801860 185 23 43 -0.23237170 186 24 43 -0.15750820 187 30 43 1.0000000 188 31 44 0.25000000 189 46 44 -0.82422480 190 50 44 0.50000000 191 66 44 1.0000000 192 32 45 0.25000000 193 47 45 -0.82422480 194 51 45 0.50000000 195 66 45 1.0000000 196 33 46 0.25000000 197 48 46 -0.82422480 198 52 46 0.50000000 199 66 46 1.0000000 200 34 47 0.25000000 201 49 47 -0.82422480 202 53 47 0.50000000 203 66 47 1.0000000 204 35 48 0.25000000 205 54 48 0.50000000 206 66 48 1.0000000 207 31 49 -0.15816260 208 32 49 -0.19477110 209 33 49 -0.23039170 210 34 49 -0.23628450 211 35 49 -0.18039000 212 40 49 -0.15816260 213 41 49 -0.19477110 214 42 49 -0.23039170 215 43 49 -0.23628450 216 44 49 -0.18039000 217 31 50 -0.97222220 218 40 50 -0.97222220 219 63 50 1.0000000 220 32 51 -0.97222220 221 41 51 -0.97222220 222 63 51 1.0000000 223 33 52 -0.97222220 224 42 52 -0.97222220 225 63 52 1.0000000 226 34 53 -0.97222220 227 43 53 -0.97222220 228 63 53 1.0000000 229 35 54 -0.97222220 230 44 54 -0.97222220 231 63 54 1.0000000 232 40 55 0.53228639E-01 233 41 55 0.75745419E-01 234 42 55 0.10610280 235 43 55 0.13338780 236 44 55 0.13153530 237 50 55 -0.10645730 238 51 55 -0.15149080 239 52 55 -0.21220560 240 53 55 -0.26677570 241 54 55 -0.26307060 242 36 56 1.8633540 243 40 56 0.47222220 244 45 56 -1.8633540 245 50 56 -0.94444440 246 65 56 1.0000000 247 37 57 1.4906830 248 41 57 0.47222220 249 45 57 -1.4906830 250 51 57 -0.94444440 251 65 57 1.0000000 252 38 58 1.1180120 253 42 58 0.47222220 254 45 58 -1.1180120 255 52 58 -0.94444440 256 65 58 1.0000000 257 39 59 0.74534160 258 43 59 0.47222220 259 45 59 -0.74534160 260 53 59 -0.94444440 261 65 59 1.0000000 262 44 60 0.47222220 263 45 60 -0.37267080 264 54 60 -0.94444440 265 65 60 1.0000000 266 36 61 -0.20699540 267 37 61 -0.23564690 268 38 61 -0.24756750 269 39 61 -0.20748730 270 45 61 1.0000000 271 46 62 1.8633540 272 50 62 0.44444440 273 55 62 -1.8633540 274 67 62 1.0000000 275 47 63 1.4906830 276 51 63 0.44444440 277 55 63 -1.4906830 278 67 63 1.0000000 279 48 64 1.1180120 280 52 64 0.44444440 281 55 64 -1.1180120 282 67 64 1.0000000 283 49 65 0.74534160 284 53 65 0.44444440 285 55 65 -0.74534160 286 67 65 1.0000000 287 54 66 0.44444440 288 55 66 -0.37267080 289 67 66 1.0000000 290 46 67 -0.14433540 291 47 67 -0.19185570 292 48 67 -0.24214980 293 49 67 -0.25411930 294 55 67 1.0000000 TEST03 Convert a sparse matrix from ST to CCS format. ST: sparse triplet, I, J, A. CCS: compressed column, I, CC, A. The ST matrix is read from the file "west_st.txt", and the CCS 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 CCS values = 294 TEST04 Convert a sparse matrix from ST to CCS format. ST: sparse triplet, I, J, A. CCS: 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 = 1 Maximum row index I_MAX = 40 Minimum col index J_MIN = 1 Maximum col index J_MAX = 40 Number of rows M = 40 Number of columns N = 40 Number of nonzeros NST = 576 Number of CCS values = 472 || ACC*X - AST*X|| = 0.241707E-11 st_to_ccs_TEST Normal end of execution. 15 September 2021 3:23:37.254 PM