On the generation of symmetric Lebesgue-like points in the triangle

We compute point sets on the triangle that have low Lebesgue constant, with sixfold symmetries and Gauss–Legendre–Lobatto distribution on the sides, up to interpolation degree 18. Such points have the best Lebesgue constants among the families of symmetric points used so far in the framework of triangular spectral elements.