157ef82e2SBarry Smith% 256fec923SBarry Smith% Performance results for PETSc on the 2d Laplacian 5 pt stencil using GMRES(30) and no 356fec923SBarry Smith% preconditioning. Run to give an indication of scalability for 457ef82e2SBarry Smith% 556fec923SBarry Smith% * matrix vector product for very sparse matrix, and 656fec923SBarry Smith% * scalability of inner products and vector updates. 756fec923SBarry Smith% 856fec923SBarry Smith% NOTE: one would NEVER use these particular solvers in practice. 956fec923SBarry Smith% 10*21789920SBarry Smith 1156fec923SBarry Smith% For a 1000 by 1000 grid, with -ksp_max_it 500 -mat_aij_no_inode 1256fec923SBarry Smith% 1356fec923SBarry Smithp1000 = [ 2 4 8 16 32]; 1456fec923SBarry Smith% 1556fec923SBarry Smith% m and s contain the times for matrix vector multiply and SLES solve for 1656fec923SBarry Smith% IBM SP2, Cray T3E and Origin2000 in that order. 1756fec923SBarry Smith% 1856fec923SBarry Smithm1000 = [ 95.54 51.87 25.74 13.63 6.88; 134.6 67.59 33.93 16.77 9.03]; 1956fec923SBarry Smiths1000 = [ 455.5 245.2 123.9 72.41 38.63; 432.3 220 111.4 56.24 30.16]; 2073f4d377SMatthew Knepleyms1000 = [ 95.54/95.54 95.54/51.87 95.54/25.74 95.54/13.63 95.54/6.88; 134.6/134.6 134.6/67.59 134.6/33.93 134.6/16.77 134.6/9.03; 1 2 4 8 16]; 2156fec923SBarry Smith% 2273f4d377SMatthew Knepleyplotpt(p1000,ms1000,'*'); title('Speedup'); 23*21789920SBarry Smith 2456fec923SBarry Smith% For a 100 by 100 grid, with -mat_aij_no_inode 2556fec923SBarry Smith% 2656fec923SBarry Smithp100 = [ 1 2 3 4 5 6 7 8]; 2756fec923SBarry Smith% 2856fec923SBarry Smith% m and s contain the times for matrix vector multiply and SLES solve for 2956fec923SBarry Smith% IBM SP2, Cray T3E and Origin2000 in that order. 3056fec923SBarry Smith% 3173f4d377SMatthew Knepleym100 = [ 1.602 1.116 .793 .638 .5297 .4633 .4181 .5808; 2.233 1.434 1.064 .8116 .6748 .6095 .5338 .4873]; 3273f4d377SMatthew Knepleys100 = [ 9.425 6.036 4.32 3.334 3.634 3.656 3.657 4.248; 9.566 5.704 4.260 3.080 3.000 2.696 2.614 2.6061]; 3357ef82e2SBarry Smith% 34