rbmatlab 0.10.01
|
00001 function res = fem_ndofs(params,grid) 00002 %function res = fem_ndofs_per_element(params,grid) 00003 % 00004 % function computing number of dofs based on 00005 % params.dimrange, params.grid and params.pdeg. 00006 % 00007 % Parameters: 00008 % grid: object of type triagrid 00009 00010 % Bernard Haasdonk 2.9.2009 00011 00012 00013 % num dofs on nodes: 00014 ndof_nodes = grid.nvertices; 00015 00016 % num dofs on edges: 00017 ndof_edges = (grid.nedges_interior + grid.nedges_boundary) * (params.pdeg-1); 00018 00019 % num dofs in interior: 00020 ndof_elements = grid.nelements * (params.pdeg-1)* (params.pdeg-2)*0.5; 00021 00022 res = grid.nelements*params.dimrange*(params.pdeg+1)*(params.pdeg+2)/2; 00023 00024 res = (ndof_nodes + ndof_edges + ndof_elements) * params.dimrange;