From nek5000-users at lists.mcs.anl.gov Mon Apr 1 19:13:27 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 1 Apr 2013 19:13:27 -0500 Subject: [Nek5000-users] Question on parameter GTHETA ( p004 ) Message-ID: Dear sir , I am learning how to use the given parameter list . My question is for the parameter p004 ( GTHETA ) : is it the gravitational acceleration g ( 9.81 m per square seconds ) ? Thank you for your help and for your time , Best Regards sincerely Giacinto M. F. Genco Master of Science in Nuclear Engineering , Purdue University Master of Science in Aerospace Engineering , Purdue University Italian " Laurea " in Mechanical Engineering , Politecnico di Bari , Bari , Italy . -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Apr 1 19:50:14 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 1 Apr 2013 19:50:14 -0500 (CDT) Subject: [Nek5000-users] Question on parameter GTHETA ( p004 ) In-Reply-To: Message-ID: Dear Giacinto, This feature is currently disabled. The best way to put gravity into the problem is to put it as a buoyancy force in either the y or z direction, depending on the orientation of your domain. An example of how this is done are in the Rayleigh Benard problems in the example suite. Best regards, Paul ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Monday, April 1, 2013 7:13:27 PM Subject: [Nek5000-users] Question on parameter GTHETA ( p004 ) Dear sir , I am learning how to use the given parameter list . My question is for the parameter p004 ( GTHETA ) : is it the gravitational acceleration g ( 9.81 m per square seconds ) ? Thank you for your help and for your time , Best Regards sincerely Giacinto M. F. Genco Master of Science in Nuclear Engineering , Purdue University Master of Science in Aerospace Engineering , Purdue University Italian " Laurea " in Mechanical Engineering , Politecnico di Bari , Bari , Italy . _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Wed Apr 3 12:24:46 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 3 Apr 2013 11:24:46 -0600 Subject: [Nek5000-users] scalar zero flux bc Message-ID: Hello Neks, I'm looking for guidance on implementing appropriate volume forcing and zero-net-flux boundary conditions for the scalar equation. I'm solving (Eq. 1) phi_t + U dot grad(phi) = - div(J) + k Laplacian(phi) where phi is the scalar, phi_t denotes time derivative, and capital letters are vector quantities. Here, J = J(phi) is modeled flux. I need zero net flux at the boundaries, which requires (-J + k grad(phi)) dot N = 0 or (Eq. 2) k grad(phi) dot N = J dot N where N is the boundary normal. For (Eq. 1), I am simply using qvol = -div(J) in userq, and I specify the constant diffusivity in uservp (udiff set for ifield.eq.2). This seems correct, but please let me know if not. In regard to the boundary condition (Eq. 2), what is the appropriate implementation in userbc and what boundary specification should be used in the .rea file ('c' or 'F')? Thanks, Mike From nek5000-users at lists.mcs.anl.gov Thu Apr 4 02:55:07 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 4 Apr 2013 02:55:07 -0500 (CDT) Subject: [Nek5000-users] scalar zero flux bc In-Reply-To: References: Message-ID: Hi Mike, You probably can get away with the following: > (Eq. 2) k grad(phi) dot N = J dot N Set the phi bc to be "f" (for "flux", evaluated in userbc). Then, in userbc, set integer e,eg e = gllel(eg) flux = scale(ix,iy,iz,e) where you precompute the array scale() in userchk, using the J from the last timestep. If you want to get fancy, you could use 2*J^n - J^{n-1}, which would give you a 2nd-order accurate BC. Note that "scale" should simply be the rhs of your Eq. 2. Nek sets this right hand side to be equal to k grad phi . nhat. I can't recall off the top of my head of q" (the surface flux) is the flux into the domain, or out of, so you should check the sign. I'm pretty certain that it is coming into the domain. The unit normal for certain points out of the domain. Please let me know if this works for you. Paul On Wed, 3 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > Hello Neks, > > I'm looking for guidance on implementing appropriate volume forcing and zero-net-flux boundary conditions for the scalar equation. > > I'm solving > > (Eq. 1) phi_t + U dot grad(phi) = - div(J) + k Laplacian(phi) > > where phi is the scalar, phi_t denotes time derivative, and capital letters are vector quantities. Here, J = J(phi) is modeled flux. > > I need zero net flux at the boundaries, which requires > > (-J + k grad(phi)) dot N = 0 > > or > > (Eq. 2) k grad(phi) dot N = J dot N > > where N is the boundary normal. > > For (Eq. 1), I am simply using > > qvol = -div(J) > > in userq, and I specify the constant diffusivity in uservp (udiff set for ifield.eq.2). This seems correct, but please let me know if not. > > In regard to the boundary condition (Eq. 2), what is the appropriate implementation in userbc and what boundary specification should be used in the .rea file ('c' or 'F')? > > Thanks, > Mike > > > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > From nek5000-users at lists.mcs.anl.gov Fri Apr 5 15:27:57 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Fri, 5 Apr 2013 16:27:57 -0400 Subject: [Nek5000-users] Nodes on the boundaries Message-ID: Hi Neks, I declared an array like: common /myarray/ tarray(lx1,ly1,lz1,lelv) and would like to set its value on the boundaries to zero. Is there any easy way of doing that?. A mask array?. Or do I need to stride through elements, faces, type of boundary conditions... Thanks. -- *Alex * -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Apr 15 10:09:16 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 15 Apr 2013 11:09:16 -0400 Subject: [Nek5000-users] Question - SIZE Parameters And Compilation Errors Message-ID: Hi All, Trying to run some initial test cases of a problems I'm working on: - Currently on single node (SMP Linux) - MPI turned off (IFMPI="false") - DOF at 5 (lx1=6, lxd=9, lelx=8) If I run with lp=1,lelt=2160, lelg=2160, for a 12x12x15 mesh, it compiles and runs fine. If I up it to lp=1,lelt=4500,lelg=4500, for a 15x15x20 mesh, it crashes the compile with: obj/drive2.o: In function `plan4_vol_': drive2.f:(.text+0x82f): relocation truncated to fit: R_X86_64_32S against symbol `scrvh_' defined in COMMON section in obj/hsmg.o drive2.f:(.text+0x837): relocation truncated to fit: R_X86_64_32S against symbol `scrvh_' defined in COMMON section in obj/hsmg.o obj/drive2.o: In function `plan3_vol_': drive2.f:(.text+0xac5): relocation truncated to fit: R_X86_64_32S against symbol `scrvh_' defined in COMMON section in obj/hsmg.o drive2.f:(.text+0xacd): relocation truncated to fit: R_X86_64_32S against symbol `scrvh_' defined in COMMON section in obj/hsmg.o obj/drive2.o: In function `plan2_vol_': drive2.f:(.text+0xeb8): relocation truncated to fit: R_X86_64_32S against symbol `scrvh_' defined in COMMON section in obj/hsmg.o drive2.f:(.text+0xec0): relocation truncated to fit: R_X86_64_32S against symbol `scrvh_' defined in COMMON section in obj/hsmg.o obj/drive2.o: In function `dofcnt_': drive2.f:(.text+0x170f): relocation truncated to fit: R_X86_64_PC32 against symbol `scrns_' defined in COMMON section in obj/convect.o obj/drive2.o: In function `settime_': drive2.f:(.text+0x7cc9): relocation truncated to fit: R_X86_64_PC32 against symbol `cprint_' defined in COMMON section in obj/navier1.o drive2.f:(.text+0x7cef): relocation truncated to fit: R_X86_64_PC32 against symbol `cprint_' defined in COMMON section in obj/navier1.o drive2.f:(.text+0x7d02): relocation truncated to fit: R_X86_64_PC32 against symbol `cprint_' defined in COMMON section in obj/navier1.o obj/plan4.o: In function `crespsp_': plan4.f:(.text+0xb23): additional relocation overflows omitted from the output collect2: ld returned 1 exit status make: *** [nek5000] Error 1 Just wondering if this is due to a true limit with the platform/word size, or if I'm screwing something up... Thanks, Murph PS - Is there any other interface to this list other than raw email? -- *Murphy Leo O'Dea, PE"Murph"* PhD Student, Mechanical Engineering /School Of Engineering And Computer Science/ /Oakland University/ mlodea at oakland.edu http://www.secs.oakland.edu/~mlodea/ -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Apr 15 10:16:36 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 15 Apr 2013 09:16:36 -0600 Subject: [Nek5000-users] Integrate Temperature Over Volume In-Reply-To: References: Message-ID: Hello Neks. I would like to integrate the temperature field over the volume of my mesh, and normalize that by the total volume. I'm hoping this can be accomplished with a line or two in the .usr file. Can you please either point me toward the correct example problem in the repo, or let me know what to enter. Thanks! --Mike From nek5000-users at lists.mcs.anl.gov Mon Apr 15 11:12:05 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 15 Apr 2013 11:12:05 -0500 (CDT) Subject: [Nek5000-users] Integrate Temperature Over Volume In-Reply-To: References: Message-ID: Hi Mike, The standard way to do this in userchk is (say): n = nx1*ny1*nz1*nelt tbar = glsc2(t,bm1,n)/voltm1 write(6,1) istep,time,tbar 1 format(i9,1p2e12.5,' tbar') Some Background: Note that tbar is computed as follows: / t_integral := | 1*t dV / = 1^T B t where the second line is the discrete inner product of the "1" vector with Temperature (t) using the diagonal mass matrix B. In nek, this is expressed as sum_i bm1_i t_i, which is evaluated using the global-scalar-product, with 2 arguments, glsc2() function shown above. glsc2() incorporates the requisite global communication to gather contributions from all processors. (Hence, it must also be called by all processors; a local variant is vlsc2().) Finally, voltm1 is just the volume of the temperature domain on mesh 1 (i.e., on the velocity/temperature mesh, no on the pressure mesh). Note that voltm1 = glsum(bm1,n) Cheers, Paul On Mon, 15 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > Hello Neks. > > I would like to integrate the temperature field over the volume of my mesh, and normalize that by the total volume. > > I'm hoping this can be accomplished with a line or two in the .usr file. > > Can you please either point me toward the correct example problem in the repo, or let me know what to enter. > > Thanks! > --Mike > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > From nek5000-users at lists.mcs.anl.gov Mon Apr 15 11:14:11 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 15 Apr 2013 11:14:11 -0500 (CDT) Subject: [Nek5000-users] Question - SIZE Parameters And Compilation Errors In-Reply-To: References: Message-ID: Hi Murph, A couple of options here: 1) Try recompiling with G = "-mcmodel=medium" in makenek. 2) Use more than one core so that lelt can remain the same but lelg can be as large as lelt*P, where P is the number of cores you'd like to use. Raw email is the standard interface as far as I know... hth, Paul On Mon, 15 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > Hi All, > > Trying to run some initial test cases of a problems I'm working on: > - Currently on single node (SMP Linux) > - MPI turned off (IFMPI="false") > - DOF at 5 (lx1=6, lxd=9, lelx=8) > > If I run with lp=1,lelt=2160, lelg=2160, for a 12x12x15 mesh, it compiles and > runs fine. > > If I up it to lp=1,lelt=4500,lelg=4500, for a 15x15x20 mesh, it crashes the > compile with: > > obj/drive2.o: In function `plan4_vol_': > drive2.f:(.text+0x82f): relocation truncated to fit: R_X86_64_32S against > symbol `scrvh_' defined in COMMON section in obj/hsmg.o > drive2.f:(.text+0x837): relocation truncated to fit: R_X86_64_32S against > symbol `scrvh_' defined in COMMON section in obj/hsmg.o > obj/drive2.o: In function `plan3_vol_': > drive2.f:(.text+0xac5): relocation truncated to fit: R_X86_64_32S against > symbol `scrvh_' defined in COMMON section in obj/hsmg.o > drive2.f:(.text+0xacd): relocation truncated to fit: R_X86_64_32S against > symbol `scrvh_' defined in COMMON section in obj/hsmg.o > obj/drive2.o: In function `plan2_vol_': > drive2.f:(.text+0xeb8): relocation truncated to fit: R_X86_64_32S against > symbol `scrvh_' defined in COMMON section in obj/hsmg.o > drive2.f:(.text+0xec0): relocation truncated to fit: R_X86_64_32S against > symbol `scrvh_' defined in COMMON section in obj/hsmg.o > obj/drive2.o: In function `dofcnt_': > drive2.f:(.text+0x170f): relocation truncated to fit: R_X86_64_PC32 against > symbol `scrns_' defined in COMMON section in obj/convect.o > obj/drive2.o: In function `settime_': > drive2.f:(.text+0x7cc9): relocation truncated to fit: R_X86_64_PC32 against > symbol `cprint_' defined in COMMON section in obj/navier1.o > drive2.f:(.text+0x7cef): relocation truncated to fit: R_X86_64_PC32 against > symbol `cprint_' defined in COMMON section in obj/navier1.o > drive2.f:(.text+0x7d02): relocation truncated to fit: R_X86_64_PC32 against > symbol `cprint_' defined in COMMON section in obj/navier1.o > obj/plan4.o: In function `crespsp_': > plan4.f:(.text+0xb23): additional relocation overflows omitted from the > output > collect2: ld returned 1 exit status > make: *** [nek5000] Error 1 > > Just wondering if this is due to a true limit with the platform/word size, or > if I'm screwing something up... > > Thanks, > Murph > > PS - Is there any other interface to this list other than raw email? > > > > -- > > *Murphy Leo O'Dea, PE"Murph"* > > PhD Student, Mechanical Engineering > > /School Of Engineering And Computer Science/ > > /Oakland University/ > > mlodea at oakland.edu > > http://www.secs.oakland.edu/~mlodea/ > > From nek5000-users at lists.mcs.anl.gov Mon Apr 15 11:14:41 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 15 Apr 2013 11:14:41 -0500 (CDT) Subject: [Nek5000-users] Integrate Temperature Over Volume In-Reply-To: Message-ID: Hi Mike, You could use n = nx1*ny1*nz1*nelt tavg = glsc2(t,bm1,n)/voltm1 Best. Aleks ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Monday, April 15, 2013 10:16:36 AM Subject: [Nek5000-users] Integrate Temperature Over Volume Hello Neks. I would like to integrate the temperature field over the volume of my mesh, and normalize that by the total volume. I'm hoping this can be accomplished with a line or two in the .usr file. Can you please either point me toward the correct example problem in the repo, or let me know what to enter. Thanks! --Mike _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Mon Apr 15 13:21:57 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 15 Apr 2013 13:21:57 -0500 Subject: [Nek5000-users] Question - SIZE Parameters And Compilation Errors In-Reply-To: References: Message-ID: Hi, You can also search the user-list email archives to see if others have had similar problems: http://lists.mcs.anl.gov/pipermail/nek5000-users/ Best, Katie On Mon, Apr 15, 2013 at 11:14 AM, wrote: > > Hi Murph, > > A couple of options here: > > 1) Try recompiling with G = "-mcmodel=medium" in makenek. > > 2) Use more than one core so that lelt can remain the same > but lelg can be as large as lelt*P, where P is the number > of cores you'd like to use. > > Raw email is the standard interface as far as I know... > > hth, > > Paul > > > On Mon, 15 Apr 2013, nek5000-users at lists.mcs.anl.**govwrote: > > Hi All, >> >> Trying to run some initial test cases of a problems I'm working on: >> - Currently on single node (SMP Linux) >> - MPI turned off (IFMPI="false") >> - DOF at 5 (lx1=6, lxd=9, lelx=8) >> >> If I run with lp=1,lelt=2160, lelg=2160, for a 12x12x15 mesh, it compiles >> and runs fine. >> >> If I up it to lp=1,lelt=4500,lelg=4500, for a 15x15x20 mesh, it crashes >> the compile with: >> >> obj/drive2.o: In function `plan4_vol_': >> drive2.f:(.text+0x82f): relocation truncated to fit: R_X86_64_32S against >> symbol `scrvh_' defined in COMMON section in obj/hsmg.o >> drive2.f:(.text+0x837): relocation truncated to fit: R_X86_64_32S against >> symbol `scrvh_' defined in COMMON section in obj/hsmg.o >> obj/drive2.o: In function `plan3_vol_': >> drive2.f:(.text+0xac5): relocation truncated to fit: R_X86_64_32S against >> symbol `scrvh_' defined in COMMON section in obj/hsmg.o >> drive2.f:(.text+0xacd): relocation truncated to fit: R_X86_64_32S against >> symbol `scrvh_' defined in COMMON section in obj/hsmg.o >> obj/drive2.o: In function `plan2_vol_': >> drive2.f:(.text+0xeb8): relocation truncated to fit: R_X86_64_32S against >> symbol `scrvh_' defined in COMMON section in obj/hsmg.o >> drive2.f:(.text+0xec0): relocation truncated to fit: R_X86_64_32S against >> symbol `scrvh_' defined in COMMON section in obj/hsmg.o >> obj/drive2.o: In function `dofcnt_': >> drive2.f:(.text+0x170f): relocation truncated to fit: R_X86_64_PC32 >> against symbol `scrns_' defined in COMMON section in obj/convect.o >> obj/drive2.o: In function `settime_': >> drive2.f:(.text+0x7cc9): relocation truncated to fit: R_X86_64_PC32 >> against symbol `cprint_' defined in COMMON section in obj/navier1.o >> drive2.f:(.text+0x7cef): relocation truncated to fit: R_X86_64_PC32 >> against symbol `cprint_' defined in COMMON section in obj/navier1.o >> drive2.f:(.text+0x7d02): relocation truncated to fit: R_X86_64_PC32 >> against symbol `cprint_' defined in COMMON section in obj/navier1.o >> obj/plan4.o: In function `crespsp_': >> plan4.f:(.text+0xb23): additional relocation overflows omitted from the >> output >> collect2: ld returned 1 exit status >> make: *** [nek5000] Error 1 >> >> Just wondering if this is due to a true limit with the platform/word >> size, or if I'm screwing something up... >> >> Thanks, >> Murph >> >> PS - Is there any other interface to this list other than raw email? >> >> >> >> -- >> >> *Murphy Leo O'Dea, PE"Murph"* >> >> PhD Student, Mechanical Engineering >> >> /School Of Engineering And Computer Science/ >> >> /Oakland University/ >> >> mlodea at oakland.edu >> >> http://www.secs.oakland.edu/~**mlodea/ >> >> >> ______________________________**_________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.**gov > https://lists.mcs.anl.gov/**mailman/listinfo/nek5000-users > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Apr 15 16:21:08 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 15 Apr 2013 16:21:08 -0500 (CDT) Subject: [Nek5000-users] postnek on a mac In-Reply-To: References: Message-ID: Hi All, I'm new to the macs and am trying to compile postnek and prenek. I'm getting the same error reported by Mark (below). Has anyone had any luck with this? Thanks! Paul On Fri, 22 Jun 2012, nek5000-users at lists.mcs.anl.gov wrote: > > Hello all, > > I have downloaded the most recent version of the software and am > trying to compile postnek on a mac running lion with gfortran and gcc. > > The command: > > ./maketools postnek > > Gives the following error: > > Warning: Deleted feature: ASSIGN statement at (1) > gcc -mcmodel=medium -c -O2 xdriver.c > xdriver.c: In function ?TypeChar?: > xdriver.c:1790: warning: comparison is always true due to limited range of data type > xdriver.c: In function ?error?: > xdriver.c:1910: warning: format not a string literal and no format arguments > xdriver.c:1910: warning: format not a string literal and no format arguments > gfortran -mcmodel=medium -c scrdmp.f > gfortran -mcmodel=medium -c coef.f > gfortran -mcmodel=medium -c postnek7.f > gfortran -mcmodel=medium -c speclib.f > gfortran -mcmodel=medium -c mxm.f > gcc -mcmodel=medium -c -DUNDERSCORE ../../nek/byte.c > gfortran -mcmodel=medium -c ssyev.f > gfortran -mcmodel=medium -c iolib.f > gfortran -mcmodel=medium -o /Users/mrp/bin/postx postnek.o postnek2.o postnek3.o postnek5.o postnek6.o tsort.o postnek8.o postnek9.o plot.o getfld.o legend.o userf.o revert.o trap.o animate.o genxyz.o screen.o g3d.o subs.o xinterface.o locglob.o postnek5a.o blas.o xdriver.o scrdmp.o coef.o postnek7.o speclib.o mxm.o byte.o ssyev.o iolib.o -L/usr/X11R6/lib -lX11 -lm > ld: in xdriver.o, in section __TEXT,__text reloc 3: length != 2 and X86_64_RELOC_GOT not supported for architecture x86_64 > collect2: ld returned 1 exit status > make[1]: *** [postx] Error 1 > make: *** [all] Error 1 > > Any ideas? > > Thanks, > > Mark > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > From nek5000-users at lists.mcs.anl.gov Mon Apr 15 16:32:03 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 15 Apr 2013 16:32:03 -0500 (CDT) Subject: [Nek5000-users] postnek on a mac In-Reply-To: Message-ID: Hello Paul, I would try recompiling clean with BIGMEM="false" in maketools and smaller upper element declarations of NELM in tools/*nek/basics.inc Best. Aleks ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Monday, April 15, 2013 4:21:08 PM Subject: Re: [Nek5000-users] postnek on a mac Hi All, I'm new to the macs and am trying to compile postnek and prenek. I'm getting the same error reported by Mark (below). Has anyone had any luck with this? Thanks! Paul On Fri, 22 Jun 2012, nek5000-users at lists.mcs.anl.gov wrote: > > Hello all, > > I have downloaded the most recent version of the software and am > trying to compile postnek on a mac running lion with gfortran and gcc. > > The command: > > ./maketools postnek > > Gives the following error: > > Warning: Deleted feature: ASSIGN statement at (1) > gcc -mcmodel=medium -c -O2 xdriver.c > xdriver.c: In function ?TypeChar?: > xdriver.c:1790: warning: comparison is always true due to limited range of data type > xdriver.c: In function ?error?: > xdriver.c:1910: warning: format not a string literal and no format arguments > xdriver.c:1910: warning: format not a string literal and no format arguments > gfortran -mcmodel=medium -c scrdmp.f > gfortran -mcmodel=medium -c coef.f > gfortran -mcmodel=medium -c postnek7.f > gfortran -mcmodel=medium -c speclib.f > gfortran -mcmodel=medium -c mxm.f > gcc -mcmodel=medium -c -DUNDERSCORE ../../nek/byte.c > gfortran -mcmodel=medium -c ssyev.f > gfortran -mcmodel=medium -c iolib.f > gfortran -mcmodel=medium -o /Users/mrp/bin/postx postnek.o postnek2.o postnek3.o postnek5.o postnek6.o tsort.o postnek8.o postnek9.o plot.o getfld.o legend.o userf.o revert.o trap.o animate.o genxyz.o screen.o g3d.o subs.o xinterface.o locglob.o postnek5a.o blas.o xdriver.o scrdmp.o coef.o postnek7.o speclib.o mxm.o byte.o ssyev.o iolib.o -L/usr/X11R6/lib -lX11 -lm > ld: in xdriver.o, in section __TEXT,__text reloc 3: length != 2 and X86_64_RELOC_GOT not supported for architecture x86_64 > collect2: ld returned 1 exit status > make[1]: *** [postx] Error 1 > make: *** [all] Error 1 > > Any ideas? > > Thanks, > > Mark > > _______________________________________________ > 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 From nek5000-users at lists.mcs.anl.gov Mon Apr 15 16:39:17 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 15 Apr 2013 16:39:17 -0500 Subject: [Nek5000-users] postnek on a mac In-Reply-To: References: Message-ID: Hi Paul, I've been running postx/prex on Snow Leopard, but I had similar issues when I first attempted to get it to compile. If I remember correctly, it has to do with which compilers you use, since Apple doesn't ship a native their own Fortran compiler. The solution that worked for me was downloading the stripped-down version of Apple's Xcode , installing Homebrew, and grabbing gfortran from Homebrew. The following link is a starter (http://kennethreitz.org/experiments/xcode-gcc-and-homebrew). This may not be the best way to solve the issue, but I know this is how I finally got things to work. Please, somebody chime in if I'm wrong. Justin On Mon, Apr 15, 2013 at 4:32 PM, wrote: > Hello Paul, > > I would try recompiling clean with > > BIGMEM="false" > > in maketools and smaller upper element declarations of NELM in tools/*nek/basics.inc > > Best. > Aleks > > > ----- Original Message ----- > From: nek5000-users at lists.mcs.anl.gov > To: nek5000-users at lists.mcs.anl.gov > Sent: Monday, April 15, 2013 4:21:08 PM > Subject: Re: [Nek5000-users] postnek on a mac > > > Hi All, > > I'm new to the macs and am trying to compile postnek and prenek. > > I'm getting the same error reported by Mark (below). > > Has anyone had any luck with this? > > Thanks! > > Paul > > > On Fri, 22 Jun 2012, nek5000-users at lists.mcs.anl.gov wrote: > >> >> Hello all, >> >> I have downloaded the most recent version of the software and am >> trying to compile postnek on a mac running lion with gfortran and gcc. >> >> The command: >> >> ./maketools postnek >> >> Gives the following error: >> >> Warning: Deleted feature: ASSIGN statement at (1) >> gcc -mcmodel=medium -c -O2 xdriver.c >> xdriver.c: In function ?TypeChar?: >> xdriver.c:1790: warning: comparison is always true due to limited range of data type >> xdriver.c: In function ?error?: >> xdriver.c:1910: warning: format not a string literal and no format arguments >> xdriver.c:1910: warning: format not a string literal and no format arguments >> gfortran -mcmodel=medium -c scrdmp.f >> gfortran -mcmodel=medium -c coef.f >> gfortran -mcmodel=medium -c postnek7.f >> gfortran -mcmodel=medium -c speclib.f >> gfortran -mcmodel=medium -c mxm.f >> gcc -mcmodel=medium -c -DUNDERSCORE ../../nek/byte.c >> gfortran -mcmodel=medium -c ssyev.f >> gfortran -mcmodel=medium -c iolib.f >> gfortran -mcmodel=medium -o /Users/mrp/bin/postx postnek.o postnek2.o postnek3.o postnek5.o postnek6.o tsort.o postnek8.o postnek9.o plot.o getfld.o legend.o userf.o revert.o trap.o animate.o genxyz.o screen.o g3d.o subs.o xinterface.o locglob.o postnek5a.o blas.o xdriver.o scrdmp.o coef.o postnek7.o speclib.o mxm.o byte.o ssyev.o iolib.o -L/usr/X11R6/lib -lX11 -lm >> ld: in xdriver.o, in section __TEXT,__text reloc 3: length != 2 and X86_64_RELOC_GOT not supported for architecture x86_64 >> collect2: ld returned 1 exit status >> make[1]: *** [postx] Error 1 >> make: *** [all] Error 1 >> >> Any ideas? >> >> Thanks, >> >> Mark >> >> _______________________________________________ >> 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 > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Mon Apr 15 20:59:16 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 15 Apr 2013 20:59:16 -0500 (CDT) Subject: [Nek5000-users] postnek on a mac In-Reply-To: References: Message-ID: Hi Aleks, Thanks -- this worked both for postnek and prenek. Paul On Mon, 15 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > Hello Paul, I would try recompiling clean with BIGMEM="false" in maketools and smaller upper element declarations of NELM in tools/*nek/basics.inc Best. Aleks ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Monday, April 15, 2013 4:21:08 PM Subject: Re: [Nek5000-users] postnek on a mac Hi All, I'm new to the macs and am trying to compile postnek and prenek. I'm getting the same error reported by Mark (below). Has anyone had any luck with this? Thanks! Paul On Fri, 22 Jun 2012, nek5000-users at lists.mcs.anl.gov wrote: > > Hello all, > > I have downloaded the most recent version of the software and am > trying to compile postnek on a mac running lion with gfortran and gcc. > > The command: > > ./maketools postnek > > Gives the following error: > > Warning: Deleted feature: ASSIGN statement at (1) > gcc -mcmodel=medium -c -O2 xdriver.c > xdriver.c: In function ?TypeChar?: > xdriver.c:1790: warning: comparison is always true due to limited range of data type > xdriver.c: In function ?error?: > xdriver.c:1910: warning: format not a string literal and no format arguments > xdriver.c:1910: warning: format not a string literal and no format arguments > gfortran -mcmodel=medium -c scrdmp.f > gfortran -mcmodel=medium -c coef.f > gfortran -mcmodel=medium -c postnek7.f > gfortran -mcmodel=medium -c speclib.f > gfortran -mcmodel=medium -c mxm.f > gcc -mcmodel=medium -c -DUNDERSCORE ../../nek/byte.c > gfortran -mcmodel=medium -c ssyev.f > gfortran -mcmodel=medium -c iolib.f > gfortran -mcmodel=medium -o /Users/mrp/bin/postx postnek.o postnek2.o postnek3.o postnek5.o postnek6.o tsort.o postnek8.o postnek9.o plot.o getfld.o legend.o userf.o revert.o trap.o animate.o genxyz.o screen.o g3d.o subs.o xinterface.o locglob.o postnek5a.o blas.o xdriver.o scrdmp.o coef.o postnek7.o speclib.o mxm.o byte.o ssyev.o iolib.o -L/usr/X11R6/lib -lX11 -lm > ld: in xdriver.o, in section __TEXT,__text reloc 3: length != 2 and X86_64_RELOC_GOT not supported for architecture x86_64 > collect2: ld returned 1 exit status > make[1]: *** [postx] Error 1 > make: *** [all] Error 1 > > Any ideas? > > Thanks, > > Mark > > _______________________________________________ > 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 _______________________________________________ Nek5000-users mailing list Nek5000-users at lists.mcs.anl.gov https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Tue Apr 16 16:37:33 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 16 Apr 2013 16:37:33 -0500 (CDT) Subject: [Nek5000-users] distance function In-Reply-To: References: Message-ID: Dear Neks, I've just added two new functions to navier5.f that might be of interest. One is a low-pass cut-off filter. Basically, ifld = 1 call cut_off_filter(u,mx,ifld) ! mx=max saved mode will project any field u() onto a polynomial of degree mx-1 (e.g., mx=2 would imply projecting u onto piecewise bilinear or trilinear in 2D or 3D). The supplemental argument ifld tells the filter to march over the number of elements associated with field "ifld" (e.g., ifld=1 implies e=1,nelv; ifld=2 implies e=1,nelt, which could be different in the case of conjugate heat transfer). The second routine is a cheap (pseudo) distance function. .For those of you who have been using this function for the past few days in the .usr file please note that there was a minor bug in that version that is now fixed in the repo. Also, the interface to the repo-version is slightly different: call cheap_dist(d,ifld) returns an approximate distance function in the array d(), which is associated with field ifld, as above (i.e. ifld=1 if you wish to march over e=1,nelv, etc.) .The reason this is a pseudo-distance function is that it actually returns the minimum tally of the distances that separate a chain of neighboring vertices. If the chain is on a line normal to the wall then this is a true distrance, otherwise it's an approximation. The choice of this distance function over more precise variants was motivated by the fact that cheap_dist() will return the proper result when you have periodic boundary conditions, whereas other approaches do not readily meet this criterion. For many applications, cheap_dist is sufficient. Please let me know if there are any issues with these routines. Cheers, Paul From nek5000-users at lists.mcs.anl.gov Thu Apr 18 20:17:04 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Thu, 18 Apr 2013 21:17:04 -0400 Subject: [Nek5000-users] Nek5000-users Digest, Vol 50, Issue 6 In-Reply-To: References: Message-ID: Paul, Thanks for your help; both of these worked... Last real machine I ran on was a Cray YMP, so I'm learning a lot about MPI. I ended up installing MPI on an SMP machine; compiled NEK for MPI with the wrappers; and ran it on multiple cores. I assume not quite as good as OpenMP, but worked great for my purposes... Thanks, Murph On 4/16/2013 1:00 PM, nek5000-users-request at lists.mcs.anl.gov wrote: > On Mon, Apr 15, 2013 at 11:14 AM, > wrote: > > > Hi Murph, > > A couple of options here: > > 1) Try recompiling with G = "-mcmodel=medium" in makenek. > > 2) Use more than one core so that lelt can remain the same > but lelg can be as large as lelt*P, where P is the number > of cores you'd like to use. > > Raw email is the standard interface as far as I know... > > hth, > > Paul > > > On Mon, 15 Apr 2013, nek5000-users at lists.mcs.anl.gov > wrote: > > Hi All, > > Trying to run some initial test cases of a problems I'm > working on: > - Currently on single node (SMP Linux) > - MPI turned off (IFMPI="false") > - DOF at 5 (lx1=6, lxd=9, lelx=8) > > If I run with lp=1,lelt=2160, lelg=2160, for a 12x12x15 mesh, > it compiles and runs fine. > > If I up it to lp=1,lelt=4500,lelg=4500, for a 15x15x20 mesh, > it crashes the compile with: > > obj/drive2.o: In function `plan4_vol_': > drive2.f:(.text+0x82f): relocation truncated to fit: > R_X86_64_32S against symbol `scrvh_' defined in COMMON section > in obj/hsmg.o > drive2.f:(.text+0x837): relocation truncated to fit: > R_X86_64_32S against symbol `scrvh_' defined in COMMON section > in obj/hsmg.o > obj/drive2.o: In function `plan3_vol_': > drive2.f:(.text+0xac5): relocation truncated to fit: > R_X86_64_32S against symbol `scrvh_' defined in COMMON section > in obj/hsmg.o > drive2.f:(.text+0xacd): relocation truncated to fit: > R_X86_64_32S against symbol `scrvh_' defined in COMMON section > in obj/hsmg.o > obj/drive2.o: In function `plan2_vol_': > drive2.f:(.text+0xeb8): relocation truncated to fit: > R_X86_64_32S against symbol `scrvh_' defined in COMMON section > in obj/hsmg.o > drive2.f:(.text+0xec0): relocation truncated to fit: > R_X86_64_32S against symbol `scrvh_' defined in COMMON section > in obj/hsmg.o > obj/drive2.o: In function `dofcnt_': > drive2.f:(.text+0x170f): relocation truncated to fit: > R_X86_64_PC32 against symbol `scrns_' defined in COMMON > section in obj/convect.o > obj/drive2.o: In function `settime_': > drive2.f:(.text+0x7cc9): relocation truncated to fit: > R_X86_64_PC32 against symbol `cprint_' defined in COMMON > section in obj/navier1.o > drive2.f:(.text+0x7cef): relocation truncated to fit: > R_X86_64_PC32 against symbol `cprint_' defined in COMMON > section in obj/navier1.o > drive2.f:(.text+0x7d02): relocation truncated to fit: > R_X86_64_PC32 against symbol `cprint_' defined in COMMON > section in obj/navier1.o > obj/plan4.o: In function `crespsp_': > plan4.f:(.text+0xb23): additional relocation overflows omitted > from the output > collect2: ld returned 1 exit status > make: *** [nek5000] Error 1 > > Just wondering if this is due to a true limit with the > platform/word size, or if I'm screwing something up... > > Thanks, > Murph > > PS - Is there any other interface to this list other than raw > email? > > > > -- > > *Murphy Leo O'Dea, PE"Murph"* > > PhD Student, Mechanical Engineering > > /School Of Engineering And Computer Science/ > > /Oakland University/ > > mlodea at oakland.edu > > > > http://www.secs.oakland.edu/~mlodea/ > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Mon Apr 22 08:30:12 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Apr 2013 08:30:12 -0500 Subject: [Nek5000-users] Boundary condition loop Message-ID: Neks, I need to loop over all element faces that are boundary faces (doesn't matter which type of boundary condition). I know one way to do this would be do e=1,nelv do f=1,2*ndim if (cbc(f,e,1).ne.'E '.and.cbc(f,e,1).ne.'P ') then ! Do computations endif enddo enddo Is there a more elegant solution for this, maybe some array that points to element/face combinations that are boundary conditions? Thanks, Josh From nek5000-users at lists.mcs.anl.gov Mon Apr 22 08:39:22 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Apr 2013 08:39:22 -0500 (CDT) Subject: [Nek5000-users] Boundary condition loop In-Reply-To: References: Message-ID: Hi Josh, Your approach is the only one I know of. Note that there are many routines around that will access the points that are on a given face pair (f,e), so you might want to check the structure of those for your "Do computations" block. Don't forget "integer e,f" :) Cheers, Paul On Mon, 22 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > Neks, > > I need to loop over all element faces that are boundary faces (doesn't matter which type of boundary condition). I know one way to do this would be > > do e=1,nelv > do f=1,2*ndim > if (cbc(f,e,1).ne.'E '.and.cbc(f,e,1).ne.'P ') then > ! Do computations > endif > enddo > enddo > > Is there a more elegant solution for this, maybe some array that points to element/face combinations that are boundary conditions? > > Thanks, > > Josh > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > From nek5000-users at lists.mcs.anl.gov Mon Apr 22 09:03:38 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Apr 2013 09:03:38 -0500 (CDT) Subject: [Nek5000-users] Boundary condition loop In-Reply-To: Message-ID: Hi Josh, one additional comment that you can loop with if statements only once to fill up your 1-D mask/pointer array which you can later loop over instead. Let me know if you need more details. Best. Aleks ----- Original Message ----- From: nek5000-users at lists.mcs.anl.gov To: nek5000-users at lists.mcs.anl.gov Sent: Monday, April 22, 2013 8:39:22 AM Subject: Re: [Nek5000-users] Boundary condition loop Hi Josh, Your approach is the only one I know of. Note that there are many routines around that will access the points that are on a given face pair (f,e), so you might want to check the structure of those for your "Do computations" block. Don't forget "integer e,f" :) Cheers, Paul On Mon, 22 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > Neks, > > I need to loop over all element faces that are boundary faces (doesn't matter which type of boundary condition). I know one way to do this would be > > do e=1,nelv > do f=1,2*ndim > if (cbc(f,e,1).ne.'E '.and.cbc(f,e,1).ne.'P ') then > ! Do computations > endif > enddo > enddo > > Is there a more elegant solution for this, maybe some array that points to element/face combinations that are boundary conditions? > > Thanks, > > Josh > _______________________________________________ > 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 From nek5000-users at lists.mcs.anl.gov Mon Apr 22 09:07:52 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Apr 2013 09:07:52 -0500 Subject: [Nek5000-users] Boundary condition loop In-Reply-To: References: Message-ID: Ah, no worries, after the 15th time of making that mistake, I always declare my integers :). Are you referring to the facexs and the facexv routines? Josh On Apr 22, 2013, at 8:39 AM, nek5000-users at lists.mcs.anl.gov wrote: > > > Hi Josh, > > Your approach is the only one I know of. > > Note that there are many routines around that will access > the points that are on a given face pair (f,e), so you might > want to check the structure of those for your "Do computations" > block. > > Don't forget "integer e,f" :) > > Cheers, Paul > > > On Mon, 22 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > >> Neks, >> >> I need to loop over all element faces that are boundary faces (doesn't matter which type of boundary condition). I know one way to do this would be >> >> do e=1,nelv >> do f=1,2*ndim >> if (cbc(f,e,1).ne.'E '.and.cbc(f,e,1).ne.'P ') then >> ! Do computations >> endif >> enddo >> enddo >> >> Is there a more elegant solution for this, maybe some array that points to element/face combinations that are boundary conditions? >> >> Thanks, >> >> Josh >> _______________________________________________ >> 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 From nek5000-users at lists.mcs.anl.gov Mon Apr 22 09:09:03 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Apr 2013 09:09:03 -0500 Subject: [Nek5000-users] Boundary condition loop In-Reply-To: References: Message-ID: Hi Aleks, I was actually thinking about doing just that, was just curious if such an array already existed. Thanks! Josh On Apr 22, 2013, at 9:03 AM, nek5000-users at lists.mcs.anl.gov wrote: > Hi Josh, > > one additional comment that you can loop with if statements only once to fill up your 1-D mask/pointer array which you can later loop over instead. > > Let me know if you need more details. > > Best. > Aleks > > > ----- Original Message ----- > From: nek5000-users at lists.mcs.anl.gov > To: nek5000-users at lists.mcs.anl.gov > Sent: Monday, April 22, 2013 8:39:22 AM > Subject: Re: [Nek5000-users] Boundary condition loop > > > > Hi Josh, > > Your approach is the only one I know of. > > Note that there are many routines around that will access > the points that are on a given face pair (f,e), so you might > want to check the structure of those for your "Do computations" > block. > > Don't forget "integer e,f" :) > > Cheers, Paul > > > On Mon, 22 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > >> Neks, >> >> I need to loop over all element faces that are boundary faces (doesn't matter which type of boundary condition). I know one way to do this would be >> >> do e=1,nelv >> do f=1,2*ndim >> if (cbc(f,e,1).ne.'E '.and.cbc(f,e,1).ne.'P ') then >> ! Do computations >> endif >> enddo >> enddo >> >> Is there a more elegant solution for this, maybe some array that points to element/face combinations that are boundary conditions? >> >> Thanks, >> >> Josh >> _______________________________________________ >> 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 > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users From nek5000-users at lists.mcs.anl.gov Mon Apr 22 09:29:27 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Apr 2013 09:29:27 -0500 Subject: [Nek5000-users] wgradm1 Message-ID: Hi Neks, I had a quick clarification question about the subroutine wgradm1. The comments suggest it is a weak gradient, but the formulation itself looks like it is computing a strong gradient (and multiplying by mass), i.e., grad(u) = W*G*D*u (W being the weights, G being the appropriate metrics, and D begin the normal derivative operator). Is this correct, or am I missing something? Thanks, Josh From nek5000-users at lists.mcs.anl.gov Mon Apr 22 09:36:34 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Mon, 22 Apr 2013 09:36:34 -0500 (CDT) Subject: [Nek5000-users] wgradm1 In-Reply-To: References: Message-ID: Hi Josh, yes - it's simply the gradient with weights, so that it can be incorporated into the rhs of a weak formulation. Paul On Mon, 22 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > Hi Neks, > > I had a quick clarification question about the subroutine wgradm1. The > comments suggest it is a weak gradient, but the formulation itself looks like > it is computing a strong gradient (and multiplying by mass), i.e., > > grad(u) = W*G*D*u > > (W being the weights, G being the appropriate metrics, and D begin the normal > derivative operator). Is this correct, or am I missing something? > > Thanks, > > Josh _______________________________________________ Nek5000-users mailing > list Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > From nek5000-users at lists.mcs.anl.gov Wed Apr 24 00:59:00 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 24 Apr 2013 11:29:00 +0530 Subject: [Nek5000-users] Benard Convection heated from below with rotation Message-ID: Hi, I am new to nek5000, I want to simulate Benard convection problem heated from below with rotation and later with magnetic field and different BCs for my research. I ran the example provided with the code. Please let me know how to proceed with my problem while using nek5000. Best Regards, venkata -------------- next part -------------- An HTML attachment was scrubbed... URL: From nek5000-users at lists.mcs.anl.gov Wed Apr 24 05:18:24 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 24 Apr 2013 12:18:24 +0200 Subject: [Nek5000-users] prenek question Message-ID: Hi Neks. I have a question about prenek. I use it to insert a cap in the mesh connecting pipe with the box, but there is an exit (call prexit) in the new version of subroutine sc_make_sphere_cap. Why is it so? Is the new version incomplete? Moreover the resolution is changed by increasing the number of radii to 7, while I'm using 3. Are you going to add additional switch to choose resolution? Regards Adam From nek5000-users at lists.mcs.anl.gov Wed Apr 24 06:05:14 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Wed, 24 Apr 2013 06:05:14 -0500 (CDT) Subject: [Nek5000-users] prenek question In-Reply-To: References: Message-ID: Hi Adam, My apologies - I didn't realize that anyone was using sc. (Glad to hear that it _is_ being used!) I think I've now restored the old functionality --- It was hardwired with a particular resolution just until it gets a bit more usage and we understand the typical use cases for the feature. Since prenek/postnek do not have standard .usr interfaces, it's often easiest just to modify the code for the application at hand, as the number of free parameters is otherwise innumerable. Please let me know if this revision works for you. Best, Paul On Wed, 24 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > Hi Neks. > > I have a question about prenek. I use it to insert a cap in the mesh > connecting pipe with the box, but there is an exit (call prexit) in the new > version of subroutine sc_make_sphere_cap. Why is it so? Is the new version > incomplete? Moreover the resolution is changed by increasing the number of > radii to 7, while I'm using 3. Are you going to add additional switch to > choose resolution? > Regards > Adam > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > From nek5000-users at lists.mcs.anl.gov Tue Apr 30 02:11:17 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 30 Apr 2013 09:11:17 +0200 Subject: [Nek5000-users] prenek question In-Reply-To: References: Message-ID: Hi It seems to work now, but I found 2 typos in iolib_no_graph.f causing pretex crashes: 1) line 153; WRITE(S,'(3G113.9)',ERR=1) changed to WRITE(S,'(3G13.9)',ERR=1) 2) line 184; WRITE(SS(NC+1:NC+5), '(I13)', changed to WRITE(SS(NC+1:NC+13),'(I13)', Thank you for help. Best regards Adam On 24/04/13 13:05, nek5000-users at lists.mcs.anl.gov wrote: > > Hi Adam, > > My apologies - I didn't realize that anyone was using sc. > (Glad to hear that it _is_ being used!) > > I think I've now restored the old functionality --- It was > hardwired with a particular resolution just until it gets > a bit more usage and we understand the typical use cases > for the feature. Since prenek/postnek do not have standard > .usr interfaces, it's often easiest just to modify the code > for the application at hand, as the number of free parameters > is otherwise innumerable. > > Please let me know if this revision works for you. > > Best, Paul > > > On Wed, 24 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > >> Hi Neks. >> >> I have a question about prenek. I use it to insert a cap in the mesh >> connecting pipe with the box, but there is an exit (call prexit) in >> the new version of subroutine sc_make_sphere_cap. Why is it so? Is >> the new version incomplete? Moreover the resolution is changed by >> increasing the number of radii to 7, while I'm using 3. Are you going >> to add additional switch to choose resolution? >> Regards >> Adam >> _______________________________________________ >> 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 From nek5000-users at lists.mcs.anl.gov Tue Apr 30 11:32:39 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 30 Apr 2013 11:32:39 -0500 (CDT) Subject: [Nek5000-users] prenek question In-Reply-To: References: Message-ID: Thanks Adam! I've updated iolib_no_graph.f Paul On Tue, 30 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: > Hi > > It seems to work now, but I found 2 typos in iolib_no_graph.f causing pretex > crashes: > 1) line 153; WRITE(S,'(3G113.9)',ERR=1) changed to WRITE(S,'(3G13.9)',ERR=1) > 2) line 184; WRITE(SS(NC+1:NC+5), '(I13)', changed to > WRITE(SS(NC+1:NC+13),'(I13)', > Thank you for help. > Best regards > Adam > > On 24/04/13 13:05, nek5000-users at lists.mcs.anl.gov wrote: >> >> Hi Adam, >> >> My apologies - I didn't realize that anyone was using sc. >> (Glad to hear that it _is_ being used!) >> >> I think I've now restored the old functionality --- It was >> hardwired with a particular resolution just until it gets >> a bit more usage and we understand the typical use cases >> for the feature. Since prenek/postnek do not have standard >> .usr interfaces, it's often easiest just to modify the code >> for the application at hand, as the number of free parameters >> is otherwise innumerable. >> >> Please let me know if this revision works for you. >> >> Best, Paul >> >> >> On Wed, 24 Apr 2013, nek5000-users at lists.mcs.anl.gov wrote: >> >>> Hi Neks. >>> >>> I have a question about prenek. I use it to insert a cap in the mesh >>> connecting pipe with the box, but there is an exit (call prexit) in the >>> new version of subroutine sc_make_sphere_cap. Why is it so? Is the new >>> version incomplete? Moreover the resolution is changed by increasing the >>> number of radii to 7, while I'm using 3. Are you going to add additional >>> switch to choose resolution? >>> Regards >>> Adam >>> _______________________________________________ >>> 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 > > _______________________________________________ > Nek5000-users mailing list > Nek5000-users at lists.mcs.anl.gov > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users > From nek5000-users at lists.mcs.anl.gov Tue Apr 30 14:42:04 2013 From: nek5000-users at lists.mcs.anl.gov (nek5000-users at lists.mcs.anl.gov) Date: Tue, 30 Apr 2013 15:42:04 -0400 Subject: [Nek5000-users] Print Statements Message-ID: Dear Nek5000 users, How can I change the frequency of print statements when I am running a simulation? Is this the IOCOMM parameter in the .rea file? Saumil Patel -------------- next part -------------- An HTML attachment was scrubbed... URL: