[petsc-dev] example broken for a while src/mat/examples/tests/ex104.c

[Hong Zhang]

On my MacAir without elemental installed,
I get clean 'make runex104'
and anticipated message from
'make runex104_elemental'
0a1,28
> [0]PETSC ERROR: --------------------- Error Message ------------------------------------
> [0]PETSC ERROR: Unknown type. Check for miss-spelling or missing external package needed for type:
> see http://www.mcs.anl.gov/petsc/documentation/installation.html#external!
> [0]PETSC ERROR: Unknown Mat type given: elemental!
> [0]PETSC ERROR: ------------------------------------------------------------------------

As Jed suggested, the error on your Mac might be result
of rounding problem in MatEqual.
We may add 'PETSC_USE_INFO' statement
in MatEqual_xxx() for such info.

Hong

On Tue, Jan 15, 2013 at 5:11 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
> On Tue, Jan 15, 2013 at 5:10 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
>>
>>
>> On Jan 15, 2013, at 5:03 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
>>
>> > On Tue, Jan 15, 2013 at 4:43 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
>> > > It works for me.
>> >
>> >    Cause you are running an install with Elemental, without elemental
>> > install it goes bad.
>> >
>> > Isn't that why the elemental tests are in TESTEXAMPLES_ELEMENTAL? The
>> > other runs fine for me without elemental, e.g.,
>>
>>   Don't for me:
>
>
> Rounding problem? How different are they?
>
>>
>>
>> Barrys-MacBook-Pro:tests barrysmith$ make runex104
>> 0a1,25
>> > [0]PETSC ERROR: --------------------- Error Message
>> > ------------------------------------
>> > [0]PETSC ERROR: Object is in wrong state!
>> > [0]PETSC ERROR: C != D!
>> > [0]PETSC ERROR:
>> > ------------------------------------------------------------------------
>> > [0]PETSC ERROR: Petsc Development HG revision:
>> > 34c1b366b0f8fabca0e1f952a026363f6ea75443  HG Date: Fri Jan 11 14:02:35 2013
>> > -0600
>> > [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: ./ex104 on a arch-gnu named Barrys-MacBook-Pro.local by
>> > barrysmith Tue Jan 15 17:09:12 2013
>> > [0]PETSC ERROR: Libraries linked from
>> > /Users/barrysmith/Src/petsc-dev/arch-gnu/lib
>> > [0]PETSC ERROR: Configure run at Fri Jan 11 21:09:18 2013
>> > [0]PETSC ERROR: Configure options --download-blacs --download-cmake
>> > --download-fftw --download-fiat --download-generator --download-hypre
>> > --download-libyaml --download-metis --download-ml --download-mpich
>> > --download-mumps --download-parmetis --download-ptscotch
>> > --download-scalapack --download-scientificpython --download-superlu
>> > --download-superlu_dist --download-triangle --download-yaml
>> > --with-afterimage --with-ams-dir=/Users/barrysmith/Src/ams-dev --with-java
>> > --with-openmp --with-pthreadclasses --with-server --with-shared-libraries
>> > PETSC_ARCH=arch-gnu --download-sowing
>> > [0]PETSC ERROR:
>> > ------------------------------------------------------------------------
>> > [0]PETSC ERROR: main() line 74 in src/mat/examples/tests/ex104.c
>> > application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0
>> > [cli_0]: aborting job:
>> > application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0
>> >
>> >
>> > ===================================================================================
>> > =   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
>> > =   EXIT CODE: 73
>> > =   CLEANING UP REMAINING PROCESSES
>> > =   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
>> >
>> > ===================================================================================
>> /Users/barrysmith/Src/petsc-dev/src/mat/examples/tests
>> Possible problem with ex104, diffs above
>>
>> >
>> > $ make PETSC_ARCH=mpich-clang ex104 runex104
>> > /home/jed/usr/mpich-clang/bin/mpicc -o ex104.o -c -fPIC -Wall
>> > -Wwrite-strings -Wno-strict-aliasing -Wno-unknown-pragmas -g3 -fno-inline
>> > -O0   -I/home/jed/petsc/include -I/home/jed/petsc/mpich-clang/include
>> > -I/usr/include -I/usr/cbind/include -I/usr/forbind/include
>> > -I/home/jed/usr/mpich-clang/include    -D__INSDIR__=src/mat/examples/tests/
>> > ex104.c
>> > /home/jed/usr/mpich-clang/bin/mpicc -fPIC -Wall -Wwrite-strings
>> > -Wno-strict-aliasing -Wno-unknown-pragmas -g3 -fno-inline -O0  -o ex104
>> > ex104.o -Wl,-rpath,/home/jed/petsc/mpich-clang/lib
>> > -L/home/jed/petsc/mpich-clang/lib  -lpetsc
>> > -Wl,-rpath,/home/jed/petsc/mpich-clang/lib -lsundials_cvode
>> > -lsundials_nvecserial -lsundials_nvecparallel -llapack -lblas -ltriangle
>> > -lX11 -lpthread -lchaco -Wl,-rpath,/usr/lib -L/usr/lib -lexoIIv2for
>> > -lexoIIv2c -lnetcdf -lm -Wl,-rpath,/home/jed/usr/mpich-clang/lib
>> > -L/home/jed/usr/mpich-clang/lib
>> > -Wl,-rpath,/usr/lib/gcc/x86_64-unknown-linux-gnu/4.7.2
>> > -L/usr/lib/gcc/x86_64-unknown-linux-gnu/4.7.2 -lmpichcxx -lstdc++ -ldl
>> > -lpmpich -lmpich -lopa -lmpl -lrt -lpthread -lgcc_s -ldl
>> > /bin/rm -f ex104.o
>> > $
>> >
>> >
>> >   /* Test MatTransposeMatMult() */
>> >   ierr =
>> > PetscObjectTypeCompare((PetscObject)A,MATELEMENTAL,&iselemental);CHKERRQ(ierr);
>> >   if (!iselemental){
>> >     ierr =
>> > MatTransposeMatMult(A,A,MAT_INITIAL_MATRIX,fill,&D);CHKERRQ(ierr); /* D =
>> > A^T*A */
>> >     ierr = MatEqual(C,D,&equal);CHKERRQ(ierr);
>> >     if (!equal) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_WRONGSTATE,"C !=
>> > D");
>> >     ierr = MatDestroy(&D);CHKERRQ(ierr);
>> >
>> >
>>
>