R与C++交互：基于指针实现cummin()函数的编译问题求解

问题背景 #

需要模仿R基础包的cummin()函数，用C++结合cpp11编写可编译的代码，同时解决现有代码的冗余问题，修复指针实现和std::min调用时的编译错误。

现有可运行但存在冗余的代码：

#include <cpp11.hpp>

using namespace cpp11;
namespace writable = cpp11::writable;

[[cpp11::register]] cpp11::doubles cummin_cpp_(cpp11::doubles x) {
  int n = x.size();

  // create a vector of 0s
  writable::doubles out(n);

  for (int i = 0; i < n; ++i) {
    out[i] = 0.0;
  }

  out[0] += x[0];

  for (int i = 1; i < n; ++i) {
    if (x[i] < out[i - 1]) {
      out[i] += x[i];
    } else {
      out[i] += out[i - 1];
    }
  }
  return out;
}

冗余点：

• 提前将输出向量填充为0
• 使用+=执行无意义的"0+目标值"操作

尝试指针实现时出现的编译错误代码：

[[cpp11::register]] cpp11::doubles cummin_cpp_(cpp11::doubles x) {
  int n = x.size();

  writable::doubles out(n);

  out[0] = x[0];

  for (int i = 1; i < n; ++i) {
    if (x[i] < out[i - 1]) {
      out[i] = x[i];
    } else {
      // code.cpp:157:18: error: no match for ‘operator*’ (operand type is ‘cpp11::writable::r_vector<double>::proxy’)
      out[i] = &(*out[i - 1]);
    }
  }
  return out;
}

错误原因：out[i-1]返回的是cpp11的代理对象（proxy），而非原始double值，无法直接对其进行取地址和解引用操作，类型不匹配导致编译失败。

使用std::min改写时的错误代码：

[[cpp11::register]] cpp11::doubles cummin_cpp_(cpp11::doubles x) {
  int n = x.size();
  writable::doubles out(n);

  out[0] = x[0];
  for (int i = 1; i < n; ++i) {
    out[i] = std::min(out[i - 1], x[i]);
  }
  return out;
}

错误信息：

error: no matching function for call to ‘min(cpp11::writable::r_vector<double>::proxy, double)’
     150 |     out[i] = std::min(out[i - 1], x[i]);

原因：std::min要求两个参数类型完全一致，而out[i-1]是代理对象，x[i]是double，类型不兼容。

解决方案：基于指针的正确实现 #

利用cpp11容器的data()方法获取底层原始指针，直接操作内存中的double值，绕过代理对象的类型问题，同时消除冗余操作：

#include <cpp11.hpp>

using namespace cpp11;
namespace writable = cpp11::writable;

[[cpp11::register]] cpp11::doubles cummin_cpp_(cpp11::doubles x) {
  int n = x.size();
  if (n == 0) {
    return writable::doubles(); // 处理空输入的边界情况
  }

  writable::doubles out(n);
  // 获取输入和输出向量的底层double指针
  const double* x_ptr = x.data();
  double* out_ptr = out.data();

  // 初始化第一个元素
  out_ptr[0] = x_ptr[0];

  // 用指针遍历计算累积最小值
  for (int i = 1; i < n; ++i) {
    out_ptr[i] = (x_ptr[i] < out_ptr[i-1]) ? x_ptr[i] : out_ptr[i-1];
  }

  return out;
}

代码说明： #

• 使用data()方法直接获取底层指针，操作原始double数据，避开cpp11的代理对象，彻底解决类型不匹配问题
• 无需提前填充输出向量为0，直接通过指针赋值，消除冗余操作
• 添加空输入的边界处理，提升代码鲁棒性
• 用三元运算符简化条件判断，代码更简洁

补充：修复std::min版本的错误 #

如果想保留std::min的写法，只需将代理对象显式转换为double类型，确保两个参数类型一致：

[[cpp11::register]] cpp11::doubles cummin_cpp_(cpp11::doubles x) {
  int n = x.size();
  if (n == 0) {
    return writable::doubles();
  }

  writable::doubles out(n);
  out[0] = x[0];

  for (int i = 1; i < n; ++i) {
    // 将代理对象转换为double，与x[i]类型匹配
    out[i] = std::min(static_cast<double>(out[i-1]), x[i]);
  }

  return out;
}

内容的提问来源于stack exchange，提问作者pachadotdev