topic Re: Proc FCMP: Program ends with no error or result in SAS Procedures

Proc FCMP: Program ends with no error or result

AmitKB — Thu, 18 Oct 2012 03:57:57 GMT

Hi all,

I was able to run form a 5000 nodes matrix in around 2 hours. But when i run for a 8000 nodes matrix it just ends in a second with no error. Is it memory issue or something else.

Appreciate you help.

Amit

Code:

proc fcmp;

array A[8000,8000] / nosymbols;

array Ident[8000,8000] / nosymbols;

array Z[8000,8000] / nosymbols;

array P[8000,8000] / nosymbols;

do i = 1 to 8000;

do j = 1 to 8000;

A[i, j] = RANUNI(12345);

if A[i, j] le .4 then A[i, j] = 0;

else A[i, j] = 1;

A[i, j] = 0.5*A[i, j];

Ident[i,j] = A[i,j];

end;

call IDENTITY (Ident);

CALL SUBTRACTMATRIX (Ident, A, Z);

CALL INV (Z, P);

rc = write_array('result', p);

run;

log:

NOTE: PROCEDURE FCMP used (Total process time):

real time 0.24 seconds

cpu time 0.01 seconds

Re: Proc FCMP: Program ends with no error or result

ArtC — Fri, 19 Oct 2012 05:49:42 GMT

The matrix being inverted seems a bit large, and memory may indeed be an issue. Have you tried increasing the memory buffers?

Although it won't help with the memory, we can speed up the processing a bit by replacing your A[I,J] IF-THEN/ELSE assignment statements.

r=ranuni(12345);

a[i,j] = 0*(r<=.4)+.5*(r>.4);

Re: Proc FCMP: Program ends with no error or result

AmitKB — Fri, 19 Oct 2012 15:19:38 GMT

Hi Art,

Appreciate your help. I have previously taken you class on Macro programming. It was really helpful. Thank you for being a great teacher.

It is memory issue. How can i increase the memory buffer as you mentioned? How does one find how much the memory buffer is?

Also I have updated my code to remove the if then else.

Regards,

Amit

Re: Proc FCMP: Program ends with no error or result

ArtC — Fri, 19 Oct 2012 16:43:25 GMT

An earlier thread on memory usage and matricies can be found at: https://communities.sas.com/message/140513#140513

The system option for memory allocation is MEMSIZE. You may also get more detailed responses in the SAS/IML forum, which specializes in issues with matrix manipulation.