On Wed, Aug 28, 2013 at 2:36 PM, Jed Brown <[email protected]> wrote:
> Satish Balay <[email protected]> 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) > This sounds better. Matt -- 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
