<br><br><div class="gmail_quote">On Sat, Aug 18, 2012 at 3:46 PM, Dirk Eddelbuettel <span dir="ltr"><<a href="mailto:edd@debian.org" target="_blank">edd@debian.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div class="im"><br>
On 18 August 2012 at 14:08, Stephen J. Barr wrote:<br>
| Hello,<br>
|<br>
| I have a quick design question. I am working on a program that uses RInside and<br>
| RcppEigen, and I want to expose a function back to R. Specifically, I have a<br>
| C++ function with the following signature:<br>
|<br>
| double nll_singleton(<br>
|     const VectorXd& theta,<br>
|     const MatrixXd& ua,<br>
|     const MatrixXd& ub,<br>
|     const MatrixXd& Xmat,<br>
|     const MatrixXd& Ymat);<br>
|<br>
| I would like to run R's optim over this function, where 'theta' is the vector<br>
| that I am optimizing over. Would I use RCpp::InternalFunction for this?<br>
<br>
</div>Per rinside_sample9.cpp, possibly.<br>
<div class="im"><br>
| nll_singleton is going to call MPI and I would rather do it this way than use<br>
| mpi over R, if this is possible.<br>
<br>
</div>Right.<br></blockquote><div>What I meant was, rather than recode everything in R and use Rmpi.</div><div><br></div><div>I guess my question was, say I use RCpp::InternalFunction with nll_singleton, can I just pass in a vector and a few matrices like:</div>

<div><br></div><div>>R.parseEval("nll_singleton(c(1,2,3), matrix(1:10, nrow=3), ...etc..")</div><div><br></div><div>?</div><div><br></div><div>I would think that I would have to define some mapping back to the Eigen matrices.</div>

<div><br></div><div>Thanks,</div><div>Stephen</div><div><br></div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Doug Bates and have both (and indepdently) used Rcpp::XPtr to pass function<br>
pointers around to C-level optimisation functions.  My use case in my "port"<br>
of DEoption to Rcpp, ie in RcppDE.<br>
<br>
Good luck,  Dirk<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Dirk Eddelbuettel | <a href="mailto:edd@debian.org">edd@debian.org</a> | <a href="http://dirk.eddelbuettel.com" target="_blank">http://dirk.eddelbuettel.com</a><br>
</font></span></blockquote></div><br>