I am using the deSolve package successfully and I think it is probably your best option if you want to do things in R. <br>The most widespread solver in C++ is the Sundials Suite as far as I know:<br><a href="https://computation.llnl.gov/casc/sundials/main.html">https://computation.llnl.gov/casc/sundials/main.html</a><br>
I think someone starting writing an interface to R, but it only covers the really basic functions.<br>Regards,<br>Jonas<br><br><div class="gmail_quote">On Wed, Oct 5, 2011 at 7:53 PM, Douglas Bates <span dir="ltr"><<a href="mailto:bates@stat.wisc.edu">bates@stat.wisc.edu</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">On Wed, Oct 5, 2011 at 12:25 PM, Richard Downe <<a href="mailto:richard-downe@uiowa.edu">richard-downe@uiowa.edu</a>> wrote:<br>
> On 10/04/2011 05:05 PM, Douglas Bates wrote:<br>
>><br>
>> I am working with pharmacokinetic/pharmacodynamic models, some of<br>
>> which can only be represented as the solutions of ordinary<br>
>> differential equations. I am looking for ODE sofware written in C++.<br>
>> The deSolve package for R has lots of wonderful functions but I may<br>
>> want to take things apart a bit and I really don't have the patience<br>
>> for wading through Fortran code to try to discover what it is doing.<br>
>><br>
>> I can, of course, google likely phrases but if anyone has experience<br>
>> with such solvers and can make any recommendations I would appreciate<br>
>> it (off-list responses are probably better than on-list)<br>
><br>
</div><div class="im">> I know it doesn't completely get you out of the fortran woods, but<br>
> LAPACK/BLAS have a very usable c++ wrapper (I haven't used it in a couple<br>
> years, but it was last I checked) called lapackpp. It's still fortran under<br>
> the hood, but the exact mechanics of how it works are fairly well publicly<br>
> documented and understood since it's been around for so long.<br>
<br>
</div>Actually I am looking for code to solve ordinary differential<br>
equations, not linear equations. We already have two packages,<br>
RcppArmadillo and RcppEigen, that interface to C++ code providing<br>
BLAS/LAPACK functionality (and much more).<br>
<div><div></div><div class="h5">_______________________________________________<br>
Rcpp-devel mailing list<br>
<a href="mailto:Rcpp-devel@lists.r-forge.r-project.org">Rcpp-devel@lists.r-forge.r-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-project.org/cgi-bin/mailman/listinfo/rcpp-devel</a><br>
</div></div></blockquote></div><br>