[Nek5000-users] r,s,t index for given point
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Fri Jul 9 11:25:41 CDT 2010
Hi,
what is the difference between "index" and "value"?
Essentially, I am looking for the information postx provides when
figuring out history points. I browsed the postx source code, located
the routine and implemented it in .usr, but I expect trouble when moving
to bigger domains. I also figured there must be something in nek.
So I took the new hpts routine (see attached my modified .usr file),
what I am getting for r,s,t for a specific point is:
1.00000000000000 -5.551115123125783E-017 -1.00000000000000
The real deal for the same point (from postx) is 4 1 1.
What am I missing? Is there a way to go from the output I have to the
r,s,t indices?
Thanks for any help,
Markus
nek5000-users at lists.mcs.anl.gov wrote:
>
> Markus,
>
> those r,s,t values (not indices) are valid - do they not correspond
> to your expectations? Off-hand one might expect that the associated
> indices would be 1,NY1/2,NZ1
>
> Paul
>
>
> On Fri, 9 Jul 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
>> Howdy Neks,
>>
>> is there a function in Nek that for a given point in x,y,z would
>> return processor number, element number and r,s,t index of the nearest
>> GLL grid point? I was desperately trying to figure out how the
>> interpolation routines (f.e. for history points) would help me, but I
>> only got garbage (-1, 0 and 1) back for the variables I would expect
>> contain r,s,t. It worked with the element number, however.
>>
>> Thanks,
>> Markus
>>
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: rsttest.usr
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20100709/74589837/attachment.ksh>
More information about the Nek5000-users
mailing list