<div class="gmail_quote">On Sun, Apr 15, 2012 at 12:55, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":1a6">  Real mathematicians conjugate the second one :-)<br></div></blockquote><div><br></div><div>I never unrderstood that choice. I like vectors to be column vectors, one forms to be row vectors, and operations to have standard fixity. How some mathematicians ended up with infix operators and postfix inner products is beyond me.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":1a6">
<br>
   The code to check is PETSc's conjugate gradient, there the conjugate does mater :-).</div></blockquote></div><div><br></div><div>Looks right to me.</div><br><div><div>       ierr = KSP_MatMult(ksp,Amat,P,W);CHKERRQ(ierr);          /*     w <- Ap         */</div>
<div>       ierr = VecXDot(P,W,&dpi);CHKERRQ(ierr);                  /*     dpi <- p'w     */</div></div><div><br></div>