c++ - 通过 Rcpp 和 bit64 R 包将最大的 int64_t 变量值从 C++ 传递到 R
问题描述
我编写了一个函数来将 2 提升到给定的幂。我想使用 64 位整数。在 R 中,bit64 包的最大和最小限制如下:
来自R
:
> bit64::lim.integer64()
integer64
[1] -9223372036854775807 9223372036854775807
这是-(2^63)
和2^63
。
但是,由于某种原因,我的Rcpp
代码只能传2^62
回 R。这是我的函数的代码,它将 2 提高到给定的幂(注意:我使用位移来实现这一点):
C++
代码:
// [[Rcpp::export]]
Rcpp::NumericVector i2_to_the_power_j ( int64_t j )
{
int64_t base = 1;
int64_t value = base << j;
// cout << "C++ value: " << value << "\n";
// Create a vector of length 1 with `value` as the sole contents
const std::vector<int64_t> v(1, value);
const size_t len = v.size();
Rcpp::NumericVector nn(len); // storage vehicle we return them in
// transfers values 'keeping bits' but changing type
// using reinterpret_cast would get us a warning
std::memcpy(&(nn[0]), &(v[0]), len * sizeof(double));
nn.attr("class") = "integer64";
return nn;
return value;
}
但是,当我在 中运行它时R
,我无法获得最大可能/限制值!
来自R
:
>library(Rcpp)
>library(bit64)
> sourceCpp("./hilbert_curve_copy.cpp")
> # I can get 2^62
> i2_to_the_power_j(62)
integer64
[1] 4611686018427387904
> # ...but I cannot get 2^63
> i2_to_the_power_j(63)
integer64
[1] <NA>
> # I cannot get 2^63, despite bit64 package claiming it can
> # handle integers of this size
> bit64::lim.integer64()
integer64
[1] -9223372036854775807 9223372036854775807
我在这里错过了什么吗?请指教,并感谢您的时间。
解决方案
我的快速猜测(被证明是正确的):最大值本身可能是标记为 NA 的值。因此,计算该值的“减一”并尝试一下。
我的快速猜测:最大值可能是标记为 NA 的值。因此,计算该值的“减一”并尝试一下
// [[Rcpp::export]]
Rcpp::NumericVector largeval ( ) {
int64_t val = 9223372036854775807LL - 1;
Rcpp::Rcout << "C++ value: " << val << "\n";
Rcpp::NumericVector dbl(1);
std::memcpy(&(dbl[0]), &val, sizeof(double));
dbl.attr("class") = "integer64";
return dbl;
}
我将其添加到您的代码中并运行它会产生:
R> largeval()
C++ value: 9223372036854775806
integer64
[1] 9223372036854775806
R>
下面的完整代码以防万一。
代码
#include <Rcpp.h>
// [[Rcpp::export]]
Rcpp::NumericVector i2_to_the_power_j ( int64_t j )
{
int64_t base = 1;
int64_t value = base << j;
// cout << "C++ value: " << value << "\n";
// Create a vector of length 1 with `value` as the sole contents
const std::vector<int64_t> v(1, value);
const size_t len = v.size();
Rcpp::NumericVector nn(len); // storage vehicle we return them in
// transfers values 'keeping bits' but changing type
// using reinterpret_cast would get us a warning
std::memcpy(&(nn[0]), &(v[0]), len * sizeof(double));
nn.attr("class") = "integer64";
return nn;
return value;
}
// [[Rcpp::export]]
Rcpp::NumericVector largeval ( ) {
int64_t val = 9223372036854775807LL - 1;
Rcpp::Rcout << "C++ value: " << val << "\n";
Rcpp::NumericVector dbl(1);
std::memcpy(&(dbl[0]), &val, sizeof(double));
dbl.attr("class") = "integer64";
return dbl;
}
/*** R
library(bit64)
# I can get 2^62
i2_to_the_power_j(62)
# ...but I cannot get 2^63
i2_to_the_power_j(63)
# I cannot get 2^63, despite bit64 package claiming it can
# handle integers of this size
bit64::lim.integer64()
largeval()
*/
推荐阅读
- java - Eclipse photon错误安装过程
- python - 如何使用python的元组替换文件中的值?
- r - 处理列名中的空格(spotfire,R)
- python - 导入 matplotlib.pyplot 缩进问题
- coldfusion - 为什么我会收到我的组件未定义的 ColdFusion 错误?
- assembly - 汇编语言:字节操作和字操作的主要区别是什么?
- gcloud - 未找到 gcloud 或权限被拒绝
- java - 使用 org.json 将 XML 转换为 JSON 时出现 NoSuchMethod 错误
- android - 如何使用 ViewModel 和 Databinding 实现验证?
- python-3.x - 对于更新 - 用于 postgres 的 psycopg2 游标