[petsc-dev] various dot products messed up for complex numbers?
    Barry Smith 
    bsmith at mcs.anl.gov
       
    Sun Apr 15 12:55:05 CDT 2012
    
    
  
On Apr 15, 2012, at 12:45 PM, Jed Brown wrote:
> On Sun, Apr 15, 2012 at 12:42, Barry Smith <bsmith at mcs.anl.gov> wrote:
> uses the complex conjugate of the FIRST argument.
> 
> So it looks like everything is major broken for the complex case.
> 
> Is this more than a documentation error in PETSc? I.e. is VecDot() _used_ in such a way that the second argument is actually conjugate, or does it use the standard convention that the first is conjugate?
   Real mathematicians conjugate the second one :-)
   The code to check is PETSc's conjugate gradient, there the conjugate does mater :-).
   Barry
    
    
More information about the petsc-dev
mailing list