<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: Times New Roman; font-size: 12pt; color: #000000'><P>Hi Paul,</P>
<P> </P>
<P>To expand on this topic, I am wondering what the convention being used to number the GL points for an element is? Could you explain this? Thinking of it terms of the r-s-t diagram for a nekton element, I know from previous examples that face 5 contains points 1 to lx1*ly1. Also, face 6 should contain points (lz1-1)*lx1*ly1 to lx1*ly1*lz1. </P>
<P> </P>
<P>What I am trying to figure out, is what points are on the remaining faces, and understanding the convention I think would be useful. Thanks for the help!</P>
<P> </P>
<P>- Michael</P>
<P><BR>----- Original Message -----<BR>From: nek5000-users@lists.mcs.anl.gov<BR>To: nek5000-users@lists.mcs.anl.gov<BR>Sent: Monday, July 12, 2010 5:13:47 AM GMT -06:00 US/Canada Central<BR>Subject: Re: [Nek5000-users] r,s,t index for given point<BR><BR>fixed in r548.<BR>Stefan<BR><BR>On Jul 12, 2010, at 10:49 AM, nek5000-users@lists.mcs.anl.gov wrote:<BR><BR>> <BR>> Hi Markus,<BR>> <BR>> lglel reads "local to global element number" so input is<BR>> local element number and output is global element number.<BR>> <BR>> I'll be seeing Stefan this week so I'm sure we can straighten<BR>> out the interpolation routines.<BR>> <BR>> You're right about the other points - the ntotv wasn't flagged<BR>> by my compiler and the .his file was taken care of by my own<BR>> scripts, so I didn't catch either of them. Sorry about that.<BR>> I've updated the source so that those issues should be<BR>> resolved.<BR>> <BR>> Paul<BR>> <BR>> <BR>> On Mon, 12 Jul 2010, nek5000-users@lists.mcs.anl.gov wrote:<BR>> <BR>>> Hi,<BR>>> <BR>>> that was not it, I'm afraid (though I thought gllel transforms from local to global).<BR>>> <BR>>> But I also realized that the hpts routine does not deliver the correct values for me. To test this, I set up a simple 3d box with inflow, outflow and symmetry all around, Re very low, initial condition=boundary condition and 4 points.<BR>>> Please find all files attached, the output for the first point is OK, but then the values seem to shift positions and don't make sense.<BR>>> <BR>>> Also, when running it for 1 step or more, I get the error message:<BR>>> "<BR>>> dump history points<BR>>> done :: dump history points<BR>>> hisfile:<BR>>> /Users/m0s1978/nektonruns/rectest/rectest.his <BR>>> ERROR: .sch file already exists. ierr= 1<BR>>> "<BR>>> <BR>>> And lastly, the newest repo won't compile, here is the error message:<BR>>> "<BR>>> Error: /Users/m0s1978/nek5trial/trunk/nek/drive2.f, line 1551: This name has already been assigned a data type. [NVTOT]<BR>>> integer*8 ntot,nvtot,nptot<BR>>> ---------------------^<BR>>> "<BR>>> <BR>>> Thanks for any help,<BR>>> Markus<BR>>> <BR>>> <BR>>> nek5000-users@lists.mcs.anl.gov wrote:<BR>>>> Hi Markus,<BR>>>> It's possible that it's returning the local (to the processor)<BR>>>> element number. You would get the global (to the computational<BR>>>> domain, defined by the .rea file) element number via:<BR>>>> <BR>>>> integer e,eg<BR>>>> <BR>>>> eg = lglel(e)<BR>>>> Paul<BR>>>> On Sun, 11 Jul 2010, nek5000-users@lists.mcs.anl.gov wrote:<BR>>>>> Hi again,<BR>>>>> Regarding finding element numbers (such as what gllel would return) for arbitrary history points, I was experimenting with the history point routine hpts a little and it seems that this function call:<BR>>>>> <BR>>>>> call findpts(inth_hpts,rcode,1,<BR>>>>> & proc,1,<BR>>>>> & elid,1,<BR>>>>> & rst,ndim,<BR>>>>> & dist,1,<BR>>>>> & pts(1,1),ndim,<BR>>>>> & pts(2,1),ndim,<BR>>>>> & pts(3,1),ndim,npoints)<BR>>>>> does not return the same element number (elid) as it can be found in the .rea file, for instance.<BR>>>>> For a handful of manually set points, I checked that the returned element number does not provide the right coordinates to enclose the desired point.<BR>>>>> Is there a way to transform the function return or am I missing something else?<BR>>>>> Thanks,<BR>>>>> Markus<BR>>>>> _______________________________________________<BR>>>>> Nek5000-users mailing list<BR>>>>> Nek5000-users@lists.mcs.anl.gov<BR>>>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<BR>>>> _______________________________________________<BR>>>> Nek5000-users mailing list<BR>>>> Nek5000-users@lists.mcs.anl.gov<BR>>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<BR>>> <BR>> _______________________________________________<BR>> Nek5000-users mailing list<BR>> Nek5000-users@lists.mcs.anl.gov<BR>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<BR><BR>_______________________________________________<BR>Nek5000-users mailing list<BR>Nek5000-users@lists.mcs.anl.gov<BR>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<BR></P></div></body></html>