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

Jed Brown jedbrown at mcs.anl.gov
Tue Apr 16 20:51:41 CDT 2013


Barry Smith <bsmith at mcs.anl.gov> writes:

>    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?

Mikolaj, can you show the code that you use for setting up this bjacobi?

>     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);

FWIW, src/snes/examples/tutorials/ex7.c also produces verbose output.  I
understand the desire to show details when someone configures each
subdomain differently, but I'd still rather have scalable output by
default.



More information about the petsc-dev mailing list