Function MultAdd

comp/l2hofespace.cpp:2648–2678

    }

    void MultAdd (double val, const BaseVector & v, BaseVector & prod) const override
    {
      auto fx = v.FV<double>();
      auto fy = prod.FV<double>();
      ParallelForRange
        (dofs->Size(), [&] (IntRange r)
         {
           constexpr size_t BS = 128;
           Matrix<> hx(BS, elbmat.Width());
           Matrix<> tmp(BS, elbmat.Height());

           for (size_t bi = r.First(); bi < r.Next(); bi+= BS)
             {
               size_t li = min2(bi+BS, r.Next());
               size_t num = li-bi;

               for (size_t i = 0; i < num; i++)
                 hx.Row(i) = fx( (*dofs)[bi+i]);

               tmp.Rows(0, num) = hx.Rows(0, num) * Trans(elbmat);

               size_t base = r.First()*tmp.Width();
               for (size_t i : Range(num*tmp.Width()))
                 tmp(i) *= rho_jac(base+i);

               hx.Rows(0, num) = tmp.Rows(0,num) * elbmat;
               for (size_t i = 0; i < num; i++)
                 fy( (*dofs)[bi+i]) += val * hx.Row(i);
             }
         });
    }


  };