Re: [Rd] package Geneland / Rgui under windows

2007-10-15 Thread Gilles GUILLOT
I forgot to say that I was the package author. I suspected a bug of R as Geneland worked fine for two years and the problem popped up with the release of R 2.6.0 And I can't see any explanation. So, any clue would help. gilles > Gilles GUILLOT wrote: >> Hi, >> I experience

[Rd] package Geneland / Rgui under windows

2007-10-15 Thread Gilles GUILLOT
Hi, I experienced a problem with the package Geneland under R 2.6.0 with windows XP professional. The commands below should simulate a dataset, then make an MCMC simulation stored in tempdir(). It works with R 2.5.1 (both GUI and command line) It works with the command line of R 2.6.0 but not wi

[Rd] RNG stuck via Fortran call

2005-11-30 Thread Gilles GUILLOT
Having not much success with my previous question I try to reformulate it: I'm simulating a Markow chain in Fortran interfaced with R. Each loop of my Fortran calls various functions of the R RNG through the wrapper given below. In a run of 100 iterations of the Markov kernel, after 20 iterati

[Rd] R binomial RNG stuck at 1 via Fortran call

2005-11-30 Thread Gilles GUILLOT
wrapper now as part as the message, sorry about previous mail #include #include void F77_SUB(rndstart)(void) { GetRNGstate(); } void F77_SUB(rndend)(void) { PutRNGstate(); } double F77_SUB(ggrnorm)(double *mu, double *sigma) { return rnorm(*mu, *sigma); } double F77_SUB(ggrexp)(double *

[Rd] R binomial RNG stuck at 1 via Fortran call

2005-11-30 Thread Gilles GUILLOT
wrapper now attached -- ÅÅÅ Gilles Guillot INRA - MIA Paris Currently working from Göteborg Stochastic Centre Eklandagatan 86 - Rum1439 Chalmers University of Technology SE 412-96 Göteborg Sweden Phone +46 31 772 3514 / Fax

[Rd] R binomial RNG stuck at 1 via Fortran call

2005-11-30 Thread Gilles GUILLOT
(1.d0,0.5d0) write(*,*) 'b=',b In certain cases, after a few iterations in the loop, b get stuck at 1 Any hint to explain that would help. Gilles -- ÅÅÅÅÅÅÅ Gilles Guillot INRA - MIA Paris Currently working fro

[Rd] R build under mandriva 10.2

2005-11-08 Thread Gilles GUILLOT
After upgrading from mandrake 10.1 to mandriva 10.2 I can't build shared archive with R 2.2.0 [EMAIL PROTECTED] src]$ R CMD SHLIB main.f sub.f wrapper.c gcc -shared -L/usr/local/lib -o main.so main.o sub.o wrapper.o -lg2c -lm -lgcc_s /usr//bin/ld: cannot find -lg2c collect2: ld returned 1 exit

Re: [Rd] pb with dyn.load - fortran code now attached

2005-10-27 Thread Gilles GUILLOT
> > is.loaded("supsmu_") > > [1] FALSE > > > is.loaded("supsmu") > > [1] TRUE > > That is a Fortran entry point, and it complies with the description > quoted. > > What is.loaded() needs depends on how the symbol would be found > and so it is not much use. OK, thanks. But still > is.loaded("fst

[Rd] pb with dyn.load - fortran code now attached

2005-10-26 Thread Gilles GUILLOT
Hi, here are a couple of strange things happening with R.2.2.0 compiled under Mandrake-10.1 At the begining of the story I have segmentation fault using package RandomFields in conjuction with some Fortran code of mine loaded with dyn.load. One of my fortran programs contains a subroutine nam

[Rd] pb with dyn.load

2005-10-26 Thread Gilles GUILLOT
Hi, here are a couple of strange things happening with R.2.2.0 compiled under Mandrake-10.1 At the begining of the story I have segmentation fault using package RandomFields in conjuction with some Fortran code of mine loaded with dyn.load. One of my fortran programs contains a subroutine nam

[Rd] calling fortran functions CHOL and DPOTRF form Fortran

2005-07-05 Thread Gilles GUILLOT
? Thanks in advance, Gilles -- _________ Gilles GUILLOT INRA -Département Mathématiques et Informatique Appliquées Unité Mixte de Recherche INRA - INAPG - ENGREF Institut National Agronomique de Paris-Grignon 16, rue Claude Bernard 752

Re: [Rd] Calling C from Fortran

2005-06-15 Thread Gilles GUILLOT
at > Fortran <-> C conventions. > > Reid Huntsinger > > -Original Message- > From: Gilles GUILLOT [mailto:[EMAIL PROTECTED] > Sent: Wednesday, June 15, 2005 3:50 AM > To: Huntsinger, Reid > Subject: Re: [Rd] Calling C from Fortran > > > Thanks Reid! >

[Rd] Calling C function from Fortran

2005-06-15 Thread Gilles GUILLOT
Hi all, the example in the R doc and the hints from Shusong Jin , Ingmar Visser and Reid Huntsinger (thanks all three) refer to the case where the function does not have arguments. I'm still looking for a proper sequence of commands to call C functions with arguemnts from R. Imagine I want

[Rd] Calling C from Fortran

2005-06-14 Thread Gilles GUILLOT
I would like to call C routines from Fortran under linux as suggested in section 5.6 of the "Writing R extensions" documentation. I'm familiar with Fortran but not with C. I understand the example provided in Fortran: subroutine testit() double precision normrnd, x call rndstart() x = normrnd()