1*2db2b1b2SPierre Jolivetoriginal A: Frobenius norm = 320.617, one norm = 273., infinity norm = 585. 2c4762a1bSJed Brown type: seqaij 38e160fb9SStefano Zampinirow 0: (0, 1.) (1, 2.) (2, 3.) (3, 4.) (4, 5.) (5, 6.) (6, 7.) (7, 8.) (8, 9.) 48e160fb9SStefano Zampinirow 1: (0, 11.) (1, 12.) (2, 13.) (3, 14.) (4, 15.) (5, 16.) (6, 17.) (7, 18.) (8, 19.) 58e160fb9SStefano Zampinirow 2: (0, 21.) (1, 22.) (2, 23.) (3, 24.) (4, 25.) (5, 26.) (6, 27.) (7, 28.) (8, 29.) 68e160fb9SStefano Zampinirow 3: (0, 31.) (1, 32.) (2, 33.) (3, 34.) (4, 35.) (5, 36.) (6, 37.) (7, 38.) (8, 39.) 78e160fb9SStefano Zampinirow 4: (0, 41.) (1, 42.) (2, 43.) (3, 44.) (4, 45.) (5, 46.) (6, 47.) (7, 48.) (8, 49.) 88e160fb9SStefano Zampinirow 5: (0, 51.) (1, 52.) (2, 53.) (3, 54.) (4, 55.) (5, 56.) (6, 57.) (7, 58.) (8, 59.) 98e160fb9SStefano Zampinirow 6: (0, 61.) (1, 62.) (2, 63.) (3, 64.) (4, 65.) (5, 66.) (6, 67.) (7, 68.) (8, 69.) 10*2db2b1b2SPierre JolivetB = A^T: Frobenius norm = 320.617, one norm = 585., infinity norm = 273. 11c4762a1bSJed Brown type: seqaij 128e160fb9SStefano Zampinirow 0: (0, 1.) (1, 11.) (2, 21.) (3, 31.) (4, 41.) (5, 51.) (6, 61.) 138e160fb9SStefano Zampinirow 1: (0, 2.) (1, 12.) (2, 22.) (3, 32.) (4, 42.) (5, 52.) (6, 62.) 148e160fb9SStefano Zampinirow 2: (0, 3.) (1, 13.) (2, 23.) (3, 33.) (4, 43.) (5, 53.) (6, 63.) 158e160fb9SStefano Zampinirow 3: (0, 4.) (1, 14.) (2, 24.) (3, 34.) (4, 44.) (5, 54.) (6, 64.) 168e160fb9SStefano Zampinirow 4: (0, 5.) (1, 15.) (2, 25.) (3, 35.) (4, 45.) (5, 55.) (6, 65.) 178e160fb9SStefano Zampinirow 5: (0, 6.) (1, 16.) (2, 26.) (3, 36.) (4, 46.) (5, 56.) (6, 66.) 188e160fb9SStefano Zampinirow 6: (0, 7.) (1, 17.) (2, 27.) (3, 37.) (4, 47.) (5, 57.) (6, 67.) 198e160fb9SStefano Zampinirow 7: (0, 8.) (1, 18.) (2, 28.) (3, 38.) (4, 48.) (5, 58.) (6, 68.) 208e160fb9SStefano Zampinirow 8: (0, 9.) (1, 19.) (2, 29.) (3, 39.) (4, 49.) (5, 59.) (6, 69.) 21c4762a1bSJed BrownMatAXPY: C = C + alpha * A, C=A, SAME_NONZERO_PATTERN 22c4762a1bSJed Brown type: seqaij 238e160fb9SStefano Zampinirow 0: (0, 2.) (1, 4.) (2, 6.) (3, 8.) (4, 10.) (5, 12.) (6, 14.) (7, 16.) (8, 18.) 248e160fb9SStefano Zampinirow 1: (0, 22.) (1, 24.) (2, 26.) (3, 28.) (4, 30.) (5, 32.) (6, 34.) (7, 36.) (8, 38.) 258e160fb9SStefano Zampinirow 2: (0, 42.) (1, 44.) (2, 46.) (3, 48.) (4, 50.) (5, 52.) (6, 54.) (7, 56.) (8, 58.) 268e160fb9SStefano Zampinirow 3: (0, 62.) (1, 64.) (2, 66.) (3, 68.) (4, 70.) (5, 72.) (6, 74.) (7, 76.) (8, 78.) 278e160fb9SStefano Zampinirow 4: (0, 82.) (1, 84.) (2, 86.) (3, 88.) (4, 90.) (5, 92.) (6, 94.) (7, 96.) (8, 98.) 288e160fb9SStefano Zampinirow 5: (0, 102.) (1, 104.) (2, 106.) (3, 108.) (4, 110.) (5, 112.) (6, 114.) (7, 116.) (8, 118.) 298e160fb9SStefano Zampinirow 6: (0, 122.) (1, 124.) (2, 126.) (3, 128.) (4, 130.) (5, 132.) (6, 134.) (7, 136.) (8, 138.) 308e160fb9SStefano Zampini 31c4762a1bSJed BrownMatAXPY: (C^T)^T = (C^T)^T + alpha * A, C=A, SAME_NONZERO_PATTERN 32c4762a1bSJed Brown type: seqaij 338e160fb9SStefano Zampinirow 0: (0, 2.) (1, 4.) (2, 6.) (3, 8.) (4, 10.) (5, 12.) (6, 14.) (7, 16.) (8, 18.) 348e160fb9SStefano Zampinirow 1: (0, 22.) (1, 24.) (2, 26.) (3, 28.) (4, 30.) (5, 32.) (6, 34.) (7, 36.) (8, 38.) 358e160fb9SStefano Zampinirow 2: (0, 42.) (1, 44.) (2, 46.) (3, 48.) (4, 50.) (5, 52.) (6, 54.) (7, 56.) (8, 58.) 368e160fb9SStefano Zampinirow 3: (0, 62.) (1, 64.) (2, 66.) (3, 68.) (4, 70.) (5, 72.) (6, 74.) (7, 76.) (8, 78.) 378e160fb9SStefano Zampinirow 4: (0, 82.) (1, 84.) (2, 86.) (3, 88.) (4, 90.) (5, 92.) (6, 94.) (7, 96.) (8, 98.) 388e160fb9SStefano Zampinirow 5: (0, 102.) (1, 104.) (2, 106.) (3, 108.) (4, 110.) (5, 112.) (6, 114.) (7, 116.) (8, 118.) 398e160fb9SStefano Zampinirow 6: (0, 122.) (1, 124.) (2, 126.) (3, 128.) (4, 130.) (5, 132.) (6, 134.) (7, 136.) (8, 138.) 40c4762a1bSJed BrownMatAXPY: C = C + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN 41c4762a1bSJed Brown type: seqaij 428e160fb9SStefano Zampinirow 0: (0, 2.) (1, 4.) (2, 6.) (3, 8.) (4, 10.) (5, 12.) (6, 14.) (7, 16.) (8, 18.) 438e160fb9SStefano Zampinirow 1: (0, 22.) (1, 24.) (2, 26.) (3, 28.) (4, 30.) (5, 32.) (6, 34.) (7, 36.) (8, 38.) 448e160fb9SStefano Zampinirow 2: (0, 42.) (1, 44.) (2, 46.) (3, 48.) (4, 50.) (5, 52.) (6, 54.) (7, 56.) (8, 58.) 458e160fb9SStefano Zampinirow 3: (0, 62.) (1, 64.) (2, 66.) (3, 68.) (4, 70.) (5, 72.) (6, 74.) (7, 76.) (8, 78.) 468e160fb9SStefano Zampinirow 4: (0, 82.) (1, 84.) (2, 86.) (3, 88.) (4, 90.) (5, 92.) (6, 94.) (7, 96.) (8, 98.) 478e160fb9SStefano Zampinirow 5: (0, 102.) (1, 104.) (2, 106.) (3, 108.) (4, 110.) (5, 112.) (6, 114.) (7, 116.) (8, 118.) 488e160fb9SStefano Zampinirow 6: (0, 122.) (1, 124.) (2, 126.) (3, 128.) (4, 130.) (5, 132.) (6, 134.) (7, 136.) (8, 138.) 49c4762a1bSJed BrownMatAXPY: (C^T)^T = (C^T)^T + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN 50c4762a1bSJed Brown type: seqaij 518e160fb9SStefano Zampinirow 0: (0, 2.) (1, 4.) (2, 6.) (3, 8.) (4, 10.) (5, 12.) (6, 14.) (7, 16.) (8, 18.) 528e160fb9SStefano Zampinirow 1: (0, 22.) (1, 24.) (2, 26.) (3, 28.) (4, 30.) (5, 32.) (6, 34.) (7, 36.) (8, 38.) 538e160fb9SStefano Zampinirow 2: (0, 42.) (1, 44.) (2, 46.) (3, 48.) (4, 50.) (5, 52.) (6, 54.) (7, 56.) (8, 58.) 548e160fb9SStefano Zampinirow 3: (0, 62.) (1, 64.) (2, 66.) (3, 68.) (4, 70.) (5, 72.) (6, 74.) (7, 76.) (8, 78.) 558e160fb9SStefano Zampinirow 4: (0, 82.) (1, 84.) (2, 86.) (3, 88.) (4, 90.) (5, 92.) (6, 94.) (7, 96.) (8, 98.) 568e160fb9SStefano Zampinirow 5: (0, 102.) (1, 104.) (2, 106.) (3, 108.) (4, 110.) (5, 112.) (6, 114.) (7, 116.) (8, 118.) 578e160fb9SStefano Zampinirow 6: (0, 122.) (1, 124.) (2, 126.) (3, 128.) (4, 130.) (5, 132.) (6, 134.) (7, 136.) (8, 138.) 588e160fb9SStefano ZampiniMatAXPY: Now without the transposition operation 598e160fb9SStefano Zampini type: seqaij 608e160fb9SStefano Zampinirow 0: (0, 2.) (1, 4.) (2, 6.) (3, 8.) (4, 10.) (5, 12.) (6, 14.) (7, 16.) (8, 18.) 618e160fb9SStefano Zampinirow 1: (0, 22.) (1, 24.) (2, 26.) (3, 28.) (4, 30.) (5, 32.) (6, 34.) (7, 36.) (8, 38.) 628e160fb9SStefano Zampinirow 2: (0, 42.) (1, 44.) (2, 46.) (3, 48.) (4, 50.) (5, 52.) (6, 54.) (7, 56.) (8, 58.) 638e160fb9SStefano Zampinirow 3: (0, 62.) (1, 64.) (2, 66.) (3, 68.) (4, 70.) (5, 72.) (6, 74.) (7, 76.) (8, 78.) 648e160fb9SStefano Zampinirow 4: (0, 82.) (1, 84.) (2, 86.) (3, 88.) (4, 90.) (5, 92.) (6, 94.) (7, 96.) (8, 98.) 658e160fb9SStefano Zampinirow 5: (0, 102.) (1, 104.) (2, 106.) (3, 108.) (4, 110.) (5, 112.) (6, 114.) (7, 116.) (8, 118.) 668e160fb9SStefano Zampinirow 6: (0, 122.) (1, 124.) (2, 126.) (3, 128.) (4, 130.) (5, 132.) (6, 134.) (7, 136.) (8, 138.) 678e160fb9SStefano Zampini 68c4762a1bSJed BrownMatAXPY: (C^H)^H = (C^H)^H + alpha * A, C=A, SAME_NONZERO_PATTERN 69c4762a1bSJed Brown type: seqaij 708e160fb9SStefano Zampinirow 0: (0, 2.) (1, 4.) (2, 6.) (3, 8.) (4, 10.) (5, 12.) (6, 14.) (7, 16.) (8, 18.) 718e160fb9SStefano Zampinirow 1: (0, 22.) (1, 24.) (2, 26.) (3, 28.) (4, 30.) (5, 32.) (6, 34.) (7, 36.) (8, 38.) 728e160fb9SStefano Zampinirow 2: (0, 42.) (1, 44.) (2, 46.) (3, 48.) (4, 50.) (5, 52.) (6, 54.) (7, 56.) (8, 58.) 738e160fb9SStefano Zampinirow 3: (0, 62.) (1, 64.) (2, 66.) (3, 68.) (4, 70.) (5, 72.) (6, 74.) (7, 76.) (8, 78.) 748e160fb9SStefano Zampinirow 4: (0, 82.) (1, 84.) (2, 86.) (3, 88.) (4, 90.) (5, 92.) (6, 94.) (7, 96.) (8, 98.) 758e160fb9SStefano Zampinirow 5: (0, 102.) (1, 104.) (2, 106.) (3, 108.) (4, 110.) (5, 112.) (6, 114.) (7, 116.) (8, 118.) 768e160fb9SStefano Zampinirow 6: (0, 122.) (1, 124.) (2, 126.) (3, 128.) (4, 130.) (5, 132.) (6, 134.) (7, 136.) (8, 138.) 77c4762a1bSJed BrownMatAXPY: C = C + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN 78c4762a1bSJed Brown type: seqaij 798e160fb9SStefano Zampinirow 0: (0, 2.) (1, 4.) (2, 6.) (3, 8.) (4, 10.) (5, 12.) (6, 14.) (7, 16.) (8, 18.) 808e160fb9SStefano Zampinirow 1: (0, 22.) (1, 24.) (2, 26.) (3, 28.) (4, 30.) (5, 32.) (6, 34.) (7, 36.) (8, 38.) 818e160fb9SStefano Zampinirow 2: (0, 42.) (1, 44.) (2, 46.) (3, 48.) (4, 50.) (5, 52.) (6, 54.) (7, 56.) (8, 58.) 828e160fb9SStefano Zampinirow 3: (0, 62.) (1, 64.) (2, 66.) (3, 68.) (4, 70.) (5, 72.) (6, 74.) (7, 76.) (8, 78.) 838e160fb9SStefano Zampinirow 4: (0, 82.) (1, 84.) (2, 86.) (3, 88.) (4, 90.) (5, 92.) (6, 94.) (7, 96.) (8, 98.) 848e160fb9SStefano Zampinirow 5: (0, 102.) (1, 104.) (2, 106.) (3, 108.) (4, 110.) (5, 112.) (6, 114.) (7, 116.) (8, 118.) 858e160fb9SStefano Zampinirow 6: (0, 122.) (1, 124.) (2, 126.) (3, 128.) (4, 130.) (5, 132.) (6, 134.) (7, 136.) (8, 138.) 86c4762a1bSJed BrownMatAXPY: (C^H)^H = (C^H)^H + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN 87c4762a1bSJed Brown type: seqaij 888e160fb9SStefano Zampinirow 0: (0, 2.) (1, 4.) (2, 6.) (3, 8.) (4, 10.) (5, 12.) (6, 14.) (7, 16.) (8, 18.) 898e160fb9SStefano Zampinirow 1: (0, 22.) (1, 24.) (2, 26.) (3, 28.) (4, 30.) (5, 32.) (6, 34.) (7, 36.) (8, 38.) 908e160fb9SStefano Zampinirow 2: (0, 42.) (1, 44.) (2, 46.) (3, 48.) (4, 50.) (5, 52.) (6, 54.) (7, 56.) (8, 58.) 918e160fb9SStefano Zampinirow 3: (0, 62.) (1, 64.) (2, 66.) (3, 68.) (4, 70.) (5, 72.) (6, 74.) (7, 76.) (8, 78.) 928e160fb9SStefano Zampinirow 4: (0, 82.) (1, 84.) (2, 86.) (3, 88.) (4, 90.) (5, 92.) (6, 94.) (7, 96.) (8, 98.) 938e160fb9SStefano Zampinirow 5: (0, 102.) (1, 104.) (2, 106.) (3, 108.) (4, 110.) (5, 112.) (6, 114.) (7, 116.) (8, 118.) 948e160fb9SStefano Zampinirow 6: (0, 122.) (1, 124.) (2, 126.) (3, 128.) (4, 130.) (5, 132.) (6, 134.) (7, 136.) (8, 138.) 958e160fb9SStefano ZampiniMatAXPY: Now without the transposition operation 968e160fb9SStefano Zampini type: seqaij 978e160fb9SStefano Zampinirow 0: (0, 2.) (1, 4.) (2, 6.) (3, 8.) (4, 10.) (5, 12.) (6, 14.) (7, 16.) (8, 18.) 988e160fb9SStefano Zampinirow 1: (0, 22.) (1, 24.) (2, 26.) (3, 28.) (4, 30.) (5, 32.) (6, 34.) (7, 36.) (8, 38.) 998e160fb9SStefano Zampinirow 2: (0, 42.) (1, 44.) (2, 46.) (3, 48.) (4, 50.) (5, 52.) (6, 54.) (7, 56.) (8, 58.) 1008e160fb9SStefano Zampinirow 3: (0, 62.) (1, 64.) (2, 66.) (3, 68.) (4, 70.) (5, 72.) (6, 74.) (7, 76.) (8, 78.) 1018e160fb9SStefano Zampinirow 4: (0, 82.) (1, 84.) (2, 86.) (3, 88.) (4, 90.) (5, 92.) (6, 94.) (7, 96.) (8, 98.) 1028e160fb9SStefano Zampinirow 5: (0, 102.) (1, 104.) (2, 106.) (3, 108.) (4, 110.) (5, 112.) (6, 114.) (7, 116.) (8, 118.) 1038e160fb9SStefano Zampinirow 6: (0, 122.) (1, 124.) (2, 126.) (3, 128.) (4, 130.) (5, 132.) (6, 134.) (7, 136.) (8, 138.) 104c4762a1bSJed Brown A: original matrix: 105c4762a1bSJed Brown type: seqaij 1068e160fb9SStefano Zampinirow 0: (0, 1.) (1, 2.) (2, 3.) (3, 4.) (4, 5.) (5, 6.) (6, 7.) (7, 8.) (8, 9.) 1078e160fb9SStefano Zampinirow 1: (0, 11.) (1, 12.) (2, 13.) (3, 14.) (4, 15.) (5, 16.) (6, 17.) (7, 18.) (8, 19.) 1088e160fb9SStefano Zampinirow 2: (0, 21.) (1, 22.) (2, 23.) (3, 24.) (4, 25.) (5, 26.) (6, 27.) (7, 28.) (8, 29.) 1098e160fb9SStefano Zampinirow 3: (0, 31.) (1, 32.) (2, 33.) (3, 34.) (4, 35.) (5, 36.) (6, 37.) (7, 38.) (8, 39.) 1108e160fb9SStefano Zampinirow 4: (0, 41.) (1, 42.) (2, 43.) (3, 44.) (4, 45.) (5, 46.) (6, 47.) (7, 48.) (8, 49.) 1118e160fb9SStefano Zampinirow 5: (0, 51.) (1, 52.) (2, 53.) (3, 54.) (4, 55.) (5, 56.) (6, 57.) (7, 58.) (8, 59.) 1128e160fb9SStefano Zampinirow 6: (0, 61.) (1, 62.) (2, 63.) (3, 64.) (4, 65.) (5, 66.) (6, 67.) (7, 68.) (8, 69.) 113c4762a1bSJed Brown B(a subset of A): 114c4762a1bSJed Brown type: seqaij 1158e160fb9SStefano Zampinirow 0: (0, 1.) (2, 3.) (4, 5.) (6, 7.) (8, 9.) 116c4762a1bSJed Brownrow 1: 1178e160fb9SStefano Zampinirow 2: (0, 21.) (2, 23.) (4, 25.) (6, 27.) (8, 29.) 118c4762a1bSJed Brownrow 3: 1198e160fb9SStefano Zampinirow 4: (0, 41.) (2, 43.) (4, 45.) (6, 47.) (8, 49.) 120c4762a1bSJed Brownrow 5: 1218e160fb9SStefano Zampinirow 6: (0, 61.) (2, 63.) (4, 65.) (6, 67.) (8, 69.) 122c4762a1bSJed BrownMatAXPY: B = B + alpha * A, SUBSET_NONZERO_PATTERN 123c4762a1bSJed Brown type: seqaij 1248e160fb9SStefano Zampinirow 0: (0, 2.) (1, 2.) (2, 6.) (3, 4.) (4, 10.) (5, 6.) (6, 14.) (7, 8.) (8, 18.) 1258e160fb9SStefano Zampinirow 1: (0, 11.) (1, 12.) (2, 13.) (3, 14.) (4, 15.) (5, 16.) (6, 17.) (7, 18.) (8, 19.) 1268e160fb9SStefano Zampinirow 2: (0, 42.) (1, 22.) (2, 46.) (3, 24.) (4, 50.) (5, 26.) (6, 54.) (7, 28.) (8, 58.) 1278e160fb9SStefano Zampinirow 3: (0, 31.) (1, 32.) (2, 33.) (3, 34.) (4, 35.) (5, 36.) (6, 37.) (7, 38.) (8, 39.) 1288e160fb9SStefano Zampinirow 4: (0, 82.) (1, 42.) (2, 86.) (3, 44.) (4, 90.) (5, 46.) (6, 94.) (7, 48.) (8, 98.) 1298e160fb9SStefano Zampinirow 5: (0, 51.) (1, 52.) (2, 53.) (3, 54.) (4, 55.) (5, 56.) (6, 57.) (7, 58.) (8, 59.) 1308e160fb9SStefano Zampinirow 6: (0, 122.) (1, 62.) (2, 126.) (3, 64.) (4, 130.) (5, 66.) (6, 134.) (7, 68.) (8, 138.) 1312f605a99SJose E. RomanMatZeroRows: 132a5225ed3SStefano Zampini type: seqaij 1338e160fb9SStefano Zampinirow 0: (0, 2.) (1, 2.) (2, 6.) (3, 4.) (4, 10.) (5, 6.) (6, 14.) (7, 8.) (8, 18.) 1348e160fb9SStefano Zampinirow 1: (0, 11.) (1, 12.) (2, 13.) (3, 14.) (4, 15.) (5, 16.) (6, 17.) (7, 18.) (8, 19.) 1358e160fb9SStefano Zampinirow 2: (0, 42.) (1, 22.) (2, 46.) (3, 24.) (4, 50.) (5, 26.) (6, 54.) (7, 28.) (8, 58.) 1368e160fb9SStefano Zampinirow 3: (0, 31.) (1, 32.) (2, 33.) (3, 34.) (4, 35.) (5, 36.) (6, 37.) (7, 38.) (8, 39.) 1378e160fb9SStefano Zampinirow 4: (0, 82.) (1, 42.) (2, 86.) (3, 44.) (4, 90.) (5, 46.) (6, 94.) (7, 48.) (8, 98.) 1388e160fb9SStefano Zampinirow 5: (0, 51.) (1, 52.) (2, 53.) (3, 54.) (4, 55.) (5, 56.) (6, 57.) (7, 58.) (8, 59.) 139a5225ed3SStefano Zampinirow 6: 1402f605a99SJose E. RomanMatShift: B = B - 2*I 1412f605a99SJose E. Roman type: seqaij 1422f605a99SJose E. Romanrow 0: (1, 2.) (2, 6.) (3, 4.) (4, 10.) (5, 6.) (6, 14.) (7, 8.) (8, 18.) 1432f605a99SJose E. Romanrow 1: (0, 11.) (1, 10.) (2, 13.) (3, 14.) (4, 15.) (5, 16.) (6, 17.) (7, 18.) (8, 19.) 1442f605a99SJose E. Romanrow 2: (0, 42.) (1, 22.) (2, 44.) (3, 24.) (4, 50.) (5, 26.) (6, 54.) (7, 28.) (8, 58.) 1452f605a99SJose E. Romanrow 3: (0, 31.) (1, 32.) (2, 33.) (3, 32.) (4, 35.) (5, 36.) (6, 37.) (7, 38.) (8, 39.) 1462f605a99SJose E. Romanrow 4: (0, 82.) (1, 42.) (2, 86.) (3, 44.) (4, 88.) (5, 46.) (6, 94.) (7, 48.) (8, 98.) 1472f605a99SJose E. Romanrow 5: (0, 51.) (1, 52.) (2, 53.) (3, 54.) (4, 55.) (5, 54.) (6, 57.) (7, 58.) (8, 59.) 1482f605a99SJose E. Romanrow 6: (6, -2.) 149