[petsc-users] Solve KSP in parallel.

[Manuel Valera]

Hello,

I'm working on solve a linear system in parallel, following ex12 of the ksp
tutorial i don't see major complication on doing so, so for a working
linear system solver with PCJACOBI and KSPGCR i did only the following
changes:

   call MatCreate(PETSC_COMM_WORLD,Ap,ierr)
!  call MatSetType(Ap,MATSEQAIJ,ierr)
  call MatSetType(Ap,MATMPIAIJ,ierr) !paralellization

  call MatSetSizes(Ap,PETSC_DECIDE,PETSC_DECIDE,nbdp,nbdp,ierr);

!  call MatSeqAIJSetPreallocationCSR(Ap,iapi,japi,app,ierr)
  call MatSetFromOptions(Ap,ierr)
!  call
MatCreateSeqAIJWithArrays(PETSC_COMM_WORLD,nbdp,nbdp,iapi,japi,app,Ap,ierr)
 call
MatCreateMPIAIJWithArrays(PETSC_COMM_WORLD,floor(real(nbdp)/sizel),PETSC_DECIDE,nbdp,nbdp,iapi,japi,app,Ap,ierr)


I grayed out the changes from sequential implementation.

So, it does not complain at runtime until it reaches KSPSolve(), with the
following error:


[1]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------
[1]PETSC ERROR: Object is in wrong state
[1]PETSC ERROR: Matrix is missing diagonal entry 0
[1]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for
trouble shooting.
[1]PETSC ERROR: Petsc Release Version 3.7.3, unknown
[1]PETSC ERROR: ./solvelinearmgPETSc


                                             � � on a arch-linux2-c-debug
named valera-HP-xw4600-Workstation by valera Mon Sep 26 13:35:15 2016
[1]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++
--with-fc=gfortran --download-fblaslapack=1 --download-mpich=1
--download-ml=1
[1]PETSC ERROR: #1 MatILUFactorSymbolic_SeqAIJ() line 1733 in
/home/valera/v5PETSc/petsc/petsc/src/mat/impls/aij/seq/aijfact.c
[1]PETSC ERROR: #2 MatILUFactorSymbolic() line 6579 in
/home/valera/v5PETSc/petsc/petsc/src/mat/interface/matrix.c
[1]PETSC ERROR: #3 PCSetUp_ILU() line 212 in
/home/valera/v5PETSc/petsc/petsc/src/ksp/pc/impls/factor/ilu/ilu.c
[1]PETSC ERROR: #4 PCSetUp() line 968 in
/home/valera/v5PETSc/petsc/petsc/src/ksp/pc/interface/precon.c
[1]PETSC ERROR: #5 KSPSetUp() line 390 in
/home/valera/v5PETSc/petsc/petsc/src/ksp/ksp/interface/itfunc.c
[1]PETSC ERROR: #6 PCSetUpOnBlocks_BJacobi_Singleblock() line 650 in
/home/valera/v5PETSc/petsc/petsc/src/ksp/pc/impls/bjacobi/bjacobi.c
[1]PETSC ERROR: #7 PCSetUpOnBlocks() line 1001 in
/home/valera/v5PETSc/petsc/petsc/src/ksp/pc/interface/precon.c
[1]PETSC ERROR: #8 KSPSetUpOnBlocks() line 220 in
/home/valera/v5PETSc/petsc/petsc/src/ksp/ksp/interface/itfunc.c
[1]PETSC ERROR: #9 KSPSolve() line 600 in
/home/valera/v5PETSc/petsc/petsc/src/ksp/ksp/interface/itfunc.c
At line 333 of file solvelinearmgPETSc.f90
Fortran runtime error: Array bound mismatch for dimension 1 of array 'sol'
(213120/106560)


This code works for -n 1 cores, but it gives this error when using more
than one core.

What am i missing?

Regards,

Manuel.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160926/123d61f7/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: solvelinearmgPETSc.f90
Type: text/x-fortran
Size: 14416 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160926/123d61f7/attachment.bin>