[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