[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