Use `NonlinearAliasSpecifier` to specify aliasing during solve

[jClugstor]

SciMLBase has a new interface for specifying which variables are aliased during a solve. I would like to put this in to use in NonlinearSolve by using the NonlinearAliasSpecifier type from SciMLBase to specify which variables are aliased. Similar to SciML/LinearSolve.jl#564, or SciML/OrdinaryDiffEq.jl#2503.

Putting this here as a tracker/reference, but also the way that NonlinearSolve chooses which variables are aliased is different than some of the other SciML packages, so I'm a bit lost on how this would work. For example the aliasing for OrdinaryDiffEq is decided in the __init function, but that's not the case for NonlinearSolve.

[ChrisRackauckas]

Was this completed?

[jClugstor]

SciML/DiffEqBase.jl#1138 puts the basic functionality in

It only lets us toggle alias_u0 though, nothing else.