March 2017 - Nek5000-users - mailman3.cels.anl.gov

Re: [Nek5000-users] Interpolation between two meshes
by nek5000-users＠lists.mcs.anl.gov 30 Mar '17

30 Mar '17

Hi Stefan, Many thanks for your reply. Please contact me via wzx1989221(a)gmail.com. Best regards, Tony ________________________________ From: nek5000-users-bounces(a)lists.mcs.anl.gov <nek5000-users-bounces(a)lists.mcs.anl.gov> on behalf of nek5000-users-request(a)lists.mcs.anl.gov <nek5000-users-request(a)lists.mcs.anl.gov> Sent: 30 March 2017 13:17 To: nek5000-users(a)lists.mcs.anl.gov Subject: Nek5000-users Digest, Vol 97, Issue 37 Today's Topics: 1. Outflow BC with negative velocity (nek5000-users(a)lists.mcs.anl.gov) 2. Re: NekNek with PN-PN (nek5000-users(a)lists.mcs.anl.gov) 3. Re: NekNek with PN-PN (nek5000-users(a)lists.mcs.anl.gov) 4. Re: Interpolation between two meshes (nek5000-users(a)lists.mcs.anl.gov) 5. Accessing the value of drag and lift coefficients in usrchk (nek5000-users(a)lists.mcs.anl.gov) Message: 4 Date: Thu, 30 Mar 2017 11:46:40 +0200 From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov <nek5000-users(a)lists.mcs.anl.gov> Subject: Re: [Nek5000-users] Interpolation between two meshes Message-ID: <mailman.4068.1490867217.2967.nek5000-users(a)lists.mcs.anl.gov> Content-Type: text/plain; charset="utf-8" Please let me know how I can contact you off list. Stefan? On 29 Mar 2017, at 14:59, "nek5000-users(a)lists.mcs.anl.gov <mailto:[email protected]> " <nek5000-users(a)lists.mcs.anl.gov <mailto:[email protected]> > wrote: 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'

1 0

Accessing the value of drag and lift coefficients in usrchk
by nek5000-users＠lists.mcs.anl.gov 30 Mar '17

30 Mar '17

Hi Nek, I was wondering whether it was possible to access the values of the drag and lift coefficient directly into usrchk? I know how to compute them but they only get printed in the logfile. I would need them for a control problem. Thanks a lot Jean-Christophe B. Loiseau Maître de conférence DynFluid, Arts et Métiers ParisTech Homepage <https://sites.google.com/site/loiseaujc/>

1 2

Re: [Nek5000-users] Interpolation between two meshes
by nek5000-users＠lists.mcs.anl.gov 30 Mar '17

30 Mar '17

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(a)lists.mcs.anl.gov <nek5000-users-bounces(a)lists.mcs.anl.gov> on behalf of nek5000-users-request(a)lists.mcs.anl.gov <nek5000-users-request(a)lists.mcs.anl.gov> Sent: 28 March 2017 22:15 To: nek5000-users(a)lists.mcs.anl.gov Subject: Nek5000-users Digest, Vol 97, Issue 35 Today's Topics: 1. Re: Interpolation between two meshes (nek5000-users(a)lists.mcs.anl.gov) 2. Re: Interpolation between two meshes (nek5000-users(a)lists.mcs.anl.gov) Message: 2 Date: Tue, 28 Mar 2017 23:14:56 +0200 From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov <nek5000-users(a)lists.mcs.anl.gov> Subject: Re: [Nek5000-users] Interpolation between two meshes Message-ID: <mailman.3928.1490735713.2967.nek5000-users(a)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:[email protected] <nek5000-users(a)lists.mcs.anl.gov> > Sent: Tuesday 28th March 2017 22:59 > To: nek5000-users(a)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

1 1

NekNek with PN-PN
by nek5000-users＠lists.mcs.anl.gov 30 Mar '17

30 Mar '17

Dear Nek users, I am coupling Navier-Stokes and Darcy equations using NekNek. I need to pass pressure from Navier-Stokes as the boundary condition for Darcy equation. Pn/Pn-2 formulations extrapolate pressure to GLL nodes and apply dsavg to compute C-0 continuous pressure field. Though the resulting pressure field is continuous it is not very smooth and exhibits oscillations near each element boundaries. Using this as the boundary condition for Darcy equation gives raise to convergence issues after some time steps. I tried to use neknek with Pn/Pn formulation which will avoid this problem. However, I endup with "ERROR: .... crs_solve: invalid handle" and the simulation breaks away after the first time step in "happy_check" subroutine. Is Neknek not implemented for Pn/Pn formulations? Thank you very much for your response, Best regards, Sudhakar Linne FLOW centre.

1 3

Outflow BC with negative velocity
by nek5000-users＠lists.mcs.anl.gov 29 Mar '17

29 Mar '17

Dear all, I've been having some trouble with the outflow BCs in an axisymmetric case with periodic velocity inflow conditions. Inflow velocity takes these values: ux=(sqrt(2.)/5.)*sin(2.*pi*0.03*time) It runs fine at Re=1000, but when set to Re=13329, which is the value I'm interested in testing, it crashes due to backflow near the boundary (with outflow BC). Is there any setting that allows the outflow to handle flow in both ways? I have already tried turb_outflow (trying to eliminate the negative flow), but it's not working properly with the axisymmetric condition. Thanks in advance for the help, MViturro

1 0

Re: [Nek5000-users] Interpolation between two meshes
by nek5000-users＠lists.mcs.anl.gov 28 Mar '17

28 Mar '17

