Add option to output Hessian, disable LP calculations in Laplace

[WardBrian]

Submission Checklist

• [x] Run unit tests: ./runTests.py src/test/unit
• [x] Run cpplint: make cpplint
• [x] Declare copyright holder and open-source license: see below

Summary

Closes #3215 Closes #3260

Intended Effect

This adds two new options to the laplace approximation service function. The first is a boolean calculate_lp, similar to #3249

The second is a structured writer which will save the result of the evaluation at the mode, namely it will produce a JSON file like

{
  "lp_mode" : 0,
  "gradient" : [ 0, 0 ],
  "Hessian" : [ [ -2.77778, 2.22222 ], [ 2.22222, -2.77778 ] ]
}

How to Verify

Tests added

Side Effects

I left the existing function for backwards compatibility. The above arguments default to true and a no-op writer, respectively.

Documentation

Added inline documentation, online doc will come with the cmdstan PR.

Copyright and Licensing

Please list the copyright holder for the work you are submitting (this will be you or your assignee, such as a university or company): Simons Foundation

By submitting this pull request, the copyright holder is agreeing to license the submitted work under the following licenses:

• Code: BSD 3-clause (https://opensource.org/licenses/BSD-3-Clause)
• Documentation: CC-BY 4.0 (https://creativecommons.org/licenses/by/4.0/)

[avehtari]

Great!