Fixed bug where rejection sampling was never used for the smoothing.
Caching of some matrices in mixed_nl_gaussian to improve performance.
Improvements to gradient_test.py
Changed number of smoothed trajectories in test_param_trans.py