<div class="gmail_quote">On Thu, Apr 5, 2012 at 04:54, Alexander Grayver <span dir="ltr"><<a href="mailto:agrayver@gfz-potsdam.de">agrayver@gfz-potsdam.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":3vi">The inconsistency is still there. I pulled petsc-dev today and came across with the same problem:<br>
<br>
src/vec/vec/impls/seq/bvec2.c<br>
239: *z = PetscRealPart(BLASdot_(&bn,xx,<u></u>&one,xx,&one));<br>
240: *z = PetscSqrtReal(*z);<br>
<br>
This doesn't work with Intel MKL. However in one of the mailing list discussions it was said there is now check for that.<br>
This check didn't work for me apparently. Any ideas why?</div></blockquote></div><br><div>As always, we need configure.log.</div><div><br></div><div>In include/petscblaslapack_uscore.h (and other places), you will find the following. So we need to know why PETSC_COMPLEX_DOT_RESULT_ARGUMENT was not found by configure.</div>
<div><br></div><div><div>#ifdef PETSC_COMPLEX_DOT_RESULT_ARGUMENT</div><div>EXTERN_C_BEGIN</div><div>extern void zdotc_(PetscScalar *,const PetscBLASInt*,const PetscScalar*,const PetscBLASInt*,const PetscScalar*,const PetscBLASInt*);</div>
<div>PETSC_STATIC_INLINE PetscScalar BLASdot_(const PetscBLASInt *n,const PetscScalar *x,const PetscBLASInt *sx,const PetscScalar *y,const PetscBLASInt *sy) {</div><div> PetscScalar tmpz;</div><div> zdotc_(&tmpz,n,x,sx,y,sy);</div>
<div> return tmpz;</div><div>}</div><div>EXTERN_C_END</div><div>#else</div><div># define BLASdot_ zdotc_</div><div>#endif</div></div>