Work-around for normal model fitting

The way it is currenty set up the MLE functions doen't seem to respect the box
constraints. For now, just return a very very large number when some values are
negative
Use default MLE algo if we aren't using box constraints