Release v1.5.1

• Add new return values soln.xmin_eval_num and soln.jacmin_eval_nums
• Allow option for parallel initial evaluations for non-random directions if npt not too large

Neither should affect the performance of the DFO-LS algorithm