Function CreateGridFunction

comp/gridfunction.cpp:1186–1204

  */

  shared_ptr<GridFunction> CreateGridFunction (shared_ptr<FESpace> space,
                                               const string & name, const Flags & flags)
  {
    /*
    shared_ptr<GridFunction> gf =
      CreateSharedVecObject<T_GridFunction, GridFunction> 
      (space->GetDimension() * int(flags.GetNumFlag("cacheblocksize",1)), 
       space->IsComplex(), space, name, flags);
    */
    shared_ptr<GridFunction> gf;
    if (space->IsComplex())
      gf = make_shared<S_GridFunction<Complex>> (space, name, flags);
    else
      gf = make_shared<S_GridFunction<double>> (space, name, flags);
    
    gf->SetCacheBlockSize(int(flags.GetNumFlag("cacheblocksize",1)));
    
    return gf;
  }