Lisandro : > BTW, do you fully agree with my previous fix? Should MatIsTranspose() > generate an error if the test is not implemented? Or should it just > return PETSC_FALSE?
I agree with you: generate an error if the test is not implemented - consistent with entire petsc style 'PETSC_FALSE' is simply one of output of this routine. Hong >> On Wed, Apr 27, 2011 at 2:04 PM, Lisandro Dalcin <dalcinl at gmail.com> >> wrote: >>> 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 >>> >> > > > > -- > 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 >
