all xi-coalescents work correctly it appears. Furthermore, I have implemented a new algorithm for generating subpartitions that scales better. With it, I managed to compute the SFS of a xi_beta coalescent for n=50 in only 14 hours and 30 min!!!