def METHOD_NAME(self, bc, index=None):
    phi = self.space.basis(bc)
    cell2dof = self.space.dof.cell2dof
    if index is None:
        bcp = np.einsum('...ij, ijk->...ik', phi, self.node[cell2dof, :])
    else:
        bcp = np.einsum('...ij, ijk->...ik', phi, self.node[cell2dof[index], :])
    bcp, _ = self.project(bcp)
    return bcp