[petsc-dev] PCApplyTranspose in BiCG

Tue Apr 10 07:19:34 CDT 2012

On Tue, Apr 10, 2012 at 04:47, Alexander Grayver <agrayver at gfz-potsdam.de>wrote:

> Barry,
>
> Well, this is not easy to apply bisection here since "good" petsc-dev had
> bugs that I patched manually (including Intel MKL BLAS incompatibility) and
> only after that I was able to run solver successfully. "Bad" petsc-dev has
> all these issues patched and names of some routines have been changed as
> well (MatCreateXXX)... However TFQMR doesn't work at all for any system I
> tried.
>
> Unfortunately, the shortest time gap between "good" and "bad" I can give
> so far is:
>
> HG revision: a4dae06b4d05ef5b78a9a9aa2cbc31**b2b158834f HG Date: Wed Feb
> 08 11:29:52 2012 -0600 --- GOOD
> HG revision: 59ca71b85d784abbe22f12c34ffcb2**ccc6dc68b0 HG Date: Tue Mar
> 20 01:05:13 2012 -0500 --- BAD
>
> I tested them like that:
> ./ex10 -f A.dat -rhs RHS.dat -ksp_type tfqmr -pc_type sor
> -ksp_monitor_true_residual -ksp_max_it 500
>
>  95 KSP preconditioned resid norm 2.117351054162e-18 true resid norm
> 1.981564630346e-15 ||r(i)||/||b|| 3.893342156134e-05   ---- GOOD
>
> 500 KSP preconditioned resid norm 3.000274090666e-13 true resid norm
> 4.577554424627e-11 ||r(i)||/||b|| 8.993895702653e-01   ---- BAD
>

$ ./ex10 -f ~/dl/A.dat -rhs ~/dl/RHS.dat -ksp_type tfqmr -pc_type sor
-ksp_monitor_true_residual
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Arguments are incompatible!
[0]PETSC ERROR: Zero diagonal on row 44!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Development HG revision:
21fb71e3b67102b07ab64b3779f68e8719c17739  HG Date: Mon Mar 19 08:24:25 2012
-0500
[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: ./ex10 on a mpich named batura by jed Tue Apr 10 07:09:45
2012
[0]PETSC ERROR: Libraries linked from /home/jed/petsc/mpich/lib
[0]PETSC ERROR: Configure run at Mon Mar 19 09:50:34 2012
[0]PETSC ERROR: Configure options --download-ams --download-blacs
--download-chaco --download-hypre --download-metis --download-ml
--download-parmetis --download-spai --download-spooles --download-sundials
--download-superlu --download-superlu_dist --download-triangle
--with-blas-lapack=/usr --with-c2html --with-cholmod-dir=/usr --with-lgrind
--with-mpi-dir=/opt/mpich2 --with-pcbddc --with-pthreadclasses
--with-shared-libraries --with-single-library=0 --with-sowing
--with-umfpack-dir=/usr --with-x -PETSC_ARCH=mpich
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: MatInvertDiagonal_SeqAIJ() line 1315 in
/home/jed/petsc/src/mat/impls/aij/seq/aij.c
[0]PETSC ERROR: MatSOR_SeqAIJ() line 1348 in
/home/jed/petsc/src/mat/impls/aij/seq/aij.c
[0]PETSC ERROR: MatSOR() line 3671 in
/home/jed/petsc/src/mat/interface/matrix.c
[0]PETSC ERROR: PCApply_SOR() line 35 in
/home/jed/petsc/src/ksp/pc/impls/sor/sor.c
[0]PETSC ERROR: PCApply() line 384 in
/home/jed/petsc/src/ksp/pc/interface/precon.c
[0]PETSC ERROR: KSPInitialResidual() line 64 in
/home/jed/petsc/src/ksp/ksp/interface/itres.c
[0]PETSC ERROR: KSPSolve_TFQMR() line 41 in
/home/jed/petsc/src/ksp/ksp/impls/tfqmr/tfqmr.c
[0]PETSC ERROR: KSPSolve() line 434 in
/home/jed/petsc/src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: main() line 327 in src/ksp/ksp/examples/tutorials/ex10.c
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

There is no space allocated at (44,44); note that counting is 0-based. This
consistency check has been in PETSc for many years, so I don't know why you
didn't trip it in earlier runs. Do you have some options in PETSC_OPTIONS
or .petscrc?

>
> If it is of any help files are here:
> http://dl.dropbox.com/u/**60982984/A.dat<http://dl.dropbox.com/u/60982984/A.dat>
> http://dl.dropbox.com/u/**60982984/RHS.dat<http://dl.dropbox.com/u/60982984/RHS.dat>
>
> Is this possible to look up changes in one file only?

$ hg log src/ksp/ksp/impls/tfqmr/

But the problem was clearly not in that file.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120410/bc56175a/attachment.html>