Method Apply

fem/diffop.cpp:107–136

  
  void DifferentialOperator ::
  Apply (const FiniteElement & fel,
         const BaseMappedIntegrationPoint & mip,
         BareSliceVector<double> x, 
         FlatVector<double> flux,
         LocalHeap & lh) const
  {
#ifndef FASTCOMPILE    
    static int cnt = 0;
    if (cnt < 3)
      {
        cnt++;
        cout << "called base class apply, type = " << typeid(*this).name() << endl;
      }
#endif
    HeapReset hr(lh);
    FlatMatrix<double,ColMajor> mat(VSDim(), fel.GetNDof(), lh);
    if (auto vsemb = GetVSEmbedding(); vsemb)
      {
        CalcMatrixVS (fel, mip, mat, lh);
        FlatVector<double> flux1(VSDim(), lh);
        flux1 = mat * x;
        flux = *vsemb * flux1;
      }
    else
      {
        CalcMatrix (fel, mip, mat, lh);        
        flux = mat * x;
      }
  }
  
  void DifferentialOperator ::
  Apply (const FiniteElement & fel,