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

1747:              if (n_fine) 
[...]
1811:                 for (j=A_diag_i[i1]+1; j < A_diag_i[i1+1]; j++)
1812:                 {
1813:                j1 = A_diag_j[j];
1814:                if (tmp_marker[j1] == i1)
1815:                {
1816:                   for (k=P_diag_i[j1]; k < P_diag_i[j1+1]; k++)
1817:                   {
1818:                  k1 = P_diag_j[k];
1819:                      alfa = A_diag_data[j]*P_diag_data[k];
1820:                      P_diag_data[tmp_array[k1]] += alfa;
1821:                      sum_C += alfa;
1822:                      sum_N += alfa;
1823:                       }
1824:                   for (k=P_offd_i[j1]; k < P_offd_i[j1+1]; k++)
1825:                   {
1826:                  k1 = P_offd_j[k];
1827:                      alfa = A_diag_data[j]*P_offd_data[k];
1828:                      P_offd_data[tmp_array_offd[k1]] += alfa;
1829:                      sum_C += alfa;
1830:                      sum_N += alfa;
1831:                       }
1832:                    }
1833:                    else
1834:                    {
1835:                       if (CF_marker[j1] != -3 && 
1836:                 (num_functions == 1 || dof_func[i1] == dof_func[j1]))
1837:                  sum_N += A_diag_data[j];
