[petsc-dev] change to MatIsTranspose() and failing test

[Hong Zhang]

 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
>