feat: add rate- and state-dependent friction model & quasi-dynamic solver.

[CusiniM]

• [x] Add Rate- and State-dependent friction model
• [x] Create a QuasiDynamicEQ solver that drives EQ simulations by:
  • uses a stress solver to update tractions on the fault (It can be any of the fracture solvers).
  • launches kernels to solve rate-and-state equations on each fault element.
• [x] If no stress solver is specified, add the possibility to add function to compute tractions on the fault (e.g., the spring-slider shear traction formula)
• [x] Add rate- and state kernels that, for now, use Euler-backward to solve the R&S equations (for now 1d on the fault).
• [x] Add spring-slider xml example and integratedTest.

For now, everything is 1D on the fault. The 2D case and the coupling to the contact solver with Lagrange multiplier will be tackled in a followup PR.

LvArray:: https://github.com/GEOS-DEV/LvArray/pull/329