Gridap and solvers from DifferentialEquations.jl

Hello,

I am using `Gridap` to solve a coupled system of differential and algebraic equations (DAE). Namely, I am dealing with the
diffusion equation in which the concentration field is coupled to some external fields that change on time scales much faster
than that of the diffusive problem. Currently, I explicitly discretize the time-propagation scheme and use an iterative scheme to propagate the solution a timestep $\Delta t$ forward in time. It should probably also be evident from the above that the problem is a multiphysics one.

For speed and convergence reasons, I would like to use a specialized DAE solver instead, such as the ones available in
[DifferentialEquations.jl](https://docs.sciml.ai/DiffEqDocs/stable/) suite. Is there a way to achieve this in `Gridap`?

I stumbled upon a [similar issue in GridapODEs](https://github.com/gridap/GridapODEs.jl/issues/4) where some attempt at a solution has been made within the module `DiffEqsWrappers`. Could I achieve the desired outcome if I [went along the same lines](https://github.com/gridap/Gridap.jl/blob/master/src/ODEs/DiffEqsWrappers/DiffEqsWrappers.jl)? I've noticed that since integration of `GridapODEs` to `Gridap` the module `DiffEqsWrappers` is no longer used:

```julia

"""

The exported names are
$(EXPORTS)
"""
module ODEs

using DocStringExtensions

include("ODETools/ODETools.jl")

include("TransientFETools/TransientFETools.jl")

# include("DiffEqsWrappers/DiffEqsWrappers.jl")

end #module

const GridapODEs = ODEs

```

Regards,
Jan