[R] Problems with sending data to FORTRAN subrutine

[Ales Ziberna]

Dear ExpeRts!



I have problem with sending data to FORTRAN subroutine on Windows XP and R 
2.1.1. When I call it using .Fortran, I get the following error:



R for Windows GUI front-end has encountered a problem and needs to close. 
We are sorry for the inconvenience.



Error signature

AppName: rgui.exe  AppVer: 2.11.50620.0           ModName: read.dll

ModVer: 0.0.0.0      Offset: 00001683



It seams that the error is somehow related to the size of the data. Here are 
the calls:

dyn.load("C:/ales/b_for/read.dll")

#The FORTRUN subroutine used to create the "read.dll" ("read.f") is attached 
(and in case something happens to the attachment also at the end of the 
mail.



n<-as.integer(134);k<-as.integer(2);M<-matrix(as.double(rnorm(n=n^2)),nrow=n,ncol=n);.Fortran("read",M=M,n=as.integer(n),clu=as.integer(sample(1:k,size=n,replace=TRUE)), 
k=as.integer(k),diag=as.integer(1),err=as.double(0.0),E=matrix(as.double(0),ncol=k,nrow=k), 
BM=matrix(as.double(0),ncol=k,nrow=k))

#This works ok



n<-as.integer(134);k<-as.integer(4);M<-matrix(as.double(rnorm(n=n^2)),nrow=n,ncol=n);.Fortran("read",M=M,n=as.integer(n),clu=as.integer(sample(1:k,size=n,replace=TRUE)), 
k=as.integer(k),diag=as.integer(1),err=as.double(0.0),E=matrix(as.double(0),ncol=k,nrow=k), 
BM=matrix(as.double(0),ncol=k,nrow=k))

#Now the k is incrised form 2 to 4 and the error occours.



n<-as.integer(40);k<-as.integer(4);M<-matrix(as.double(rnorm(n=n^2)),nrow=n,ncol=n);.Fortran("read",M=M,n=as.integer(n),clu=as.integer(sample(1:k,size=n,replace=TRUE)), 
k=as.integer(k),diag=as.integer(1),err=as.double(0.0),E=matrix(as.double(0),ncol=k,nrow=k), 
BM=matrix(as.double(0),ncol=k,nrow=k))

#If I leave the k at 4 and reduce the n to 40, then it works.



Any suggestions are welcomed. Thank you in advance!



Best regards,

Ales Ziberna





The commands used to generate "read.dll".

g77 -c read.f

R CMD SHLIB read.o



The FORTRAN subroutine ("read.f"):

        subroutine read(M,n,clu,k,diag,err,E,BM)

        INTEGER n, clu, k, i, j, ii, nA, nAD

        DOUBLE PRECISION M, E, BM, A, AD, vecA, vecAD, err, mean, temp, ss

        LOGICAL diag

        DIMENSION M(n,n), clu(n), E(k,k), BM(k,k), A(k,k,n*n), AD(k,n), 
nA(k,k), nAD(k), vecA(n*n), vecAD(n)





        end