[Nek5000-users] question on subroutine sem2lex

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Mon Jun 2 19:57:30 CDT 2014


Hi,Aleks:
    Thank you for your suggestion.
Xianbei







At 2014-06-03 00:13:58, nek5000-users at lists.mcs.anl.gov wrote:

Hi Xianbei,

I would suggest checking out a stand-alone interpolation tool

nek5_svn/trunk/tools/int_tp/README

Best.
Aleks


From:nek5000-users-bounces at lists.mcs.anl.gov [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: Monday, June 02, 2014 8:05 AM
To: Nek5000
Subject: Re: [Nek5000-users] question on subroutine sem2lex


Hi,Paul:
    You are right. The code seems just out put the u in another ix,iy,iz,le way in spectral data, not what we expected as ix,iy,iz, where ix(max)=(nelx-1)*lx1+1. So is there any better way to output the data ready for the energy spectrum and two-point correlation already exists?

Xianbei

 







At 2014-06-02 08:22:55, nek5000-users at lists.mcs.anl.gov wrote:



Hi Xianbei,


It should be eg.   


However, note that it's really not as simple as this....You'll need to understand
where your data resides, etc.


I think this code example is quite dated and definitely not tested as it is not
part of the standard source.


Paul


From:nek5000-users-bounces at lists.mcs.anl.gov [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: Monday, June 02, 2014 12:46 AM
To: Nek5000
Subject: [Nek5000-users] question on subroutine sem2lex


Hi,all:
    In the parallel version of sem2lex, there is a part of code wrote as following:
if (gllnid(eg).eq.nid) then  ! element eg is on this processor
           eg = eg+1    
            le = gllel(ge)            ! le is local element number
            do iz=1,nz1
            do iy=1,ny1
            do ix=1,nx1

However, ge doesn't appear anywhere except here,(le=gllel(ge)), so is this should be "le=gllel(eg)"? Or just my misunderstanding?

Xianbei




来自网易手机号码邮箱了解更多



来自网易手机号码邮箱了解更多
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140603/92bc36f5/attachment.html>


More information about the Nek5000-users mailing list