found and fixed a bug with set_P_mp that was producing disagreement with set_P. The multi-processing algorithms now all agree with their single CPU counterparts
added model_average routine
added a test of the model_average routine with fixed_bandwidth models