def get_mapping(self, qp_coors, weights, poly_space=None, mode='surface'): """ Get the mapping for given quadrature points, weights, and polynomial space. Returns ------- cmap : CMapping instance The surface mapping. """ poly_space = get_default(poly_space, self.poly_space) bf_g = self.get_base(qp_coors, diff=True) if nm.allclose(bf_g, 0.0): raise ValueError('zero base function gradient!') cmap = CMapping(self.n_el, qp_coors.shape[0], self.dim, poly_space.n_nod, mode=mode) cmap.describe(self.coors, self.conn, bf_g, None, weights) return cmap
def get_mapping(self, qp_coors, weights, poly_space=None, ori=None): """ Get the mapping for given quadrature points, weights, and polynomial space. Returns ------- cmap : CMapping instance The volume mapping. """ poly_space = get_default(poly_space, self.poly_space) bf_g = self.get_base(qp_coors, diff=True) ebf_g = poly_space.eval_base(qp_coors, diff=True, ori=ori, force_axis=True) flag = ori is not None cmap = CMapping(self.n_el, qp_coors.shape[0], self.dim, poly_space.n_nod, mode='volume', flag=flag) cmap.describe(self.coors, self.conn, bf_g, ebf_g, weights) return cmap