during the multiplication step. Also I think it's strange that the sum (and sum_sq) functions use the mean values instead of just using the sum function of the miller/flex array. Doesn't this induce some overhead? With regards, Dipl. Phys. Jan M. Simons Institute of Crystallography RWTH Aachen University