<div dir="ltr">And is this just a simple cube problem in elasticity as it looks like from your input file name?</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Jan 30, 2014 at 4:40 AM, Thomas Gross <span dir="ltr"><<a href="mailto:tgross@ilsb.tuwien.ac.at" target="_blank">tgross@ilsb.tuwien.ac.at</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Please find enclosed the output for GAMG using  -mg_levels_ksp_max_it 1:<br>
<div class="im"><br>
Prometheus:<br>
-ksp_type cg -pc_type prometheus -log_summary -ksp_monitor -ksp_view -options_left > Prometheus_Large_Log_New.txt<br>
<br>
Gamg:<br>
</div>-ksp_type cg -pc_type gamg -pc_gamg_type agg -pc_gamg_agg_nsmooths 1 -log_summary -ksp_monitor -ksp_view -options_left -mg_levels_ksp_type richardson -mg_levels_pc_type sor -mg_levels_ksp_max_it 1 >GAMG_Rich_max_it.txt<br>

<br>
Best regards,<br>
Thomas<br>
<br>
<br><br>
<br>
On Jan 30, 2014, at 2:14 AM, Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br>
<br>
> Thomas Gross <<a href="mailto:tgross@ilsb.tuwien.ac.at">tgross@ilsb.tuwien.ac.at</a>> writes:<br>
><br>
>> I have now compiled Petsc 3.2.0 with Prometheus and Petsc 3.4.3 with GAMG using the same compiler. Furthermore I am using the same MPI implementation for both runs. Still Prometheus is considerably faster (see attached log files).<br>

>><br>
>> Prometheus:<br>
>> -ksp_type cg -pc_type prometheus -log_summary -ksp_monitor -ksp_view -options_left > Prometheus_Large_Log_New.txt<br>
>><br>
>> Gamg:<br>
>> -ksp_type cg -pc_type gamg -pc_gamg_type agg -pc_gamg_agg_nsmooths 1 -log_summary -ksp_monitor -ksp_view -options_left -mg_levels_ksp_max_it 1 > GAMG_Large_max_it_log_New.txt<br>
>> -ksp_type cg -pc_type gamg -pc_gamg_type agg -pc_gamg_agg_nsmooths 1 -log_summary -ksp_monitor -ksp_view -options_left mg_levels_ksp_type richardson -mg_levels_pc_type sor > GAMG_Large_Rich_Log_New.txt<br>
><br>
> Can you add -mg_levels_ksp_max_it 1 to this last configuration?  I see a<br>
> lot of time in MatSOR.  We can actually speed that up slightly with some<br>
> extra caching, which may be worthwhile.<br>
><br>
> We're seeing a bit lower performance in MatMult with GAMG, perhaps<br>
> because we are not using block formats specialized for elasticity.<br>
><br>
> Mark, what else is different?  What does Prometheus do differently in<br>
> setup (not the bottleneck here, but I'm curious).<br>
<br>
<br></blockquote></div><br></div>