Method Update

comp/facetfespace.cpp:765–994

  }
  
  void FacetFESpace :: Update()
  {
    FESpace::Update();
    
    if(print) 
      *testout << " FacetFEspace with order " << order << " rel_order " << rel_order << " var_order " << var_order << endl; 

    bool first_update = GetTimeStamp() < ma->GetTimeStamp();
    if (first_update) timestamp = NGS_Object::GetNextTimeStamp();
    
    if (low_order_space)
      low_order_space -> Update();

    nel = ma->GetNE();
    nfa = ma->GetNFacets(); 

    if (first_update)
      {
	int p = 0; 
	if(!var_order) p = order; 
    
	order_facet.SetSize(nfa);
	order_facet = p;

	fine_facet.SetSize(nfa);
	fine_facet = false; 
    
	Array<int> fanums;

	for (Ngs_Element el : ma->Elements<BND>())
	  if (DefinedOn(el))
	    fine_facet[el.Facets()] = true;
    
	for (int i = 0; i < nel; i++)
	  {
	    ElementId ei(VOL, i);
	    if (!DefinedOn(ei)) continue;
	    ELEMENT_TYPE eltype=ma->GetElType(ei); 
	    const POINT3D * points = ElementTopology :: GetVertices (eltype);	


	    if (ma->GetDimension() == 1)
	      {
		auto fanums = ma->GetElVertices(ei);
		for (auto f : fanums)
		  fine_facet[f] = true;
              }
            
	    else if (ma->GetDimension() == 2)
	      {
		auto fanums = ma->GetElEdges(ei);
		for (auto f : fanums)
		  fine_facet[f] = true;
            
		if(var_order)
		  {
		    IVec<3> el_orders = ma->GetElOrders(i);