1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
|
### Coercion and Methods for Symmetric Packed Matrices
dsp2dsy <- function(from) .Call(dspMatrix_as_dsyMatrix, from)
dsp2C <- function(from) dsy2C(.Call(dspMatrix_as_dsyMatrix, from))
setAs("dspMatrix", "dsyMatrix", dsp2dsy)
## setAs("dspMatrix", "dsCMatrix", dsp2C)
setAs("dspMatrix", "CsparseMatrix", dsp2C)
setAs("dspMatrix", "sparseMatrix", dsp2C)
## dge <--> dsp via dsy
.dense2sp <- function(from) .dsy2dsp(.dense2sy(from))
setAs("dgeMatrix", "dspMatrix", .dense2sp)
setAs("matrix", "dspMatrix",
function(from) .dense2sp(..2dge(from)))
## S3-matrix <--> dsp via dsy
setAs("dspMatrix", "matrix", function(from) .dsy2mat(dsp2dsy(from)))
setMethod("rcond", signature(x = "dspMatrix", norm = "character"),
function(x, norm, ...)
.Call(dspMatrix_rcond, x, norm),
valueClass = "numeric")
setMethod("rcond", signature(x = "dspMatrix", norm = "missing"),
function(x, norm, ...)
.Call(dspMatrix_rcond, x, "O"),
valueClass = "numeric")
setMethod("BunchKaufman", signature(x = "dspMatrix"),
function(x, ...) .Call(dspMatrix_trf, x))
## Should define multiplication from the right
setMethod("solve", signature(a = "dspMatrix", b = "missing"),
function(a, b, ...) .Call(dspMatrix_solve, a),
valueClass = "dspMatrix")
setMethod("solve", signature(a = "dspMatrix", b = "matrix"),
function(a, b, ...) .Call(dspMatrix_matrix_solve, a, b),
valueClass = "dgeMatrix")
setMethod("solve", signature(a = "dspMatrix", b = "ddenseMatrix"),
function(a, b, ...) .Call(dspMatrix_matrix_solve, a, b),
valueClass = "dgeMatrix")
##setMethod("solve", signature(a = "dspMatrix", b = "numeric"),
## function(a, b, ...)
## .Call(dspMatrix_matrix_solve, a, as.matrix(b)),
## valueClass = "dgeMatrix")
## No longer needed
## setMethod("solve", signature(a = "dspMatrix", b = "integer"),
## function(a, b, ...) {
## storage.mode(b) <- "double"
## .Call(dspMatrix_matrix_solve, a, as.matrix(b))
## }, valueClass = "dgeMatrix")
setMethod("norm", signature(x = "dspMatrix", type = "character"),
function(x, type, ...)
if(identical("2", type)) norm2(x) else .Call(dspMatrix_norm, x, type),
valueClass = "numeric")
setMethod("norm", signature(x = "dspMatrix", type = "missing"),
function(x, type, ...) .Call(dspMatrix_norm, x, "O"),
valueClass = "numeric")
setMethod("t", signature(x = "dspMatrix"),
function(x) .dsy2dsp(t(dsp2dsy(x))),
valueClass = "dspMatrix")
setMethod("diag", signature(x = "dspMatrix"),
function(x, nrow, ncol) .Call(dspMatrix_getDiag, x))
setMethod("diag<-", signature(x = "dspMatrix"),
function(x, value) .Call(dspMatrix_setDiag, x, value))
|
