Consistently use ycur for user RHS routines

[drreynolds]

This PR standardizes ARKODE calls to the user-supplied right-hand-side functions so that wherever possible, these RHS calls are provided the user-supplied solution vector that was passed to ARKodeEvolve.

The notable exceptions to this rule are:

1. the Hermite temporal interpolation module,
2. the provided preconditioners ARKBANDPRE and ARKBBDPRE,
3. banded or dense linear solvers with automatically-approximated Jacobian matrices,
4. iterative linear solvers with automatically-approximated Jacobian-times-vector products,
5. temporal root-finding,
6. discrete adjoint modules in ARKStep or ERKStep,
7. the SPRKStep stepper, and
8. LSRKStep's use of the automated dominant eigenvalue estimation module.

However, for most ARKODE integrators is it now possible for users to consistently receive back their own state vector for RHS calls, allowing them to leverage auxiliary data stored therein (e.g., derived quantities or MPI communication buffers).

[drreynolds]

This is the fourth in a sequence of PRs for ARKODE to support our FES collaborators. Please review this one after #815. When the stack is complete, we'll do a final PR into develop

[drreynolds]

This PR should be ready to review -- all CI issues resulting from the new release have been resolved.