Method Update

comp/hdivhofespace.cpp:430–621

  }
  
  void HDivHighOrderFESpace :: Update()
  {
    FESpace::Update();
    // In discontinuous spaces the order on edges and faces  
    // are also set to the maximal order of adjacent elements  
    // and not to the element order (Motivation: Equilibrated_EE) 

    // SZ hack since is not supported(tested) yet 
    rel_curl_order= rel_order; 
    curl_order = order; 

    bool first_update = GetTimeStamp() < ma->GetTimeStamp();
    if (first_update) timestamp = NGS_Object::GetNextTimeStamp();
    
    if (low_order_space)
      low_order_space -> Update();
    
    // int nv = ma->GetNV();
    size_t nel = ma->GetNE();
    size_t nfa = ma->GetNFacets();
    size_t dim = ma->GetDimension();

    if (first_update)
      {
	order_facet.SetSize(nfa);
	order_inner.SetSize(nel); 
	order_inner_curl.SetSize(nel); 
	fine_facet.SetSize(nfa); 
	boundary_facet.SetSize(nfa); 
   
	boundary_facet = false;
	/*
	  for (int i = 0; i < ma->GetNSE(); i++)
	  {
	  Array<int> elfacets; 
	  ma->GetSElFacets (i,elfacets); 
	  boundary_facet[elfacets[0]] = true;
	  }
	*/
	for (auto el : ma->Elements(BND))
	  boundary_facet[el.Facets()] = true;

	// cout << " order hdiv " << order << endl; 
	// cout << " curl_order hdiv " << curl_order << endl; 

	int p = 0, pc = 0; 
    
	if(!var_order)
	  {
	    p = order; 
	    pc = curl_order; 
	  } 
    
	order_facet = pc;
	order_inner = p;
	order_inner_curl = pc;
	fine_facet = 0; //!!!!