[petsc-dev] Error during KSPDestroy

Alexander Grayver agrayver at gfz-potsdam.de
Sun May 6 09:16:57 CDT 2012


On 06.05.2012 15:34, Matthew Knepley wrote:
> On Sun, May 6, 2012 at 9:24 AM, Alexander Grayver 
> <agrayver at gfz-potsdam.de <mailto:agrayver at gfz-potsdam.de>> wrote:
>
>     Hm, valgrind gives a lot of output like that (see full log in
>     previous message):
>
>
> Can you run this with --download-f-blas-lapack? This sounds much more 
> like an MKL bug.

I did:
--download-scalapack --download-blacs --download-blas-lapack 
--with-precision=double --with-scalar-type=complex

The error is still there. I checked "ldd solveTest", mkl is not used for 
sure. This is not an MKL problem I guess:

==13600== Invalid read of size 8
==13600==    at 0x58636AF: dlasq3_ (in /usr/local/lib/liblapack.so.3.2.2)
==13600==    by 0x5862C84: dlasq2_ (in /usr/local/lib/liblapack.so.3.2.2)
==13600==    by 0x5861F2C: dlasq1_ (in /usr/local/lib/liblapack.so.3.2.2)
==13600==    by 0x571A479: zbdsqr_ (in /usr/local/lib/liblapack.so.3.2.2)
==13600==    by 0x57466A7: zgesvd_ (in /usr/local/lib/liblapack.so.3.2.2)
==13600==    by 0x694687: KSPComputeExtremeSingularValues_GMRES 
(gmreig.c:46)
==13600==    by 0x69C62A: KSPComputeExtremeSingularValues (itfunc.c:47)
==13600==    by 0x44E02C: main (solveTest.c:62)
==13600==  Address 0x10826b90 is 16 bytes before a block of size 832 alloc'd
==13600==    at 0x4C25D66: memalign (vg_replace_malloc.c:694)
==13600==    by 0x4B5ACB: PetscMallocAlign (mal.c:30)
==13600==    by 0x686181: KSPSetUp_GMRES (gmres.c:73)
==13600==    by 0x69E6FE: KSPSetUp (itfunc.c:239)
==13600==    by 0x6A090C: KSPSolve (itfunc.c:402)
==13600==    by 0x44E009: main (solveTest.c:61)

The weird thing is that the it gives correct result, so zgesvd works fine.

And also running this program with 10 iterations in valgrind doesn't 
produce error. The low above is with 100 iterations.
Without valgrind the error is always there.

-- 
Regards,
Alexander

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/a7b674d7/attachment.html>


More information about the petsc-dev mailing list