Method MultTrans

bem/fmmoperator.hpp:174–196

    }

    void MultTrans(const BaseVector & x, BaseVector & y) const override
    {
        static Timer tall("ngbem fmm apply Trans "+KERNEL::Name()); RegionTimer reg(tall);

        auto shape = KERNEL::Shape();
        auto matx = x.FV<typename KERNEL::value_type>().AsMatrix(xpts.Size(), shape[0]);
        auto maty = y.FV<typename KERNEL::value_type>().AsMatrix(ypts.Size(), shape[1]);

        maty = 0;
        auto singmp = kernel.target.CreateMultipoleExpansion (cy, ry, fmm_params);
        ParallelFor (ypts.Size(), [&](int i){
          kernel.target.AddSource(*singmp, ypts[i], ynv[i], matx.Row(i));
        });
        singmp->CalcMP();
        auto regmp = kernel.source.CreateLocalExpansion (cx, rx, fmm_params);
        ParallelFor (xpts.Size(), [&](int i){
          regmp->AddTarget(xpts[i]);
        });
        regmp->CalcMP(singmp);
        ParallelFor (xpts.Size(), [&](int i) {
          kernel.source.EvaluateMP(*regmp, xpts[i], xnv[i], maty.Row(i));
        });
    }

    BaseMatrix::OperatorInfo GetOperatorInfo () const override
    {