When fitting, call the llh function again to make sure parameters are set
Make sure to copy diagonal
Add support for returning expected mean
Average points that are close together when adding new observation