/home/kcamus/qaas_runs/169-443-9681/intel/AMG/build/AMG/AMG/parcsr_ls/par_lr_interp.c: 1494 - 1655
--------------------------------------------------------------------------------

1494:            for (jj = S_diag_i[i]; jj < S_diag_i[i+1]; jj++)
[...]
1503:              if (CF_marker[i1] >= 0)
1504:              {
1505:                if (P_marker[i1] < jj_begin_row)
1506:                {
1507:                  P_marker[i1] = jj_counter;
1508:                  P_diag_j[jj_counter]    = fine_to_coarse[i1];
1509:                  P_diag_data[jj_counter] = zero;
1510:                  jj_counter++;
1511:                }
1512:              }
1513:              else  if (CF_marker[i1] != -3)
1514:              {
1515:                P_marker[i1] = strong_f_marker;
1516:                for (kk = S_diag_i[i1]; kk < S_diag_i[i1+1]; kk++)
1517:                {
1518:                  k1 = S_diag_j[kk];
1519:                  if (CF_marker[k1] >= 0)
1520:                  {
1521:                    if(P_marker[k1] < jj_begin_row)
1522:                    {
1523:                      P_marker[k1] = jj_counter;
1524:                      P_diag_j[jj_counter] = fine_to_coarse[k1];
1525:                      P_diag_data[jj_counter] = zero;
1526:                      jj_counter++;
1527:                    }
1528:                  }
1529:                }
1530:                if(num_procs > 1)
1531:                {
1532:                  for (kk = S_offd_i[i1]; kk < S_offd_i[i1+1]; kk++)
1533:                  {
1534:                    if(col_offd_S_to_A)
1535:                     k1 = col_offd_S_to_A[S_offd_j[kk]];
1536:                    else
1537:                     k1 = S_offd_j[kk];
1538:                    if(CF_marker_offd[k1] >= 0)
1539:                    {
1540:                      if(P_marker_offd[k1] < jj_begin_row_offd)
1541:                      {
1542:                        P_marker_offd[k1] = jj_counter_offd;
1543:                        P_offd_j[jj_counter_offd] = k1;
1544:                        P_offd_data[jj_counter_offd] = zero;
1545:                        jj_counter_offd++;
[...]
1655:                    for(jj1 = A_offd_i[i1]; jj1 < A_offd_i[i1+1]; jj1++)
