Re: Help: Unable to allocate memory for work

I ran the following on our local IBM ... it ran fine.
It took a *very* long time [ :-( , not sure why]
but it encountered no problems.

============
begin
   jdm = 30
   nsvd = 3

 do N=50,51
   idm = 64*N
   x = random_normal( -5, 5, (/idm,jdm/) )
   y = random_normal(-20,50, (/idm,jdm/) )

   ncols=idm
   ncolz=ncols

   homlft = new((/nsvd,ncols/),float)
   hetlft = new((/nsvd,ncols/),float)
   homrgt = new((/nsvd,ncolz/),float)
   hetrgt = new((/nsvd,ncolz/),float)

   svd_y = svdstd(x,y,nsvd,homlft,hetlft,homrgt,hetrgt)
   print(svd_y)

   delete(x) delete an recreate for 51
   delete(y)
   delete(homlft)
   delete(hetlft)
   delete(homrgt)
   delete(hetrgt)
   delete(svd_y)
 end do

end
======

====================================> N=50

Variable: svd_y
Type: float
Total Size: 12 bytes
            3 values
Number of Dimensions: 1
Dimensions and sizes: [3]
Coordinates:
Number Of Attributes: 5
  bk : <ARRAY>
  ak : <ARRAY>
  lapack_err : 0
  condn : -0.0006778455
  fnorm : -0.006238905
(0) 4.291725
(1) 4.211056
(2) 4.128187
(0) svdstd

====================================> N=51

Variable: svd_y
Type: float
Total Size: 12 bytes
            3 values
Number of Dimensions: 1
Dimensions and sizes: [3]
Coordinates:
Number Of Attributes: 5
  bk : <ARRAY>
  ak : <ARRAY>
  lapack_err : 0
  condn : -0.002549693
  fnorm : -0.0003648346
(0) 4.359809
(1) 4.188076
(2) 4.127023

>
>I still got the same error.
>% ulimit -a
>core file size (blocks, -c) unlimited
>data seg size (kbytes, -d) unlimited
>file size (blocks, -f) unlimited
>max memory size (kbytes, -m) unlimited
>open files (-n) unlimited
>pipe size (512 bytes, -p) 64
>stack size (kbytes, -s) unlimited
>cpu time (seconds, -t) unlimited
>max user processes (-u) 1024
>virtual memory (kbytes, -v) unlimited
>
>And this time, I use a IBM p690 with 64G mem, but still got same error
>:Unable to allocate memory for work.
>I test a script with left and right array with same size of 64*50*30,it
>works.
>I check the system info. mem use < 1% .
>But change idm to 64*51, it can not work with the error.
>the script is below:
>~~~~~~~~~~~~~~~~~~~~
>idm=64*50 ; 51
>jdm=30
>x=new((/idm,jdm/),"float")
>y=new((/idm,jdm/),"float")
>... ...
>nsvd = 3
>ncols=idm
>ncolz=ncols
> homlft = new((/nsvd,ncols/),float)
> hetlft = new((/nsvd,ncols/),float)
> homrgt = new((/nsvd,ncolz/),float)
> hetrgt = new((/nsvd,ncolz/),float)
> svd_y = svdstd(x,y,nsvd,homlft,hetlft,homrgt,hetrgt)
> system("rm testout.nc")
>end
>~~~~~~~~~~~~~~~~~~~~~~~~
>
>
>How should I do?
>
>>From: Kate Hedstrom <kate_at_arsc.edu>
>>To: nothingtolose terry <nothingtolose_999_at_hotmail.com>
>>CC: ncl-talk_at_ucar.edu
>>Subject: Re: Help: Unable to allocate memory for work
>>Date: Mon, 13 Mar 2006 11:29:18 -0900
>>
>>On Mon, Mar 13, 2006 at 03:55:39PM +0000, nothingtolose terry wrote:
>> > Hi, I write a short script for svdstd, left and right field array size
>is
>> > about
>> > 8000*30 and 8200*30 float numbers.
>> > I get error: fatal:svdstd: Unable to allocate memory for work array
>> >
>> > How should I do?
>> >
>> > BTW, I run it at IBM P690 machine with 32 G Mem.
>> >
>> > Thank you very much.
>>
>>Try running "ulimit -a" to see what your personal limits are. There may
>>be 32 GB, but they only let you use part of it.
>>
>>f2n1 8% ulimit -a
>>time(seconds) 7200
>>file(blocks) unlimited
>>data(kbytes) 2097152
>>stack(kbytes) unlimited
>>memory(kbytes) 2097152
>>coredump(blocks) 20000
>>nofiles(descriptors) 2000
>>
>>The numbers to note here are datasize and memoryuse, both currently set
>>to 2 GB. In this case, I'm actually on the front end to a P690 and
>>would have to submit a batch job to use the P690 proper - and the
>>script would have to adjust the limits to be larger than the default:
>>
>>ulimit -dS unlimited
>>ulimit -mS unlimited
>>
>>Kate
>>--
>>Kate Hedström Arctic Region Supercomputing Center
>>kate_at_arsc.edu University of Alaska, Fairbanks
>>Phone: 907-450-8678
>
>_________________________________________________________________
>ÓëÁª»úµÄÅóÓѽøÐн»Á÷£¬ÇëʹÓà MSN Messenger: http://messenger.msn.com/cn
>
>_______________________________________________
>ncl-talk mailing list
>ncl-talk_at_ucar.edu
>http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Mar 15 2006 - 09:33:53 MST