[petsc-users] Error with MatView on multiple processors

[Matthew Knepley]

On Fri, May 27, 2016 at 6:34 PM, Paul Urbanczyk <gomer at stanford.edu> wrote:

> Hello,
>
> I'm having trouble with the MatView function drawing the matrix
> structure(s) when I execute my code on multiple processors.
>
> When I run on a single processor, the code runs fine, and the graphics
> window displays cleanly.
>

Lets start with an example. I do this

  cd $PETSC_DIR/src/snes/examples/tutorials
  make ex19

Make sure it runs

  ./ex19 -snes_monitor

Make sure it runs in parallel (maybe you need
$PETSC_DIR/$PETSC_ARCH/bin/mpiexec)

  mpiexec -n 2 ./ex19 -snes_monitor

Make sure it can draw

  mpiexec -n 2 ./ex19 -snes_monitor -mat_view draw -draw_pause 1

This runs fine for me. Can you try it?

  Thanks,

     Matt

When I run with multiple processors, I get error messages (see below).
>
> The matrices are constructed with DMCreateMatrix(da, &A_matrix).
>
> I then set the values with
> MatSetValuesStencil(A_matrix,1,&row,2,col_A,value_A,INSERT_VALUES).
>
> Finally, I call MatAssemblyBegin(A_matrix,MAT_FINAL_ASSEMBLY) and
> MatAssemblyEnd(A_matrix,MAT_FINAL_ASSEMBLY).
>
> I also test that the matrices are assembled with MatAssembled(A_matrix,
> &is_assembled_bool), and it appears they are successfully assembled.
>
> Any help/advice is greatly appreciated.
>
> Thanks in advance!
>
> -Paul Urbanczyk
>
> [0]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
> [0]PETSC ERROR: Invalid argument
> [0]PETSC ERROR: Wrong type of object: Parameter # 1
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
> for trouble shooting.
> [0]PETSC ERROR: Petsc Release Version 3.7.1, unknown
> [0]PETSC ERROR: ./urbanSCFD on a arch-linux2-c-debug named prometheus by
> gomer Fri May 27 16:29:01 2016
> [0]PETSC ERROR: Configure options --with-cc=mpicc --with-cxx=mpicxx
> --with-fc=mpif90
> [0]PETSC ERROR: #1 AOApplicationToPetsc() line 267 in
> /home/gomer/local/petsc/src/vec/is/ao/interface/ao.c
> [0]PETSC ERROR: #2 MatView_MPI_DA() line 557 in
> /home/gomer/local/petsc/src/dm/impls/da/fdda.c
> [0]PETSC ERROR: #3 MatView() line 901 in
> /home/gomer/local/petsc/src/mat/interface/matrix.c
> [1]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
> [1]PETSC ERROR: Invalid argument
> [1]PETSC ERROR: Wrong type of object: Parameter # 1
> [1]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
> for trouble shooting.
> [1]PETSC ERROR: Petsc Release Version 3.7.1, unknown
> [1]PETSC ERROR: ./urbanSCFD on a arch-linux2-c-debug named prometheus by
> gomer Fri May 27 16:29:01 2016
> [1]PETSC ERROR: Configure options --with-cc=mpicc --with-cxx=mpicxx
> --with-fc=mpif90
> [1]PETSC ERROR: #1 AOApplicationToPetsc() line 267 in
> /home/gomer/local/petsc/src/vec/is/ao/interface/ao.c
> [1]PETSC ERROR: #2 MatView_MPI_DA() line 557 in
> /home/gomer/local/petsc/src/dm/impls/da/fdda.c
> [1]PETSC ERROR: #3 MatView() line 901 in
> /home/gomer/local/petsc/src/mat/interface/matrix.c
> [0]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
> [0]PETSC ERROR: Invalid argument
> [1]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
> [1]PETSC ERROR: Invalid argument
> [1]PETSC ERROR: Wrong type of object: Parameter # 1
> [1]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
> for trouble shooting.
> [1]PETSC ERROR: Petsc Release Version 3.7.1, unknown
> [1]PETSC ERROR: ./urbanSCFD on a arch-linux2-c-debug named prometheus by
> gomer Fri May 27 16:29:01 2016
> [1]PETSC ERROR: Configure options --with-cc=mpicc --with-cxx=mpicxx
> --with-fc=mpif90
> [1]PETSC ERROR: #4 AOApplicationToPetsc() line 267 in
> /home/gomer/local/petsc/src/vec/is/ao/interface/ao.c
> [1]PETSC ERROR: #5 MatView_MPI_DA() line 557 in
> /home/gomer/local/petsc/src/dm/impls/da/fdda.c
> [1]PETSC ERROR: #6 MatView() line 901 in
> /home/gomer/local/petsc/src/mat/interface/matrix.c
> [0]PETSC ERROR: Wrong type of object: Parameter # 1
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
> for trouble shooting.
> [0]PETSC ERROR: Petsc Release Version 3.7.1, unknown
> [0]PETSC ERROR: ./urbanSCFD on a arch-linux2-c-debug named prometheus by
> gomer Fri May 27 16:29:01 2016
> [0]PETSC ERROR: Configure options --with-cc=mpicc --with-cxx=mpicxx
> --with-fc=mpif90
> [0]PETSC ERROR: #4 AOApplicationToPetsc() line 267 in
> /home/gomer/local/petsc/src/vec/is/ao/interface/ao.c
> [0]PETSC ERROR: #5 MatView_MPI_DA() line 557 in
> /home/gomer/local/petsc/src/dm/impls/da/fdda.c
> [0]PETSC ERROR: #6 MatView() line 901 in
> /home/gomer/local/petsc/src/mat/interface/matrix.c
>
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160527/66a91b5c/attachment-0001.html>