/home/kcamus/qaas_runs/169-443-9681/intel/AMG/build/AMG/AMG/parcsr_ls/par_multi_interp.c: 917 - 1133
--------------------------------------------------------------------------------

917:           for (i=0; i < n_coarse; i++)
[...]
1072:           for (i=thread_start; i < thread_stop; i++)
1073:           {
1074:              i1 = pass_array[i];
1075:              for (j=S_diag_i[i1]; j < S_diag_i[i1+1]; j++)
1076:              {
1077:                 j1 = S_diag_j[j];
1078:                 if (assigned[j1] == pass-1)
1079:                 {
1080:                    j_start = P_diag_start[j1];
1081:                    j_end = j_start+P_diag_i[j1+1];
1082:                    for (k=j_start; k < j_end; k++)
1083:                    {
1084:                       k1 = P_diag_pass[pass-1][k];
1085:                       if (P_marker[k1] != -i1-1)
1086:                       {
1087:                           P_diag_pass[pass][cnt_nz++] = k1;
1088:                           P_marker[k1] = -i1-1;
1089:                       }
1090:                    }
1091:                    j_start = P_offd_start[j1];
1092:                    j_end = j_start+P_offd_i[j1+1];
1093:                    for (k=j_start; k < j_end; k++)
1094:                    {
1095:                       k1 = P_offd_pass[pass-1][k];
1096:                       if (P_marker_offd[k1] != -i1-1)
1097:                       {
1098:                          P_offd_pass[pass][cnt_nz_offd++] = k1;
1099:                          P_marker_offd[k1] = -i1-1;
1100:                       }
1101:                    }
1102:                 }
1103:              }
1104:              for (j=S_offd_i[i1]; j < S_offd_i[i1+1]; j++)
1105:              {
1106:                 j1 = S_offd_j[j];
1107:                 if (assigned_offd[j1] == pass-1)
1108:                 {
1109:                    j_start = Pext_start[j1];
1110:                    j_end = j_start+Pext_i[j1+1];
1111:                    for (k=j_start; k < j_end; k++)
1112:                    {
1113:                       k1 = Pext_pass[pass][k];
1114:                       if (k1 < 0)
1115:                       {
1116:                          if (P_marker[-k1-1] != -i1-1)
1117:                          {
1118:                             P_diag_pass[pass][cnt_nz++] = -k1-1;
1119:                             P_marker[-k1-1] = -i1-1;
1120:                          }
1121:                       }
1122:                       else if (P_marker_offd[k1] != -i1-1)
1123:                       {
1124:                          P_offd_pass[pass][cnt_nz_offd++] = k1;
1125:                          P_marker_offd[k1] = -i1-1;
[...]
1133:           if ( (n_coarse_offd) || (new_num_cols_offd  == local_index+1) )
