[petsc-users] PCView only for one block ?

[Mikolaj]

Thank you Matt.

I have done what you suggested i.e.,

call PCBJacobiGetSubKSP(pc,nlocal,PETSC_NULL,subksp,ierr)
call KSPView(subksp(1),PETSC_VIEWER_STDOUT_WORLD,ierr)

and it works perfectly.

Mikolaj.

> On Tue, Apr 16, 2013 at 4:05 AM, Mikolaj <mikolaj.szydlarski at cea.fr 
> <mailto:mikolaj.szydlarski at cea.fr>> wrote:
>
>     Hello,
>
>     I am using a block preconditioners and I wonder if it is possible
>     to print PC data structure for sub-solvers limited only to the one
>     block. By default PCView 'pollute' std output with the information
>     about every single sub-setup from each MPI process ( e.g. below )
>     and I haven't found so far a solution for printing it only from
>     one processor.
>
>
> Pull out the subpc you want and call PCView. There is nothing from the 
> command line.
>
>    Matt
>
>     call PCView(pc,PETSC_VIEWER_STDOUT_WORLD,ierr) after setup of
>     sub-solvers:
>
>     PC Object: 16 MPI processes
>       type: bjacobi
>         block Jacobi: number of blocks = 16
>         Local solve info for each block is in the following KSP and PC
>     objects:
>       [0] number of local blocks = 1, first local block number = 0
>         [0] local block number 0
>
>     And then 16 x
>
>     {
>
>         KSP Object:    (sub_)     1 MPI processes
>           type: preonly
>           maximum iterations=10000, initial guess is zero
>           tolerances:  relative=1e-05, absolute=1e-50, divergence=10000
>           left preconditioning
>           using NONE norm type for convergence test
>         PC Object:    (sub_)     1 MPI processes
>           type: ilu
>             ILU: out-of-place factorization
>             1 level of fill
>             tolerance for zero pivot 2.22045e-14
>
>     }
>
>     Thank you in advance for any hint,
>
>     best regards,
>
>     Mikolaj
>
>     p.s. My petsc version: petsc-3.3-p5 and piece of code for printing:
>
>
>
>
>
>
>
>
> -- 
> 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/20130416/760c570b/attachment.html>