def __init__(self, mesh_el_group, order, index, *, unit_nodes):
basis = mp.orthonormal_basis_for_space(mp.QN(mesh_el_group.dim, order),
mp.Hypercube(mesh_el_group.dim))
super().__init__(mesh_el_group,
order,
index,
basis=basis,
unit_nodes=unit_nodes)
def space(self):
return mp.QN(self.dim, self.order)
def from_mesh_interp_matrix(self):
meg = self.mesh_el_group
meg_space = mp.QN(meg.dim, meg.order)
return mp.resampling_matrix(
mp.basis_for_space(meg_space, self._shape).functions,
self.unit_nodes, meg.unit_nodes)