def METHOD_NAME(self, point, index=np.s_[:], p=None, scaled=True):
    """
    p >= 0
    """
    smgphi = self.smspace.METHOD_NAME(point, index=index, p=p)
    ldof = self.number_of_local_dofs(p=p, doftype='cell')
    shape = point.shape[:-1]+(ldof, 2, 2) 
    gphi = np.zeros(shape, dtype=self.ftype)
    gphi[..., :ldof//2, : , 0] = smgphi
    gphi[..., -ldof//2:, : , 1] = smgphi
    return gphi