added extra parameter to put penalty on estimated ellipsoids