Alexander, Sorry, never mind. I found the code in the earlier email.,
barry On May 6, 2012, at 12:49 PM, Barry Smith wrote: > > Alexander, > > No attached test program. (only valgrind log). > > Barry > > On May 6, 2012, at 8:21 AM, Alexander Grayver wrote: > >> 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 >> >> <valgrind.zip> >