<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<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> nek5000-users-bounces@lists.mcs.anl.gov <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users-request@lists.mcs.anl.gov <nek5000-users-request@lists.mcs.anl.gov><br>
<b>Sent:</b> 28 March 2017 22:15<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<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>
(nek5000-users@lists.mcs.anl.gov)<br>
2. Re: Interpolation between two meshes<br>
(nek5000-users@lists.mcs.anl.gov)<br>
<br>
Message: 2<br>
Date: Tue, 28 Mar 2017 23:14:56 +0200<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
Message-ID:<br>
<mailman.3928.1490735713.2967.nek5000-users@lists.mcs.anl.gov><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:nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
> Sent: Tuesday 28th March 2017 22:59<br>
> To: nek5000-users@lists.mcs.anl.gov<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>
</body>
</html>