[petsc-users] Advice Being Sought

[David Scott]

On 03/05/2013 15:16, David Scott wrote:

> I have inserted
>      call PetscViewerBinaryOpen(PETSC_COMM_WORLD, 'DifficultSystem', &
>        FILE_MODE_WRITE, viewer, ierr)
>      call KSPGetOperators(ksp, A, PETSC_NULL_OBJECT, &
>        PETSC_NULL_INTEGER, ierr)
>      call MatView(A, viewer, ierr)
>      call KSPGetRhs(ksp, b, ierr)
>      call VecView(b, viewer, ierr)
>      call PetscViewerDestroy(viewer, ierr)
> in my code after KSPSolve.
>
> You may download the result for the 29th time step from
>    http://www2.epcc.ed.ac.uk/~dscott6/
>
> I can generate data for other time steps if necessary.
>
> David

I forgot to say that the program ran on 128 MPI processes with the 
following arguments
   -pc_type gamg -pc_gamg_agg_nsmooths 1 -ksp_rtol 0.0001

David
-- 
Dr. D. M. Scott
Applications Consultant
Edinburgh Parallel Computing Centre
Tel. 0131 650 5921

The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.