Add option crb.solve-fem-monolithic=true to compare RB approximation using
affine parameter decomposition to FEM approximation using monolithic
formulation.
Need the model implements computeMonolithicFormulation() to use it.
But this function is optional and returns the matrix and rhs to use.
Up the model to use option crb.solve-fem-monolithic=true