[petsc-users] Question on SLEPc + computing SVD with a "matrix free" matrix

[Jose E. Roman]

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 <petsc-users at mcs.anl.gov> 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