Method CreateSmoothingBlocks

comp/fespace.cpp:1604–2072

  shared_ptr<Table<int>> FESpace :: CreateSmoothingBlocks (const Flags & flags) const
  {
    bool eliminate_internal =
      flags.GetDefineFlag("eliminate_internal") ||
      flags.GetDefineFlag("condense");
      
    auto freedofs = GetFreeDofs(eliminate_internal);

    if (flags.GetDefineFlag("subassembled"))
      {
        freedofs = make_shared<BitArray>(*freedofs);
        for (auto i : Range(GetNDof()))
          if (GetDofCouplingType(i) != WIREBASKET_DOF)
            freedofs->Clear(i);
      }

    if (flags.AnyFlagDefined("additional_dirichlet_constraints"))
      {
        Region reg = std::any_cast<Region>(flags.GetAnyFlag("additional_dirichlet_constraints"));
        BitArray dofs = GetDofs(reg);
        dofs.Invert();
        freedofs = make_shared<BitArray>(*freedofs);
        freedofs->And(dofs);
      }

    FilteredTableCreator creator(freedofs.get());

    /*
    for ( ; !creator.Done(); creator++)
      {
	for (size_t i = 0; i < nd; i++)
	  if (freedofs->Test(i))
	    creator.Add (i, i);
      }
    */


    // new version for ngs24:

    Array<string> blocktypes { flags.GetStringListFlag("blocktype") };
    std::map<string, optional<string>> filtered_blocktypes;
    
    if (flags.GetStringFlag("blocktype")=="vertex" ||
        flags.GetStringFlag("blocktype")=="edge" ||
        flags.GetStringFlag("blocktype")=="face" ||
        flags.GetStringFlag("blocktype")=="facet" ||
        flags.GetStringFlag("blocktype")=="element" ||        
        flags.GetStringFlag("blocktype")=="vertexpatch" ||
        flags.GetStringFlag("blocktype")=="edgepatch" ||                
        flags.GetStringFlag("blocktype")=="vertexedge")
      blocktypes += flags.GetStringFlag("blocktype");
    
    for (auto fulltype : blocktypes)
      {
        auto pos = fulltype.find(':');
        string type, filter;
        if (pos != std::string::npos)
          {