Can you also call nekgsync() before and after findpts. Just trying to understand if it really hangs there. -----Original message----- > From:[email protected] <nek5000-users(a)lists.mcs.anl.gov> > Sent: Tuesday 28th March 2017 22:59 > To: nek5000-users(a)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 > > -----Original message----- > > From:[email protected] <nek5000-users(a)lists.mcs.anl.gov> > > Sent: Tuesday 28th March 2017 22:46 > > To: nek5000-users(a)lists.mcs.anl.gov > > Subject: Re: [Nek5000-users] Interpolation between two meshes > > > > Hi Stefan, > > > > In case the attachments didn't go through, please see below for the .usr file and log file. > > I only used 'userchk' and 'usrdat' in the .usr file for interpolation. The other subroutines are empty. > > (I added several print lines in 'gfldr.f' to see where the interpolation halted.) > > > > Best regards, > > Tony > > > > .USR FILE: > > c----------------------------------------------------------------------- > > subroutine userchk() > > include 'SIZE' > > include 'TOTAL' > > include 'ZPER' > > include 'CTIMER' > > include 'RESTART' > > > > character*132 sourcefld > > > > sourcefld='source.f00001' > > call gfldr(sourcefld) > > call outpost(vx,vy,vz,pr,t,'rst') > > call exitt() > > > > return > > end > > c----------------------------------------------------------------------- > > subroutine usrdat() ! This routine to modify element vertices > > include 'SIZE' > > include 'TOTAL' > > > > parameter (ielx=150) > > parameter (iely=22) > > parameter (ielz=142) > > integer i,e,eg,ex,ey,ez > > real cony > > > > cony = 2.36d0 > > > > do e=1,nelv > > eg = lglel(e) > > if(eg<=ielx*iely*ielz) then > > call get_exyz(ex,ey,ez,eg,ielx,iely,ielz) > > > > do i=1,8 > > if(i<=2 .or. i==5 .or. i==6) then > > if((yc(i,e)>0.0e0.and.yc(i,e)<1.0e0).OR. > > & (yc(i,e)>1.0e0.and.yc(i,e)<2.0e0)) then > > yc(i,e) = 2.0e0*Real(ey-1)/Real(iely) > > end if > > else if(i==3 .or. i==4 .or. i>=7) then > > if((yc(i,e)>0.0e0.and.yc(i,e)<1.0e0).OR. > > & (yc(i,e)>1.0e0.and.yc(i,e)<2.0e0)) then > > yc(i,e) = 2.0e0*Real(ey)/Real(iely) > > end if > > end if > > end do > > > > end if > > end do > > > > do e=1,nelv > > eg = lglel(e) > > if(eg==1) print*,yc(1,e),yc(2,e),yc(5,e),yc(6,e) > > if(eg==1) print*,yc(3,e),yc(4,e),yc(7,e),yc(8,e) > > do i=1,8 > > if(yc(i,e)>0.0e0 .AND. yc(i,e)<1.0e0) then > > yc(i,e) = yc(i,e)-1.0e0 > > yc(i,e) = tanh(cony*yc(i,e))/tanh(cony) > > yc(i,e) = yc(i,e)1.0e0 > > else if(yc(i,e)>1.0e0 .AND. yc(i,e)<2.0e0) then > > yc(i,e) = yc(i,e)-1.0e0 > > yc(i,e) = tanh(cony*yc(i,e))/tanh(cony) > > yc(i,e) = yc(i,e)1.0e0 > > end if > > end do > > end do > > > > return > > end > > c----------------------------------------------------------------------- > > > > > > LOG FILE: > > > > /----------------------------------------------------------\\ > > | _ __ ______ __ __ ______ ____ ____ ____ | > > | / | / // ____// //_/ / ____/ / __ \\/ __ \\/ __ \\ | > > | / |/ // __/ / ,< /___ \\ / / / // / / // / / / | > > | / /| // /___ / /| | ____/ / / /_/ // /_/ // /_/ / | > > | /_/ |_//_____//_/ |_|/_____/ \\___/ \\___/ \\___/ | > > | | > > |----------------------------------------------------------| > > | | > > | NEK5000: Open Source Spectral Element Solver | > > | COPYRIGHT (c) 2008-2017 UCHICAGO ARGONNE, LLC | > > | Version: 17.0.0 | > > | Web: http://nek5000.mcs.anl.gov | > > | | > > \\----------------------------------------------------------/ > > > > > > Number of processors: 720 > > REAL wdsize : 8 > > INTEGER wdsize : 4 > > Timer accuracy : 9.54E-08 > > Reading /work/e01/e01/tony/geometry_nonosci_cavity/interp.rea > > > > mapping elements to processors > > Reading /work/e01/e01/tony/geometry_nonosci_cavity/interp.map > > element load imbalance: 1 652 653 > > done :: mapping elements to processors > > > > 0 objects found > > 118 Parameters from file: > > 1 1.00000 p001 DENSITY > > 2 -3150.00 p002 VISCOS > > 12 -0.11750E-01 p012 DT > > 15 500.000 p015 IOSTEP > > 21 0.100000E-07 p021 DIVERGENCE > > 22 0.100000E-07 p022 HELMHOLTZ > > 26 2.00000 p026 COURANT/NTAU > > 27 3.00000 p027 TORDER > > 59 0.00000 p059 !=0 --> full Jac. eval. for each el. > > 63 8.00000 p063 =8 --> force 8-byte output > > 65 1.00000 p065 #iofiles (eg, 0 or 64); <0 --> sep. dirs > > 66 6.00000 p066 output : <0=ascii, else binary > > 67 6.00000 p067 restart: <0=ascii, else binary > > 93 20.0000 p093 Number of previous pressure solns saved > > 94 5.00000 p094 start projecting velocity after p94 step > > 95 5.00000 p095 start projecting pressure after p95 step > > 99 3.00000 p099 dealiasing: <0--> off/3--> old/4--> new > > 103 0.10000E-01 p103 weight of stabilizing filter (.01) > > > > nelgt/nelgv/lelt: 469500 469500 653 > > lx1 /lx2 /lx3 : 8 8 8 > > done :: read .rea file 45.222 sec > > > > setup mesh topology > > Right-handed check complete for 469500 elements. OK. > > setvert3d: 8 60855000 162267000 60855000 60855000 > > call usrsetvert > > done :: usrsetvert > > > > gs_setup: 10325502 unique labels shared > > pairwise times (avg, min, max): 0.00044986 0.00033021 0.000500798 > > crystal router : 0.00145524 0.00143301 0.0014781 > > used all_to_all method: pairwise > > handle bytes (avg, min, max): 2.64212e06 2463020 2796252 > > buffer bytes (avg, min, max): 491593 352288 632640 > > setupds time 5.3504E-01 seconds 0 8 60855000 469500 > > 8 max multiplicity > > done :: setup mesh topology > > > > 0.000000000000000E000 0.000000000000000E000 0.000000000000000E000 > > 0.000000000000000E000 > > 9.090909090909091E-002 9.090909090909091E-002 9.090909090909091E-002 > > 9.090909090909091E-002 > > call usrdat > > done :: usrdat > > > > generate geometry data > > NOTE: All elements deformed , param(59) ^=0 > > done :: generate geometry data > > > > call usrdat2 > > done :: usrdat2 > > > > regenerate geometry data 1 > > NOTE: All elements deformed , param(59) ^=0 > > done :: regenerate geometry data 1 > > > > verify mesh topology > > 0.000000000000000E000 6.00000000000000 Xrange > > -0.120000000000000 2.00000000000000 Yrange > > 0.000000000000000E000 49.7000000000000 Zrange > > done :: verify mesh topology > > > > IFTRAN = T > > IFFLOW = T > > IFHEAT = F > > IFSPLIT = T > > IFLOMACH = F > > IFUSERVP = F > > IFUSERMV = F > > IFSTRS = F > > IFCHAR = T > > IFCYCLIC = F > > IFAXIS = F > > IFMVBD = F > > IFMELT = F > > IFMODEL = F > > IFKEPS = F > > IFMOAB = F > > IFNEKNEK = F > > IFSYNC = T > > > > IFVCOR = F > > IFINTQ = F > > IFCWUZ = F > > IFSWALL = F > > IFGEOM = F > > IFSURT = F > > IFWCNO = F > > > > IFTMSH for field 1 = F > > IFADVC for field 1 = T > > IFNONL for field 1 = F > > > > Dealiasing enabled, lxd= 12 > > verify mesh topology > > 0.000000000000000E000 6.00000000000000 Xrange > > -0.120000000000000 2.00000000000000 Yrange > > 0.000000000000000E000 49.7000000000000 Zrange > > done :: verify mesh topology > > > > call usrdat3 > > done :: usrdat3 > > > > set initial conditions > > nekuic (1) for ifld 1 > > call nekuic for vel > > xyz min 0.0000 -0.12000 0.0000 > > uvwpt min 0.0000 0.0000 0.0000 0.0000 0.0000 > > xyz max 6.0000 2.0000 49.700 > > uvwpt max 0.0000 0.0000 0.0000 0.0000 0.0000 > > done :: set initial conditions > > > > gridpoints unique/tot: 162267000 240384000 > > dofs: 159972750 162104250 > > > > Initialization successfully completed 46.588 sec > > > > nsteps=0 -> skip time loop > > running solver in post processing mode > > > > call userchk > > call gfldr > > byte swap: F 6.543210 -2.9312772E35 > > sourcefld X 0.000000000000000E000 6.00000000000000 > > sourcefld Y -0.120000000000000 2.00000000000000 > > sourcefld Z 0.000000000000000E000 49.7000000000000 > > CHECK interp prepare > > CHECK interp read > > > > ----------- > > From: nek5000-users-bounces(a)lists.mcs.anl.gov <nek5000-users-bounces(a)lists.mcs.anl.gov> on behalf of nek5000-users-request(a)lists.mcs.anl.gov <nek5000-users-request(a)lists.mcs.anl.gov> > > Sent: 28 March 2017 20:51 > > To: nek5000-users(a)lists.mcs.anl.gov > > Subject: Nek5000-users Digest, Vol 97, Issue 33 > > > > Message: 3 > > Date: Tue, 28 Mar 2017 19:51:04 0000 > > From: nek5000-users(a)lists.mcs.anl.gov > > To: "nek5000-users(a)lists.mcs.anl.gov" > > <nek5000-users(a)lists.mcs.anl.gov> > > Subject: Re: [Nek5000-users] Interpolation between two meshes > > Message-ID: > > <mailman.3915.1490730701.2967.nek5000-users(a)lists.mcs.anl.gov> > > Content-Type: text/plain; charset="iso-8859-1" > > > > Hi Stefan, > > > > > > Many thanks for your reply. I'm running my simulations on ARCHER (Cray XC30 supercomputer) and I'm using the Intel compiler. > > > > > > The current simulation has 469500 elements with polynomial order of 7. The source field I used for interpolation has 383700 elements with the same polynomial order. The size of the source field is about 10GB. I'm not sure if it can be easily transferred. But > > I have attached my .usr and log files herewith (also the gfldr.f with few added checking lines). I ran the interpolation using 720 processors. > > > > > > Please have a look and let me know if you need further information about the case. Thank you very much for you help. > > > > > > Best regards, > > > > Tony > > _______________________________________________ > > Nek5000-users mailing list > > Nek5000-users(a)lists.mcs.anl.gov > > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > _______________________________________________ > Nek5000-users mailing list > Nek5000-users(a)lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users >

