<div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Dec 12, 2021 at 3:19 PM Karthikeyan Chockalingam - STFC UKRI <<a href="mailto:karthikeyan.chockalingam@stfc.ac.uk" target="_blank">karthikeyan.chockalingam@stfc.ac.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">





<div lang="EN-GB">
<div>
<p class="MsoNormal"><span>Thank for your response that was helpful. I have a couple of questions:<u></u><u></u></span></p>
<p class="MsoNormal"><span><u></u> <u></u></span></p>
<p style="margin-left:54pt">
<u></u><span><span>(i)<span style="font:7pt "Times New Roman"">                 
</span></span></span><u></u><span>How can I control the level of refinement? I tried to pass the flag “-ex56_dm_refine 0” but that didn’t stop the refinement from
</span>8 giving 32 cubes.</p></div></div></blockquote><div><br></div><div>I answered this question recently but ex56 clobbers ex56_dm_refine in the convergence loop. I have an MR that prints a warning if you provide a ex56_dm_refine. </div><div><br></div><div>* snes/ex56 runs a convergence study and confusingly sets the options manually, thus erasing your -ex56_dm_refine.<div><br></div><div>* To refine, use -max_conv_its N <3>, this sets the number of steps of refinement. That is, the length of the convergence study</div><div><br></div><div>* You can adjust where it starts from with -cells i,j,k <1,1,1></div><div><div>You do want to set this if you have multiple MPI processes so that the size of this mesh is the number of processes. That way it starts with one cell per process and refines from there.<br></div><div><br></div></div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div lang="EN-GB"><div><p style="margin-left:54pt"><span><u></u><u></u></span></p>
<p style="margin-left:54pt">
<u></u><span><span>(ii)<span style="font:7pt "Times New Roman"">               
</span></span></span><u></u>What does -cell 2,2,1 correspond to? </p></div></div></blockquote><div><br></div><div>The initial mesh or mesh_0. The convergence test uniformly refines this mesh. So if you want to refine this twice you could use -cells 8,8,4</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div lang="EN-GB"><div><p style="margin-left:54pt">How can I determine the total number of dofs? </p></div></div></blockquote><div>Unfortunately, that is not printed but you can calculate from the initial cell grid, the order of the element and the refinement in each iteration of the convergence tests.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div lang="EN-GB"><div><p style="margin-left:54pt">So that I can perform a scaling study by changing the input of the flag -cells.<span><u></u><u></u></span></p>
<p class="MsoNormal"><span><br></span></p></div></div></blockquote><div><br></div><div>You can and the convergence test gives you data for a strong speedup study in one run. Each solve is put in its own "stage" of the output and you want to look at KSPSolve lines in the log_view output. <br></div></div></div>