/beegfs/hackathon/users/eoseret/qaas_runs_test/178-668-9699/intel/AMG/build/AMG/AMG/parcsr_ls/ams.c: 3662 - 3684
--------------------------------------------------------------------------------

3662: #pragma omp parallel for private(i,ii,jj,res) HYPRE_SMP_SCHEDULE
3663: #endif
3664:       for (i = 0; i < n; i++)
3665:       {
3666:          /*-----------------------------------------------------------
3667:           * If diagonal is nonzero, relax point i; otherwise, skip it.
3668:           *-----------------------------------------------------------*/
3669:          if (A_diag_data[A_diag_i[i]] != zero)
3670:          {
3671:             res = f_data[i];
3672:             for (jj = A_diag_i[i]; jj < A_diag_i[i+1]; jj++)
3673:             {
3674:                ii = A_diag_j[jj];
3675:                res -= A_diag_data[jj] * Vtemp_data[ii];
3676:             }
3677:             for (jj = A_offd_i[i]; jj < A_offd_i[i+1]; jj++)
3678:             {
3679:                ii = A_offd_j[jj];
3680:                res -= A_offd_data[jj] * Vext_data[ii];
3681:             }
3682:             u_data[i] += (relax_weight*res)/l1_norms[i];
3683:          }
3684:       }
