xplorer0nr

11-14-2005, 05:41 PM

Hello all... Hope someone knows a fix these two questions. Thanks in advance...

Q1. I am having trouble creating the matrix y(n,n) below.

I got a stack overflow error whenever the prgram hits the line declaring the matrix y. the matrix a and its size are known to the subroutine and i checked this. I am using Compac Visual Fortran V 6.5. Question is: how can i fix this? I would appreciate any suggestions to this.

Subroutine ludcmp(a,indx,d)

double precision,DIMENSION(:,:),INTENT(INOUT)::a

double precision,DIMENSION(size(a),size(a))::y

Q2. The following line fails to run for me. This is the line that fails under the same Subroutine ludcmp. Again with stack overflow error. the function outerprod() is well defined but has a line similar to the one in Q1 above. I am not sure if its the parallel change in the indices or something else...

Any ideas how to fix this?

do j=1,n

a(j+1:n,j+1:n)=a(j+1:n,j+1:n)-outerprod(a(j+1:n,j),a(j,j+1:n))

enddo

Thank you very much...

/Mongi

Q1. I am having trouble creating the matrix y(n,n) below.

I got a stack overflow error whenever the prgram hits the line declaring the matrix y. the matrix a and its size are known to the subroutine and i checked this. I am using Compac Visual Fortran V 6.5. Question is: how can i fix this? I would appreciate any suggestions to this.

Subroutine ludcmp(a,indx,d)

double precision,DIMENSION(:,:),INTENT(INOUT)::a

double precision,DIMENSION(size(a),size(a))::y

Q2. The following line fails to run for me. This is the line that fails under the same Subroutine ludcmp. Again with stack overflow error. the function outerprod() is well defined but has a line similar to the one in Q1 above. I am not sure if its the parallel change in the indices or something else...

Any ideas how to fix this?

do j=1,n

a(j+1:n,j+1:n)=a(j+1:n,j+1:n)-outerprod(a(j+1:n,j),a(j,j+1:n))

enddo

Thank you very much...

/Mongi