[Octave-bug-tracker] [bug #32120] complex times complex failure

[Michael C. Grant]

Follow-up Comment #17, bug #32120 (project octave):

OK, here's a patch that seems to do the trick.

Here's what it does: in ax_blas.m4, it creates a new define USE_CBLAS_COMPLEX
that is set to 1 when the "-framework vecBlas" version of the BLAS is
detected. This define is then used in CMatrix.cc and fCMatrix.cc to replace
the calls to xzdotu, etc. to corresponding CBLAS calls.

I've tested this and found it does indeed produce the correct results for both
single and double precision.

It seems to me that it might be a good idea to generalize this: instead of
defining USE_CBLAS_COMPLEX for the specific case of vecLib, you could test for
the existence of these functions in whatever BLAS library was successfully
called, and then define USE_CBLAS_COMPLEX whenever they are found.

By the way, USE_BLASWRAP and the corresponding code in blaswrap.c should
probably be dropped. It seems to me this is much less fragile solution.

(file #30922)
    _______________________________________________________

Additional Item Attachment:

File name: cblas.diff                     Size:4 KB


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?32120>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/