1 0

Re: [Nek5000-users] Interpolation between two meshes
by nek5000-users＠lists.mcs.anl.gov 28 Mar '17

28 Mar '17

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 -----Original message----- > From:[email protected] <nek5000-users(a)lists.mcs.anl.gov> > Sent: Tuesday 28th March 2017 22:46 > To: nek5000-users(a)lists.mcs.anl.gov > Subject: Re: [Nek5000-users] Interpolation between two meshes > > Hi Stefan, > > In case the attachments didn't go through, please see below for the .usr file and log file. > I only used 'userchk' and 'usrdat' in the .usr file for interpolation. The other subroutines are empty. > (I added several print lines in 'gfldr.f' to see where the interpolation halted.) > > Best regards, > Tony > > .USR FILE: > c----------------------------------------------------------------------- > subroutine userchk() > include 'SIZE' > include 'TOTAL' > include 'ZPER' > include 'CTIMER' > include 'RESTART' > > character*132 sourcefld > > sourcefld='source.f00001' > call gfldr(sourcefld) > call outpost(vx,vy,vz,pr,t,'rst') > call exitt() > > return > end > c----------------------------------------------------------------------- > subroutine usrdat() ! This routine to modify element vertices > include 'SIZE' > include 'TOTAL' > > parameter (ielx=150) > parameter (iely=22) > parameter (ielz=142) > integer i,e,eg,ex,ey,ez > real cony > > cony = 2.36d0 > > do e=1,nelv > eg = lglel(e) > if(eg<=ielx*iely*ielz) then > call get_exyz(ex,ey,ez,eg,ielx,iely,ielz) > > do i=1,8 > if(i<=2 .or. i==5 .or. i==6) then > if((yc(i,e)>0.0e0.and.yc(i,e)<1.0e0).OR. > & (yc(i,e)>1.0e0.and.yc(i,e)<2.0e0)) then > yc(i,e) = 2.0e0*Real(ey-1)/Real(iely) > end if > else if(i==3 .or. i==4 .or. i>=7) then > if((yc(i,e)>0.0e0.and.yc(i,e)<1.0e0).OR. > & (yc(i,e)>1.0e0.and.yc(i,e)<2.0e0)) then > yc(i,e) = 2.0e0*Real(ey)/Real(iely) > end if > end if > end do > > end if > end do > > do e=1,nelv > eg = lglel(e) > if(eg==1) print*,yc(1,e),yc(2,e),yc(5,e),yc(6,e) > if(eg==1) print*,yc(3,e),yc(4,e),yc(7,e),yc(8,e) > do i=1,8 > if(yc(i,e)>0.0e0 .AND. yc(i,e)<1.0e0) then > yc(i,e) = yc(i,e)-1.0e0 > yc(i,e) = tanh(cony*yc(i,e))/tanh(cony) > yc(i,e) = yc(i,e)1.0e0 > else if(yc(i,e)>1.0e0 .AND. yc(i,e)<2.0e0) then > yc(i,e) = yc(i,e)-1.0e0 > yc(i,e) = tanh(cony*yc(i,e))/tanh(cony) > yc(i,e) = yc(i,e)1.0e0 > end if > end do > end do > > return > end > c----------------------------------------------------------------------- > > > LOG FILE: > > /----------------------------------------------------------\\ > | _ __ ______ __ __ ______ ____ ____ ____ | > | / | / // ____// //_/ / ____/ / __ \\/ __ \\/ __ \\ | > | / |/ // __/ / ,< /___ \\ / / / // / / // / / / | > | / /| // /___ / /| | ____/ / / /_/ // /_/ // /_/ / | > | /_/ |_//_____//_/ |_|/_____/ \\___/ \\___/ \\___/ | > | | > |----------------------------------------------------------| > | | > | NEK5000: Open Source Spectral Element Solver | > | COPYRIGHT (c) 2008-2017 UCHICAGO ARGONNE, LLC | > | Version: 17.0.0 | > | Web: http://nek5000.mcs.anl.gov | > | | > \\----------------------------------------------------------/ > > > Number of processors: 720 > REAL wdsize : 8 > INTEGER wdsize : 4 > Timer accuracy : 9.54E-08 > Reading /work/e01/e01/tony/geometry_nonosci_cavity/interp.rea > > mapping elements to processors > Reading /work/e01/e01/tony/geometry_nonosci_cavity/interp.map > element load imbalance: 1 652 653 > done :: mapping elements to processors > > 0 objects found > 118 Parameters from file: > 1 1.00000 p001 DENSITY > 2 -3150.00 p002 VISCOS > 12 -0.11750E-01 p012 DT > 15 500.000 p015 IOSTEP > 21 0.100000E-07 p021 DIVERGENCE > 22 0.100000E-07 p022 HELMHOLTZ > 26 2.00000 p026 COURANT/NTAU > 27 3.00000 p027 TORDER > 59 0.00000 p059 !=0 --> full Jac. eval. for each el. > 63 8.00000 p063 =8 --> force 8-byte output > 65 1.00000 p065 #iofiles (eg, 0 or 64); <0 --> sep. dirs > 66 6.00000 p066 output : <0=ascii, else binary > 67 6.00000 p067 restart: <0=ascii, else binary > 93 20.0000 p093 Number of previous pressure solns saved > 94 5.00000 p094 start projecting velocity after p94 step > 95 5.00000 p095 start projecting pressure after p95 step > 99 3.00000 p099 dealiasing: <0--> off/3--> old/4--> new > 103 0.10000E-01 p103 weight of stabilizing filter (.01) > > nelgt/nelgv/lelt: 469500 469500 653 > lx1 /lx2 /lx3 : 8 8 8 > done :: read .rea file 45.222 sec > > setup mesh topology > Right-handed check complete for 469500 elements. OK. > setvert3d: 8 60855000 162267000 60855000 60855000 > call usrsetvert > done :: usrsetvert > > gs_setup: 10325502 unique labels shared > pairwise times (avg, min, max): 0.00044986 0.00033021 0.000500798 > crystal router : 0.00145524 0.00143301 0.0014781 > used all_to_all method: pairwise > handle bytes (avg, min, max): 2.64212e06 2463020 2796252 > buffer bytes (avg, min, max): 491593 352288 632640 > setupds time 5.3504E-01 seconds 0 8 60855000 469500 > 8 max multiplicity > done :: setup mesh topology > > 0.000000000000000E000 0.000000000000000E000 0.000000000000000E000 > 0.000000000000000E000 > 9.090909090909091E-002 9.090909090909091E-002 9.090909090909091E-002 > 9.090909090909091E-002 > call usrdat > done :: usrdat > > generate geometry data > NOTE: All elements deformed , param(59) ^=0 > done :: generate geometry data > > call usrdat2 > done :: usrdat2 > > regenerate geometry data 1 > NOTE: All elements deformed , param(59) ^=0 > done :: regenerate geometry data 1 > > verify mesh topology > 0.000000000000000E000 6.00000000000000 Xrange > -0.120000000000000 2.00000000000000 Yrange > 0.000000000000000E000 49.7000000000000 Zrange > done :: verify mesh topology > > IFTRAN = T > IFFLOW = T > IFHEAT = F > IFSPLIT = T > IFLOMACH = F > IFUSERVP = F > IFUSERMV = F > IFSTRS = F > IFCHAR = T > IFCYCLIC = F > IFAXIS = F > IFMVBD = F > IFMELT = F > IFMODEL = F > IFKEPS = F > IFMOAB = F > IFNEKNEK = F > IFSYNC = T > > IFVCOR = F > IFINTQ = F > IFCWUZ = F > IFSWALL = F > IFGEOM = F > IFSURT = F > IFWCNO = F > > IFTMSH for field 1 = F > IFADVC for field 1 = T > IFNONL for field 1 = F > > Dealiasing enabled, lxd= 12 > verify mesh topology > 0.000000000000000E000 6.00000000000000 Xrange > -0.120000000000000 2.00000000000000 Yrange > 0.000000000000000E000 49.7000000000000 Zrange > done :: verify mesh topology > > call usrdat3 > done :: usrdat3 > > set initial conditions > nekuic (1) for ifld 1 > call nekuic for vel > xyz min 0.0000 -0.12000 0.0000 > uvwpt min 0.0000 0.0000 0.0000 0.0000 0.0000 > xyz max 6.0000 2.0000 49.700 > uvwpt max 0.0000 0.0000 0.0000 0.0000 0.0000 > done :: set initial conditions > > gridpoints unique/tot: 162267000 240384000 > dofs: 159972750 162104250 > > Initialization successfully completed 46.588 sec > > nsteps=0 -> skip time loop > running solver in post processing mode > > call userchk > call gfldr > byte swap: F 6.543210 -2.9312772E35 > sourcefld X 0.000000000000000E000 6.00000000000000 > sourcefld Y -0.120000000000000 2.00000000000000 > sourcefld Z 0.000000000000000E000 49.7000000000000 > CHECK interp prepare > CHECK interp read > > ----------- > From: nek5000-users-bounces(a)lists.mcs.anl.gov <nek5000-users-bounces(a)lists.mcs.anl.gov> on behalf of nek5000-users-request(a)lists.mcs.anl.gov <nek5000-users-request(a)lists.mcs.anl.gov> > Sent: 28 March 2017 20:51 > To: nek5000-users(a)lists.mcs.anl.gov > Subject: Nek5000-users Digest, Vol 97, Issue 33 > > Message: 3 > Date: Tue, 28 Mar 2017 19:51:04 0000 > From: nek5000-users(a)lists.mcs.anl.gov > To: "nek5000-users(a)lists.mcs.anl.gov" > <nek5000-users(a)lists.mcs.anl.gov> > Subject: Re: [Nek5000-users] Interpolation between two meshes > Message-ID: > <mailman.3915.1490730701.2967.nek5000-users(a)lists.mcs.anl.gov> > Content-Type: text/plain; charset="iso-8859-1" > > Hi Stefan, > > > Many thanks for your reply. I'm running my simulations on ARCHER (Cray XC30 supercomputer) and I'm using the Intel compiler. > > > The current simulation has 469500 elements with polynomial order of 7. The source field I used for interpolation has 383700 elements with the same polynomial order. The size of the source field is about 10GB. I'm not sure if it can be easily transferred. But > I have attached my .usr and log files herewith (also the gfldr.f with few added checking lines). I ran the interpolation using 720 processors. > > > Please have a look and let me know if you need further information about the case. Thank you very much for you help. > > > Best regards, > > Tony > _______________________________________________ > Nek5000-users mailing list > Nek5000-users(a)lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

