feat: SCEC benchmarks BP6 and BP7

[CusiniM]

Collection of changes and scripts needed for the solution of the SCEC benchmarks.

• Add support for backslip to quasi-dynamic earthquake solver
• Fix initialization of the system, ensuring that the rate-and-state equations are compatible with stresses

[CusiniM]

@victorapm for this physics solver, which uses this MGR strategy, the boomer AMG setup crashes in parallel during setup. Without the filter everything seems to work fine.

The case that I have been using to debug it is: https://github.com/GEOS-DEV/GEOS/blob/develop/inputFiles/lagrangianContactMechanics/LagrangeContactBubbleStab_FixedSlip_smoke.xml and you can reduce the size to :
nx = 4 and ny = 2 and still reproduce the issue. A partition of -x 2 is sufficient to trigger the crash.

[CusiniM]

this is the stack trace of the crash:

Number of element for each fracture state: stick:          380 | new slip:            0 | slip:             0 | open:             0
    Attempt:  0, ConfigurationIter:  0, NewtonIter:  0
        ( Rsolid ) = ( 0.00e+00 )        ( Rt  ) = (    1.549432e+01  )        ( R ) = ( 1.55e+01 )
        MGR preconditioner: numComponentsPerField = [3, 3, 3]
Received signal 11: Segmentation fault: 11

** StackTrace of 0 frames **
Frame 0: _sigtramp 
=====

Received signal 15: Terminated: 15

** StackTrace of 24 frames **
Frame 0: _sigtramp 
Frame 1: opal_progress 
Frame 2: mca_pml_ob1_iprobe 
Frame 3: MPI_Iprobe 
Frame 4: hypre_MPI_Iprobe 
Frame 5: hypre_DataExchangeList 
Frame 6: hypre_ParCSRCommPkgCreateApart_core 
Frame 7: hypre_ParCSRCommPkgCreateApart 
Frame 8: hypre_MatvecCommPkgCreate 
Frame 9: hypre_BoomerAMGBuildModExtPEInterpHost 
Frame 10: hypre_BoomerAMGBuildModExtPEInterp 
Frame 11: hypre_BoomerAMGSetup 
Frame 12: hypre_MGRSetup 
Frame 13: geos::HyprePreconditioner::setup(geos::HypreMatrix const&) 
Frame 14: geos::HypreSolver::setup(geos::HypreMatrix const&) 
Frame 15: geos::PhysicsSolverBase::solveLinearSystem(geos::DofManager const&, geos::HypreMatrix&, geos::HypreVector&, geos::HypreVector&) 
Frame 16: geos::PhysicsSolverBase::solveNonlinearSystem(double const&, double const&, int, geos::DomainPartition&) 
Frame 17: geos::PhysicsSolverBase::nonlinearImplicitStep(double const&, double const&, int, geos::DomainPartition&) 
--------------------------------------------------------------------------
prterun noticed that process rank 1 with PID 5091 on node onix exited on
signal 6 (Abort trap: 6).

[victorapm]

Thanks for reporting. Let me take a look