Re: [Rd] Building R for windows tools link

2005-05-12 Thread Gregor GORJANC
Prof Brian Ripley wrote:
 On Wed, 11 May 2005, Gorjanc Gregor wrote:
 
 Hello!

 I think that link to page of Duncan Murdoch (bellow) should be given

 http://www.murdoch-sutherland.com/Rtools/

 at

 http://developer.r-project.org/.
 
 
 Are you saying that it is given incorrectly?  I can't see it.
 That site is for the development of R, not for users building R, so it
 does not seem appropriate to link from there.
No, it isn't given incorrectly. I was surprised that it isn't there. I
agree to some point with you, but building R is probably also part of the
development.

 It's easy to launch Google and find mentioned page, but it took me quite
 some time to do this, because I was sure, that it is noted somewhere
 on {www,CRAN}.r-project.org. Maybe I didn't read carefully ...
 
 It _is_ in the R-admin manual, the place building R for Windows is
 discussed (as from 2.1.0).  Apart from in that manual I am unaware of it
 being on the main R web sites.
OK, it looks like I'll have to update a bookmark in my head.

-- 
Lep pozdrav / With regards,
Gregor Gorjanc

--
University of Ljubljana
Biotechnical FacultyURI: http://www.bfro.uni-lj.si/MR/ggorjan
Zootechnical Department mail: gregor.gorjanc at bfro.uni-lj.si
Groblje 3   tel: +386 (0)1 72 17 861
SI-1230 Domzale fax: +386 (0)1 72 17 888
Slovenia, Europe
--
One must learn by doing the thing; for though you think you know it,
 you have no certainty until you try. Sophocles ~ 450 B.C.

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] Patch to address (PR#7853) -- tested briefly, seems to

2005-05-12 Thread ripley
Thank you for the patch.

To clarify: this is not a bug.  ?.C says

  The mapping of the types of R arguments to C or Fortran arguments
  in '.C' or '.Fortran' is

R  C   Fortran
integerint *   integer
numericdouble *double precision
- or - float * real
complexRcomplex *  double complex
logicalint *   integer
character  char ** [see below]
list   SEXP *  not allowed
other  SEXPnot allowed

and `raw' is in the `other' category.  So this would be a change in 
behaviour and needs some changes to the documentation.  We will do that 
for 2.2.0.

Two hints (more for R-devel than just you).

1) Please provide a patch against the current sources, not 2.0.0.  (it 
worked, with a bit of fiddling).

2) Please _attach_ patches, as your mailer wrapped this one.


