[petsc-dev] the complex dot business

[Matthew Knepley]

On Mon, Feb 27, 2012 at 10:55 AM, Barry Smith <bsmith at mcs.anl.gov> wrote:

>
>   I'm afraid the sudden appearance of the blas dot for complex numbers
> problem is my fault.
>
>   Traditionally since 199x we never used the complex dot product from blas
> because of the return complex type problem. VecDot_Seq() used #ifdef for
> complex and a simple loop.
>
>   Since some blas norm suck I off-hand told Shri to use dot for norm and
> wasn't thinking the complex case.
>
>   Since then, of course, the complex dot problem has been coming up.
>
>   Maybe it is best to continue to never use complex blas dot and just
> change the VecNorm_ routines to use norm for complex and not dot.? Just
> throw away this silly problem and time sink?


I just pushed a fix for this. The pthread stuff will still have to be fixed
since I can't build it.

   Matt


>
>    Barry
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120227/30ed833a/attachment.html>