<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
On 06.05.2012 15:34, Matthew Knepley wrote:
<blockquote
cite="mid:CAMYG4G=T_XL_rxriG6yxXZni+-CwbnxWPDKQY_NcH=vWwu86bQ@mail.gmail.com"
type="cite">On Sun, May 6, 2012 at 9:24 AM, Alexander Grayver <span
dir="ltr"><<a moz-do-not-send="true"
href="mailto:agrayver@gfz-potsdam.de" target="_blank">agrayver@gfz-potsdam.de</a>></span>
wrote:<br>
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
0.8ex; border-left: 1px solid rgb(204, 204, 204);
padding-left: 1ex;">
<div bgcolor="#ffffff" text="#000000"> Hm, valgrind gives a
lot of output like that (see full log in previous message):</div>
</blockquote>
<div><br>
</div>
<div>Can you run this with --download-f-blas-lapack? This sounds
much more like an MKL bug.</div>
</div>
</blockquote>
<br>
I did:<br>
--download-scalapack --download-blacs --download-blas-lapack
--with-precision=double --with-scalar-type=complex<br>
<br>
The error is still there. I checked "ldd solveTest", mkl is not used
for sure. This is not an MKL problem I guess:<br>
<br>
==13600== Invalid read of size 8<br>
==13600== at 0x58636AF: dlasq3_ (in
/usr/local/lib/liblapack.so.3.2.2)<br>
==13600== by 0x5862C84: dlasq2_ (in
/usr/local/lib/liblapack.so.3.2.2)<br>
==13600== by 0x5861F2C: dlasq1_ (in
/usr/local/lib/liblapack.so.3.2.2)<br>
==13600== by 0x571A479: zbdsqr_ (in
/usr/local/lib/liblapack.so.3.2.2)<br>
==13600== by 0x57466A7: zgesvd_ (in
/usr/local/lib/liblapack.so.3.2.2)<br>
==13600== by 0x694687: KSPComputeExtremeSingularValues_GMRES
(gmreig.c:46)<br>
==13600== by 0x69C62A: KSPComputeExtremeSingularValues
(itfunc.c:47)<br>
==13600== by 0x44E02C: main (solveTest.c:62)<br>
==13600== Address 0x10826b90 is 16 bytes before a block of size 832
alloc'd<br>
==13600== at 0x4C25D66: memalign (vg_replace_malloc.c:694)<br>
==13600== by 0x4B5ACB: PetscMallocAlign (mal.c:30)<br>
==13600== by 0x686181: KSPSetUp_GMRES (gmres.c:73)<br>
==13600== by 0x69E6FE: KSPSetUp (itfunc.c:239)<br>
==13600== by 0x6A090C: KSPSolve (itfunc.c:402)<br>
==13600== by 0x44E009: main (solveTest.c:61)<br>
<br>
The weird thing is that the it gives correct result, so zgesvd works
fine. <br>
<br>
And also running this program with 10 iterations in valgrind doesn't
produce error. The low above is with 100 iterations.<br>
Without valgrind the error is always there. <br>
<br>
<pre class="moz-signature" cols="72">--
Regards,
Alexander</pre>
</body>
</html>