Method Update

comp/facetsurffespace.cpp:350–414

  void FacetSurfaceFESpace :: Update()
  {
    FESpace::Update();
    
    if(print) 
      *testout << " FacetSurfaceFEspace with order " << order << " rel_order " << rel_order << " var_order " << var_order << endl; 

    nel = ma->GetNSE();
    nfa = ma->GetNEdges(); 
    first_edge_dof.SetSize(nfa+1);
    first_edge_dof = 0;
    
    // first_edge_dof.SetSize(nfa+1); 
    // first_edge_dof = nfa;

    if(ma->GetDimension() == 3)
      {
        /*
        for(int i = 0; i < nfa; i++)
          {
            first_edge_dof[i] = ndof;
            ndof += order+1;
          }
        first_edge_dof[nfa] = ndof;              
        */
        for (auto el : ma->Elements(BND))
          for (auto edge : el.Edges())
            first_edge_dof[edge] = order+1;
      }
    else if(ma->GetDimension() == 2)
      {
        for (auto el : ma->Elements(BND))
          for (auto vertex : el.Vertices())
            first_edge_dof[vertex] = 1;
      }
    else
      {
	throw Exception("Only implemented for 3d and 2d!");
      }

    size_t ndof = 0;
    for (size_t i = 0; i < nfa; i++)
      {
        size_t tmp = first_edge_dof[i];
        first_edge_dof[i] = ndof;
        ndof += tmp;
      }
    first_edge_dof[nfa] = ndof;    

    
    SetNDof (ndof);
    /*
    while (ma->GetNLevels() > ndlevel.Size())
      ndlevel.Append (ndof);
    ndlevel.Last() = ndof;
    */
    UpdateCouplingDofArray();