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>