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

[Barry Smith]

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:

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);
> 
>