[petsc-dev] Error during KSPDestroy

[Alexander Grayver]

On 06.05.2012 15:07, Matthew Knepley wrote:
>
>>         Hello,
>>
>>         I use KSP and random rhs to compute largest singular value:
>>
>>
>>     1) Is this the whole program? If not, this can be caused by
>>     memory corruption somewhere else. This is what I suspect.
>
>     Matt,
>
>     I can reproduce error using attached test programm and this matrix
>     (7 mb):
>     http://dl.dropbox.com/u/60982984/A.dat
>
>
> I run it fine with the latest petsc-dev:
>
>   1.405802e+00
>
> Can you valgrind it on your machine?
>

I did:
valgrind --tool=memcheck -q --num-callers=20 --log-file=valgrind.log.%p 
/solveTest -ksp_monitor_true_residual -log_summary -mat_type aij 
-ksp_rtol 1.0e-10 -malloc off

The error is better constrained:

==20287== Invalid read of size 8
==20287==    at 0x7874B4C: opal_os_path (in 
/opt/mpi/intel/openmpi-1.4.2/lib/libopen-pal.so.0.0.0)
==20287==    by 0x75F2E27: orte_session_dir_finalize (in 
/opt/mpi/intel/openmpi-1.4.2/lib/libopen-rte.so.0.0.0)
==20287==    by 0x76012E8: orte_errmgr_base_error_abort (in 
/opt/mpi/intel/openmpi-1.4.2/lib/libopen-rte.so.0.0.0)
==20287==    by 0x73396E9: ompi_mpi_abort (in 
/opt/mpi/intel/openmpi-1.4.2/lib/libmpi.so.0.0.2)
==20287==    by 0x734F36E: PMPI_Abort (in 
/opt/mpi/intel/openmpi-1.4.2/lib/libmpi.so.0.0.2)
==20287==    by 0x7499AB: PetscDefaultSignalHandler (signal.c:169)
==20287==    by 0x749267: PetscSignalHandler_Private (signal.c:53)
==20287==    by 0x924B9DF: ??? (in /lib64/libc-2.11.1.so)
==20287==    by 0x535D9E: VecDestroyVecs (vector.c:653)
==20287==    by 0x68B61D: KSPReset_GMRES (gmres.c:258)
==20287==    by 0x6A9D39: KSPReset (itfunc.c:733)
==20287==    by 0x6AA839: KSPDestroy (itfunc.c:780)
==20287==    by 0x4509F8: main (solveTest.c:66)
==20287==  Address 0xbde4860 is 0 bytes inside a block of size 2 alloc'd
==20287==    at 0x4C26B9B: malloc (vg_replace_malloc.c:263)
==20287==    by 0x92876DF: vasprintf (in /lib64/libc-2.11.1.so)
==20287==    by 0x9266C67: asprintf (in /lib64/libc-2.11.1.so)
==20287==    by 0x75F1701: orte_util_convert_vpid_to_string (in 
/opt/mpi/intel/openmpi-1.4.2/lib/libopen-rte.so.0.0.0)
==20287==    by 0x75F2D4A: orte_session_dir_finalize (in 
/opt/mpi/intel/openmpi-1.4.2/lib/libopen-rte.so.0.0.0)
==20287==    by 0x76012E8: orte_errmgr_base_error_abort (in 
/opt/mpi/intel/openmpi-1.4.2/lib/libopen-rte.so.0.0.0)
==20287==    by 0x73396E9: ompi_mpi_abort (in 
/opt/mpi/intel/openmpi-1.4.2/lib/libmpi.so.0.0.2)
==20287==    by 0x734F36E: PMPI_Abort (in 
/opt/mpi/intel/openmpi-1.4.2/lib/libmpi.so.0.0.2)
==20287==    by 0x7499AB: PetscDefaultSignalHandler (signal.c:169)
==20287==    by 0x749267: PetscSignalHandler_Private (signal.c:53)
==20287==    by 0x924B9DF: ??? (in /lib64/libc-2.11.1.so)
==20287==    by 0x535D9E: VecDestroyVecs (vector.c:653)
==20287==    by 0x68B61D: KSPReset_GMRES (gmres.c:258)
==20287==    by 0x6A9D39: KSPReset (itfunc.c:733)
==20287==    by 0x6AA839: KSPDestroy (itfunc.c:780)
==20287==    by 0x4509F8: main (solveTest.c:66)

Full log is attached.

Important.
If I comment this line:
KSPComputeExtremeSingularValues(ksp, &maxx, &minx);

It works.

-- 

Regards,
Alexander

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/7d19a9e4/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: valgrind.zip
Type: application/x-zip-compressed
Size: 4572 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120506/7d19a9e4/attachment.bin>