Collapse and expand auxilliary variables in nonlinear gradients etc

[johanlofberg]

Only work in non-auxilliary variables in the solvers, use auxilliary temporaily only as function values and gradients are computed