Function PrecomputeCacheCF

fem/coefficient.cpp:7821–7852

void PrecomputeCacheCF (const Array<CoefficientFunction*> & cachecfs, BaseMappedIntegrationRule & mir,
                        LocalHeap & lh)
{
  // static Timer t("Precompute CacheCF");
  // RegionTracer reg(TaskManager::GetThreadId(), t);
  
  auto & trafo = mir.GetTransformation();
  ProxyUserData * ud = static_cast<ProxyUserData*> (trafo.userdata);
  if (!ud) throw Exception ("don't have a userdata");
  
  new (&ud->caches) FlatArray<pair<const CoefficientFunction*, void*>> (cachecfs.Size(), lh);

  for (int i = 0; i < cachecfs.Size(); i++)
    {
      ud->caches[i].first = cachecfs[i];
      // or complex ...
      if (cachecfs[i]->IsComplex())
        {
          auto * mat = new (lh) FlatMatrix<Complex>(mir.Size(), cachecfs[i]->Dimension(), lh);
          static_cast<CacheCoefficientFunction*>(cachecfs[i])->Func() -> Evaluate (mir, *mat);
          ud->caches[i].second = mat;
        }
      else
        {
          auto * mat = new (lh) FlatMatrix<double>(mir.Size(), cachecfs[i]->Dimension(), lh);
          static_cast<CacheCoefficientFunction*>(cachecfs[i])->Func() -> Evaluate (mir, *mat);
          ud->caches[i].second = mat;
        }
        
    }

}