implement memory friendly fitting on the simplex in cython
general simplex fitting to non-Euclidean stabilizer
add partial volume C routines