[petsc-dev] Double BLAS in PETSc complex

Matthew Knepley knepley at gmail.com
Wed Aug 28 14:53:02 CDT 2013

On Wed, Aug 28, 2013 at 2:36 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:

> 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:

This sounds better.


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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130828/df961342/attachment.html>

More information about the petsc-dev mailing list