From ihaka@stat.auckland.ac.nz  Mon Aug 11 22:14:42 1997
From: ihaka@stat.auckland.ac.nz (Ross Ihaka)
Date: Tue, 12 Aug 1997 09:14:42 +1200 (NZST)
Subject: R-alpha: eigen
In-Reply-To: <Pine.SUN.3.95.970811094129.432B-100000@zen>
References: <199708111336.PAA00446@aragorn.ci.tuwien.ac.at>
 <Pine.SUN.3.95.970811094129.432B-100000@zen>
Message-ID: <199708112114.JAA09123@stat1.stat.auckland.ac.nz>

Thomas Lumley writes:
 > 
 > eigen is behaving very strangely in my copy of patch 3 (solaris 2.5.1). 
 > The calls to .Fortran are returning absolute rubbish and are also mutating
 > things that shouldn't be mutated. 

I think that this is my fault for making a "quick fix" without really
testing it.  The problem (as pointed out by Martin) is in dotcode.c.

The function naoktrim should read

static SEXP naoktrim(SEXP s, int * len, int *naok, int *dup)
{
        SEXP value;

        if(s == R_NilValue) {
                value = R_NilValue;
                *naok = 0;
                *dup = 1;	/* <---- 1 not 0 */
                *len = 0;
        }

Sorry about that.
	Ross
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-


R-alpha: eigen
Ross Ihaka