Hi Nikhil,
The RB method is really not intended to be used with such a large number of
Aq matrices, since the computational cost of RB grows quite fast with the
number of terms in the affine expansion (e.g. IIRC the error bound used in
the greedy depends on the number of terms to the power of 4).
Hi Nikhil,
Typically, you would grab the "raw" underlying Mat object and do the PETSc
call yourself. This would be accessible from the libMesh::PetscMatrix
object. Typically the ImplicitSystem has the system matrix, but I'm not
familiar with where these matrices would be cached on the RB side of
t
