[Nek5000-users] findpts returns ix iy iz iel

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

Dear Jian,


rst will not generally map to the discrete point set (i,j,k) within an element.


One strategy that would be reasonably fast would be to bin the intervals

on either side of each Gauss Lobatto Legendre (GLL) point on the 1D interval [-1,1]

and then to search for the GLL point that is closest to the r value returned from

findpts.  Then repeat for s and t.


In fact, you could generate a mapping function u(r) that would map the GLL points to a

uniform grid [1,lx1] then you could interrogate u(r) and and round the output via:


i = u(r) + 0.5


Then repeat for s-->j, and t--> k.


hth


Paul


________________________________
From: Nek5000-users <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
Sent: Saturday, July 28, 2018 6:02:39 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] findpts returns ix iy iz iel

Hi Nek,

What I want is the ix, iy, iz and iel at some points pts. So the findpts would be called with rst returned. But as far as I can see, rst are the coordinates of such element. May I has the grid point in that element, i.e. ix [1,nx1], iy [1,ny1], iz [1,nz1]?

Kind regards,

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