Method Update

comp/l2hofespace.cpp:1018–1060

  void L2HighOrderFESpace :: Update()
  {
    FESpace::Update();
    if(low_order_space) low_order_space -> Update();

    // nel = ma->GetNE();

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

    if (first_update)
      {
	order_inner.SetSize(ma->GetNE());

	order_inner = IVec<3>(order);

	if(var_order)
	  for (auto i : Range(ma->GetNE()))
            order_inner[i] = ma->GetElOrders(i)+IVec<3>(rel_order);
        
        for (auto i : Range(ma->GetNE()))
	  {
	    ElementId ei(VOL,i);
	    order_inner[i] = order_inner[i] + IVec<3> (et_bonus_order[ma->GetElType(ei)]);
	    order_inner[i] = Max(order_inner[i], IVec<3>(0));
	    if (!DefinedOn (ei))
	      order_inner[i] = all_dofs_together ? -1 : 0;
	  }
	if(print)
	  *testout << " order_inner (l2ho) " << order_inner << endl;
      }

    UpdateDofTables();
    /*
    while (ma->GetNLevels() > ndlevel.Size())
      ndlevel.Append (ndof);
    ndlevel.Last() = ndof;
    */
    // if(low_order_space) prol->Update(*this);
    if (prol) prol->Update(*this);

    UpdateCouplingDofArray();
  }

  void L2HighOrderFESpace :: UpdateCouplingDofArray()
  {