Hi,<br>Here is the sample code:<br>// [[Rcpp::export]]<br>List betahat(Function ker, double t0, NumericMatrix Xr, NumericMatrix yr, NumericVector tr, double h, int m) {<br> int n = Xr.nrow(), p = Xr.ncol();<br>
arma::mat X(Xr.begin(), n, p, false);<br> arma::mat y(yr.begin(), n, 1, false);<br> arma::colvec t(tr.begin(), tr.size(), false);<br> arma::mat T = X;<br> T.each_col() %= (t-t0)/h;<br> arma::vec K =as<arma::vec>(ker(tr-t0,h))/<div id=":18y">
m;<br>
double L1 = arma::accu(K%X.col(0)%X.col(0));<br> double L2 = arma::accu(K%X.col(0)%X.col(1));<br> double L3 = arma::accu(K%X.col(1)%X.col(1));<br> double L4 = arma::accu(K%X.col(0)%T.col(0));<br> double L5 = arma::accu(K%X.col(1)%T.col(0));<br>
double L6 = arma::accu(K%X.col(1)%T.col(1));<br> double L7 = arma::accu(K%T.col(0)%T.col(0));<br> double L8 = arma::accu(K%T.col(0)%T.col(1));<br> double L9 = arma::accu(K%T.col(1)%T.col(1));<br> double R1 = arma::accu(K%X.col(0)%y);<br>
double R2 = arma::accu(K%X.col(1)%y);<br> double R3 = arma::accu(K%T.col(0)%y);<br> double R4 = arma::accu(K%T.col(1)%y);<br> arma::mat L(2*p,2*p);<br> L(0,0)=L1;L(0,1)=L2;L(0,2)=L4;L(0,3)=L5;<br> L(1,0)=L2;L(1,1)=L3;L(1,2)=L5;L(1,3)=L6;<br>
L(2,0)=L4;L(2,1)=L5;L(2,2)=L7;L(2,3)=L8;<br> L(3,0)=L5;L(3,1)=L6;L(3,2)=L8;L(3,3)=L9;<br> arma::mat R(2*p,1);<br> R(0,0)=R1;R(1,0)=R2;R(2,0)=R3;R(3,0)=R4;<br> arma::vec betahat = arma::solve(L,R);<div class="yj6qo ajU">
<div id=":19c" class="ajR" tabindex="0"><img class="ajT" src="https://mail.google.com/mail/u/0/images/cleardot.gif"></div></div><div class="im HOEnZb adL"><br> arma::colvec betahat0(betahat.begin(),betahat.size()/2,false);<br>
return List::create(Named("betahat") = betahat0);<br>}<br><br>I will call this function repeatedly, and at some point, it went wrong with the following error:<br>Error in betahat(ker, x, X, y, t, h, m) :<br> promise already under evaluation: recursive default argument reference or earlier problems?<br>
Calls: system.time ... apply -> FUN -> betahat -> .External -> cpp_exception<br>Execution halted<br><br>I have no idea what's this error. I am just wondering whether I need to free some memory in this code. Thanks.<br>
</div></div><br clear="all"><div><div>Best wishes!</div><div> </div><div>Honglang Wang</div><div> </div><div>Office C402 Wells Hall</div><div>Department of Statistics and Probability</div><div>Michigan State University</div>
<div>1579 I Spartan Village, East Lansing, MI 48823</div><div><a href="mailto:wangho16@msu.edu" target="_blank">wangho16@msu.edu</a></div></div><br>
<br><br><div class="gmail_quote">On Thu, Dec 13, 2012 at 4:57 PM, <span dir="ltr"><<a href="mailto:romain@r-enthusiasts.com" target="_blank">romain@r-enthusiasts.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Hello,<br>
<br>
What memory do you want to free ? Can you give an example of code where there is memory you'd like to "free" ?<br>
<br>
Romain<br>
<br>
Le 2012-12-13 22:52, Honglang Wang a écrit :<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Dear All,<br>
How to free memory in Rcpp? What's the command? Thanks.<br>
<br>
Best wishes!<br>
<br>
Honglang Wang<br>
<br>
Office C402 Wells Hall<br>
Department of Statistics and Probability<br>
Michigan State University<br>
1579 I Spartan Village, East Lansing, MI 48823<br>
</div></div><a href="mailto:wangho16@msu.edu" target="_blank">wangho16@msu.edu</a> [1]<br>
<br>
<br>
Links:<br>
------<br>
[1] mailto:<a href="mailto:wangho16@msu.edu" target="_blank">wangho16@msu.edu</a><br>
</blockquote>
<br>
______________________________<u></u>_________________<br>
Rcpp-devel mailing list<br>
<a href="mailto:Rcpp-devel@lists.r-forge.r-project.org" target="_blank">Rcpp-devel@lists.r-forge.r-<u></u>project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel" target="_blank">https://lists.r-forge.r-<u></u>project.org/cgi-bin/mailman/<u></u>listinfo/rcpp-devel</a></blockquote></div><br>