/home/kcamus/qaas_runs/169-443-9681/intel/AMG/build/AMG/AMG/parcsr_mv/par_csr_matop.c: 865 - 989
--------------------------------------------------------------------------------

865:     for (i1 = ns; i1 < ne; i1++)
[...]
874:       if ( allsquare ) 
875:       {
876:          B_marker[i1] = jj_count_diag;
877:          C_diag_data[jj_count_diag] = zero;
878:          C_diag_j[jj_count_diag] = i1;
879:          jj_count_diag++;
[...]
886:       if (num_cols_offd_A)
887:       {
888:          for (jj2 = A_offd_i[i1]; jj2 < A_offd_i[i1+1]; jj2++)
889:          {
890:             i2 = A_offd_j[jj2];
891:             a_entry = A_offd_data[jj2];
[...]
897:             for (jj3 = B_ext_offd_i[i2]; jj3 < B_ext_offd_i[i2+1]; jj3++)
898:             {
899:                i3 = num_cols_diag_B+B_ext_offd_j[jj3];
[...]
907:                if (B_marker[i3] < jj_row_begin_offd)
908:                {
909:                   B_marker[i3] = jj_count_offd;
910:                   C_offd_data[jj_count_offd] = a_entry*B_ext_offd_data[jj3];
911:                   C_offd_j[jj_count_offd] = i3-num_cols_diag_B;
912:                   jj_count_offd++;
913:                }
914:                else
915:                   C_offd_data[B_marker[i3]] += a_entry*B_ext_offd_data[jj3];
916:             }
917:             for (jj3 = B_ext_diag_i[i2]; jj3 < B_ext_diag_i[i2+1]; jj3++)
918:             {
919:                i3 = B_ext_diag_j[jj3];
920:                if (B_marker[i3] < jj_row_begin_diag)
921:                {
922:                   B_marker[i3] = jj_count_diag;
923:                   C_diag_data[jj_count_diag] = a_entry*B_ext_diag_data[jj3];
924:                   C_diag_j[jj_count_diag] = i3;
925:                   jj_count_diag++;
926:                }
927:                else
928:                   C_diag_data[B_marker[i3]] += a_entry*B_ext_diag_data[jj3];
[...]
937:       for (jj2 = A_diag_i[i1]; jj2 < A_diag_i[i1+1]; jj2++)
938:       {
939:          i2 = A_diag_j[jj2];
940:          a_entry = A_diag_data[jj2];
[...]
946:          for (jj3 = B_diag_i[i2]; jj3 < B_diag_i[i2+1]; jj3++)
947:          {
948:             i3 = B_diag_j[jj3];
[...]
956:             if (B_marker[i3] < jj_row_begin_diag)
957:             {
958:                B_marker[i3] = jj_count_diag;
959:                C_diag_data[jj_count_diag] = a_entry*B_diag_data[jj3];
960:                C_diag_j[jj_count_diag] = i3;
961:                jj_count_diag++;
962:             }
963:             else
964:             {
965:                C_diag_data[B_marker[i3]] += a_entry*B_diag_data[jj3];
966:             }
967:          }
968:          if (num_cols_offd_B)
969:          {
970:             for (jj3 = B_offd_i[i2]; jj3 < B_offd_i[i2+1]; jj3++)
971:             {
972:                i3 = num_cols_diag_B+map_B_to_C[B_offd_j[jj3]];
[...]
980:                if (B_marker[i3] < jj_row_begin_offd)
981:                {
982:                   B_marker[i3] = jj_count_offd;
983:                   C_offd_data[jj_count_offd] = a_entry*B_offd_data[jj3];
984:                   C_offd_j[jj_count_offd] = i3-num_cols_diag_B;
985:                   jj_count_offd++;
986:                }
987:                else
988:                {
989:                   C_offd_data[B_marker[i3]] += a_entry*B_offd_data[jj3];
