Just use the GNU Scientific Library (GSL) for RNG. (But also include my attempt at copying over the relevant bits.)

Add Mersenne Twister RNG library/header. Generalize to stream<< any std::map<T> as CSV-like string.

Print mutset_t (to stream using operator<<) as CSV-like string.

Start of Wright-Fisher dynamics using STL. Mutations stored as std::set, species counts in std::map.