Method Apply

bem/diffopwithfactor.hpp:102–117

    void Apply (const FiniteElement & fel,
                const BaseMappedIntegrationPoint & mip,
                BareSliceVector<Complex> x, 
                FlatVector<Complex> flux,
                LocalHeap & lh) const override
    {
      auto dims = factor->Dimensions();

      HeapReset hr(lh);
      FlatVector<Complex> tmpflux(dims[1], lh);
      FlatVector<Complex> factorx(dims[0]*dims[1], lh);
      
      diffop -> Apply (fel, mip, x, tmpflux, lh);
      factor -> Evaluate (mip, factorx);
      flux = factorx.AsMatrix(dims[0], dims[1]) * tmpflux;
    }

    
    void Apply (const FiniteElement & fel,