[petsc-users] Non-hermitian KSP

[Bikash Kanungo]

Hi,


I'm trying to solve a KSP using GMRES in which the matrix (say, A) is
non-hermitian. I'm using a shell matrix for A.

In this context I have the following few questions:


   1. I guess one needs to provide appropriate user-defined function
   override for MatMultHermitianTranspose instead of MatMultTranspose (i.e.,
   MATOP_MULT_TRANSPOSE). However, I am unable to find any
   MATOP_MULT_HERMITIAN_TRANSPOSE flag. So if GMRES expects
   MatlMultHermitianTranspose then should I override the MatMultTranspose in a
   way to actually perform MatMultHermitianTranspose? Or am I mistaken in my
   assumption that GMRES requires MatMultHermitianTranspose instead of
   MatMultTranspose?

   2. For the same problem, I want to provide my own preconditioner matrix
   (say, P) as a shell a matrix. What operations (e.g., MATOP_MULT,
   MATOP_MULT_TRANSPOSE etc.) should I override through user-defined
   functions?


Thanks,
Bikash

-- 
Bikash S. Kanungo
PhD Student
Computational Materials Physics Group
Mechanical Engineering
University of Michigan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170220/6113bfd5/attachment.html>