[petsc-dev] the complex dot business

Barry Smith bsmith at mcs.anl.gov
Mon Feb 27 10:55:21 CST 2012


   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?

    Barry




More information about the petsc-dev mailing list