Method Evaluate

comp/gridfunction.cpp:1653–1747

  }

  void GridFunctionCoefficientFunction ::   
  Evaluate (const SIMD_BaseMappedIntegrationRule & ir,
            BareSliceMatrix<SIMD<double>> bvalues) const
  {
    // if (gf -> GetLevelUpdated() < gf->GetFESpace()->GetMeshAccess()->GetNLevels())
    if (gf -> GetLevelUpdated() < gf->GetMeshAccess()->GetNLevels())
      {
        bvalues.AddSize(Dimension(), ir.Size()) = 0.0;
        return;
      }
    
    ProxyUserData * ud = (ProxyUserData*)ir.GetTransformation().userdata;
    if (ud)
      {
        if (ud->HasMemory(this) && ud->Computed(this))
          {
            bvalues.AddSize(Dimension(), ir.Size()) = ud->GetAMemory(this);
            return;
          }
      }
    

    
    LocalHeapMem<100000> lh2("GridFunctionCoefficientFunction - Evaluate 3c");
    // static Timer timer ("GFCoeffFunc::Eval-vec", NoTracing);
    // RegionTimer reg (timer);
    auto values = bvalues.AddSize(Dimension(), ir.Size());
    const ElementTransformation & trafo = ir.GetTransformation();
    
    int elnr = trafo.GetElementNr();
    VorB vb = trafo.VB();
    ElementId ei(vb, elnr);

    if (!trafo.BelongsToMesh ((void*)(fes->GetMeshAccess().get())))
      {
        throw ExceptionNOSIMD ("SIMD - evaluation not available for different meshes");
        // for (int i = 0; i < ir.Size(); i++)
        // Evaluate (ir[i], values.Row(i));
        return;
      }
    
    if (!fes->DefinedOn(vb,trafo.GetElementIndex())) 
      { 
        values = 0.0; 
        return;
      }
    
    const FiniteElement & fel = fes->GetFE (ei, lh2);
    int dim = fes->GetDimension();

    ArrayMem<int, 50> dnums;
    fes->GetDofNrs (ei, dnums);
    
    VectorMem<50> elu(dnums.Size()*dim);

    gf->GetElementVector (comp, dnums, elu);
    fes->TransformVec (ei, elu, TRANSFORM_SOL);
    /*