Resolvent operator with FF/PETSc (MatMatSolve?)

OK, I understand. I am not so familiar with tgv before.

Thanks very much for your quick reply!

Hi Chris,

I went through your code. Can you pinpoint the file where you calculate the inverse of R=(iωMq+J)āˆ’1. I am not that much experienced with PETSc and I know this is too much to ask, but if you can help it will be helpful.

thanks

See here. Note that we never compute the resolvent operator explicitly. We only compute its action on a vector.