[petsc-users] MatRARt

[Hong Zhang]

John :

> Should MatRARt work for SeqAIJ matrices?  I don't understand what is
> wrong with this test code:
>

The code runs well on Mac and Linux using petsc-3.3 or petsc-dev. I cannot
repeat the error.
However, as you noticed, MatRARt calls coloring, which is experimental,
not well-tested, and the code needs to be cleaned :-(

Can you test MatRARt on
petsc-3.3/src/mat/examples/tests/ex94 with attached matrix data file by

./ex94 -f0 medium -f1 medium

Do it crash?

>From flops analysis and experiments, MatRARt() is not as efficient as
MatPtAP,
likely costs twice as MatPtAP.

Hong





> #include <petsc.h>
>
> #undef __FUNCT__
> #define __FUNCT__ "main"
> int main(int argc,char **args)
> {
>   Mat            A,R,C;
>   PetscInt       ai[5] = {0,2,4,6,8};
>   PetscInt       aj[8] = {0,3,1,2,1,2,0,3};
>   PetscReal      av[8] = {1,2,3,4,5,6,7,8};
>   PetscInt       ri[3] = {0,1,2};
>   PetscInt       rj[2] = {1,2};
>   PetscReal      rv[2] = {1,1};
>   PetscErrorCode ierr;
>
>   PetscInitialize(&argc,&args,(char *)0,(char *)0);
>
>   ierr =
> MatCreateSeqAIJWithArrays(PETSC_COMM_WORLD,4,4,ai,aj,av,&A);CHKERRQ(ierr);
>   ierr =
> MatCreateSeqAIJWithArrays(PETSC_COMM_WORLD,2,4,ri,rj,rv,&R);CHKERRQ(ierr);
>
>   //this is ok:
>   //ierr = MatMatMult(R,A,MAT_INITIAL_MATRIX,1.0,&C);CHKERRQ(ierr);
>
>   //this is not:
>   ierr = MatRARt(A,R,MAT_INITIAL_MATRIX,1.0,&C);CHKERRQ(ierr);
>
>   ierr = PetscFinalize();
>   return 0;
> }
>
>
> When I run it, it says:
>
> [0]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type!
> [0]PETSC ERROR: MatGetColumnIJ() not supported for matrix type seqaij!
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 2, Fri Jul 13
> 15:42:00 CDT 2012
> [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: ./rart_simple on a linux-deb named lagrange.tomato by
> jfe Wed Jul 18 13:07:33 2012
> [0]PETSC ERROR: Libraries linked from
> /home/jfe/local/petsc-3.3-p2/linux-debug/lib
> [0]PETSC ERROR: Configure run at Wed Jul 18 13:00:59 2012
> [0]PETSC ERROR: Configure options --with-x=0
> --download-f-blas-lapack=1 --with-mpi=1 --with-mpi-shared=1
> --with-mpi=1 --download-mpich=1 --with-debugging=1
> --with-gnu-compilers=yes --with-shared-libraries=1 --with-c++-support
> --with-clanguage=C
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: MatTransposeColoringCreate_SeqAIJ() line 1294 in
> /home/jfe/local/petsc-3.3-p2/src/mat/impls/aij/seq/matmatmult.c
> [0]PETSC ERROR: MatTransposeColoringCreate() line 9318 in
> /home/jfe/local/petsc-3.3-p2/src/mat/interface/matrix.c
> [0]PETSC ERROR: MatRARtSymbolic_SeqAIJ_SeqAIJ() line 342 in
> /home/jfe/local/petsc-3.3-p2/src/mat/impls/aij/seq/matrart.c
> [0]PETSC ERROR: MatRARt_SeqAIJ_SeqAIJ() line 541 in
> /home/jfe/local/petsc-3.3-p2/src/mat/impls/aij/seq/matrart.c
> [0]PETSC ERROR: MatRARt() line 8405 in
> /home/jfe/local/petsc-3.3-p2/src/mat/interface/matrix.c
> [0]PETSC ERROR: main() line 25 in
> src/ksp/ksp/examples/tutorials/rart_simple.c
>
> What am I doing wrong?
>
> Thanks,
> John
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120718/1b6177b2/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: medium
Type: application/octet-stream
Size: 29136 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120718/1b6177b2/attachment-0001.obj>