On Thu, 12 May 2005 [EMAIL PROTECTED] wrote:

 --- dotcode.c2004-09-05 02:31:57.0 -0700
 +++ /usr/local/src/R-2.0.0/src/main/dotcode.c2005-05-11
 12:24:11.0 -0700
 @@ -190,6 +190,7 @@
 static void *RObjToCPtr(SEXP s, int naok, int dup, int narg, int Fort,
 const char *name, R_toCConverter **converter,
   int targetType)
 {
 +unsigned char *rawptr;
 int *iptr;
 float *sptr;
 double *rptr;
 @@ -228,6 +229,16 @@
 }

 switch(TYPEOF(s)) {
 +case RAWSXP:
 +n = LENGTH(s);
 +rawptr = RAW(s);
 +if (dup) {
 +rawptr = (unsigned char *) R_alloc(n, sizeof(unsigned char));
 +for (i = 0; i  n; i++)
 +rawptr[i] = RAW(s)[i];
 +}
 +return (void *) rawptr;
 +break;
 case LGLSXP:
 case INTSXP:
 n = LENGTH(s);
 @@ -329,6 +340,7 @@

 static SEXP CPtrToRObj(void *p, SEXP arg, int Fort,
 R_NativePrimitiveArgType type)
 {
 +unsigned char *rawptr;
 int *iptr, n=length(arg);
 float *sptr;
 double *rptr;
 @@ -339,6 +351,12 @@
 SEXP s, t;

 switch(type) {
 +case RAWSXP:
 +s = allocVector(type, n);
 +rawptr = (unsigned char *)p;
 +for (i = 0; i  n; i++)
 +RAW(s)[i] = rawptr[i];
 +break;
 case LGLSXP:
 case INTSXP:
 s = allocVector(type, n);

 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel



-- 
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] Building R for windows tools link

2005-05-12 Thread Prof Brian Ripley
On Thu, 12 May 2005, Gregor GORJANC wrote:
Prof Brian Ripley wrote:
On Wed, 11 May 2005, Gorjanc Gregor wrote:
Hello!
I think that link to page of Duncan Murdoch (bellow) should be given
http://www.murdoch-sutherland.com/Rtools/
at
http://developer.r-project.org/.

Are you saying that it is given incorrectly?  I can't see it.
That site is for the development of R, not for users building R, so it
does not seem appropriate to link from there.
No, it isn't given incorrectly. I was surprised that it isn't there. I
agree to some point with you, but building R is probably also part of the
development.
As the front page says
`This site is intended as an intermediate repository for more or less 
finalized ideas and plans for the R statistical system.'

What exactly is unclear about that?  How to build current R is not part of 
the development of R.

--
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


[Rd] Windows freezes plotting large line plots (PR#7856)

2005-05-12 Thread timo . becker
Full_Name: Timo Becker
Version: 2.1.0
OS: Windows XP
Submission from: (NULL) (193.170.87.36)


Hi!

The following command causes Windows XP to freeze so that I can only pull the
plug:
 plot(seq(1, 44100), rnorm(44100), type=l)

This does not happen with
 plot(seq(1, 44100), rnorm(44100))
and
 plot(seq(1,100), rnorm(100), type=l)

platform i386-pc-mingw32
arch i386   
os   mingw32
system   i386, mingw32  
status  
major2  
minor1.0
year 2005   
month04 
day  18 
language R

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] Building R for windows tools link

2005-05-12 Thread Gregor GORJANC
Prof Brian Ripley wrote:
 On Thu, 12 May 2005, Gregor GORJANC wrote:
 
 Prof Brian Ripley wrote:

 On Wed, 11 May 2005, Gorjanc Gregor wrote:

 Hello!

 I think that link to page of Duncan Murdoch (bellow) should be given

 http://www.murdoch-sutherland.com/Rtools/

 at

 http://developer.r-project.org/.



 Are you saying that it is given incorrectly?  I can't see it.
 That site is for the development of R, not for users building R, so it
 does not seem appropriate to link from there.

 No, it isn't given incorrectly. I was surprised that it isn't there. I
 agree to some point with you, but building R is probably also part of the
 development.
 
 
 As the front page says
 
 `This site is intended as an intermediate repository for more or less
 finalized ideas and plans for the R statistical system.'
 
 What exactly is unclear about that?  How to build current R is not part
 of the development of R.

Nothing is unclear. I still think that building is part of the development.
I'm sure that R-devils, uups I meant R-devels, don't want to develop
something that is not buildable. Anyway, lets end with this disscussion. I
will remember that mentioned link is in R-admin manual.

-- 
Lep pozdrav / With regards,
Gregor Gorjanc

--
University of Ljubljana
Biotechnical FacultyURI: http://www.bfro.uni-lj.si/MR/ggorjan
Zootechnical Department mail: gregor.gorjanc at bfro.uni-lj.si
Groblje 3   tel: +386 (0)1 72 17 861
SI-1230 Domzale fax: +386 (0)1 72 17 888
Slovenia, Europe
--
One must learn by doing the thing; for though you think you know it,
 you have no certainty until you try. Sophocles ~ 450 B.C.

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] Building R for windows tools link

2005-05-12 Thread Peter Dalgaard
Gregor GORJANC [EMAIL PROTECTED] writes:

  As the front page says
  
  `This site is intended as an intermediate repository for more or less
  finalized ideas and plans for the R statistical system.'
  
  What exactly is unclear about that?  How to build current R is not part
  of the development of R.
 
 Nothing is unclear. I still think that building is part of the development.
 I'm sure that R-devils, uups I meant R-devels, don't want to develop
 something that is not buildable. Anyway, lets end with this disscussion. I
 will remember that mentioned link is in R-admin manual.

Just to be clear: The point is that the developer site was never
intended as being complete and polished, and nobody wants to commit to
making it so. It was set up for things like preliminary white papers
which were too large to handle as part of regular email traffic within
R Core, plus some items that we kept needing to remind ourselves
(and some others) about. It is open to the public, but it is not a
service. 

-- 
   O__   Peter Dalgaard Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics 2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark  Ph: (+45) 35327918
~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


RE: [Rd] bug in modulus operator %% (PR#7852)

2005-05-12 Thread Prof Brian Ripley
I've now found a Windows system that does this.  This is also Windows XP, 
fully patched, and with the same rw2010.  So it may be chip-specific: the 
one that works is a P4 and the one that does not is a latest Pentium M.

I am not sure that the guarantee on the help page has been supported for a 
while.  I've altered the code so it is more likely to be.

BTW,
options(digits=20)
1 %% 0.001
[1] 0.0009792
shows why the original is not a bug in R.
On Thu, 12 May 2005 [EMAIL PROTECTED] wrote:
On Wed, 11 May 2005 [EMAIL PROTECTED] wrote:
Yes, you are correct. I had only checked one of my platforms. Linux
works as you suggest. But for me on Windows,
x - 1
y - 0.2
x %/% y
[1] 5  ## I get a 4 in Linux
I get 5 on Windows, but
(x %% y) + y * (x %/% y)
[1] 1
so is there a problem particular to your Windows runtime?

version
_ =20
platform i386-pc-mingw32
arch i386  =20
os   mingw32   =20
system   i386, mingw32 =20
status =20
major2 =20
minor1.0   =20
year 2005  =20
month04=20
day  18=20
language R =20
-Original Message-
From: Peter Dalgaard [mailto:[EMAIL PROTECTED]
Sent: Wednesday, May 11, 2005 4:14 PM
To: McGehee, Robert
Cc: [EMAIL PROTECTED]; Peter Dalgaard; [EMAIL PROTECTED];
[EMAIL PROTECTED]; r-devel@stat.math.ethz.ch
Subject: Re: [Rd] bug in modulus operator %% (PR#7852)
McGehee, Robert [EMAIL PROTECTED] writes:
Yes, but from ?%%:
It is guaranteed that 'x =3D=3D (x %% y) + y * (x %/% y)' (up to =
rounding
error) ...
=20
(R 2.1.0)
x - 1
y - 0.2
x %% y
[1] 0.2
(x %% y) + y * (x %/% y)
[1] 1.2
=20
Certainly 1 does not equal 1.2 as the documentation would suggest, and
these seem like large enough numbers to not be effected by rounding
errors or lack of precision.
Now that looks a bit odd, but it isn't universal:
x - 1
y - 0.2
x %% y
[1] 0.2
x %/% y
[1] 4
(x %% y) + y * (x %/% y)
[1] 1
So what platform was that happening on?

--
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] bug, feature of mistery?

2005-05-12 Thread Uwe Ligges
Giuseppe Ragusa wrote:
Uwe,

- Where is *lexical* scoping involved?
Abuse of notation. I intended to say scoping.
  - Are you really calling you code from a clean workspace?
Yes it is clean.
- Why don't use pass g through optim() to f? Please do so, because it 
might be a scoping problem.
Tried that. Still no luck. z- does not get assigned even if i pass g
through optim

- The .C call in psi() should not matter unless you are doing strange 
things in the part you omitted (...).

I do not think is the call to C. And I did not omit, psi() is the last
call. And again, same code works fine on the other machine. 
Can you send me a reproducible example (off list), please (inlcuding 
example data)?

Uwe


On Wed, 11 May 2005, Uwe Ligges wrote:

Giuseppe Ragusa wrote:

I have two machines a linux_amd64_x86 (gentoo_amd64) and a linux_x86. Both 
run
R-2.1.0. I have a very long program (hopefully will become a package)
that works perfectly on the linux_amd_x64. Great means no error, no
problems and results that, where the analytic solution exists,
coincide with it. I have problem making the code run on the x64 machine. I
am baffled. The same code on the same version of R, different arch,
is behaving differently. 

After hour of debugging, I traced down what is triggering the error on
the x86 machine.
the code snippet is the following similar to:
R
g - h(d)
f - function( lambda )
{
z - g %*% lambda
sum( psi(z) )
...
}
optim( init.value, f )
R
The function f() is using lexical scoping to get obtain g. The
function psi() is a call to a wrapper function that call a (.C) C
function doing some simple calculation on z. 
- Where is *lexical* scoping involved?
- Are you really calling you code from a clean workspace?
- Why don't use pass g through optim() to f? Please do so, because it 
might be a scoping problem.

- The .C call in psi() should not matter unless you are doing strange 
things in the part you omitted (...).

Uwe Ligges


What's the problem? When f() is called, g is there, lambda is there,
but the assignment  z - g %*% lambda results in a matrix of NaN. This
happens from the second time f() is called, i.e. the first time f() is
called from optim() after the C call has been made. The error is then
that passing a NaN vector to .C results in halted execution. 

If I debug f() during the call to optim, I can without problem assign
z the correct value, but during the execution z is matrix(NaN, nr,
1). 

At this point I can think of the following:
1) the external C code has errors 

It is not a programming error, because when called from console it
 returns the right results. Also, remember, the program work on my
 other machine (the 64 bit);
2) R error (I do not think so)
3) Compiling error
Can be the gcc is messing thing around?
on the 32bit machine
gcc version 3.3.5  (Gentoo Linux 3.3.5-r1, ssp-3.3.2-3, pie-8.7.7.1)
on the 64-bit machine
gcc version 3.4.3 20041125 (Gentoo Linux 3.4.3-r1, ssp-3.4.3-0, pie-8.7.7)
Any help, suggestions, thoughts?
Thank you.
Giuseppe Ragusa
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] Windows freezes plotting large line plots (PR#7856)

2005-05-12 Thread Ted Harding
On 12-May-05 [EMAIL PROTECTED] wrote:
 This is already addressed in R-patched: we do ask you to check
 before filing a bug report so as not to waste our time.
  (It was also discussed on R-help: see the archives.)
 
 It is a bug *in Windows* and not one in R: nothing an application
 can do should be that fatal (and this is a perfectly legitimate
 Windows GDI call).  I trust that you will file a bug report in the
 correct place.
 
From the CHANGES file:
 
The support for mitred lines encountered a Windows problem
with a dramatic slowdown if more than 1000 or so segments
were plotted in a polyline.  We now plot in bunches of 1000,
which is less accurate but  works around the Windows problem.

Clearly a neccessary change, even though it is less accurate.
But does this change apply only to Windows compilations (i.e.
conditional compilation), or will we all have to put up with it?

Ted.



E-Mail: (Ted Harding) [EMAIL PROTECTED]
Fax-to-email: +44 (0)870 094 0861
Date: 12-May-05   Time: 08:26:43
-- XFMail --

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] Windows freezes plotting large line plots (PR#7856)

2005-05-12 Thread Prof Brian Ripley
On Thu, 12 May 2005 [EMAIL PROTECTED] wrote:
On 12-May-05 [EMAIL PROTECTED] wrote:
This is already addressed in R-patched: we do ask you to check
before filing a bug report so as not to waste our time.
 (It was also discussed on R-help: see the archives.)
It is a bug *in Windows* and not one in R: nothing an application
can do should be that fatal (and this is a perfectly legitimate
Windows GDI call).  I trust that you will file a bug report in the
correct place.
From the CHANGES file:
   The support for mitred lines encountered a Windows problem
   with a dramatic slowdown if more than 1000 or so segments
   were plotted in a polyline.  We now plot in bunches of 1000,
   which is less accurate but  works around the Windows problem.
Clearly a neccessary change, even though it is less accurate.
But does this change apply only to Windows compilations (i.e.
conditional compilation), or will we all have to put up with it?
It applies only to the windows() device (including win.metafile etc).
--
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


RE: [Rd] bug in modulus operator %% (PR#7852)

2005-05-12 Thread McGehee, Robert
Could very easily be chipset, but I'm actually running a Pentium 4
myself. Pentium 4 1.9GHz, on a Windows XP Professional Version 2002
Service Professional, Service Pack 1. 

Also of interest,
 1 %% 0.2
[1] 0.2
 1 %% 0.25
[1] 0



-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] 
Sent: Thursday, May 12, 2005 3:39 AM
To: r-devel@stat.math.ethz.ch
Subject: RE: [Rd] bug in modulus operator %% (PR#7852)


I've now found a Windows system that does this.  This is also Windows
XP, 
fully patched, and with the same rw2010.  So it may be chip-specific:
the 
one that works is a P4 and the one that does not is a latest Pentium M.

I am not sure that the guarantee on the help page has been supported for
a 
while.  I've altered the code so it is more likely to be.

BTW,

 options(digits=20)
 1 %% 0.001
[1] 0.0009792

shows why the original is not a bug in R.


On Thu, 12 May 2005 [EMAIL PROTECTED] wrote:

 On Wed, 11 May 2005 [EMAIL PROTECTED] wrote:

 Yes, you are correct. I had only checked one of my platforms. Linux
 works as you suggest. But for me on Windows,

 x - 1
 y - 0.2
 x %/% y
 [1] 5  ## I get a 4 in Linux

 I get 5 on Windows, but

 (x %% y) + y * (x %/% y)
 [1] 1

 so is there a problem particular to your Windows runtime?



 version
 _ =20
 platform i386-pc-mingw32
 arch i386  =20
 os   mingw32   =20
 system   i386, mingw32 =20
 status =20
 major2 =20
 minor1.0   =20
 year 2005  =20
 month04=20
 day  18=20
 language R =20


 -Original Message-
 From: Peter Dalgaard [mailto:[EMAIL PROTECTED]
 Sent: Wednesday, May 11, 2005 4:14 PM
 To: McGehee, Robert
 Cc: [EMAIL PROTECTED]; Peter Dalgaard;
[EMAIL PROTECTED];
 [EMAIL PROTECTED]; r-devel@stat.math.ethz.ch
 Subject: Re: [Rd] bug in modulus operator %% (PR#7852)


 McGehee, Robert [EMAIL PROTECTED] writes:

 Yes, but from ?%%:
 It is guaranteed that 'x =3D=3D (x %% y) + y * (x %/% y)' (up to =
 rounding
 error) ...
 =20
 (R 2.1.0)
 x - 1
 y - 0.2
 x %% y
 [1] 0.2
 (x %% y) + y * (x %/% y)
 [1] 1.2
 =20
 Certainly 1 does not equal 1.2 as the documentation would suggest,
and
 these seem like large enough numbers to not be effected by rounding
 errors or lack of precision.

 Now that looks a bit odd, but it isn't universal:

 x - 1
 y - 0.2
 x %% y
 [1] 0.2
 x %/% y
 [1] 4
 (x %% y) + y * (x %/% y)
 [1] 1

 So what platform was that happening on?


-- 
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] bug in modulus operator %% (PR#7852)

2005-05-12 Thread Kjetil Brinchmann Halvorsen
Prof Brian Ripley wrote:
I've now found a Windows system that does this.  This is also Windows 
XP, fully patched, and with the same rw2010.  So it may be 
chip-specific: the one that works is a P4 and the one that does not is 
a latest Pentium M.

I am not sure that the guarantee on the help page has been supported 
for a while.  I've altered the code so it is more likely to be.

BTW,
options(digits=20)
1 %% 0.001
[1] 0.0009792
shows why the original is not a bug in R.
Why is that not a bug in R? On my machine (windows XP, rw2010 from CRAN) 
I get:

test - function(x, y) (x %% y) + y * ( x %/% y ) # should be x
 test(1, 0.001)
[1] 1.001
 test(1, 0.1)
[1] 1.1
 test(1, 1)
[1] 1
 test(1, 0.01)
[1] 1.01
and this differences (well, not the third one) cannot be said to be 
rounding
error.

Kjetil
On Thu, 12 May 2005 [EMAIL PROTECTED] wrote:
On Wed, 11 May 2005 [EMAIL PROTECTED] wrote:
Yes, you are correct. I had only checked one of my platforms. Linux
works as you suggest. But for me on Windows,
x - 1
y - 0.2
x %/% y
[1] 5  ## I get a 4 in Linux

I get 5 on Windows, but
(x %% y) + y * (x %/% y)
[1] 1
so is there a problem particular to your Windows runtime?

version
_ =20
platform i386-pc-mingw32
arch i386  =20
os   mingw32   =20
system   i386, mingw32 =20
status =20
major2 =20
minor1.0   =20
year 2005  =20
month04=20
day  18=20
language R =20
-Original Message-
From: Peter Dalgaard [mailto:[EMAIL PROTECTED]
Sent: Wednesday, May 11, 2005 4:14 PM
To: McGehee, Robert
Cc: [EMAIL PROTECTED]; Peter Dalgaard; [EMAIL PROTECTED];
[EMAIL PROTECTED]; r-devel@stat.math.ethz.ch
Subject: Re: [Rd] bug in modulus operator %% (PR#7852)
McGehee, Robert [EMAIL PROTECTED] writes:
Yes, but from ?%%:
It is guaranteed that 'x =3D=3D (x %% y) + y * (x %/% y)' (up to =
rounding
error) ...
=20
(R 2.1.0)
x - 1
y - 0.2
x %% y
[1] 0.2
(x %% y) + y * (x %/% y)
[1] 1.2
=20
Certainly 1 does not equal 1.2 as the documentation would suggest, and
these seem like large enough numbers to not be effected by rounding
errors or lack of precision.

Now that looks a bit odd, but it isn't universal:
x - 1
y - 0.2
x %% y
[1] 0.2
x %/% y
[1] 4
(x %% y) + y * (x %/% y)
[1] 1
So what platform was that happening on?



--
Kjetil Halvorsen.
Peace is the most effective weapon of mass construction.
  --  Mahdi Elmandjra

--
No virus found in this outgoing message.
Checked by AVG Anti-Virus.
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


[Rd] segfault on large number of open brackets (PR#7859)

2005-05-12 Thread m . e . bruche
Full_Name: Max
Version: R 2.0.1 (2004-11-15)
OS: Gentoo Linux
Submission from: (NULL) (158.143.49.181)


I leaned on the ( key by accident, and it looks as if R segfaults on a large
number (88 or more in my case) open brackets:

Script started on Thu May 12 15:18:04 2005
$ R

R : Copyright 2004, The R Foundation for Statistical Computing
Version 2.0.1  (2004-11-15), ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for a HTML browser interface to help.
Type 'q()' to quit R.

 ((
(
+ (
Segmentation fault
$ uname -a
Linux max 2.6.11-gentoo-r6 #2 Sat May 7 19:24:52 GMT 2005 i686 Intel(R)
Pentium(R) 4 CPU 2.40GHz GenuineIntel GNU/Linux
$ exit

Script done on Thu May 12 15:18:16 2005

Regards,

Max

R compiled with
gcc 3.3.5
glibc-2.3.4

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] bug in modulus operator %% (PR#7852)

2005-05-12 Thread Prof Brian Ripley
On Thu, 12 May 2005, Kjetil Brinchmann Halvorsen wrote:
Prof Brian Ripley wrote:
I've now found a Windows system that does this.  This is also Windows XP, 
fully patched, and with the same rw2010.  So it may be chip-specific: the 
one that works is a P4 and the one that does not is a latest Pentium M.

I am not sure that the guarantee on the help page has been supported for a 
while.  I've altered the code so it is more likely to be.

BTW,
options(digits=20)
1 %% 0.001
[1] 0.0009792
shows why the original is not a bug in R.
Why is that not a bug in R? On my machine (windows XP, rw2010 from CRAN) I
Because 0 = y %% x  y, so that adding any multiple of 0.001 takes it out 
of range.

get:
test - function(x, y) (x %% y) + y * ( x %/% y ) # should be x
test(1, 0.001)
[1] 1.001
test(1, 0.1)
[1] 1.1
test(1, 1)
[1] 1
test(1, 0.01)
[1] 1.01
and this differences (well, not the third one) cannot be said to be rounding
error.
I think they can.  The problem is that x %/% y is affected by rounding 
error in the representation of 0.01, as in

1 %/% 0.01
[1] 99
on my machine.  That is due to rounding error.
--
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


RE: [Rd] segfault on large number of open brackets (PR#7859)

2005-05-12 Thread Liaw, Andy
I believe that's been fixed in R-2.1.0.  Please check.

I got the same segfault with R-2.0.1 on SuSE Linux x86_64, but on both that
machine and my WinXP latop, I get syntax error.  E.g.,

R : Copyright 2005, The R Foundation for Statistical Computing
Version 2.1.0 Patched (2005-05-12), ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

  Natural language support but running in an English locale

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for a HTML browser interface to help.
Type 'q()' to quit R.

 (((
+ (((
Error: syntax error

Andy


 From: [EMAIL PROTECTED]
 
 Full_Name: Max
 Version: R 2.0.1 (2004-11-15)
 OS: Gentoo Linux
 Submission from: (NULL) (158.143.49.181)
 
 
 I leaned on the ( key by accident, and it looks as if R 
 segfaults on a large
 number (88 or more in my case) open brackets:
 
 Script started on Thu May 12 15:18:04 2005
 $ R
 
 R : Copyright 2004, The R Foundation for Statistical Computing
 Version 2.0.1  (2004-11-15), ISBN 3-900051-07-0
 
 R is free software and comes with ABSOLUTELY NO WARRANTY.
 You are welcome to redistribute it under certain conditions.
 Type 'license()' or 'licence()' for distribution details.
 
 R is a collaborative project with many contributors.
 Type 'contributors()' for more information and
 'citation()' on how to cite R or R packages in publications.
 
 Type 'demo()' for some demos, 'help()' for on-line help, or
 'help.start()' for a HTML browser interface to help.
 Type 'q()' to quit R.
 
  
 ((
 
 (
 + (
 Segmentation fault
 $ uname -a
 Linux max 2.6.11-gentoo-r6 #2 Sat May 7 19:24:52 GMT 2005 
 i686 Intel(R)
 Pentium(R) 4 CPU 2.40GHz GenuineIntel GNU/Linux
 $ exit
 
 Script done on Thu May 12 15:18:16 2005
 
 Regards,
 
 Max
 
 R compiled with
 gcc 3.3.5
 glibc-2.3.4
 
 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel
 
 


__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


RE: [Rd] segfault on large number of open brackets (PR#7859)

2005-05-12 Thread andy_liaw
I believe that's been fixed in R-2.1.0.  Please check.

I got the same segfault with R-2.0.1 on SuSE Linux x86_64, but on both that
machine and my WinXP latop, I get syntax error.  E.g.,

R : Copyright 2005, The R Foundation for Statistical Computing
Version 2.1.0 Patched (2005-05-12), ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

  Natural language support but running in an English locale

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for a HTML browser interface to help.
Type 'q()' to quit R.

 (((
+ (((
Error: syntax error

Andy


 From: [EMAIL PROTECTED]
 
 Full_Name: Max
 Version: R 2.0.1 (2004-11-15)
 OS: Gentoo Linux
 Submission from: (NULL) (158.143.49.181)
 
 
 I leaned on the ( key by accident, and it looks as if R 
 segfaults on a large
 number (88 or more in my case) open brackets:
 
 Script started on Thu May 12 15:18:04 2005
 $ R
 
 R : Copyright 2004, The R Foundation for Statistical Computing
 Version 2.0.1  (2004-11-15), ISBN 3-900051-07-0
 
 R is free software and comes with ABSOLUTELY NO WARRANTY.
 You are welcome to redistribute it under certain conditions.
 Type 'license()' or 'licence()' for distribution details.
 
 R is a collaborative project with many contributors.
 Type 'contributors()' for more information and
 'citation()' on how to cite R or R packages in publications.
 
 Type 'demo()' for some demos, 'help()' for on-line help, or
 'help.start()' for a HTML browser interface to help.
 Type 'q()' to quit R.
 
  
 ((
 
 (
 + (
 Segmentation fault
 $ uname -a
 Linux max 2.6.11-gentoo-r6 #2 Sat May 7 19:24:52 GMT 2005 
 i686 Intel(R)
 Pentium(R) 4 CPU 2.40GHz GenuineIntel GNU/Linux
 $ exit
 
 Script done on Thu May 12 15:18:16 2005
 
 Regards,
 
 Max
 
 R compiled with
 gcc 3.3.5
 glibc-2.3.4
 
 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel
 
 


__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] segfault on large number of open brackets (PR#7859)

2005-05-12 Thread ligges
Liaw, Andy wrote:

 I believe that's been fixed in R-2.1.0.  Please check.
 
 I got the same segfault with R-2.0.1 on SuSE Linux x86_64, but on both that
 machine and my WinXP latop, I get syntax error.  E.g.,
 
 R : Copyright 2005, The R Foundation for Statistical Computing
 Version 2.1.0 Patched (2005-05-12), ISBN 3-900051-07-0
 
 R is free software and comes with ABSOLUTELY NO WARRANTY.
 You are welcome to redistribute it under certain conditions.
 Type 'license()' or 'licence()' for distribution details.
 
   Natural language support but running in an English locale
 
 R is a collaborative project with many contributors.
 Type 'contributors()' for more information and
 'citation()' on how to cite R or R packages in publications.
 
 Type 'demo()' for some demos, 'help()' for on-line help, or
 'help.start()' for a HTML browser interface to help.
 Type 'q()' to quit R.
 
 
(((
 
 + (((
 Error: syntax error


With R-2.1.0 (release!) I get syntax error on Windows, no error (at 
least not after the first couple of hundred openings) on Linux, but the 
described *segfault* on Solaris 5.7 (UltraSparc).
[I am too lazy to try out R-patched on the slowish Solaris machine yet]

Uwe Ligges


 Andy
 
 
 
From: [EMAIL PROTECTED]

Full_Name: Max
Version: R 2.0.1 (2004-11-15)
OS: Gentoo Linux
Submission from: (NULL) (158.143.49.181)


I leaned on the ( key by accident, and it looks as if R 
segfaults on a large
number (88 or more in my case) open brackets:

Script started on Thu May 12 15:18:04 2005
$ R

R : Copyright 2004, The R Foundation for Statistical Computing
Version 2.0.1  (2004-11-15), ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for a HTML browser interface to help.
Type 'q()' to quit R.


((

(
+ (
Segmentation fault
$ uname -a
Linux max 2.6.11-gentoo-r6 #2 Sat May 7 19:24:52 GMT 2005 
i686 Intel(R)
Pentium(R) 4 CPU 2.40GHz GenuineIntel GNU/Linux
$ exit

Script done on Thu May 12 15:18:16 2005

Regards,

Max

R compiled with
gcc 3.3.5
glibc-2.3.4

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel



 
 
 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] segfault on large number of open brackets (PR#7859)

2005-05-12 Thread ripley
Yes, it is coincidence.

I have found the error and have a fix: more detailed message to follow.

Brian

On Thu, 12 May 2005 [EMAIL PROTECTED] wrote:

 Liaw, Andy wrote:

 I believe that's been fixed in R-2.1.0.  Please check.

 I got the same segfault with R-2.0.1 on SuSE Linux x86_64, but on both that
 machine and my WinXP latop, I get syntax error.  E.g.,

 R : Copyright 2005, The R Foundation for Statistical Computing
 Version 2.1.0 Patched (2005-05-12), ISBN 3-900051-07-0

 R is free software and comes with ABSOLUTELY NO WARRANTY.
 You are welcome to redistribute it under certain conditions.
 Type 'license()' or 'licence()' for distribution details.

   Natural language support but running in an English locale

 R is a collaborative project with many contributors.
 Type 'contributors()' for more information and
 'citation()' on how to cite R or R packages in publications.

 Type 'demo()' for some demos, 'help()' for on-line help, or
 'help.start()' for a HTML browser interface to help.
 Type 'q()' to quit R.


 (((

 + (((
 Error: syntax error


 With R-2.1.0 (release!) I get syntax error on Windows, no error (at
 least not after the first couple of hundred openings) on Linux, but the
 described *segfault* on Solaris 5.7 (UltraSparc).
 [I am too lazy to try out R-patched on the slowish Solaris machine yet]

 Uwe Ligges


 Andy



 From: [EMAIL PROTECTED]

 Full_Name: Max
 Version: R 2.0.1 (2004-11-15)
 OS: Gentoo Linux
 Submission from: (NULL) (158.143.49.181)


 I leaned on the ( key by accident, and it looks as if R
 segfaults on a large
 number (88 or more in my case) open brackets:

 Script started on Thu May 12 15:18:04 2005
 $ R

 R : Copyright 2004, The R Foundation for Statistical Computing
 Version 2.0.1  (2004-11-15), ISBN 3-900051-07-0

 R is free software and comes with ABSOLUTELY NO WARRANTY.
 You are welcome to redistribute it under certain conditions.
 Type 'license()' or 'licence()' for distribution details.

 R is a collaborative project with many contributors.
 Type 'contributors()' for more information and
 'citation()' on how to cite R or R packages in publications.

 Type 'demo()' for some demos, 'help()' for on-line help, or
 'help.start()' for a HTML browser interface to help.
 Type 'q()' to quit R.


 ((
 
 (
 + (
 Segmentation fault
 $ uname -a
 Linux max 2.6.11-gentoo-r6 #2 Sat May 7 19:24:52 GMT 2005
 i686 Intel(R)
 Pentium(R) 4 CPU 2.40GHz GenuineIntel GNU/Linux
 $ exit

 Script done on Thu May 12 15:18:16 2005

 Regards,

 Max

 R compiled with
 gcc 3.3.5
 glibc-2.3.4

 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel





 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel

 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel



-- 
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


RE: [Rd] segfault on large number of open brackets (PR#7859)

2005-05-12 Thread ripley
I suspect that is coincidence, as I can see no source change that would 
account for it.

2.0.1 and R-devel work for me and 2.1.0 segfaults with

999 KeepSource = *LOGICAL(GetOption(install(keep.source), 
R_NilValue));
(gdb) bt
#0  0x080ca229 in ParseInit () at gram.y:999
#1  0x080cb305 in token () at gram.y:1530

I am not sure that is right, but I am pretty sure that

static char contextstack[50], *contextp;

just above is an unchecked nesting limit, as in

static void IfPush(void)
{
 if (*contextp==LBRACE ||
*contextp=='['||
*contextp=='('||
*contextp == 'i')
*++contextp = 'i';
}

However, although I thought this was the problem, it seems only to apply 
to nested if()s.  So I do not currently have a solution.


On Thu, 12 May 2005 [EMAIL PROTECTED] wrote:

 I believe that's been fixed in R-2.1.0.  Please check.

 I got the same segfault with R-2.0.1 on SuSE Linux x86_64, but on both that
 machine and my WinXP latop, I get syntax error.  E.g.,

 R : Copyright 2005, The R Foundation for Statistical Computing
 Version 2.1.0 Patched (2005-05-12), ISBN 3-900051-07-0

 R is free software and comes with ABSOLUTELY NO WARRANTY.
 You are welcome to redistribute it under certain conditions.
 Type 'license()' or 'licence()' for distribution details.

  Natural language support but running in an English locale

 R is a collaborative project with many contributors.
 Type 'contributors()' for more information and
 'citation()' on how to cite R or R packages in publications.

 Type 'demo()' for some demos, 'help()' for on-line help, or
 'help.start()' for a HTML browser interface to help.
 Type 'q()' to quit R.

 (((
 + (((
 Error: syntax error

 Andy


 From: [EMAIL PROTECTED]

 Full_Name: Max
 Version: R 2.0.1 (2004-11-15)
 OS: Gentoo Linux
 Submission from: (NULL) (158.143.49.181)


 I leaned on the ( key by accident, and it looks as if R
 segfaults on a large
 number (88 or more in my case) open brackets:

 Script started on Thu May 12 15:18:04 2005
 $ R

 R : Copyright 2004, The R Foundation for Statistical Computing
 Version 2.0.1  (2004-11-15), ISBN 3-900051-07-0

 R is free software and comes with ABSOLUTELY NO WARRANTY.
 You are welcome to redistribute it under certain conditions.
 Type 'license()' or 'licence()' for distribution details.

 R is a collaborative project with many contributors.
 Type 'contributors()' for more information and
 'citation()' on how to cite R or R packages in publications.

 Type 'demo()' for some demos, 'help()' for on-line help, or
 'help.start()' for a HTML browser interface to help.
 Type 'q()' to quit R.


 ((
 
 (
 + (
 Segmentation fault
 $ uname -a
 Linux max 2.6.11-gentoo-r6 #2 Sat May 7 19:24:52 GMT 2005
 i686 Intel(R)
 Pentium(R) 4 CPU 2.40GHz GenuineIntel GNU/Linux
 $ exit

 Script done on Thu May 12 15:18:16 2005

 Regards,

 Max

 R compiled with
 gcc 3.3.5
 glibc-2.3.4

 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel




 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel



-- 
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] Implementing R on IBM p690 cluster Jump

2005-05-12 Thread Brian D Ripley
As far as I know R 2.0.1 was built on AIX.  A quick diff

gannet% diff ~/R/R-2.0.1/src/modules/X11/Makefile
~/R/R-2.1.0/src/modules/X11/Makefile
54c54
   $(SHLIB_LINK) -o $@ $(R_X11_la_LDFLAGS) $(R_X11_la_OBJECTS)
$(R_X11_la_LIBADD) $(LIBS)
---
   $(SHLIB_LINK) -o $@ $(R_X11_la_LDFLAGS) $(R_X11_la_OBJECTS)
$(R_X11_la_LIBADD)
128c128,129
   ../../../src/include/Rmodules/RX11.h
---
   ../../../src/include/Rmodules/RX11.h \
   ../../../src/include/Rconnections.h

suggests that AIX (alone) required $(LIBS) there.  As it is harmful on
some other platforms (the intention is to pick up entry points in
R.bin/libR rather than duplicate copies of static libraries) I suspect we
need an AIX-specific workaround (in particular the libintl library is not
going to be in LIBS).

This probably applies to other modules and even shared objects in
packages.

Is there anyone with an AIX machine who is able to work through this and
provide appropriate diffs?


On Thu, 12 May 2005, Bill Northcott wrote:

 On 11/05/2005, at 11:26 PM, Kurt Hornik wrote:
  There seem to be problems in the autoconf stuff.  cos and sin are
  being located in libm but the link line you show has no '-lm' which
  accounts for 6 of the 7 missing symbols.  Also gettext seems to have
  been located by the configure script, but there is no option on the
  link line to include it in the link. Also search configure.in for
  strings like powerpc-ibm-aix to see what stuff is being done
  specifically for AIX.  It may be out of date in relation to your
  system.
 
 
 
  If it was my problem, I would have look at those bits of the
  configure
  script and the sections of configure.in that generated them.  Then I
  would compare the AC_SUBST variables with the makefile.in that
  generated the link.  It seems there is some non- communication there.

  Yes most of the undefined symbols seem to come from libm.  But they
  would also be in libm on all other platforms which have such, and
  we do
  not include it there.  On Linux I get
 
  gcc-3.4 -shared -L/usr/local/lib -o R_X11.so  dataentry.lo
  devX11.lo rotated.lo rbitmap.lo  -lSM -lICE -L/usr/X11R6/lib -lX11
  -ljpeg -lpng -lz -L../../../lib -lR

 What happens on Linux is just not really relevant to what happens on
 AIX, although it might explain why the errors have not been picked
 up.  A number of platforms do not have a separate libm.  Instead they
 have a symlink to some general C library.  So typically libm is
 linked by default on these platforms.  It looks as if that is not the
 case on AIX, and that a specific -lm flag is needed.
 
  and this works.  As Brian already wrote, we are building a *module*
  here, so I assume that we are not getting the flags for building
  modules
  right.  Assuming that libtool can, it might be helpful to look at the
  output of BUILDDIR/libtool --config.
 
 The error messages are about missing symbols.

 Bill

 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel


-- 
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] Implementing R on IBM p690 cluster Jump

2005-05-12 Thread Kurt Hornik
 Brian D Ripley writes:

 As far as I know R 2.0.1 was built on AIX.  A quick diff
 gannet% diff ~/R/R-2.0.1/src/modules/X11/Makefile
 ~/R/R-2.1.0/src/modules/X11/Makefile
 54c54
$(SHLIB_LINK) -o $@ $(R_X11_la_LDFLAGS) $(R_X11_la_OBJECTS)
 $(R_X11_la_LIBADD) $(LIBS)
 ---
 $(SHLIB_LINK) -o $@ $(R_X11_la_LDFLAGS) $(R_X11_la_OBJECTS)
 $(R_X11_la_LIBADD)
 128c128,129
../../../src/include/Rmodules/RX11.h
 ---
 ../../../src/include/Rmodules/RX11.h \
 ../../../src/include/Rconnections.h

 suggests that AIX (alone) required $(LIBS) there.  As it is harmful on
 some other platforms (the intention is to pick up entry points in
 R.bin/libR rather than duplicate copies of static libraries) I suspect we
 need an AIX-specific workaround (in particular the libintl library is not
 going to be in LIBS).

 This probably applies to other modules and even shared objects in
 packages.

Yes.  It would still be good if someone on AIX could provide us with the
output of BUILDDIR/libtool --config.

-k

 Is there anyone with an AIX machine who is able to work through this
 and provide appropriate diffs?


 On Thu, 12 May 2005, Bill Northcott wrote:

 On 11/05/2005, at 11:26 PM, Kurt Hornik wrote:
  There seem to be problems in the autoconf stuff.  cos and sin are
  being located in libm but the link line you show has no '-lm' which
  accounts for 6 of the 7 missing symbols.  Also gettext seems to have
  been located by the configure script, but there is no option on the
  link line to include it in the link. Also search configure.in for
  strings like powerpc-ibm-aix to see what stuff is being done
  specifically for AIX.  It may be out of date in relation to your
  system.
 
 
 
  If it was my problem, I would have look at those bits of the
  configure
  script and the sections of configure.in that generated them.  Then I
  would compare the AC_SUBST variables with the makefile.in that
  generated the link.  It seems there is some non- communication there.
 
  Yes most of the undefined symbols seem to come from libm.  But they
  would also be in libm on all other platforms which have such, and
  we do
  not include it there.  On Linux I get
 
  gcc-3.4 -shared -L/usr/local/lib -o R_X11.so  dataentry.lo
  devX11.lo rotated.lo rbitmap.lo  -lSM -lICE -L/usr/X11R6/lib -lX11
  -ljpeg -lpng -lz -L../../../lib -lR
 
 What happens on Linux is just not really relevant to what happens on
 AIX, although it might explain why the errors have not been picked
 up.  A number of platforms do not have a separate libm.  Instead they
 have a symlink to some general C library.  So typically libm is
 linked by default on these platforms.  It looks as if that is not the
 case on AIX, and that a specific -lm flag is needed.
 
  and this works.  As Brian already wrote, we are building a *module*
  here, so I assume that we are not getting the flags for building
  modules
  right.  Assuming that libtool can, it might be helpful to look at the
  output of BUILDDIR/libtool --config.
 
 The error messages are about missing symbols.
 
 Bill
 
 __
 R-devel@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-devel
 

 -- 
 Brian D. Ripley,  [EMAIL PROTECTED]
 Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
 University of Oxford, Tel:  +44 1865 272861 (self)
 1 South Parks Road, +44 1865 272866 (PA)
 Oxford OX1 3TG, UKFax:  +44 1865 272595

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


RE: [Rd] segfault on large number of open brackets (PR#7859)

2005-05-12 Thread Prof Brian Ripley
On Thu, 12 May 2005 [EMAIL PROTECTED] wrote:
I suspect that is coincidence, as I can see no source change that would
account for it.
2.0.1 and R-devel work for me and 2.1.0 segfaults with
999 KeepSource = *LOGICAL(GetOption(install(keep.source), 
R_NilValue));
(gdb) bt
#0  0x080ca229 in ParseInit () at gram.y:999
#1  0x080cb305 in token () at gram.y:1530
I am not sure that is right, but I am pretty sure that
static char contextstack[50], *contextp;
just above is an unchecked nesting limit, as in
static void IfPush(void)
{
if (*contextp==LBRACE ||
*contextp=='['||
*contextp=='('||
*contextp == 'i')
*++contextp = 'i';
}
However, although I thought this was the problem, it seems only to apply
to nested if()s.  So I do not currently have a solution.
I was almost right, just a different instance of *++contextp.
Fixed in R-patched, to give
(((
Error: contextstack overflow
--
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


[Rd] Bug in axis labels (PR#7860)

2005-05-12 Thread murdoch
Bob O'hara wrote:
 I'm a bit reluctant to call anything a bug: I know it's usually my 
 incompetence instead.  In this case, I can't see what else it is, 
 although it may be a bug in Windows.
 
 The problem comes from trying to create a .png of a figure in Windows 
 XP, with R2.1.0.  On the screen it looks OK, but in the .png the text 
 for the x label is smaller than for the y label.  The problem seems to 
 be created by using axes=F in plot().  Here's some code to demonstrate 
 the problem:
 
 XX=1:5;  YY=1:5
 png(thing.png)
 par(mfrow=c(2,1))
  plot(XX, YY, axes=F, xlab=X label, ylab=Y label)
  plot(XX, YY, xlab=X label, ylab=Y label)
 dev.off()
 
 I get a similar problem with bmp() and jpeg(), but not postscript() or 
 pdf().  If I draw the figure in the window in R and then save it as a 
 .png from the menu, the png looks fine.
 
 Everything works fine in R1.5.1 (I know, I know!) on Debian Linux, which 
 is why I'm reporting this as a Windows problem.

I see the problem in the first figure, and it's still there in a recent 
r-devel build.  My test was also on Windows. Since png() uses the 
Windows graphics driver, it does look like an R bug in the Windows 
graphics driver, but I'd like to hear from someone on a different 
platform...

Duncan Murdoch

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] Bug in axis labels (PR#7860)

2005-05-12 Thread Peter Dalgaard
[EMAIL PROTECTED] writes:

  Everything works fine in R1.5.1 (I know, I know!) on Debian Linux, which 
  is why I'm reporting this as a Windows problem.
 
 I see the problem in the first figure, and it's still there in a recent 
 r-devel build.  My test was also on Windows. Since png() uses the 
 Windows graphics driver, it does look like an R bug in the Windows 
 graphics driver, but I'd like to hear from someone on a different 
 platform...

I don't see it with 2.1.0 on Linux (FC3).

-- 
   O__   Peter Dalgaard Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics 2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark  Ph: (+45) 35327918
~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


[Rd] R make errors

2005-05-12 Thread Ambika Sundaresan
I tried to install the latest version of R (2-1.0)in anAIX machine. I was able 
to tar, and configure with no
errors.  However, when I said make, I got the following errors:
gcc -Wl,-bM:SRE -Wl,-H512 -Wl,-T512 -Wl,-bnoentry -Wl,-bexpall
-Wl,-bI:../../../etc/R.exp -L/usr/local/lib -o R_X11.so  dataentry.lo devX11.lo
rotated.lo rbitmap.lo -lSM -lICE -lX11
ld: 0711-317 ERROR: Undefined symbol: .log10
ld: 0711-317 ERROR: Undefined symbol: .floor
ld: 0711-317 ERROR: Undefined symbol: .libintl_gettext
ld: 0711-317 ERROR: Undefined symbol: .pow
ld: 0711-317 ERROR: Undefined symbol: .sin
ld: 0711-317 ERROR: Undefined symbol: .cos
ld: 0711-317 ERROR: Undefined symbol: .tan
ld: 0711-345 Use the -bloadmap or -bnoquiet option to obtain more information.
collect2: ld returned 8 exit status
make: 1254-004 The error code from the last command is 1.
Stop.
make: 1254-004 The error code from the last command is 2.
I understand I have to edit the Makefile.  Can someone tell me which 
Makefile(folder) and also which flags to edit.  I greatly appreciate your help. 
 Thanks
Ambika Sundaresan
Grad. Student
UCSD.
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


RE: [Rd] segfault on large number of open brackets (PR#7859)

2005-05-12 Thread Thomas Lumley
And this also fixes PR#7022 (a highly nested polynomial)
-thomas
On Thu, 12 May 2005, Prof Brian Ripley wrote:
On Thu, 12 May 2005 [EMAIL PROTECTED] wrote:
I suspect that is coincidence, as I can see no source change that would
account for it.
2.0.1 and R-devel work for me and 2.1.0 segfaults with
999 KeepSource = *LOGICAL(GetOption(install(keep.source), 
R_NilValue));
(gdb) bt
#0  0x080ca229 in ParseInit () at gram.y:999
#1  0x080cb305 in token () at gram.y:1530

I am not sure that is right, but I am pretty sure that
static char contextstack[50], *contextp;
just above is an unchecked nesting limit, as in
static void IfPush(void)
{
if (*contextp==LBRACE ||
*contextp=='['||
*contextp=='('||
*contextp == 'i')
*++contextp = 'i';
}
However, although I thought this was the problem, it seems only to apply
to nested if()s.  So I do not currently have a solution.
I was almost right, just a different instance of *++contextp.
Fixed in R-patched, to give
(((
Error: contextstack overflow
--
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
Thomas Lumley   Assoc. Professor, Biostatistics
[EMAIL PROTECTED]   University of Washington, Seattle
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] R make errors

2005-05-12 Thread Prof Brian Ripley
Please see the thread entitled:
Implementing R on  IBM p690 cluster Jump
in this list's archives.  We don't have a solution currently, but we do 
need some help.

On Thu, 12 May 2005, Ambika Sundaresan wrote:
I tried to install the latest version of R (2-1.0)in anAIX machine. I was 
able to tar, and configure with no
errors.  However, when I said make, I got the following errors:

   gcc -Wl,-bM:SRE -Wl,-H512 -Wl,-T512 -Wl,-bnoentry -Wl,-bexpall
-Wl,-bI:../../../etc/R.exp -L/usr/local/lib -o R_X11.so  dataentry.lo 
devX11.lo
rotated.lo rbitmap.lo -lSM -lICE -lX11
ld: 0711-317 ERROR: Undefined symbol: .log10
ld: 0711-317 ERROR: Undefined symbol: .floor
ld: 0711-317 ERROR: Undefined symbol: .libintl_gettext
ld: 0711-317 ERROR: Undefined symbol: .pow
ld: 0711-317 ERROR: Undefined symbol: .sin
ld: 0711-317 ERROR: Undefined symbol: .cos
ld: 0711-317 ERROR: Undefined symbol: .tan
ld: 0711-345 Use the -bloadmap or -bnoquiet option to obtain more 
information.
collect2: ld returned 8 exit status
make: 1254-004 The error code from the last command is 1.

Stop.
make: 1254-004 The error code from the last command is 2.
I understand I have to edit the Makefile.  Can someone tell me which 
Makefile(folder) and also which flags to edit.  I greatly appreciate your 
help.  Thanks

Ambika Sundaresan
Grad. Student
UCSD.
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel

--
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595
__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


[Rd] RFC822/M400 Mail Network -- Delivery Report (PR#7861)

2005-05-12 Thread root
--Delivery-Status-Report

Not delivered to: [EMAIL PROTECTED]
user agent unavailable

--Delivery-Status-Report
Content-Type: message/delivery-status

Original-Envelope-Id: in*vsnl*rfc987;4284390b4708000mimey2k
X400-Content-Identifier: 050513105010+053
Reporting-MTA: x400; /ADMD=VSNL/C=IN
DSN-Gateway: smtp; terminator1.vsnl.net.in

Final-Recipient: rfc822;
S=chen/G=inter/PRMD=gems/ADMD=vsnl/[EMAIL PROTECTED]
Action: failed
Diagnostic-Code: x400; 1 4
Status: 5.1.1
Last-Attempted-Date: 13 May 2005 10:36 +0530
X-Supplementary: Hello World
--Delivery-Status-Report--

__
R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel