std::pair<pt, pt> leaf(pt p, std::ostream *os = nullptr) const { pt p1 = geodesique(p, os); p.reverse(); pt p2 = geodesique(p, os); return (p1 < p2) ? ptpair{p1, p2} : ptpair{p2, p1}; }