[R-pkg-devel] not resolved from current namespace error

[Tom Wainwright]

Not sure, but your problem might be answered in the .Fortran() help page:

All Fortran compilers known to be usable to compile R map symbol names to
> lower case, and so does .Fortran.
>

I've been caught by that before, and found that using all lowercase names
for Fortran routines in R is safest.

Tom Wainwright

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
The contents of this message are mine personally and do not
necessarily reflect any position of the Government or the
National Oceanic and Atmospheric Administration.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


On Tue, Aug 11, 2015 at 10:32 AM, Ignacio Martinez <ignacio82 at gmail.com>
wrote:

> I'm trying to create a package that uses a MPI Fortran module. I have a
> working
> version <https://github.com/ignacio82/MyPi> of that package that uses a
> Fortran module without MPI.
>
> When I run the function `FMPIpi(DARTS = 5000, ROUNDS = 100, cores=2)` I get
> the following errors:
>
> > FMPIpi(DARTS = 5000, ROUNDS = 100, cores=2)
> 2 slaves are spawned successfully. 0 failed.
> master (rank 0, comm 1) of size 3 is running on: 2d60fd60575b
> slave1 (rank 1, comm 1) of size 3 is running on: 2d60fd60575b
> slave2 (rank 2, comm 1) of size 3 is running on: 2d60fd60575b
> Error in .Fortran("MPIpi", avepi = as.numeric(1), DARTS =
> as.integer(DARTS),  :
>  "mpipi" not resolved from current namespace (MyPi)
>
> It looks like something is wrong in my NAMESPACE. This is what I have
> there:
>
> export(Pibenchmark)
> export(Fpi)
> export(FMPIpi)
> export(Rpi)
> useDynLib(MyPi)
> exportPattern("^[[:alpha:]]+")
> This is my Fortran module:
>
> Module Fortranpi
> USE MPI
> IMPLICIT NONE
> contains
> subroutine dboard(darts, dartsscore)
>  integer, intent(in)                    :: darts
>  double precision, intent(out)          :: dartsscore
>  double precision                       :: x_coord, y_coord
>  integer                                :: score, n
>
> score = 0
> do n = 1, darts
>  call random_number(x_coord)
>  call random_number(y_coord)
>
>  if ((x_coord**2 + y_coord**2) <= 1.0d0) then
>  score = score + 1
>  end if
> end do
>
> dartsscore = 4.0d0*score/darts
>
> end subroutine dboard
>
> subroutine pi(avepi, DARTS, ROUNDS) bind(C, name="pi_")
>  use, intrinsic                         :: iso_c_binding, only : c_double,
> c_int
>  real(c_double), intent(out)            ::  avepi
>  integer(c_int), intent(in)             ::  DARTS, ROUNDS
>  integer                                ::  MASTER, rank, i, n
>  integer, allocatable                   ::  seed(:)
>  double precision                       ::  pi_est, homepi, pirecv, pisum
>
> ! we set it to zero in the sequential run
> rank = 0
> ! initialize the random number generator
> ! we make sure the seed is different for each task
> call random_seed()
> call random_seed(size = n)
> allocate(seed(n))
> seed = 12 + rank*11
> call random_seed(put=seed(1:n))
> deallocate(seed)
>
> avepi = 0
> do i = 0, ROUNDS-1
>  call dboard(darts, pi_est)
>  ! calculate the average value of pi over all iterations
>  avepi = ((avepi*i) + pi_est)/(i + 1)
> end do
> end subroutine pi
>
>
> subroutine MPIpi(avepi, DARTS, ROUNDS) bind(C, name="MPIpi_")
> use, intrinsic                         :: iso_c_binding, only : c_double,
> c_int
> real(c_double), intent(out)            :: avepi
> integer(c_int), intent(in)             :: DARTS, ROUNDS
> integer                                :: i, n, mynpts, ierr, numprocs,
> proc_num
> integer, allocatable                   :: seed(:)
> double precision                       :: pi_est, y, sumpi
>
>  call mpi_init(ierr)
>  call mpi_comm_size(MPI_COMM_WORLD, numprocs, ierr)
>  call mpi_comm_rank(MPI_COMM_WORLD, proc_num, ierr)
>
>  if (numprocs .eq. 0) then
> mynpts = ROUNDS - (numprocs-1)*(ROUNDS/numprocs)
>  else
> mynpts = ROUNDS/numprocs
>  endif
>
>  ! initialize the random number generator
>  ! we make sure the seed is different for each task
>  call random_seed()
>  call random_seed(size = n)
>  allocate(seed(n))
>  seed = 12 + proc_num*11
>  call random_seed(put=seed(1:n))
>  deallocate(seed)
>
>  y=0.0d0
> do i = 1, mynpts
> call dboard(darts, pi_est)
> y = y + pi_est
>  end do
>
>  call mpi_reduce(y, sumpi, 1, mpi_double_precision, mpi_sum, 0, &
>  mpi_comm_world, ierr)
>  if (proc_num==0) avepi = sumpi/ROUNDS
>  call mpi_finalize(ierr)
> end subroutine MPIpi
>
> end module Fortranpi
> and this is my R function
>
> #'@export
> FMPIpi <- function(DARTS, ROUNDS, cores) {
>  Rmpi::mpi.spawn.Rslaves(nslaves=cores)
>  retvals <- .Fortran("MPIpi", avepi = as.numeric(1), DARTS =
>  as.integer(DARTS), ROUNDS =  as.integer(ROUNDS))
>  return(retvals$avepi)
> }
>
> What am I doing wrong?
>
> Thanks a lot!
>
> Ignacio
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-package-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-package-devel
>

	[[alternative HTML version deleted]]