Default callback function for printing optimization progress in LBFGS

As the disp and iprint options are deprecated with the new scipy version, it currently doesnt print the optimization progress and is hard to track for a long optimization run.

We already have an option for callback function, and we could add some default printing function to track the progress.