Method Update

comp/fespace.cpp:438–515

  }
  
  void FESpace :: Update()
  {
    if (print)
      {
 	*testout << "Update FESpace, type = " << typeid(*this).name() << endl;
	*testout << "name = " << name << endl;
      }

    /*
    for (int i = 0; i < specialelements.Size(); i++)
      delete specialelements[i]; 
    specialelements.SetSize(0);
    */
    
    ma->UpdateBuffers();  // is free if netgen-mesh did not change
    int dim = ma->GetDimension();
    
    dirichlet_vertex.SetSize (ma->GetNV());
    dirichlet_edge.SetSize (ma->GetNEdges());
    if (dim == 3)
      dirichlet_face.SetSize (ma->GetNFaces());
    
    dirichlet_vertex = false;
    dirichlet_edge = false;
    dirichlet_face = false;

    /*
    {
      if (dirichlet_boundaries.Size())
        for (Ngs_Element ngel : ma->Elements(BND))  // .OmpSplit())
          if (dirichlet_boundaries[ngel.GetIndex()])
            {
              dirichlet_vertex[ngel.Vertices()] = true;
              if (dim >= 2)
                dirichlet_edge[ngel.Edges()] = true;
              if (dim == 3)
                dirichlet_face[ngel.Faces()[0]] = true;
            }
    }
    */
    for (auto vb : { BND, BBND, BBBND })
      {
        auto & dc = dirichlet_constraints[int(vb)];
        if (dc.Size())
          for (Ngs_Element ngel : ma->Elements(vb)) 
            if (dc[ngel.GetIndex()])
              {
                dirichlet_vertex[ngel.Vertices()] = true;
                // if (dim >= 2)
                if (dim-vb >= 1)
                  dirichlet_edge[ngel.Edges()] = true;
                // if (dim == 3)
                if (dim-vb >= 2)
                  dirichlet_face[ngel.Faces()[0]] = true;
              }
      }