[petsc-dev] Fwd: [petsc-users] PCView only for one block ?

[Barry Smith]

   Someone has been mucking with PCView_BJacobi? Under normal circumstances where the same PC type is used on each block it only prints a single block. Not sure why it suddenly prints something for each block?


    if (jac->same_local_solves) {
      ierr = PetscViewerASCIIPrintf(viewer,"  Local solve is same for all blocks, in the following KSP and PC objects:\n");CHKERRQ(ierr);
      if (jac->ksp && !jac->psubcomm) {
        ierr = PetscViewerGetSingleton(viewer,&sviewer);CHKERRQ(ierr);
        if (!rank) {
          ierr = PetscViewerASCIIPushTab(viewer);CHKERRQ(ierr);
          ierr = KSPView(jac->ksp[0],sviewer);CHKERRQ(ierr);
          ierr = PetscViewerASCIIPopTab(viewer);CHKERRQ(ierr);
        }
        ierr = PetscViewerRestoreSingleton(viewer,&sviewer);CHKERRQ(ierr);
      } else if (jac->psubcomm && !jac->psubcomm->color) {
        ierr = PetscViewerASCIIGetStdout(mpjac->psubcomm->comm,&sviewer);CHKERRQ(ierr);
        ierr = PetscViewerASCIIPushTab(viewer);CHKERRQ(ierr);
        ierr = KSPView(*(jac->ksp),sviewer);CHKERRQ(ierr);
        ierr = PetscViewerASCIIPopTab(viewer);CHKERRQ(ierr);
      }
    } else {



Begin forwarded message:

> From: Jed Brown <jedbrown at mcs.anl.gov>
> Subject: Re: [petsc-users] PCView only for one block ?
> Date: April 16, 2013 9:36:53 AM CDT
> To: Matthew Knepley <knepley at gmail.com>, PETSc users list <petsc-users at mcs.anl.gov>
> Reply-To: PETSc users list <petsc-users at mcs.anl.gov>
> 
> Matthew Knepley <knepley at gmail.com> writes:
> 
>> Pull out the subpc you want and call PCView. There is nothing from the
>> command line.
> 
> I would be in favor of changing the default command line behavior.  The
> current behavior makes -ksp_view output hopeless for large runs.  I
> think -<solver>_view should always produce scalable output by default.