[go: up one dir, main page]

Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Allow implicit conversion to Rcpp::RObject #114

Open
DavisVaughan opened this issue May 6, 2019 · 4 comments
Open

Allow implicit conversion to Rcpp::RObject #114

DavisVaughan opened this issue May 6, 2019 · 4 comments

Comments

@DavisVaughan
Copy link
Contributor
DavisVaughan commented May 6, 2019

Currently, something like this is allowed with no problems:

SEXP rray__exp_impl(const xt::rarray<rlogical>& x) {
  xt::rarray<double> res = xt::exp(x);
  return res;
}

The (potentially safer) approach would be to return an Rcpp::RObject like this:

Rcpp::RObject rray__exp_impl(const xt::rarray<rlogical>& x) {
  xt::rarray<double> res = xt::exp(x);
  return res;
}

But this gives me a "no viable conversion" error as it does not know how to convert to RObject from xt::rarray<double>. I can do something like this, but it is just a bit annoying:

Rcpp::RObject rray__exp_impl(const xt::rarray<rlogical>& x) {
  xt::rarray<double> res = xt::exp(x);
  return Rcpp::as<Rcpp::RObject>(res);
}

Do you know if its possible to add this implicit conversion?

PS: If you not aware, RObjects are wrappers around SEXPs that Rcpp added that help automatically manage the underlying storage properties (they are also the base class that all core Rcpp classes inherit from). So you don't have to call PROTECT() and UNPROTECT() like you would have to do on the raw R SEXP objects

@wolfv
Copy link
Member
wolfv commented May 7, 2019

Hey, yes, it is definitely possible to add an implicit conversion like that. You'd just need to add a function like this:

rcontainer ... 
{
    operator Rcpp::RObject () {
         return Rcpp::as<Rcpp::RObject>(*this);
    }
}

However, I am unsure where to extend xtensor-r to support this. There is also a chance that Rcpp has an internal mechanism to do this (like they have for the SEXP casting, if you take a look at rcpp_extensions header). Maybe we can ask them?

@DavisVaughan
Copy link
Contributor Author
DavisVaughan commented May 7, 2019

I know that

  • as<T> is for SEXP -> C++ type
  • wrap() is for C++ type -> SEXP

as summarized nicely in the first code block here
http://dirk.eddelbuettel.com/code/rcpp/Rcpp-extending.pdf

But I guess this is technically rarray<double> -> RObject which are 2 C++ types, so maybe it shouldn't go through that mechanism.

I think this last example (repeated below) might only be working because rcontainer inherits from PreserveStorage which I think has the capabilities to be cast like this because it "looks like" a SEXP. I'm not really sure.

Rcpp::RObject rray__exp_impl(const xt::rarray<rlogical>& x) {
  xt::rarray<double> res = xt::exp(x);
  return Rcpp::as<Rcpp::RObject>(res);
}

I can ask someone if you aren't sure either

@wolfv
Copy link
Member
wolfv commented May 7, 2019

This works because you're explicitly asking for this conversion.

C++ allows one implicit conversion only. So implicitly you can convert to a SEXP like this: SEXP msxp = x; (where x is an rarray).

Now RObject has (most likely) a constructor from SEXP (which would require two conversions for implicit conversion) (rarray -> SEXP -> RObject). If you ask for it explicitly it would also work (RObject( (SEXP) x))

Does explicitly asking for RObject(x) work?

Btw. I think rarray ... use the exact same mechanisms as RObject so in theory at least it shouldn't matter that your object stays an rarray. Out of curiosity, what benefits do you get from the RObject cast?

@DavisVaughan
Copy link
Contributor Author
DavisVaughan commented May 7, 2019

Yea, both of these work:

// [[Rcpp::depends(xtensor)]]
// [[Rcpp::plugins(cpp14)]]

#include <xtensor-r/rarray.hpp>

// [[Rcpp::export(rng = false)]]
Rcpp::RObject rray__exp_explicit(const xt::rarray<rlogical>& x) {
  xt::rarray<double> res = xt::exp(x);
  return Rcpp::RObject((SEXP) res);
}

// [[Rcpp::export(rng = false)]]
Rcpp::RObject rray__exp_direct(const xt::rarray<rlogical>& x) {
  xt::rarray<double> res = xt::exp(x);
  return Rcpp::RObject(res);
}

As to why I'd do this, it's for the fact that RObject abstracts over the fact that I return a double, vs logical, vs integer, so the return type can be a bit more flexible.

A few functions, like the identity function +x in base R, return:

  • doubles if given doubles
  • integers if given integers
  • integers if given logicals

So I wanted to be able to do:

template <typename T>
xt::rarray<T> rray__identity_impl(const xt::rarray<T>& x) {
  return xt::operator+(x);
}

But that won't work for logicals because the return type should be different than the input type. So I can do something like:

template <typename T>
Rcpp::RObject rray__identity_impl(const xt::rarray<T>& x) {
  return xt::operator+(x);
}

Rcpp::RObject rray__identity_impl(const xt::rarray<rlogical>& x) {
  xt::rarray<int> res = xt::operator+(x);
  return Rcpp::as<Rcpp::RObject>(res);
}

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants