/home/kcamus/qaas_runs/169-443-9681/intel/AMG/build/AMG/AMG/parcsr_ls/par_strength.c: 1714 - 1797
--------------------------------------------------------------------------------

1714:          for (ic = ic_begin; ic < ic_end; ic++)
[...]
1720:              HYPRE_Int i1 = coarse_to_fine[ic];
1721:        
1722:              HYPRE_Int jj_row_begin_diag = num_nonzeros_diag;
1723:              HYPRE_Int jj_row_begin_offd = num_nonzeros_offd;
1724: 
1725:              C_diag_i[ic] = num_nonzeros_diag;
1726:              if (num_cols_offd_C)
1727:              {
1728:                 C_offd_i[ic] = num_nonzeros_offd;
1729:              }
1730: 
1731:              for (jj1 = S_diag_i[i1]; jj1 < S_diag_i[i1+1]; jj1++)
1732:              {
1733:                  i2 = S_diag_j[jj1];
1734:                  if (CF_marker[i2] > 0)
1735:                  {
1736:                     index = fine_to_coarse[i2];
1737:                     if (S_marker[index] < jj_row_begin_diag)
1738:                     {
1739:                        S_marker[index] = num_nonzeros_diag;
1740:                        num_nonzeros_diag++;
1741:                     }
1742:                  }
1743:                  for (jj2 = S_diag_i[i2]; jj2 < S_diag_i[i2+1]; jj2++)
1744:                  {
1745:                     i3 = S_diag_j[jj2];
1746:                     if (CF_marker[i3] > 0)
1747:                     {
1748:                        index = fine_to_coarse[i3];
1749:                        if (index != ic && S_marker[index] < jj_row_begin_diag)
1750:                        {
1751:                           S_marker[index] = num_nonzeros_diag;
1752:                           num_nonzeros_diag++;
1753:                        }
1754:                     }
1755:                  }
1756:                  for (jj2 = S_offd_i[i2]; jj2 < S_offd_i[i2+1]; jj2++)
1757:                  {
1758:                     i3 = S_offd_j[jj2];
1759:                     if (CF_marker_offd[i3] > 0)
1760:                     {
1761:                        index = map_S_to_C[i3];
1762:                        if (S_marker_offd[index] < jj_row_begin_offd)
1763:                        {
1764:                           S_marker_offd[index] = num_nonzeros_offd;
1765:                           num_nonzeros_offd++;
1766:                        }
1767:                     }
1768:                  }
1769:              }
1770:              for (jj1 = S_offd_i[i1]; jj1 < S_offd_i[i1+1]; jj1++)
1771:              {
1772:                  i2 = S_offd_j[jj1];
1773:                  if (CF_marker_offd[i2] > 0)
1774:                  {
1775:                     index = map_S_to_C[i2];
1776:                     if (S_marker_offd[index] < jj_row_begin_offd)
1777:                     {
1778:                        S_marker_offd[index] = num_nonzeros_offd;
1779:                        num_nonzeros_offd++;
1780:                     }
1781:                  }
1782:                  for (jj2 = S_ext_diag_i[i2]; jj2 < S_ext_diag_i[i2+1]; jj2++)
1783:                  {
1784:                     i3 = S_ext_diag_j[jj2];
1785:                     if (i3 != ic && S_marker[i3] < jj_row_begin_diag)
1786:                     {
1787:                        S_marker[i3] = num_nonzeros_diag;
1788:                        num_nonzeros_diag++;
1789:                     }
1790:                  }
1791:                  for (jj2 = S_ext_offd_i[i2]; jj2 < S_ext_offd_i[i2+1]; jj2++)
1792:                  {
1793:                     i3 = S_ext_offd_j[jj2];
1794:                     if (S_marker_offd[i3] < jj_row_begin_offd)
1795:                     {
1796:                        S_marker_offd[i3] = num_nonzeros_offd;
1797:                        num_nonzeros_offd++;