1 0

Re: [Nek5000-users] Interpolation between two meshes
by nek5000-users＠lists.mcs.anl.gov 28 Mar '17

28 Mar '17

Hi Stefan, In case the attachments didn't go through, please see below for the .usr file and log file. I only used 'userchk' and 'usrdat' in the .usr file for interpolation. The other subroutines are empty. (I added several print lines in 'gfldr.f' to see where the interpolation halted.) Best regards, Tony .USR FILE: c----------------------------------------------------------------------- subroutine userchk() include 'SIZE' include 'TOTAL' include 'ZPER' include 'CTIMER' include 'RESTART' character*132 sourcefld sourcefld='source.f00001' call gfldr(sourcefld) call outpost(vx,vy,vz,pr,t,'rst') call exitt() return end c----------------------------------------------------------------------- subroutine usrdat() ! This routine to modify element vertices include 'SIZE' include 'TOTAL' parameter (ielx=150) parameter (iely=22) parameter (ielz=142) integer i,e,eg,ex,ey,ez real cony cony = 2.36d0 do e=1,nelv eg = lglel(e) if(eg<=ielx*iely*ielz) then call get_exyz(ex,ey,ez,eg,ielx,iely,ielz) do i=1,8 if(i<=2 .or. i==5 .or. i==6) then if((yc(i,e)>0.0e0.and.yc(i,e)<1.0e0).OR. & (yc(i,e)>1.0e0.and.yc(i,e)<2.0e0)) then yc(i,e) = 2.0e0*Real(ey-1)/Real(iely) end if else if(i==3 .or. i==4 .or. i>=7) then if((yc(i,e)>0.0e0.and.yc(i,e)<1.0e0).OR. & (yc(i,e)>1.0e0.and.yc(i,e)<2.0e0)) then yc(i,e) = 2.0e0*Real(ey)/Real(iely) end if end if end do end if end do do e=1,nelv eg = lglel(e) if(eg==1) print*,yc(1,e),yc(2,e),yc(5,e),yc(6,e) if(eg==1) print*,yc(3,e),yc(4,e),yc(7,e),yc(8,e) do i=1,8 if(yc(i,e)>0.0e0 .AND. yc(i,e)<1.0e0) then yc(i,e) = yc(i,e)-1.0e0 yc(i,e) = tanh(cony*yc(i,e))/tanh(cony) yc(i,e) = yc(i,e)+1.0e0 else if(yc(i,e)>1.0e0 .AND. yc(i,e)<2.0e0) then yc(i,e) = yc(i,e)-1.0e0 yc(i,e) = tanh(cony*yc(i,e))/tanh(cony) yc(i,e) = yc(i,e)+1.0e0 end if end do end do return end c----------------------------------------------------------------------- LOG FILE: /----------------------------------------------------------\\ | _ __ ______ __ __ ______ ____ ____ ____ | | / | / // ____// //_/ / ____/ / __ \\/ __ \\/ __ \\ | | / |/ // __/ / ,< /___ \\ / / / // / / // / / / | | / /| // /___ / /| | ____/ / / /_/ // /_/ // /_/ / | | /_/ |_//_____//_/ |_|/_____/ \\___/ \\___/ \\___/ | | | |----------------------------------------------------------| | | | NEK5000: Open Source Spectral Element Solver | | COPYRIGHT (c) 2008-2017 UCHICAGO ARGONNE, LLC | | Version: 17.0.0 | | Web: http://nek5000.mcs.anl.gov | | | \\----------------------------------------------------------/ Number of processors: 720 REAL wdsize : 8 INTEGER wdsize : 4 Timer accuracy : 9.54E-08 Reading /work/e01/e01/tony/geometry_nonosci_cavity/interp.rea mapping elements to processors Reading /work/e01/e01/tony/geometry_nonosci_cavity/interp.map element load imbalance: 1 652 653 done :: mapping elements to processors 0 objects found 118 Parameters from file: 1 1.00000 p001 DENSITY 2 -3150.00 p002 VISCOS 12 -0.11750E-01 p012 DT 15 500.000 p015 IOSTEP 21 0.100000E-07 p021 DIVERGENCE 22 0.100000E-07 p022 HELMHOLTZ 26 2.00000 p026 COURANT/NTAU 27 3.00000 p027 TORDER 59 0.00000 p059 !=0 --> full Jac. eval. for each el. 63 8.00000 p063 =8 --> force 8-byte output 65 1.00000 p065 #iofiles (eg, 0 or 64); <0 --> sep. dirs 66 6.00000 p066 output : <0=ascii, else binary 67 6.00000 p067 restart: <0=ascii, else binary 93 20.0000 p093 Number of previous pressure solns saved 94 5.00000 p094 start projecting velocity after p94 step 95 5.00000 p095 start projecting pressure after p95 step 99 3.00000 p099 dealiasing: <0--> off/3--> old/4--> new 103 0.10000E-01 p103 weight of stabilizing filter (.01) nelgt/nelgv/lelt: 469500 469500 653 lx1 /lx2 /lx3 : 8 8 8 done :: read .rea file 45.222 sec setup mesh topology Right-handed check complete for 469500 elements. OK. setvert3d: 8 60855000 162267000 60855000 60855000 call usrsetvert done :: usrsetvert gs_setup: 10325502 unique labels shared pairwise times (avg, min, max): 0.00044986 0.00033021 0.000500798 crystal router : 0.00145524 0.00143301 0.0014781 used all_to_all method: pairwise handle bytes (avg, min, max): 2.64212e+06 2463020 2796252 buffer bytes (avg, min, max): 491593 352288 632640 setupds time 5.3504E-01 seconds 0 8 60855000 469500 8 max multiplicity done :: setup mesh topology 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000 0.000000000000000E+000 9.090909090909091E-002 9.090909090909091E-002 9.090909090909091E-002 9.090909090909091E-002 call usrdat done :: usrdat generate geometry data NOTE: All elements deformed , param(59) ^=0 done :: generate geometry data call usrdat2 done :: usrdat2 regenerate geometry data 1 NOTE: All elements deformed , param(59) ^=0 done :: regenerate geometry data 1 verify mesh topology 0.000000000000000E+000 6.00000000000000 Xrange -0.120000000000000 2.00000000000000 Yrange 0.000000000000000E+000 49.7000000000000 Zrange done :: verify mesh topology IFTRAN = T IFFLOW = T IFHEAT = F IFSPLIT = T IFLOMACH = F IFUSERVP = F IFUSERMV = F IFSTRS = F IFCHAR = T IFCYCLIC = F IFAXIS = F IFMVBD = F IFMELT = F IFMODEL = F IFKEPS = F IFMOAB = F IFNEKNEK = F IFSYNC = T IFVCOR = F IFINTQ = F IFCWUZ = F IFSWALL = F IFGEOM = F IFSURT = F IFWCNO = F IFTMSH for field 1 = F IFADVC for field 1 = T IFNONL for field 1 = F Dealiasing enabled, lxd= 12 verify mesh topology 0.000000000000000E+000 6.00000000000000 Xrange -0.120000000000000 2.00000000000000 Yrange 0.000000000000000E+000 49.7000000000000 Zrange done :: verify mesh topology call usrdat3 done :: usrdat3 set initial conditions nekuic (1) for ifld 1 call nekuic for vel xyz min 0.0000 -0.12000 0.0000 uvwpt min 0.0000 0.0000 0.0000 0.0000 0.0000 xyz max 6.0000 2.0000 49.700 uvwpt max 0.0000 0.0000 0.0000 0.0000 0.0000 done :: set initial conditions gridpoints unique/tot: 162267000 240384000 dofs: 159972750 162104250 Initialization successfully completed 46.588 sec nsteps=0 -> skip time loop running solver in post processing mode call userchk call gfldr byte swap: F 6.543210 -2.9312772E+35 sourcefld X 0.000000000000000E+000 6.00000000000000 sourcefld Y -0.120000000000000 2.00000000000000 sourcefld Z 0.000000000000000E+000 49.7000000000000 CHECK interp prepare CHECK interp read ________________________________ From: nek5000-users-bounces(a)lists.mcs.anl.gov <nek5000-users-bounces(a)lists.mcs.anl.gov> on behalf of nek5000-users-request(a)lists.mcs.anl.gov <nek5000-users-request(a)lists.mcs.anl.gov> Sent: 28 March 2017 20:51 To: nek5000-users(a)lists.mcs.anl.gov Subject: Nek5000-users Digest, Vol 97, Issue 33 Message: 3 Date: Tue, 28 Mar 2017 19:51:04 +0000 From: nek5000-users(a)lists.mcs.anl.gov To: "nek5000-users(a)lists.mcs.anl.gov" <nek5000-users(a)lists.mcs.anl.gov> Subject: Re: [Nek5000-users] Interpolation between two meshes Message-ID: <mailman.3915.1490730701.2967.nek5000-users(a)lists.mcs.anl.gov> Content-Type: text/plain; charset="iso-8859-1" Hi Stefan, Many thanks for your reply. I'm running my simulations on ARCHER (Cray XC30 supercomputer) and I'm using the Intel compiler. The current simulation has 469500 elements with polynomial order of 7. The source field I used for interpolation has 383700 elements with the same polynomial order. The size of the source field is about 10GB. I'm not sure if it can be easily transferred. But I have attached my .usr and log files herewith (also the gfldr.f with few added checking lines). I ran the interpolation using 720 processors. Please have a look and let me know if you need further information about the case. Thank you very much for you help. Best regards, Tony

