<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><br>
</div>
<div>Thanks Sandeep.</div>
<div><br>
</div>
Your approach should work too.
<div><br>
</div>
<div>Can you elaborate more on the difference ?</div>
<div><br>
</div>
<div>It may be because there is a dsavg() in comp_vort3... but I think I turned that off.</div>
<div><br>
</div>
<div>Comp_vort3 might be a bit faster, but we normally don't care about performance for</div>
<div>this routine, so I think your way is cleaner and thus better.</div>
<div><br>
</div>
<div>I'm concerned however that they are different (which likely points to a problem in the</div>
<div>new comp_vort3 routine).</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
</div>
<div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF786475" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> nek5000-users-bounces@lists.mcs.anl.gov [nek5000-users-bounces@lists.mcs.anl.gov] on behalf of nek5000-users@lists.mcs.anl.gov [nek5000-users@lists.mcs.anl.gov]<br>
<b>Sent:</b> Monday, March 27, 2017 11:16 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] comp_vort3, nelt, conj_ht<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div class="gmail_default" style="font-size:small">Thank you Paul. Now comp_vort3 doesnt give zeros in the solid domain.</div>
<div class="gmail_default" style="font-size:small">What is the difference between comp_vort3 and the calculation of curl using </div>
<div class="gmail_default"> call gradm1(dAxdx,dAxdy,dAxdz,t(1,1,1,1,1)) </div>
<div class="gmail_default"> call gradm1(dAydx,dAydy,dAydz,t(1,1,1,1,2)) </div>
<div class="gmail_default"> call gradm1(dAzdx,dAzdy,dAzdz,t(1,1,1,1,3))</div>
<div class="gmail_default"> call sub3(bfx1,dAzdy,dAydz,m) </div>
<div class="gmail_default"> call sub3(bfy1,dAxdz,dAzdx,m) </div>
<div class="gmail_default"> call sub3(bfz1,dAydx,dAxdy,m) </div>
<div class="gmail_extra">
<div class="gmail_default" style="font-size:small">​When I compare the glsc3 of a curls computed by comp_vort3 and gradm1, there is a difference of 0.5% to 1%.</div>
<div class="gmail_default" style="font-size:small">Which procedure suits well, if I want to use u x curl in the userq</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">Sincerely,</div>
<div class="gmail_default" style="font-size:small">Sandeep</div>
<br>
</div>
<div class="gmail_extra">
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<br>
<br>
<br>
Hi Sandeep,<br>
<br>
This issue should now be resolved.<br>
<br>
comp_vort3() will now compute the curl of the supplied vector field over all elements, e=1,...,nelt.<br>
<br>
Moreover, prepost() has been modified to not put zeros in the velocity field when e > nelv.<br>
<br>
So, you should be able to do something like:<br>
<br>
call comp_vort3(curlv,w1,w2,u,v,w)<br>
<br>
call outpost(curlv(1,1),curlv(1,2),<wbr>culrv(1,3),pr,t,' ')<br>
<br>
and have the right thing happen.<br>
<br>
Paul<br>
______________________________<wbr>__<br>
From: <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.mc<wbr>s.anl.gov</a> [<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.m<wbr>cs.anl.gov</a>] on behalf
of <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a> [<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>]<br>
Sent: Thursday, March 23, 2017 2:57 PM<br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
Subject: [Nek5000-users] comp_vort3, nelt, conj_ht<br>
<br>
Hello Neks,<br>
I modified conj_ht example and defined three scalars. NEL>NELV. The domain of scalars is bigger than the fluid.<br>
I create a vector by taking the curl of three scalars, using<br>
common /myvort/ b1(lx1,ly1,lz1,lelt,3),<br>
call comp_vort3(b1,work1,work2,t(1,<wbr>1,1,1,1),<br>
$ t(1,1,1,1,2),t(1,1,1,1,3))<br>
<br>
b1 is computed correctly in the fluid domain, but is solid domain it showing zero.<br>
Can someone tell me how to modify the comp_vort3 subroutine, to correctly compute b1 over the entire domain.<br>
I tried to modify call dudxyz(work1,w,rym1,sym1,tym1,<wbr>jacm1,2,2), changing the imsh=2 (the last second parameter in the above call), but the the result remains the same.<br>
I use the three components of computed b1 in userq<br>
Thank you in advance.<br>
<br>
Sincerely,<br>
Sandeep<br>
<br>
<br>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>