[petsc-users] About MatTransposeMatMult

Thu Sep 26 13:41:14 CDT 2013

Hello Hong,

I am attaching the whole error messages. I am using Petsc3.4.1. Thank you.

Joon

[0]PETSC ERROR: --------------------- Error Message ------------------------------------
[0]PETSC ERROR: Arguments are incompatible!
[0]PETSC ERROR: MatTransposeMatMult requires A, seqaij, to be compatible with B, seqdense!
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.4.1, Jun, 10, 2013
[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: ./tenfac on a linux named rossmann-b000.rcac.purdue.edu by choi240 Thu Sep 26 14:39:01 2013
[0]PETSC ERROR: Libraries linked from /apps/rhel5/petsc-3.4.1/64/impi-4.1.0.030_intel-13.0.1.117_ind64/linux/lib
[0]PETSC ERROR: Configure run at Sun Jun 23 10:20:51 2013
[0]PETSC ERROR: Configure options --with-cc=mpiicc --with-cxx=mpiicpc --with-fc=mpiifort --download-sowing --with-scalar-type=real --with-shared-libraries=1 --with-pic=1 --with-clanguage=C++ --with-fortran --with-fortran-kernels=1 --with-64-bit-indices=1 --with-debugging=0 --with-blas-lapack-dir=/opt/intel/composer_xe_2013.1.117/mkl/lib/intel64 --COPTFLAGS=-O3 --CXXOPTFLAGS=-O3 --FOPTFLAGS=-O3 --download-hdf5=no --download-metis=no --download-parmetis=no --download-superlu_dist=no --download-mumps=no --download-scalapack=yes --download-blacs=no --download-hypre=no --download-spooles=no
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: MatTransposeMatMult() line 8850 in /apps/rhel5/petsc-3.4.1/64/impi-4.1.0.030_intel-13.0.1.117_ind64/src/mat/interface/matrix.c
[0]PETSC ERROR: cal_M1() line 36 in tenfac.cpp
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind[0]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors
[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run
[0]PETSC ERROR: to get more information on the crash.
[0]PETSC ERROR: --------------------- Error Message ------------------------------------
[0]PETSC ERROR: Signal received!
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.4.1, Jun, 10, 2013
[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: ./tenfac on a linux named rossmann-b000.rcac.purdue.edu by choi240 Thu Sep 26 14:39:01 2013
[0]PETSC ERROR: Libraries linked from /apps/rhel5/petsc-3.4.1/64/impi-4.1.0.030_intel-13.0.1.117_ind64/linux/lib
[0]PETSC ERROR: Configure run at Sun Jun 23 10:20:51 2013
[0]PETSC ERROR: Configure options --with-cc=mpiicc --with-cxx=mpiicpc --with-fc=mpiifort --download-sowing --with-scalar-type=real --with-shared-libraries=1 --with-pic=1 --with-clanguage=C++ --with-fortran --with-fortran-kernels=1 --with-64-bit-indices=1 --with-debugging=0 --with-blas-lapack-dir=/opt/intel/composer_xe_2013.1.117/mkl/lib/intel64 --COPTFLAGS=-O3 --CXXOPTFLAGS=-O3 --FOPTFLAGS=-O3 --download-hdf5=no --download-metis=no --download-parmetis=no --download-superlu_dist=no --download-mumps=no --download-scalapack=yes --download-blacs=no --download-hypre=no --download-spooles=no
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: User provided function() line 0 in unknown directory unknown file
application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0

----- Original Message -----
From: "Hong Zhang" <hzhang at mcs.anl.gov>
To: "Karl Rupp" <rupp at mcs.anl.gov>
Cc: "Joon Hee Choi" <choi240 at purdue.edu>, "PETSc users list" <petsc-users at mcs.anl.gov>
Sent: Thursday, September 26, 2013 10:50:54 AM
Subject: Re: [petsc-users] About MatTransposeMatMult

Joon, 
Please send ENTIRE error message when reporting problems, from which we would know which petsc version and specific line in the petsc library that produced error. 

Hong 

On Thu, Sep 26, 2013 at 4:55 AM, Karl Rupp < rupp at mcs.anl.gov > wrote: 

Hi, 

I am trying to calculate M = X^T * B. X is a seqaij and B is a seqdense matrix. 

This is not yet allowed in the latest 3.4. release. 

However, I used MatTransposeMatMult() but I got an error. Can someone why I got the error and how I should modify my code? I am attaching my simple code, error message, and result. Thank you very much. 

Do you use the master branch from our Bitbucket repository? It contains the commit where Hong added the respective functionality: 
https://bitbucket.org/petsc/ petsc/commits/ 8949adfd119cb1d44c9fb5cf3dc01f 4b6a02252d 

Best regards, 
Karli 