In 2009, Richard Schwartz proved that any obtuse triangle whose largest angle is $\leq100^{\circ}$ has a stable periodic billiard orbit. My question then, is:
How can I reproduce Schwartz's result using a simulation?
More specifically,
How can I develop a numerical method converging to a $(P_0,V_0)$ (initial position, initial direction) giving a Schwartz periodic orbit?
I've created a functioning simulation which produces paths in a triangular billiard table, given the initial position and direction of the ball:
However, I have no idea how to even begin computationally reproducing Schwartz's result. A brute-force approach of supplying the ball with every initial condition (position, direction) for every triangle with angles $100\leq \alpha \leq 180$ is simply infeasible for an infinite search space. Given this, I would very much appreciate any insights from the Math Stack Exchange community for how to go about reproducing this result. And of course, I can make any modifications to my program as necessary.
Best Answer
The explanation is given in chapter 2 of the cited paper. Every orbit corresponds to an unfolding, which in turn can be described by a "word" formed by digits $123$, a digit for each "bounce" in the orbit, where $1$ represents the shortest side, $2$ the next-shortest side and $3$ the longest side. The first digit in the word represents the starting side. Subsequent digits represent a reflection of the last constructed triangle about one of its sides, giving the next triangle in the unfolding.
If the word satisfies a condition given in Lemma 2.8, then the starting side of the first triangle is parallel to the corresponding side of the last triangle. In this case they show that any straight line joining a starting point $F$ on the starting side with the corresponding point $F'$ on the last side represents a periodic orbit, provided segment $FF'$ lies all within the unfolding.
You can see below an example of closed orbit, taken from the paper (Chapter 2, Figure 1) and corresponding to $W=12323131232313$. Any straight line joining a starting point $F$ on side $A_1B_1$ with the corresponding point $F'$ on side $A_8B_8$, all lying within the unfolding, represents a periodic orbit (dashed line).
EDIT.
One can see below another example, corresponding to $W=213213$. All starting positions $F$ lead to a periodic orbit of period $6$.
Of course one must check that segment $FF'$ be all inside the unfolding, otherwise we get a wrong path (see figure below). Checking is not difficult: line $FF'$ must cross the sides of the triangles in the order given by $W$, as in figure above. If that is not the case (in figure below the order is $213123$) then $FF'$ goes outside the unfolding.
But for a particular position of $F$ the orbit degenerates into twice an orbit of period $3$, see figure below. In that case we could just use the first three triangles, corresponding to $W=213$.
Hence periodic orbits are also possible when the condition of Lemma 2.8 is not satisfied. But they exist only for some positions of starting point $F$.