As.matrix for unconstrained parameters

Is there a way to get unconstrained parameters with warmup without the permutations in rstan?

Like extract(fit, unconstrained = TRUE, permuted = FALSE, inc_warmup = TRUE)?

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.