OK, I understand. I am not so familiar with tgv before.
Thanks very much for your quick reply!
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.