[Nek5000-users] Interpolation between two meshes

[nek5000-users at lists.mcs.anl.gov]

Please let me know how I can contact you off list.

Stefan 

On 29 Mar 2017, at 14:59, "nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> " <nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> > wrote:

Hi Stefan,



Thank you very much for your quick reply.



1. Yes, my source field is in double precision since I need to keep the accuracy.

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 'CHECK interp findpts before' and it never gets updated again.

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.

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.



I guess it might be something to do with the memory. What do you think? Is there any solution to this?


Best regards,
Tony

        if(nid==0) print*,'CHECK interp findpts before'
        call nekgsync()
        call findpts(inth_gfldr,
     &               rcode,1,
     &               proc,1,
     &               elid,1,
     &               rst,ndim,
     &               dist,1,
     &               xm1,1,
     &               ym1,1,
     &               zm1,1,ntot)
        call nekgsync()
        if(nid==0) print*,'CHECK interp findpts after'

--------------------------------
From: nek5000-users-bounces at lists.mcs.anl.gov <mailto:nek5000-users-bounces at lists.mcs.anl.gov> <nek5000-users-bounces at lists.mcs.anl.gov <mailto:nek5000-users-bounces at lists.mcs.anl.gov> > on behalf of nek5000-users-request at lists.mcs.anl.gov <mailto:nek5000-users-request at lists.mcs.anl.gov> <nek5000-users-request at lists.mcs.anl.gov <mailto:nek5000-users-request at lists.mcs.anl.gov> >
Sent: 28 March 2017 22:15
To: nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> 
Subject: Nek5000-users Digest, Vol 97, Issue 35
 
Today's Topics:

   1. Re: Interpolation between two meshes
      (nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> )
   2. Re: Interpolation between two meshes
      (nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> )

Message: 2
Date: Tue, 28 Mar 2017 23:14:56 +0200
From: nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> 
To: nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov>       <nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> >
Subject: Re: [Nek5000-users] Interpolation between two meshes
Message-ID:
        <mailman.3928.1490735713.2967.nek5000-users at lists.mcs.anl.gov <mailto:mailman.3928.1490735713.2967.nek5000-users at lists.mcs.anl.gov> >
Content-Type: text/plain; charset=utf-8

Can you also call nekgsync() before and after findpts. Just trying to understand if it really hangs there.
 
-----Original message-----
> From:nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> <nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> >
> Sent: Tuesday 28th March 2017 22:59
> To: nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> 
> Subject: Re: [Nek5000-users] Interpolation between two meshes
> 
> We have tested glfdr() using different cases/machines and it worked fine. With the information below I cannot say what's going wrong.
> 
> - Is your source fld in double precision?
> - Do have any additional information where it hangs or crashes?
> - Can you please test gfldr() say with turbChannel (NekExamples).
> 
> Stefan


_______________________________________________

Nek5000-users mailing list

Nek5000-users at lists.mcs.anl.gov <mailto:Nek5000-users at lists.mcs.anl.gov> 

https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170330/35431d47/attachment.html>