Question

我正在将R函数转换为Rcpp函数。一切都很好，但我在使用标准的最大和最小功能方面遇到了困难。以下代码：

#include <math.h>
#include <RcppArmadillo.h>
#include <algorithm>
#include <iostream>

// [[Rcpp::depends(RcppArmadillo)]]

using namespace Rcpp;
using namespace arma;
using namespace std;

double f4RandomC(double x, double a, double b) {
  double out, temp;

  temp=(log( (a*(1-x)) / (a+x) )) /log(b) ;
  out= std::min(1,temp );
  return out;
}

返回错误＆＃34;没有matchinf函数用于调用min（int，＆amp; double）。如果可能的话，我会想要使用std :: library min函数

Answer 1

只需将std::min(1,temp)更改为std::min(1.0,temp)：

#include <cmath>
#include <Rcpp.h>
// [[Rcpp::export]]
double f4RandomC(double x, double a, double b) {
  double out, temp;

  temp =(log( (a*(1-x)) / (a+x) )) /log(b) ;
  out = std::min(1.0,temp );
  return out;
}

我假设这与std::min

的模板定义有关

template <class T> const T& min (const T& a, const T& b);

仅根据一种类型（T）定义，而您传递了两种数据类型（int和double）。

或者由于您只是比较两个值，因此可以通过将std::min替换为三元运算符（?:）来更简洁地执行此操作：

double f4RandomC(double x, double a, double b) {
  double temp;

  temp =(log( (a*(1-x)) / (a+x) )) /log(b) ;
  return temp < 1 ? temp : 1;
}

我猜测operator<的类型推导比std::min更灵活。

std::min的其他两个选项：

// [[Rcpp::export]]
double f4RandomC2(double x, double a, double b) {
  double out, temp;
  int z = 1;
  temp =(log( (a*(1-x)) / (a+x) )) /log(b) ;
  out = std::min( static_cast<double>(z),temp );
  return out;
}

// [[Rcpp::export]]
double f4RandomC3(double x, double a, double b) {
  double out, temp;
  int z = 1;
  temp =(log( (a*(1-x)) / (a+x) )) /log(b) ;
  out = std::min<double>( z,temp );
  return out;
}

虽然在这种情况下，将1更改为1.0比仅（不必要地）定义int z更容易，只是稍后将其转换为双倍。

通过阅读函数/类定义（与大多数编程语言一样），您可以学到很多东西 - cplusplus.com和cppreference.com是相当标准的资源 - 而且它往往会使编译器错误看起来更加神秘

Rcpp程序中的最小值和最大值

1 个答案: