sbaij & hermitian mat-vec multiplication

Hong Zhang hzhang at mcs.anl.gov
Fri Oct 30 08:45:39 CDT 2009


Christos,

> Hi guys. Following up on a previous user request, I've been looking  
> at supporting hermitian mat-vec multiplication for the SBAIJ  
> implementation. I haven't played with petsc before, so please  
> forgive me for not understanding the internals. So, SBAIJ uses the  
> transpose of the upper band where access to the lower band is  
> needed. Am I right to believe that what is needed to be done is to  
> have the SBAIJ data retriever (such as MatGetValues_MPISBAIJ)  
> modified such that on sub-diagonal accesses the conjugate of the  
> original values are returned instead (possibly by updating the  
> contents of the "v" argument once filled in)?

Yes.
Be careful in parallel implementation though,

Hong



More information about the petsc-dev mailing list