On Jan 15, 2013, at 4:35 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
> It works for me.
Cause you are running an install with Elemental, without elemental install
it goes bad.
/* 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);
}
>
> $ make ex104 runex104 runex104_elemental runex104_elemental_2
> /opt/mpich/bin/mpicc -o ex104.o -c -fPIC -Wall -Wwrite-strings
> -Wno-strict-aliasing -Wno-unknown-pragmas -g3 -fno-inline -O0 -fopenmp
> -I/home/jed/petsc/include -I/home/jed/petsc/mpich/include
> -I/home/jed/usr/clique-mpich/include -I/usr/include -I/usr/cbind/include
> -I/usr/forbind/include -I/opt/mpich/include
> -D__INSDIR__=src/mat/examples/tests/ ex104.c
> /opt/mpich/bin/mpicc -fopenmp -fPIC -Wall -Wwrite-strings
> -Wno-strict-aliasing -Wno-unknown-pragmas -g3 -fno-inline -O0 -fopenmp -o
> ex104 ex104.o -Wl,-rpath,/home/jed/petsc/mpich/lib
> -L/home/jed/petsc/mpich/lib -lpetscmat -lpetscvec -lpetscsys
> -Wl,-rpath,/home/jed/petsc/mpich/lib -L/home/jed/petsc/mpich/lib -lspooles
> -Wl,-rpath,/home/jed/usr/clique-mpich/lib -L/home/jed/usr/clique-mpich/lib
> -lclique -lparmetis-addons -lmetis-addons -lelemental -lplcg -lpmrrr
> -lsuperlu_4.3 -lml -Wl,-rpath,/opt/mpich/lib -L/opt/mpich/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++ -lspai
> -lsundials_cvode -lsundials_nvecserial -lsundials_nvecparallel -lHYPRE
> -lmpichcxx -lstdc++ -lsuperlu_dist_3.2 -Wl,-rpath,/usr/lib -L/usr/lib
> -lcholmod -lamd -lcolamd -lcamd -lccolamd -lumfpack -lcholmod -lcamd
> -lccolamd -lcolamd -lamd -lsuitesparseconfig -llapack -lblas -ltriangle -lX11
> -lpthread -lparmetis -lmetis -lchaco -lamspub -lexoIIv2for -lexoIIv2c
> -lnetcdf -lhdf5_fortran -lhdf5 -lhdf5hl_fortran -lhdf5_hl -lz -lm -lmpichf90
> -lgfortran -lm -lgfortran -lm -lquadmath -lm -lmpichcxx -lstdc++ -ldl -lmpich
> -lopa -lmpl -lrt -lpthread -lgcc_s -ldl
> /bin/rm -f ex104.o
> $
>
>
> On Tue, Jan 15, 2013 at 4:30 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
>
> This example has been broken for a while, did it ever work?
>
> =========================================
> 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 16:28:10 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
> >
>
>
