| 70 | } |
| 71 | |
| 72 | void JKM_FESpace::Update() |
| 73 | { |
| 74 | FESpace::Update(); |
| 75 | first_edge_dof.SetSize(ma->GetNEdges()+1); |
| 76 | first_element_dof.SetSize(ma->GetNE()+1); |
| 77 | |
| 78 | size_t ndof = 0; |
| 79 | for (size_t i = 0; i < ma->GetNEdges(); i++) |
| 80 | { |
| 81 | first_edge_dof[i] = ndof; |
| 82 | ndof += 2; |
| 83 | } |
| 84 | first_edge_dof[ma->GetNEdges()] = ndof; |
| 85 | |
| 86 | for (size_t i = 0; i < ma->GetNE(); i++) |
| 87 | { |
| 88 | first_element_dof[i] = ndof; |
| 89 | ndof += 9; |
| 90 | } |
| 91 | first_element_dof[ma->GetNE()] = ndof; |
| 92 | |
| 93 | SetNDof (ndof); |
| 94 | } |
| 95 | |
| 96 | void JKM_FESpace::FinalizeUpdate() |
| 97 | { |