On Fri, Dec 19, 2008 at 4:39 PM, Rafael Santos Coelho <span dir="ltr"><<a href="mailto:rafaelsantoscoelho@gmail.com">rafaelsantoscoelho@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hey, Barry,<br><br>here's what you asked:</blockquote><div><br>You are taking 25K + iterations of GMRES(30). This I think is probably pretty sensitive to the<br>exact nature of the FD approximation. Did you try varying the h parameter?<br>
<br> Matt<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br><div>--------------------------------------------------------------------------------------------------------<br>
JFNK<br>
KSP SOLVER: GMRES<br>
MESH SIZE : 512 x 512 unknows<br>
NUMBER OF PROCESSORS : 24</div><br>Linear solve converged due to CONVERGED_RTOL iterations 25042<br>Linear solve converged due to CONVERGED_RTOL iterations 33804<br>Linear solve converged due to CONVERGED_RTOL iterations 33047<br>
Linear solve converged due to CONVERGED_RTOL iterations 21219<br>SNES Object:<br> type: ls<br> line search variant: SNESLineSearchCubic<br> alpha=0.0001, maxstep=1e+08, steptol=1e-12<br> maximum iterations=100, maximum function evaluations=1000000<br>
tolerances: relative=1e-08, absolute=1e-50, solution=1e-08<br> total number of linear solver iterations=113112<br> total number of function evaluations=116893<br> KSP Object:<br> type: gmres<br> GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<br>
GMRES: happy breakdown tolerance 1e-30<br> maximum iterations=1000000, initial guess is zero<br> tolerances: relative=1e-05, absolute=1e-50, divergence=10000<br> left preconditioning<br> PC Object:<br> type: none<br>
linear system matrix = precond matrix:<br> Matrix Object:<br> type=mffd, rows=262144, cols=262144<br> matrix-free approximation:<br> err=1e-07 (relative error in function evaluation)<br> Using wp compute h routine<br>
Computes normA<br> Does not compute normU<br>Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE<br><br> Max Max/Min Avg Total <br>
Time (sec): 8.080e+02 1.00009 8.080e+02<br>Objects: 6.600e+01 1.00000 6.600e+01<br>Flops: 1.044e+11 1.01176 1.036e+11 2.485e+12<br>
Flops/sec: 1.292e+08 1.01184 1.282e+08 3.076e+09<br>Memory: 5.076e+06 1.01530 1.207e+08<br>MPI Messages: 4.676e+05 2.00000 3.702e+05 8.884e+06<br>
MPI Message Lengths: 4.002e+08 2.00939 8.623e+02 7.661e+09<br>MPI Reductions: 9.901e+03 1.00000<br><br>Object Type Creations Destructions Memory Descendants' Mem.<br><br>--- Event Stage 0: Main Stage<br>
<br>SNES 1 1 124 0<br>Krylov Solver 1 1 16880 0<br>Preconditioner 1 1 0 0<br>
Distributed array 1 1 46568 0<br>Index Set 6 6 135976 0<br>Vec 46 46 3901252 0<br>
Vec Scatter 3 3 0 0<br>IS L to G Mapping 1 1 45092 0<br>MatMFFD 1 1 0 0<br>
Matrix 4 4 1011036 0<br>Viewer 1 1 0 0<br><br>OptionTable: -ksp_converged_reason<br>OptionTable: -ksp_max_it 1000000<br>
OptionTable: -ksp_type gmres<br>OptionTable: -log_summary<br>OptionTable: -option_table<br>OptionTable: -par 4.0<br>OptionTable: -pc_type none<br>OptionTable: -snes_converged_reason<br>OptionTable: -snes_max_funcs 1000000<br>
OptionTable: -snes_max_it 100<br>OptionTable: -snes_mf<br>OptionTable: -snes_view<br>OptionTable: -xdiv 512<br>OptionTable: -ydiv 512<br><div>--------------------------------------------------------------------------------------------------------<br>
<br><div>--------------------------------------------------------------------------------------------------------</div>NK<br>
KSP SOLVER: GMRES<br>
MESH SIZE : 512 x 512 unknows<br>
NUMBER OF PROCESSORS : 24<br></div><br>Linear solve converged due to CONVERGED_RTOL iterations 25038<br>Linear solve converged due to CONVERGED_RTOL iterations 25995<br>Linear solve converged due to CONVERGED_RTOL iterations 26769<br>
SNES Object:<br> type: ls<br> line search variant: SNESLineSearchCubic<br> alpha=0.0001, maxstep=1e+08, steptol=1e-12<br> maximum iterations=100, maximum function evaluations=1000000<br> tolerances: relative=1e-08, absolute=1e-50, solution=1e-08<br>
total number of linear solver iterations=77802<br> total number of function evaluations=4<br> KSP Object:<br> type: gmres<br> GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<br>
GMRES: happy breakdown tolerance 1e-30<br> maximum iterations=1000000, initial guess is zero<br> tolerances: relative=1e-05, absolute=1e-50, divergence=10000<br> left preconditioning<br> PC Object:<br> type: none<br>
linear system matrix = precond matrix:<br> Matrix Object:<br> type=mpiaij, rows=262144, cols=262144<br> total: nonzeros=1308672, allocated nonzeros=1308672<br> not using I-node (on process 0) routines<br>
Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE<br><br> Max Max/Min Avg Total <br>Time (sec): 5.329e+02 1.00008 5.329e+02<br>
Objects: 6.300e+01 1.00000 6.300e+01<br>Flops: 6.553e+10 1.01230 6.501e+10 1.560e+12<br>Flops/sec: 1.230e+08 1.01230 1.220e+08 2.928e+09<br>
Memory: 4.989e+06 1.01589 1.186e+08<br>MPI Messages: 3.216e+05 1.99999 2.546e+05 6.111e+06<br>MPI Message Lengths: 2.753e+08 2.00939 8.623e+02 5.269e+09<br>
MPI Reductions: 6.704e+03 1.00000<br><br>Object Type Creations Destructions Memory Descendants' Mem.<br><br>--- Event Stage 0: Main Stage<br><br>SNES 1 1 124 0<br>
Krylov Solver 1 1 16880 0<br>Preconditioner 1 1 0 0<br>Distributed array 1 1 46568 0<br>
Index Set 6 6 135976 0<br>Vec 45 45 3812684 0<br>Vec Scatter 3 3 0 0<br>
IS L to G Mapping 1 1 45092 0<br>Matrix 3 3 1011036 0<br>Viewer 1 1 0 0<br>
<br>OptionTable: -ksp_converged_reason<br>OptionTable: -ksp_max_it 1000000<br>OptionTable: -ksp_type gmres<br>OptionTable: -log_summary<br>OptionTable: -option_table<br>OptionTable: -par 4.0<br>OptionTable: -pc_type none<br>
OptionTable: -snes_converged_reason<br>OptionTable: -snes_max_funcs 1000000<br>OptionTable: -snes_max_it 100<br>OptionTable: -snes_view<br>OptionTable: -xdiv 512<br>OptionTable: -ydiv 512<br><div>
--------------------------------------------------------------------------------------------------------</div><br>
</blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>