Updating quasi-Newton matrices with limited storage

Updating quasi-Newton matrices with limited storage