---------------
-- Phase plots:
---------------

-- a plot of two dimensions:
phasePlot2 :: LA.Vector Double
           -> LA.Matrix Double
           -> [Species]
           -> (Species,Species)
           -> IO ()
phasePlot2 ts soln ss ss'
    = plotPhase
      (filter (\(s,_)-> (s == (specName (fst ss'))) || s == (specName (snd ss')))
       (zip (map specName ss) (map LA.toList (LA.toColumns soln))))
phasePlot2 :: LA.Vector Double
           -> LA.Matrix Double
           -> [Species]
           -> (Species,Species)
           -> IO ()
phasePlot2 ts soln ss ss'
    = plotPhase
      (filter (\(s,_)-> (s == (specName (fst ss'))) || s == (specName (snd ss')))
       (zip (map specName ss) (map LA.toList (LA.toColumns soln))))
phasePlot2 ts soln ss ss'
    = plotPhase
      (filter (\(s,_)-> (s == (specName (fst ss'))) || s == (specName (snd ss')))
       (zip (map specName ss) (map LA.toList (LA.toColumns soln))))