<div dir="ltr">Hi Ami,<div><br></div><div>As long as lhis times the number of processors is greater than your total number of points (as the comment in the SIZE file suggests), then the SIZE file shouldn't cause any issues.</div>
<div><br></div><div>How many total points are you trying to interpolate? If your <a href="http://hpts.in">hpts.in</a> file is not too large, could you post that (as well as the relevant portion of your .usr file where you call the hpts routine)</div>
<div><br></div><div>- Josh</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Dec 13, 2013 at 1:07 PM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Hi,<br><br></div>I have a problem with hpts() function to interpolating data. I was using hpts() function without any problem in userchk and by input it was giving me the interpolation, but after some days I don't what change I did, whatever I interpolate now it just gives me multiplication of certain points(like in one case is 107 points and in other it is 90), I think there is some problems with setting parameters on SIZE, below is my SIZE file for one of my cases which I can just get interpolation with multiplication of 107 points (if I have less than 107 points in <a href="http://hpts.in" target="_blank">hpts.in</a>, it just gives me nothing in htps.out!!):<br>
<br>C Dimension file to be included <br>C<br>C HCUBE array dimensions<br>C<br> parameter (ldim=3)<br> parameter (lx1=6,ly1=lx1,lz1=lx1,lelt=410,lelv=lelt)<br> parameter (lxd=10,lyd=lxd,lzd=lxd)<br>
parameter (lelx=1,lely=1,lelz=1)<br> <br> parameter (lzl=3 + 2*(ldim-3))<br> <br> parameter (lx2=lx1-2)<br> parameter (ly2=ly1-2)<br> parameter (lz2=lz1-2)<br> parameter (lx3=lx1)<br> parameter (ly3=ly1)<br>
parameter (lz3=lz1)<br><br> parameter (lp = 512)<br> parameter (lelg = 16000)<br>c<br>c parameter (lpelv=lelv,lpelt=lelt,lpert=3) ! perturbation<br>c parameter (lpx1=lx1,lpy1=ly1,lpz1=lz1) ! array sizes<br>
c parameter (lpx2=lx2,lpy2=ly2,lpz2=lz2)<br>c<br> parameter (lpelv=1,lpelt=1,lpert=1) ! perturbation<br> parameter (lpx1=1,lpy1=1,lpz1=1) ! array sizes<br> parameter (lpx2=1,lpy2=1,lpz2=1)<br>
c<br>c parameter (lbelv=lelv,lbelt=lelt) ! MHD<br>c parameter (lbx1=lx1,lby1=ly1,lbz1=lz1) ! array sizes<br>c parameter (lbx2=lx2,lby2=ly2,lbz2=lz2)<br>c<br> parameter (lbelv=1,lbelt=1) ! MHD<br>
parameter (lbx1=1,lby1=1,lbz1=1) ! array sizes<br> parameter (lbx2=1,lby2=1,lbz2=1)<br> <br>C LX1M=LX1 when there are moving meshes; =1 otherwise<br> parameter (lx1m=1,ly1m=1,lz1m=1)<br> parameter (ldimt= 2) ! 3 passive scalars + T<br>
parameter (ldimt1=ldimt+1)<br> parameter (ldimt3=ldimt+3)<br>c<br>c Note: In the new code, LELGEC should be about sqrt(LELG)<br>c<br> PARAMETER (LELGEC = 1)<br> PARAMETER (LXYZ2 = 1)<br> PARAMETER (LXZ21 = 1)<br>
<br> PARAMETER (LMAXV=LX1*LY1*LZ1*LELV)<br> PARAMETER (LMAXT=LX1*LY1*LZ1*LELT)<br> PARAMETER (LMAXP=LX2*LY2*LZ2*LELV)<br> PARAMETER (LXZ=LX1*LZ1)<br> PARAMETER (LORDER=3)<br> PARAMETER (MAXOBJ=4,MAXMBR=LELT*6)<br>
PARAMETER (lhis=200) ! # of pts a proc reads from <a href="http://hpts.in" target="_blank">hpts.in</a><br> ! Note: lhis*np > npoints in <a href="http://hpts.in" target="_blank">hpts.in</a><br>
C<br>C Common Block Dimensions<br>
C<br> PARAMETER (LCTMP0 =2*LX1*LY1*LZ1*LELT)<br> PARAMETER (LCTMP1 =4*LX1*LY1*LZ1*LELT)<br>C<br>C The parameter LVEC controls whether an additional 42 field arrays<br>C are required for Steady State Solutions. If you are not using<br>
C Steady State, it is recommended that LVEC=1.<br>C<br> PARAMETER (LVEC=1)<br>C<br>C Uzawa projection array dimensions<br>C<br> parameter (mxprev = 20)<br> parameter (lgmres = 30)<br>C<br>C Split projection array dimensions<br>
C<br> parameter(lmvec = 1)<br> parameter(lsvec = 1)<br> parameter(lstore=lmvec*lsvec)<br>c<br>c NONCONFORMING STUFF<br>c<br> parameter (maxmor = lelt)<br>C<br>C Array dimensions<br>C<br> COMMON/DIMN/NELV,NELT,NX1,NY1,NZ1,NX2,NY2,NZ2<br>
$,NX3,NY3,NZ3,NDIM,NFIELD,NPERT,NID<br> $,NXD,NYD,NZD<br><br>c automatically added by makenek<br> parameter(lxo = lx1) ! max output grid size (lxo>=lx1)<br><br>c automatically added by makenek<br> parameter(lpart = 1 ) ! max number of particles<br>
<br>c automatically added by makenek<br> integer ax1,ay1,az1,ax2,ay2,az2<br> parameter (ax1=lx1,ay1=ly1,az1=lz1,ax2=lx2,ay2=ly2,az2=lz2) ! running averages<br><br>c automatically added by makenek<br> parameter (lxs=1,lys=lxs,lzs=(lxs-1)*(ldim-2)+1) !New Pressure Preconditioner<br>
<br>c automatically added by makenek<br> parameter (lfdm=0) ! == 1 for fast diagonalization method<br><br></div>Is anyone has idea what is the problem in this case?!<br><br></div>Thanks,<br>Ami<br><div><div><div><br>
</div></div></div></div>
<br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Josh Camp<br><br>"All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke
</div>