After wasting as considerable amount of time debugging issues related to uninitialized variables passed to MatIsTranspose, I've pushed this:
changeset: 18869:3bdc748ba073 tag: tip user: Lisandro Dalcin <dalcinl at gmail.com> date: Wed Apr 27 15:59:56 2011 -0300 summary: MatIsTranspose(): error if any matrix does not support the test And now, the best below fails, becase SeqBAIJ does not support transpose test. [dalcinl at trantor tests]$ cd $PETSC_DIR/src/ksp/ksp/examples/tests [dalcinl at trantor tests]$ make runex32_aij_sbaij_icc make: [runex32_aij_sbaij_icc] Error 56 (ignored) 1,19c1,7 < [0]PETSC ERROR: --------------------- Error Message ------------------------------------ < [0]PETSC ERROR: No support for this operation for this object type! < [0]PETSC ERROR: Matrix of type <seqbaij> does not support checking for transpose! < [0]PETSC ERROR: ------------------------------------------------------------------------ < [0]PETSC ERROR: Petsc Development HG revision: a02bb9bcaf33cf16868e2771c3c8306d52235c69 HG Date: Fri Apr 22 13:15:40 2011 -0500 < [0]PETSC ERROR: See docs/changes/index.html for recent updates. < [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. < [0]PETSC ERROR: See docs/index.html for manual pages. < [0]PETSC ERROR: ------------------------------------------------------------------------ < [0]PETSC ERROR: ./ex32 on a arch-linu named trantor by dalcinl Wed Apr 27 16:00:08 2011 < [0]PETSC ERROR: Libraries linked from /home/devel/petsc/dev/arch-linux2-c-debug/lib < [0]PETSC ERROR: Configure run at Sat Apr 23 17:40:22 2011 < [0]PETSC ERROR: Configure options --with-shared-libraries=1 PETSC_ARCH=arch-linux2-c-debug < [0]PETSC ERROR: ------------------------------------------------------------------------ < [0]PETSC ERROR: MatIsTranspose() line 4383 in /home/devel/petsc/dev/src/mat/interface/matrix.c < [0]PETSC ERROR: main() line 69 in src/ksp/ksp/examples/tests/ex32.c < application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0 < rank 0 in job 3534 trantor_56062 caused collective abort of all ranks < exit status of rank 0: return code 56 --- > 0 KSP Residual norm 0.210648 > 1 KSP Residual norm 0.07826 > 2 KSP Residual norm 0.0115565 > 3 KSP Residual norm 0.00113989 > 4 KSP Residual norm 0.00013248 > 5 KSP Residual norm 9.04476e-06 > 6 KSP Residual norm 4.52149e-07 Possible problem with ex32_aij_sbaij_icc, diffs above -- Lisandro Dalcin --------------- CIMEC (INTEC/CONICET-UNL) Predio CONICET-Santa Fe Colectora RN 168 Km 472, Paraje El Pozo 3000 Santa Fe, Argentina Tel: +54-342-4511594 (ext 1011) Tel/Fax: +54-342-4511169
