As.matrix for unconstrained parameters

Not really. You can call unconstrain_pars once per iteration or you can specify diagnostic_file and read in the unconstrained parameters from the CSV.