1 0

Re: [Nek5000-users] Interpolation between two meshes
by nek5000-users＠lists.mcs.anl.gov 28 Mar '17

28 Mar '17

Hi Stefan, Many thanks for your reply. I'm running my simulations on ARCHER (Cray XC30 supercomputer) and I'm using the Intel compiler. The current simulation has 469500 elements with polynomial order of 7. The source field I used for interpolation has 383700 elements with the same polynomial order. The size of the source field is about 10GB. I'm not sure if it can be easily transferred. But I have attached my .usr and log files herewith (also the gfldr.f with few added checking lines). I ran the interpolation using 720 processors. Please have a look and let me know if you need further information about the case. Thank you very much for you help. Best regards, Tony ________________________________ From: nek5000-users-bounces(a)lists.mcs.anl.gov <nek5000-users-bounces(a)lists.mcs.anl.gov> on behalf of nek5000-users-request(a)lists.mcs.anl.gov <nek5000-users-request(a)lists.mcs.anl.gov> Sent: 28 March 2017 13:19 To: nek5000-users(a)lists.mcs.anl.gov Subject: Nek5000-users Digest, Vol 97, Issue 31 Send Nek5000-users mailing list submissions to nek5000-users(a)lists.mcs.anl.gov To subscribe or unsubscribe via the World Wide Web, visit https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users or, via email, send a message with subject or body 'help' to nek5000-users-request(a)lists.mcs.anl.gov You can reach the person managing the list at nek5000-users-owner(a)lists.mcs.anl.gov When replying, please edit your Subject line so it is more specific than "Re: Contents of Nek5000-users digest..." Today's Topics: 1. Re: Interpolation between two meshes (nek5000-users(a)lists.mcs.anl.gov) 2. Re: Copy velocity slice from one boundary to another (nek5000-users(a)lists.mcs.anl.gov) ---------------------------------------------------------------------- Message: 1 Date: Tue, 28 Mar 2017 09:57:33 +0200 From: nek5000-users(a)lists.mcs.anl.gov To: nek5000-users(a)lists.mcs.anl.gov <nek5000-users(a)lists.mcs.anl.gov> Subject: Re: [Nek5000-users] Interpolation between two meshes Message-ID: <mailman.3869.1490687873.2967.nek5000-users(a)lists.mcs.anl.gov> Content-Type: text/plain; charset=utf-8 Can you please provide more information about your case (machine, compiler) and include the .usr and log file. If the case is reasonably small, I can try to run it myself. -----Original message----- > From:[email protected] <nek5000-users(a)lists.mcs.anl.gov> > Sent: Tuesday 28th March 2017 4:03 > To: nek5000-users(a)lists.mcs.anl.gov > Subject: Re: [Nek5000-users] Interpolation between two meshes > > Hi Stefan, > > Many thanks for your reply. I have tried the new 'gfldr' subroutine in two ways. > > First, I copied the whole 'gfldr.f' into my usr file and ran the interpolation with my legacy nek master code > (tuned few lines to match the old byte_mpi.f format). The source field was read in successfully and I also confirmed the coordinates from the source field. However, the interpolation halted again at 'findpts' in 'gfldr_intp' subroutine. > > Second, I ran the interpolation directly with the latest master nek code. The same thing happened and the interpolation halted at 'findpts' and never moved on. > > In addition, I also tried reducing polynomial order (lx1) as Mattias suggested, but with no success. > > I really don't know what causes the problem when interpolating between two large meshes. Do you have any workaround? > > Please help me out. Thank you very much. > > Best regards, > Tony > > > ----------- > From: nek5000-users-bounces(a)lists.mcs.anl.gov <nek5000-users-bounces(a)lists.mcs.anl.gov> on behalf of nek5000-users-request(a)lists.mcs.anl.gov <nek5000-users-request(a)lists.mcs.anl.gov> > Sent: 27 March 2017 18:00 > To: nek5000-users(a)lists.mcs.anl.gov > Subject: Nek5000-users Digest, Vol 97, Issue 29 > Send Nek5000-users mailing list submissions to > ??????? nek5000-users(a)lists.mcs.anl.gov > > To subscribe or unsubscribe via the World Wide Web, visit > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users> > or, via email, send a message with subject or body 'help' to > ??????? nek5000-users-request(a)lists.mcs.anl.gov > > You can reach the person managing the list at > ??????? nek5000-users-owner(a)lists.mcs.anl.gov > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Nek5000-users digest..." > > > Today's Topics: > > ?? 1. Re: Interpolation between two meshes > ????? (nek5000-users(a)lists.mcs.anl.gov) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 27 Mar 2017 18:21:52 0200 > From: nek5000-users(a)lists.mcs.anl.gov > To: nek5000-users(a)lists.mcs.anl.gov????? <nek5000-users(a)lists.mcs.anl.gov> > Subject: Re: [Nek5000-users] Interpolation between two meshes > Message-ID: > ??????? <mailman.3273.1490631729.2967.nek5000-users(a)lists.mcs.anl.gov> > Content-Type: text/plain; charset=utf-8 > > Just give it a whirl. > It's important for us to see if the latest version fixes your problem. I have some confidence that it will ;) > ? > -----Original message----- > > From:[email protected] <nek5000-users(a)lists.mcs.anl.gov> > > Sent: Monday 27th March 2017 18:19 > > To: nek5000-users(a)lists.mcs.anl.gov > > Subject: Re: [Nek5000-users] Interpolation between two meshes > > > > Hi Stefan and Mattias, > > > > Thank you both for your kind advice. > > > > Stefan, I had a look at the new 'gfldr' subroutine. It seems that the spectral interpolation procedure is the same as before, however, the latest subroutine breaks down the number of fields to interpolate in one go. I still don't see why the old g2gi hangs > >? at 'findpts'. Is it mainly because of memory size? (I have been using my own copy of the legacy nek, so I'm quite keen on making it work with g2gi.) > > > > Mattias, I will give it a go with your suggestion. By reducing the polynomial order, did you mean the order of the old field or the order of the new geometry field, or both? > > > > Thank you again for both your help. > > > > Best regards, > > Tony > > > > > > ----------- > > From: nek5000-users-bounces(a)lists.mcs.anl.gov <nek5000-users-bounces(a)lists.mcs.anl.gov> on behalf of nek5000-users-request(a)lists.mcs.anl.gov <nek5000-users-request(a)lists.mcs.anl.gov> > > Sent: 27 March 2017 12:57 > > To: nek5000-users(a)lists.mcs.anl.gov > > Subject: Nek5000-users Digest, Vol 97, Issue 27 > > Send Nek5000-users mailing list submissions to > > ??????? nek5000-users(a)lists.mcs.anl.gov > > > > To subscribe or unsubscribe via the World Wide Web, visit > > > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users> <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>> > > or, via email, send a message with subject or body 'help' to > > ??????? nek5000-users-request(a)lists.mcs.anl.gov > > > > You can reach the person managing the list at > > ??????? nek5000-users-owner(a)lists.mcs.anl.gov > > > > When replying, please edit your Subject line so it is more specific > > than "Re: Contents of Nek5000-users digest..." > > > > > > Today's Topics: > > > > ?? 1. Interpolation between two meshes (nek5000-users(a)lists.mcs.anl.gov) > > ?? 2. Re: Interpolation between two meshes > > ????? (nek5000-users(a)lists.mcs.anl.gov) > > ?? 3. Re: Copy velocity slice from one boundary to another > > ????? (nek5000-users(a)lists.mcs.anl.gov) > > ?? 4. Robin Boundary Condition -- Navier's Slip Length > > ????? (nek5000-users(a)lists.mcs.anl.gov) > > ?? 5. Re: Interpolation between two meshes > > ????? (nek5000-users(a)lists.mcs.anl.gov) > > > > > > ---------------------------------------------------------------------- > > > > Message: 1 > > Date: Mon, 27 Mar 2017 04:37:46 0000 > > From: nek5000-users(a)lists.mcs.anl.gov > > To: Nek5000 <nek5000-users(a)lists.mcs.anl.gov> > > Subject: [Nek5000-users] Interpolation between two meshes > > Message-ID: > > ??????? <mailman.3211.1490589471.2967.nek5000-users(a)lists.mcs.anl.gov> > > Content-Type: text/plain; charset="iso-8859-1" > > > > Dear Neks, > > > > I have been using 'g2gi' to do interpolation between different meshes. Recently, I have noticed that when I interpolate between two big meshes (say 300000 elements in each mesh), the interpolation hangs at 'call findpts' and never moves on. I tried increasing > >? number of processors and changing 'lpart', but none of these worked. > > > > I'm wondering what could cause the interpolation to hang at 'findpts'. Has anyone encountered similar issue? > > > > Hope anyone could give me some advice. Thank you very much in advance. > > > > Best regards, > > Tony > >

1 0

Re: [Nek5000-users] Project a solution onto a polynomial space
by nek5000-users＠lists.mcs.anl.gov 28 Mar '17

28 Mar '17

Dear all, Thank you for your responses. I've taken some time to think about your comments and I have some more questions: Paul: Are you suggesting to represent the moments on the 2*M grid and not project it back onto the original space? I am assuming this means that the moments would be in the space of polynomials of order 2*N. Am I understanding correctly? Aleks: You talk about truncating to polyn. order < lx1. Would this have the same effect as projecting onto the original space? Also, is the "cut_off_filter" subroutine the correct way to implement this? Thanks, Juan Diego

1 0