[Nek5000-users] problem with interpolation by hpts

[nek5000-users at lists.mcs.anl.gov]

Hi Ami,

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.

How many total points are you trying to interpolate?  If your hpts.in 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)

- Josh


On Fri, Dec 13, 2013 at 1:07 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Hi,
>
> 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 hpts.in, it just gives me
> nothing in htps.out!!):
>
> C     Dimension file to be included
> C
> C     HCUBE array dimensions
> C
>       parameter (ldim=3)
>       parameter (lx1=6,ly1=lx1,lz1=lx1,lelt=410,lelv=lelt)
>       parameter (lxd=10,lyd=lxd,lzd=lxd)
>       parameter (lelx=1,lely=1,lelz=1)
>
>       parameter (lzl=3 + 2*(ldim-3))
>
>       parameter (lx2=lx1-2)
>       parameter (ly2=ly1-2)
>       parameter (lz2=lz1-2)
>       parameter (lx3=lx1)
>       parameter (ly3=ly1)
>       parameter (lz3=lz1)
>
>       parameter (lp = 512)
>       parameter (lelg = 16000)
> c
> c     parameter (lpelv=lelv,lpelt=lelt,lpert=3)  ! perturbation
> c     parameter (lpx1=lx1,lpy1=ly1,lpz1=lz1)     ! array sizes
> c     parameter (lpx2=lx2,lpy2=ly2,lpz2=lz2)
> c
>       parameter (lpelv=1,lpelt=1,lpert=1)        ! perturbation
>       parameter (lpx1=1,lpy1=1,lpz1=1)           ! array sizes
>       parameter (lpx2=1,lpy2=1,lpz2=1)
> c
> c     parameter (lbelv=lelv,lbelt=lelt)          ! MHD
> c     parameter (lbx1=lx1,lby1=ly1,lbz1=lz1)     ! array sizes
> c     parameter (lbx2=lx2,lby2=ly2,lbz2=lz2)
> c
>       parameter (lbelv=1,lbelt=1)                ! MHD
>       parameter (lbx1=1,lby1=1,lbz1=1)           ! array sizes
>       parameter (lbx2=1,lby2=1,lbz2=1)
>
> C     LX1M=LX1 when there are moving meshes; =1 otherwise
>       parameter (lx1m=1,ly1m=1,lz1m=1)
>       parameter (ldimt= 2)                       ! 3 passive scalars + T
>       parameter (ldimt1=ldimt+1)
>       parameter (ldimt3=ldimt+3)
> c
> c     Note:  In the new code, LELGEC should be about sqrt(LELG)
> c
>       PARAMETER (LELGEC = 1)
>       PARAMETER (LXYZ2  = 1)
>       PARAMETER (LXZ21  = 1)
>
>       PARAMETER (LMAXV=LX1*LY1*LZ1*LELV)
>       PARAMETER (LMAXT=LX1*LY1*LZ1*LELT)
>       PARAMETER (LMAXP=LX2*LY2*LZ2*LELV)
>       PARAMETER (LXZ=LX1*LZ1)
>       PARAMETER (LORDER=3)
>       PARAMETER (MAXOBJ=4,MAXMBR=LELT*6)
>       PARAMETER (lhis=200)         ! # of pts a proc reads from hpts.in
>                                    ! Note: lhis*np > npoints in hpts.in
> C
> C     Common Block Dimensions
> C
>       PARAMETER (LCTMP0 =2*LX1*LY1*LZ1*LELT)
>       PARAMETER (LCTMP1 =4*LX1*LY1*LZ1*LELT)
> C
> C     The parameter LVEC controls whether an additional 42 field arrays
> C     are required for Steady State Solutions.  If you are not using
> C     Steady State, it is recommended that LVEC=1.
> C
>       PARAMETER (LVEC=1)
> C
> C     Uzawa projection array dimensions
> C
>       parameter (mxprev = 20)
>       parameter (lgmres = 30)
> C
> C     Split projection array dimensions
> C
>       parameter(lmvec = 1)
>       parameter(lsvec = 1)
>       parameter(lstore=lmvec*lsvec)
> c
> c     NONCONFORMING STUFF
> c
>       parameter (maxmor = lelt)
> C
> C     Array dimensions
> C
>       COMMON/DIMN/NELV,NELT,NX1,NY1,NZ1,NX2,NY2,NZ2
>      $,NX3,NY3,NZ3,NDIM,NFIELD,NPERT,NID
>      $,NXD,NYD,NZD
>
> c automatically added by makenek
>       parameter(lxo   = lx1) ! max output grid size (lxo>=lx1)
>
> c automatically added by makenek
>       parameter(lpart = 1  ) ! max number of particles
>
> c automatically added by makenek
>       integer ax1,ay1,az1,ax2,ay2,az2
>       parameter (ax1=lx1,ay1=ly1,az1=lz1,ax2=lx2,ay2=ly2,az2=lz2) !
> running averages
>
> c automatically added by makenek
>       parameter (lxs=1,lys=lxs,lzs=(lxs-1)*(ldim-2)+1) !New Pressure
> Preconditioner
>
> c automatically added by makenek
>       parameter (lfdm=0)  ! == 1 for fast diagonalization method
>
> Is anyone has idea what is the problem in this case?!
>
> Thanks,
> Ami
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>


-- 
Josh Camp

"All that is necessary for the triumph of evil is that good men do nothing"
-- Edmund Burke
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20131217/2987554f/attachment.html>