What are strategies for generating the following types of pictures:
Here's what's going on here. Take a toda flow in 3 variables. The equations of motion are:
$$\frac{d}{dt}a_k=2(b_k^2-b_{k-1}^2),$$ $$\frac{d}{dt}b_k=b_k(a_{k+1}-a_k),$$
where $a_i,b_i>0$, which can be represented as a commutator involving the matrix $L(a,b)$, which is a tridiagonal matrix with $a_i$ on the diagonal and $b_i$ on the superdiagonals. The hexagon above incorporates the following phenomenon of the Toda flow of $L(a,b)$:
where the edges and vertices are labeled by the structure of the $L$ matrix there. My question is, what are the possible ways of giving this hexagon a sense of distance (thus being able to draw Toda flow trajectories). Essentially this should somehow depend on the relative size of $a_i,b_i$'s, how close $a_i$ is to a particular eigenvalue, etc. For example, the $b_i$ are all exponentially decreasing, so it seems like one needs to look at ratios of $b_i$ to conclude which side of the hexagon we are closer to. This still seems very arbitrary and results in singularities at the boundary (which in higher dimensions become an issue since Toda flows perfectly well cross over such boundaries).
I know that the $n$ dimensional Toda flow is diffeomorphic to the $n-1$ sphere via looking at the vector $(u_{11},u_{12},\ldots,u_{1n})$ of normalized first-components of the corresponding eigenvectors of $L$. Is it as simple as projecting these onto the embedded permutahedron? Here's the picture for $n=4$,
What is a good coordinate system to draw these trajectories on the hexagon/permutahedron?
Reference: ODE's and the Symmetric Eigenvalue Problem