Faster support function of template reach-set

[mforets]

Currently querying the support function of a TemplateReachSet is suboptimal because it recomputes instead of using the stored field sf. The difficulty is that the template directions dirs do not necessarily relate to variables, although in practice, we often set dirs to be a canonical vector such as [0, 1, 0, .., 0].