<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto"><div></div><div>Please let me know how I can contact you off list.</div><div><br></div><div>Stefan </div><div><br>On 29 Mar 2017, at 14:59, "<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>" <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br><br></div><blockquote type="cite"><div>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Hi Stefan,</p>
<p><br>
</p>
<p>Thank you very much for your quick reply.</p>
<p><br>
</p>
<p>1. Yes, my source field is in double precision since I need to keep the accuracy.</p>
<p>2. I added several print lines in gfldr.f to check where the interpolation hangs (see below). I also put 'call nekgsync()' before and after findpts as you suggested. I can confirm that findpts is exactly where the interpolation hangs. The last output in
the logfile is '<span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px">CHECK interp findpts before</span>' and it never gets updated again.<br>
</p>
<p>3. I've just tested gfldr() with turbChannel example. Everything worked fine and the source field was interpolated onto the turbChannel mesh (512 elements). So I moved on increasing the number of elements of turbChannel case.</p>
<p>First, I used 234300 elements in the turbChannel example (half of my current simulation size) and the interpolation was done successfully. Then, I used 468600 elements for the turbChannel example (similar to my current simulation size), the interpolation
hangs again at findpts.</p>
<p><br>
</p>
<p>I guess it might be something to do with the memory. What do you think? Is there any solution to this?</p>
<br>
Best regards,<br>
Tony<br>
<br>
<div> if(nid==0) print*,'CHECK interp findpts before'<br>
call nekgsync()<br>
call findpts(inth_gfldr,<br>
& rcode,1,<br>
& proc,1,<br>
& elid,1,<br>
& rst,ndim,<br>
& dist,1,<br>
& xm1,1,<br>
& ym1,1,<br>
& zm1,1,ntot)<br>
call nekgsync()<br>
if(nid==0) print*,'CHECK interp findpts after'</div>
<br>
<div style="color:rgb(19,19,18)">
<div>
<hr style="display:inline-block; width:98%" tabindex="-1">
<div id="x_divRplyFwdMsg" dir="ltr"><font color="#000000" face="Calibri, sans-serif" style="font-size:11pt"><b>From:</b> <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">nek5000-users-bounces@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">nek5000-users-bounces@lists.mcs.anl.gov</a>> on behalf of <a href="mailto:nek5000-users-request@lists.mcs.anl.gov">nek5000-users-request@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users-request@lists.mcs.anl.gov">nek5000-users-request@lists.mcs.anl.gov</a>><br>
<b>Sent:</b> 28 March 2017 22:15<br>
<b>To:</b> <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
<b>Subject:</b> Nek5000-users Digest, Vol 97, Issue 35</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt"></span></font><br>
<font size="2"><span style="font-size:10pt">
<div class="PlainText">Today's Topics:<br>
<br>
1. Re: Interpolation between two meshes<br>
(<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>)<br>
2. Re: Interpolation between two meshes<br>
(<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>)<br>
<br>
Message: 2<br>
Date: Tue, 28 Mar 2017 23:14:56 +0200<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>><br>
Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
Message-ID:<br>
<<a href="mailto:mailman.3928.1490735713.2967.nek5000-users@lists.mcs.anl.gov">mailman.3928.1490735713.2967.nek5000-users@lists.mcs.anl.gov</a>><br>
Content-Type: text/plain; charset=utf-8<br>
<br>
Can you also call nekgsync() before and after findpts. Just trying to understand if it really hangs there.
<br>
<br>
-----Original message-----<br>
> From:<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>><br>
> Sent: Tuesday 28th March 2017 22:59<br>
> To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
> Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
> <br>
> We have tested glfdr() using different cases/machines and it worked fine. With the information below I cannot say what's going wrong.
<br>
> <br>
> - Is your source fld in double precision?<br>
> - Do have any additional information where it hangs or crashes?<br>
> - Can you please test gfldr() say with turbChannel (NekExamples).<br>
> <br>
> Stefan<br>
<br>
</div>
</span></font></div>
</div>
<pre>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</pre></div></blockquote></body></html>