def METHOD_NAME(self, point, index=np.s_[:], p=None):
    smphi = 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)
    phi = np.zeros(shape, dtype=self.ftype)
    phi[..., :ldof//2, 0] = smphi
    phi[..., -ldof//2:, 1] = smphi
    return phi