【问题标题】:Converting R split() function to C++将 R split() 函数转换为 C++
【发布时间】:2021-08-20 07:52:57
【问题描述】:

考虑 R 中的可重现示例:

test <- c(1:12)
> test
 [1]  1  2  3  4  5  6  7  8  9 10 11 12

预期结果:

test.list <- split(test, gl(2, 3))
> test.list
$`1`
[1] 1 2 3 7 8 9

$`2`
[1]  4  5  6 10 11 12

我正在尝试用 C++ 编写等效代码,以生成并返回由 test.list 生成的两个向量。请注意,我在 C++ 中处于尴尬的新手阶段。

【问题讨论】:

标签: c++ r split rcpp


【解决方案1】:

我们可以使用@jignatius 的好答案,并将其设为 R 可调用函数。为简单起见,我将其保留在NumericVector;我们这里有很多答案,显示根据运行时有效负载在 NumericVectorIntegerVector 之间切换。

代码

#include <Rcpp.h>

// [[Rcpp::export]]
Rcpp::List mysplit(Rcpp::NumericVector nums, int n, int size) {
    std::vector<std::vector<double>> result(n);
    int i = 0;
    auto beg = nums.cbegin();
    auto end = nums.cend();

    while (beg != nums.cend()) {
        //get end iterator safely
        auto next = std::distance(beg, end) >= size ? beg + size : end;
        //insert into result
        result[i].insert(result[i].end(), beg, next);
        //advance iterator
        beg = next;
        i = (i + 1) % n;
    }

    Rcpp::List ll;
    for (const auto&v : result)
        ll.push_back(v);

    return ll;
}

/*** R
testvec <- 1:12
mysplit(testvec, 2, 3)
*/

输出

> Rcpp::sourceCpp("~/git/stackoverflow/68858728/answer.cpp")

> testvec <- 1:12

> mysplit(testvec, 2, 3)
[[1]]
[1] 1 2 3 7 8 9

[[2]]
[1]  4  5  6 10 11 12

> 

原来的问题有一个小错误,我们不需要调用gl();只需要这两个标量。

【讨论】:

  • 显然你是Rcpp的作者,它可以让你在R中使用C++。非常有用。
  • 我是 Rcpp 背后团队的一员,是的。事实上,它有它的用途。更多 here on GitHub例如 here on SO.
【解决方案2】:

试试这个,它会创建一个向量向量,其中包含交替块中的源元素:

#include <iostream>
#include <vector>

template<typename T>
std::vector<std::vector<T>> split(std::vector<T> nums, int n, int size)
{
    std::vector<std::vector<T>> result(n);
    int i = 0;
    auto beg = nums.cbegin();
    auto end = nums.cend();

    while (beg != nums.cend()) {
        //get end iterator safely
        auto next = std::distance(beg, end) >= size ? beg + size : end;
        //insert into result
        result[i].insert(result[i].end(), beg, next);
        //advance iterator
        beg = next;
        i = (i + 1) % n;
    }

    return result;
}

int main()
{
    std::vector<int> vnums = { 1, 2, 3,  4,  5,  6,  7,  8,  9, 10, 11, 12 };
    auto vectors = split(vnums, 2, 3);

    for (const auto& v : vectors)
    {
        for (auto num : v) {
            std::cout << num << " ";
        }
        std::cout << std::endl;
    }
}

Demo

【讨论】:

  • 很好,虽然不能从 R 中使用。我会在一个变体中添加它,感谢您。
  • 感谢这两个解决方案,但特别感谢您提供了一个可从 R 中使用的解决方案,而我正是需要的。
猜你喜欢
  • 2012-07-04
  • 1970-01-01
  • 2011-12-24
  • 2020-11-02
  • 2022-06-11
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
相关资源
最近更新 更多