Flowchart of the algorithm followed for the computation of LCS. We note that it requires integration of two dynamical systems. First, the flow of the original dynamical system on a grid must be computed, in order to estimate its Jacobian. After that, we take as initial condition (i.c.) local extrema of the singular values and integrate the field given by the corresponding right-singular vectors.