Basic unit testing framework.

Check for file-opening / existence errors (vs segfault). write_mod for general Hamiltonian. --help flag.

Add write_mod option to write population at regular time intervals (two-site/two-allele only).

Locate segfault for small N (population size = 0). Seed as additional run parameter. Concatenate source files?

Special two-site + two-allele (two-locus + binary) Hamiltonian with all parameters passed as command line arguments.

More output for non-epitope runs. Assuage compiler warnings of unused parameters and signed/unsigned int comparison.

Wright-Fisher code by Tom Butler (via John Barton).

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.