← Back to branch summary

~showard314/ubuntu/karmic/r-base/remove_start_comments

« back to all changes in this revision

Viewing changes to src/library/base/R/eigen.R

  • Committer: Bazaar Package Importer
  • Author(s): Dirk Eddelbuettel
  • Date: 2009-04-17 06:56:48 UTC
  • mfrom: (1.3.1 upstream) (6.1.2 sid)
  • Revision ID: james.westby@ubuntu.com-20090417065648-9mgm2udb2eanluot
* New upstream version released this morning

* debian/rules: Turn optimisation back to -O3 on alpha (as elsewhere) 
  with thanks to Kurt Roeckx for applying a fix to gcc

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/library/base/R/eigen.R
17
17
 
18
18
isSymmetric <- function(object, ...) UseMethod("isSymmetric")
19
19
 
20
 
isSymmetric.matrix <- function(object, tol = 100*.Machine$double.eps, ...) {
 
20
isSymmetric.matrix <- function(object, tol = 100*.Machine$double.eps, ...)
 
21
{
21
22
    if(!is.matrix(object)) return(FALSE) ## we test for  symmetric *matrix*
22
23
    ## cheap pretest: is it square?
23
24
    d <- dim(object)
24
 
    if(d[1] != d[2]) return(FALSE)
 
25
    if(d[1L] != d[2L]) return(FALSE)
25
26
    test <-
26
27
        if(is.complex(object))
27
28
            all.equal.numeric(object, Conj(t(object)), tolerance = tol, ...)
82
83
                          dbl.n,
83
84
                          dbl.n,
84
85
                          double(2*n),
85
 
                          ierr = integer(1),
 
86
                          ierr = integer(1L),
86
87
                          PACKAGE="base")
87
88
            if (z$ierr)
88
89
                stop(gettextf("'ch' returned code %d in 'eigen'", z$ierr),
101
102
                          vectors = x,
102
103
                          dbl.n,
103
104
                          dbl.n,
104
 
                          ierr = integer(1),
 
105
                          ierr = integer(1L),
105
106
                          PACKAGE="base")
106
107
            if (z$ierr)
107
108
                stop(gettextf("'rs' returned code %d in 'eigen'", z$ierr),
126
127
                          dbl.n,
127
128
                          dbl.n,
128
129
                          dbl.n,
129
 
                          ierr = integer(1),
 
130
                          ierr = integer(1L),
130
131
                          PACKAGE="base")
131
132
            if (z$ierr)
132
133
                stop(gettextf("'cg' returned code %d in 'eigen'", z$ierr),
147
148
                          vectors = x,
148
149
                          integer(n),
149
150
                          dbl.n,
150
 
                          ierr = integer(1),
 
151
                          ierr = integer(1L),
151
152
                          PACKAGE="base")
152
153
            if (z$ierr)
153
154
                stop(gettextf("'rg' returned code %d in 'eigen'", z$ierr),
154
155
                     domain = NA)
155
 
            ind <- z$ivalues > 0
 
156
            ind <- z$ivalues > 0L
156
157
            if(any(ind)) {#- have complex (conjugated) values
157
158
                ind <- seq.int(n)[ind]
158
159
                z$values <- complex(real=z$values,imaginary=z$ivalues)