[petsc-dev] Double BLAS in PETSc complex

Jed Brown jedbrown at mcs.anl.gov
Wed Aug 28 14:36:04 CDT 2013


Satish Balay <balay at mcs.anl.gov> writes:

> On Wed, 28 Aug 2013, Matthew Knepley wrote:
>
>> The last error I have is needing the double version of BLAS when we are
>> using complex. How do we handle that?
>
> Hm - I'm not sure.. Perhaps we need to add in a new prefix for
> Blas(PetscReal) [and corresponding macros?] to petscblaslapack.h?
> The #ifdef login in there will get more complex...

What if we _always_ defined real and complex prototypes, so the only
#ifdef logic would be in precision.  Then we define the scalar versions
something like this:

#define BLAS_Scalar_GEMM  SCALAR_SWITCH(BLAS_Real_GEMM,BLAS_Complex_GEMM)
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130828/a24cc11b/attachment.sig>


More information about the petsc-dev mailing list