[Nek5000-users] Interpolation between two meshes

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

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 <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users-request at lists.mcs.anl.gov <nek5000-users-request at lists.mcs.anl.gov>
Sent: 28 March 2017 22:15
To: 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)
   2. Re: Interpolation between two meshes
      (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
To: nek5000-users at lists.mcs.anl.gov      <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>
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 <nek5000-users at lists.mcs.anl.gov>
> Sent: Tuesday 28th March 2017 22:59
> To: 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

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