From kjell.konis at epfl.ch  Wed Jul  1 10:30:22 2009
From: kjell.konis at epfl.ch (Kjell Konis)
Date: Wed, 1 Jul 2009 10:30:22 +0200
Subject: [Rd] Matrix with random number
In-Reply-To: <30454.53644.qm@web59716.mail.ac4.yahoo.com>
References: <30454.53644.qm@web59716.mail.ac4.yahoo.com>
Message-ID: <0C5135DD-C2AB-4E6B-A8A6-5813B7FA17E6@epfl.ch>

Hi Fabio,

Your function myrbeta returns void so assigning the output isn't going  
to work.  Instead you need to call it like a FORTRAN subroutine.   
Also, I added arguments for the parameters of the beta and moved the  
fseedi and fseedo calls outside of the loop.

This is a pretty basic FORTRAN programming question and there are lots  
of books that can help you learn FORTRAN.


     subroutine mat(x,l,c,pa,pb)
     integer l,c
     double precision x(l,c), a
     integer i,j
     call fseedi()
     do j = 1, c
        do i = 1, l
          call myrbeta(x(i,j),pa,pb)
         enddo
       enddo
     call fseedo()
     end


In R call it like this:

storage.mode(x) <- "double"
.Fortran("mat", x = x, as.integer(l), as.integer(c), as.double(1),  
as.double(2))


Cheers,
Kjell



On 30 juin 09, at 20:02, Fabio Mathias wrote:

> Thanks Mr. Barry Rowlingson
>
> However, the matrix appears to zeros!
>
> Notice the code below! Please!
>
> Code in fortran
>
>     subroutine mat(x,l,c,a)
>     integer l,c
>     double precision x(l,c), a
>     integer i,j
>      do j = 1, c
>        do i = 1, l
>         call fseedi()
>          x(i,j) = myrbeta(a,1,2)
>         call fseedo()
>            enddo
>           enddo
>     end
>
> In R:
>
> dyn.load("func.so")
> x <- matrix(0,5,6)
> l <- nrow(x)
> c <-
> ncol(x)
> a <- 0
>
> ..Fortran("mat", x = x, l, c, as.double(a))
>
> Results:
>
> $x
>      [,1] [,2] [,3] [,4] [,5] [,6]
> [1,]    0    0    0    0    0    0
> [2,]    0    0    0    0    0    0
> [3,]    0    0    0    0    0    0
> [4,]    0    0    0    0    0    0
> [5,]    0    0    0    0    0    0
>
> [[2]]
> [1] 5
>
> [[3]]
> [1] 6
>
> [[4]]
> [1] 1
>
>
> Thanks!!!
>
>
>              F?bio Mathias Corr?a                       UFLA
>
>
>
>
>       
> ____________________________________________________________________________________
> [[elided Yahoo spam]]
>
> 	[[alternative HTML version deleted]]
>
> <ATT00001.txt>



[Rd] Matrix with random number
Kjell Konis