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
|
### =========================================================================
### Manipulation of array selections
### -------------------------------------------------------------------------
###
### Like base::arrayInd() but faster and accepts a matrix for 'dim' (with 1
### row per element in 'Lindex').
Lindex2Mindex <- function(Lindex, dim, use.names=FALSE)
{
if (!isTRUEorFALSE(use.names))
stop("'use.names' must be TRUE or FALSE")
## 'dim' can be a matrix so it's important to use storage.mode()
## instead of as.integer().
if (storage.mode(dim) == "double")
storage.mode(dim) <- "integer"
## 'Lindex' and 'dim' will be fully checked at the C level.
.Call2("C_Lindex2Mindex", Lindex, dim, use.names, PACKAGE="DelayedArray")
}
Mindex2Lindex <- function(Mindex, dim, use.names=FALSE, as.integer=FALSE)
{
if (!isTRUEorFALSE(use.names))
stop("'use.names' must be TRUE or FALSE")
if (!isTRUEorFALSE(as.integer))
stop("'as.integer' must be TRUE or FALSE")
## 'dim' and/or 'Mindex' can be matrices so it's important to use
## storage.mode() instead of as.integer(). Also, unlike as.integer(),
## this preserves the names/dimnames.
if (storage.mode(dim) == "double")
storage.mode(dim) <- "integer"
if (storage.mode(Mindex) == "double")
storage.mode(Mindex) <- "integer"
## 'Mindex' and 'dim' will be fully checked at the C level.
.Call2("C_Mindex2Lindex", Mindex, dim, use.names,
as.integer, PACKAGE="DelayedArray")
}
### 'aind' must be an integer matrix or vector (a vector is treated like
### a 1-row matrix).
### Return a numeric vector with one element per row in 'aind'.
linearInd <- function(aind, dim)
{
.Defunct("Mindex2Lindex")
Mindex2Lindex(aind, dim, use.names=TRUE)
}
|
