Hello,
Thanks for your response.
@David
Reducing the number of Aq matrices is not really ideal for my application.
I would like to try and see if I can reduce the memory requirement in some
other way.
@Paul
Indeed PetscMatrix is used for Reduced Basis also. I went through the
PetscMatrix class d
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
