<div dir="ltr"><div><div><div>Dear Paul,<br><br></div>Thank you for the answer. I will use gradm1 then.<br><br></div>Best,<br></div>Nicolas<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Oct 3, 2017 at 11:25 PM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div>

<div id="m_-1009409361474000332divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p><br>
</p>
<p>Hi Nicolas,</p>
<p><br>
</p>
<p>You can use gradm1 to get your strong divergence.</p>
<p><br>
</p>
<p>Laplacian is a bit trickier since the 2nd derivative really is well-defined.</p>
<p>Some people compute the gradient and then apply dsavg() to each component</p>
<p>prior to taking the divergence of that field.</p>
<p><br>
</p>
<p>Right now there are no tools to do exactly the thing you want, so using gradm1()</p>
<p>is as good as any option.</p>
<p><br>
</p>
<p>Note that gradm1() is as efficient (or more so) than some of the other tools</p>
<p>to compute gradient in Nek.  In order to get (say) du/dx, you have to compute</p>
<p>du/dr , du/ds , du/dt, and then get du/dx from the chain rule.   du/dy, du/dz</p>
<p>also come from the chain rule, almost for free, even though you don't need them.</p>
<p>So you can simplify everything just through repeated use of gradm1.</p>
<p><br>
</p>
<p>Best, Paul</p>
<p><br>
</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_-1009409361474000332divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Nek5000-users <<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.<wbr>mcs.anl.gov</a>> on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<wbr>gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<wbr>gov</a>><br>
<b>Sent:</b> Tuesday, October 3, 2017 12:12:59 PM<br>
<b>To:</b> <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
<b>Subject:</b> [Nek5000-users] Strong divergence and Laplacian operator</font>
<div> </div>
</div><div><div class="h5">
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>Hi all,<br>
<br>
</div>
I need to compute some operators in the strong form directly for a project, i.e. without the mass matrix included in the operator.<br>
<br>
</div>
For the convection and gradient operators, I am using the functions convop and gradm1 respectively. However, I would also need to compute strong divergence and Laplacian and I did not find the corresponding functions (opdiv and axhelm are computing the weak
 form if I understand correctly). Do such functions exist? Or what modifications should I bring to the existing routines to change them to strong form?<br>
<br>
</div>
Best regards,<br>
</div>
Nicolas<br>
</div>
</div>
</div></div></div>

<br>______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>