uniquedouble: Map positive 64-bit integers onto unique doubles

In obreschkow/cooltools: Practical Tools for Scientific Computations and Visualizations

Map positive 64-bit integers onto unique doubles

Description

Uses a reversible one-to-one mapping to converts non-negative 64-bit integers, produced via the bit64 package, to double-precision floating-point values for faster handling of 64-bit integer indices, e.g. comparison of indices, count number of unique indices, etc. The double values are completely different from the corresponding integer values.

Usage

uniquedouble(x, inverse = FALSE)

Arguments

x	If inverse=FALSE (default), this is a scalar/vector/array of non-negative integers, typically 64-bit integers. If inverse=TRUE, it is a scalar/vector/array of doubles.
inverse	Logical flag. If TRUE, doubles are converted back to their original 64-bit integers.

Value

Scalar/vector/array of double values, if inverse=FALSE, or integer64 values, if inverse=TRUE.

Author(s)

Danail Obreschkow

Examples

# produce example vector of non-negative 64-bit integers
n = 1e3
x = c(bit64::as.integer64(0:n),
      bit64::as.integer64('9218868437227405311')+c(-n:n),
      bit64::lim.integer64()[2]-c(0:n))

# map this vector onto unique doubles
y = uniquedouble(x)

# check if the double-vector can be inverted back to the integer-vector
cat(sprintf("Check inversion = %s.\n", all(x==uniquedouble(y, inverse=TRUE))))

# measure the time taken to compare all the int64 values
tick("Compare original integer64 values. ")
comparison = TRUE
for(i in seq_along(x)) comparison = comparison & x[i]==x[i]
tock(comparison, fmt="Time taken = %.3fs. Test result = %s.\n")

# measure the time taken to compare all the corresponding double values
tick("Compare corresponding double values. ")
comparison = TRUE
for(i in seq_along(y)) comparison = comparison & y[i]==y[i]
tock(comparison, fmt="Time taken = %.3fs. Test result = %s.\n")

obreschkow/cooltools documentation built on Nov. 16, 2024, 2:46 a.m.