<div dir="ltr">On Wed, Aug 28, 2013 at 2:36 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5">Satish Balay <<a href="mailto:balay@mcs.anl.gov">balay@mcs.anl.gov</a>> writes:<br>

<br>
> On Wed, 28 Aug 2013, Matthew Knepley wrote:<br>
><br>
>> The last error I have is needing the double version of BLAS when we are<br>
>> using complex. How do we handle that?<br>
><br>
> Hm - I'm not sure.. Perhaps we need to add in a new prefix for<br>
> Blas(PetscReal) [and corresponding macros?] to petscblaslapack.h?<br>
> The #ifdef login in there will get more complex...<br>
<br>
</div></div>What if we _always_ defined real and complex prototypes, so the only<br>
#ifdef logic would be in precision.  Then we define the scalar versions<br>
something like this:<br>
<br>
#define BLAS_Scalar_GEMM  SCALAR_SWITCH(BLAS_Real_GEMM,BLAS_Complex_GEMM)<br>
</blockquote></div><br>This sounds better.</div><div class="gmail_extra"><br></div><div class="gmail_extra">   Matt<br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>