Yes, you have to set it with SVDSetImplicitTranspose(), but then the matrix-free matrix should implement both "A.v" and "A'.v" operations. Jose
> El 24 jun 2020, a las 23:25, Ernesto Prudencio via petsc-users > <[email protected]> escribió: > > Hi, > > Is it possible to compute a SVD for a “matrix free” matrix? > > At first, it seems it would be ok with a MatCreateMFFD() and > MatMFFDSetFunction(), because one could then provide the routine that > computes “A . v” for any given v, which is an operation needed by SLEPc. > However, one would also need to set up the SVD object in SLEPc with an > implicit transpose. Would that be possible? > > Thanks in advance, > > Ernesto. > > Schlumberger-Private
