COMP: check in intermediate work on compute dispersion
it compiles, but it isn't complete.