revised user-specified initial inverse Hessian Code

This is exactly the code discussed in the previous PR, but built off the most recent master branch.
``nothing`` isn't a ``Matrix``, so this method is not callable.