If you have an explicit sparse matrix you are storing you can run with -ksp_view_binary and email to us the resulting file called binaryoutput at petsc-maint at mcs.anl.gov if the matrix is smaller than 20 megabytes. If it is larger you can anonymously ftp it to the ftp site called ftp.mcs.anl.gov and put it in the directory incoming then send email to petsc-maint at mcs.anl.gov and tell us where the file is.
Barry On Oct 27, 2011, at 9:45 AM, Klaij, Christiaan wrote: >> What happens at restart to the residual norm? That is does the >> residual norm change dramatically on the first iteration after >> the restart? This is an indication of loss of orthogonality. >> >> You can put in a restart of 1000 and run for 1000 >> iterations. What happens then? >> >> Barry > > 1000 is a bit much but I did 100 without restart and 100 with > default restart (every 30 its). I don't see significant change on > the first iteration after the restart. > > > dr. ir. Christiaan Klaij > CFD Researcher > Research & Development > E mailto:C.Klaij at marin.nl > T +31 317 49 33 44 > > MARIN > 2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands > T +31 317 49 39 11, F +31 317 49 32 45, I www.marin.nl >
