<div class="gmail_quote">On Fri, Jan 20, 2012 at 16:03, Nystrom, William D <span dir="ltr"><<a href="mailto:wdn@lanl.gov">wdn@lanl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":4ye">I'm getting ready to use petsc-dev as part of our test suite for doing<br>
stability testing of some new hardware platforms with gpus.  I'm<br>
interested in testing both with and without using gpus.  So far, I have<br>
not extensively explored the various example problems that come with<br>
PETSc.  Mainly, I have run src/ksp/ksp/examples/tutorials/ex2.c<br></div></blockquote><div><br></div><div>This is a horrible example because the parallel decomposition is 1D. You can run essentially the same model with a better decomposition (and inside a nonlinear problem) with src/snes/examples/tutorials/ex5.c (which can also do geometric multigrid, but you don't have to use it).</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":4ye">
which is a 2d Laplace's equation problem.  I'm looking for recommendations<br>
of PETSc example problems that I should try, especially ones that might<br>
be more challenging numerically than ex2.c or more realistic of real<br>
world problems.  Does anyone reading this email have a favorite PETSc<br>
example problem that you would suggest I consider and try?</div></blockquote></div><div><br></div><div>Here are some options and profiling output (on a BG/P) for a 3D ice flow problem (src/snes/examples/tutorials/ex48.c):</div>
<br><div><a href="https://github.com/jedbrown/tme-ice/tree/master/shaheen/b">https://github.com/jedbrown/tme-ice/tree/master/shaheen/b</a></div><div><br></div><div><br></div><div>src/ksp/ksp/examples/tutorials/ex43.c (2D Stokes, preconditioned using fieldsplit inside multigrid)</div>
<div><br></div><div>mpiexec -n 8 ./ex43 -stokes_ksp_type gcr -stokes_ksp_gcr_restart 60 -stokes_ksp_norm_type unpreconditioned -stokes_ksp_rtol 1e-8 -c_str 3 -sinker_eta0 1.0 -sinker_eta1 100 -sinker_dx 0.4 -sinker_dy 0.3 -mx 1024 -my 1024 -stokes_ksp_monitor -o sc1024.dat -stokes_pc_type mg -stokes_mg_levels_pc_type fieldsplit -stokes_pc_mg_galerkin -stokes_mg_levels_pc_fieldsplit_block_size 3 -stokes_mg_levels_pc_fieldsplit_0_fields 0,1 -stokes_mg_levels_pc_fieldsplit_1_fields 2 -stokes_mg_levels_fieldsplit_0_pc_type sor -stokes_mg_levels_ksp_type chebychev -stokes_mg_levels_ksp_max_it 1 -stokes_mg_levels_ksp_chebychev_estimate_eigenvalues 0,0.2,0,1.1 -stokes_pc_mg_levels 5 -log_summary</div>
<div><br></div><div><br></div><div>src/ksp/ksp/examples/tutorials/ex56.c (3D elasticity, process count should be 8^n)</div><div>mpiexec -n 64 ./ex56 -ne 255 -alpha 1.e-3 -ksp_monitor_short -ksp_type cg -pc_gamg_type sa  -ksp_norm_type unpreconditioned -pc_gamg_verbose -pc_gamg_repartition -pc_gamg_process_eq_limit 10 -pc_gamg_coarse_eq_limit 20 -pc_gamg_mat_partitioning_type ptscotch -log_summary -pc_gamg_mat_partitioning_ptscotch_strategy scalability</div>