Method Mult

comp/hcurlamg.cpp:954–992

  template<typename SCAL>
  void HCurlAMG_Matrix<SCAL> :: Mult(const BaseVector& f, BaseVector& u) const
  {
    static Timer timer("HCurlAMG::Mult");
    RegionTimer regt(timer);
    static Timer timer_n("Node correction");
    static Timer timer_c("Coarse correction");

    u = 0.;
    smoother->Smooth(u, f, param.smoothing_steps);
    
    auto residuum = f.CreateVector();

    if (gradient)
      {
        residuum = f - (*mat) * u;
        auto op = gradient*node_h1*trans_gradient;
        auto size = op->Height();
        u.Range(size) += *op * residuum.Range(size);
      }
    
    {
      RegionTimer rt(timer_c);
      residuum = f - (*mat) * u;

      auto op = prolongation * coarse_precond * restriction;
      auto size = op->Height();
      u.Range(size) += *op * residuum.Range(size);
    }

    if (gradient)
      {
        residuum = f - (*mat) * u;
        auto op = gradient*node_h1*trans_gradient;
        auto size = op->Height();
        u.Range(size) += *op * residuum.Range(size);
      }
    
    smoother->SmoothBack(u, f, param.smoothing_steps);    
  }


  /* ************************* HCurlAMG preconditioner ************************** */