[Nek5000-users] problem with interpolation by hpts

Tue Dec 17 09:47:23 CST 2013

H

On Tue, Dec 17, 2013 at 8:04 AM, <nek5000-users at lists.mcs.anl.gov> wrote:

> 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.infile 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
>
> _______________________________________________
> 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: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20131217/5eade0a4/attachment-0001.html>