On Oct 22, 2008, at 9:15 AM, zhifeng sheng wrote: > 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. >
KSPCG and use KSPCGSetType(ksp,KSP_CG_HERMITIAN ) > By the way, if I really need conjugate transpose function for > complex matrices, must I implement it myself? yes > > (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 >>>>> >>>>> > >>>> >>> >> >
