1. Add OpenMP parallelism to gradient calculation. 2. Add code for computing approximate Hessian, central-limit uncertainties. 3. Fix a bug where we were specifying the transition probabilities backward in the theta vector. Not really a bug, but it meant that theta was actually not the equilibrium probabilities, as I had indicated.