Dear all suppose I have a complex matrix (Hermitian positive definite) to solve, which KSP solver(s) can support solving it? PS: for this moment, I don't need to take into account that the matrix is hermitian.
By the way, if I really need conjugate transpose function for complex matrices, must I implement it myself? (I mean, is there any function that I can make use of?) Thanks Best regards Zhifeng Barry Smith wrote: > > There is only support for CG with Hermitian transpose, ksp_cg_type > symmetric or hermitian, KSPCGSetType() > the others only support complex, no Hermitian transpose. > > Barry > > On Oct 22, 2008, at 2:40 AM, zhifeng sheng wrote: > >> Hi, >> >> you mean the conjugate transpose for complex matrix is not supported? >> >> then how can you implement the iterative solvers for complex >> matrices? because, some iterative solvers need it. >> >> Thanks >> Best regards >> >> >> >> Hong Zhang wrote: >>> >>> Zhifeng, >>> >>> We do not have support for matrix operations on Hermitian matrix yet. >>> Hong >>> >>> On Mon, 20 Oct 2008, zhifeng sheng wrote: >>> >>>> Dear all >>>> >>>> I am using this MatMatMultTranspose function for complex matrices, >>>> but it seems to be doing something weird. >>>> >>>> for instance, if I have complex matrix A, and I compute A^T*A with >>>> this function, it does not generate a Hermitian >>>> <http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatIsHermitian.html> >>>> >>>> matrix. >>>> >>>> I am thinking that maybe the function take the transpose of A >>>> instead of the conjugate transpose .... >>>> >>>> Do you know how I can get an A^H*A instead of A^T*A for complex >>>> matrices? >>>> >>>> Thanks a lot >>>> Best regards >>>> Zhifeng >>>> <http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatIsHermitian.html> >>>> >>>> >>> >> >
