On Mon, Feb 27, 2012 at 10:55 AM, Barry Smith <[email protected]> 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