[WIP] Move iterative solver interface and internals to Thyra

[maxfirmbach]

Description and Context

This PR aims to move the iterative solver internals to Thyra with the goal of:

• Removing the inheritance of LinalgPrecondOperator from Epetra_Operator
• Removing the necessity to call epetra_operator() during solve with overwritten methods like Apply() for our block sparse matrix implementation, which should cut down the need of inheritance of BlockSparseMatrix from Epetra_Operator
• Remove the need to store the preconditioner object as Epetra_Operator
• Moves the LinalgPrecondOperator directly to the projection based preconditioner solver internals
• Enables the use of individual linear solver status tests on specific multi-physics blocks with GMRES
• Simplifies some code paths

I put this in draft PR mode as I already know that this does not work with AMGnxn and with the monolithic multigrid from MueLu. Currently this PR serves to see how many tests can pass with the changes made.

Related Issues and Pull Requests

Related to #1400