Constrained Covariance

[spinkney]

Try adding a small epsilon to the diagonal. See for more info Cholesky decomposition failed for precision/covariance matrix - #5 by anon75146577

...
matrix[N_eta, N_eta] zeta_Lcov = add_diag(cholesky_decompose(zeta_cov), 1e-6);