Method Mult

linalg/commutingAMG.cpp:285–315

  
  void AMG_H1 :: Mult (const BaseVector & x, BaseVector & y) const
  {
    static Timer timer("H1-AMG::Mult");
    RegionTimer reg (timer);

    if (inv)
      {
	y = (*inv) * x;
	return;
      }

    auto hv = pmat->CreateColVector();
    auto wc = coarsemat->CreateColVector();
    auto dc = coarsemat->CreateColVector();

    y = 0;
    jacobi->GSSmooth (y, x);

    if (recAMG)
      {
	hv = x - (*pmat) * y;
	dc = Transpose (*prol) * hv;

	if (recAMG)
	  recAMG -> Mult (dc, wc);

	y += (*prol) * wc;
      }

    jacobi->GSSmoothBack (y, x);
  }