Numerically check gradient computations at a few values
Compute obj.fn. and gradient for softmax, and check gradients numerically
Updated convergence timing