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.
