R-alpha: bug in svd

Andreas Weingessel Andreas.Weingessel@ci.tuwien.ac.at
Fri, 17 Oct 1997 16:42:49 +1000


I use R Version 0.60 Alpha (September 18, 1997) on a Linux Pentium
(Debian 1.3) and on a Sparc-Sun-Solaris 2.5.=20


R> svd(matrix(1:16,4,4))
=09yields on both machines
Error: error  4  in dsvdc

R> svd(matrix(1:20,4,5))
=09gives a result on the Linux computer
$d
[1]  0  0  0 NA

$u
     [,1] [,2] [,3] [,4]
[1,]    1    0    0    0
[2,]    0    1    0    0
[3,]    0    0    1    0
[4,]    0    0    0    1

$v
     [,1] [,2] [,3] [,4]
[1,]    1    0    0    0
[2,]    0    1    0    0
[3,]    0    0   NA   NA
[4,]    0    0   NA   NA
[5,]    0    0   NA   NA

but causes a Segmentation Fault on the solaris machine.

=09Andreas Weingessel

***********************************************************************=
*
*                          Andreas Weingessel                          =
*
***********************************************************************=
*
* Institut f=FCr Statistik      *                Tel: (+43 1) 58801 454=
1 *
* Technische Universit=E4t Wien *                Fax: (+43 1)  504 14 9=
8 *
* Wiedner Hauptstr. 8-10/1071 *     Andreas.Weingessel@ci.tuwien.ac.at =
*
* A-1040 Wien, Austria        * http://www.ci.tuwien.ac.at/~weingessel =
*
***********************************************************************=
*



=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
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
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-