<div dir="ltr">On Mon, Sep 30, 2013 at 4:42 PM, Douglas Bates <span dir="ltr"><<a href="mailto:bates@stat.wisc.edu" target="_blank">bates@stat.wisc.edu</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div class="im">On Mon, Sep 30, 2013 at 4:16 PM, Dirk Eddelbuettel <span dir="ltr"><<a href="mailto:edd@debian.org" target="_blank">edd@debian.org</a>></span> wrote:<br>
</div><div class="gmail_extra"><div class="gmail_quote"><div class="im">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><br>
Hi Martyn,<br>
<div><br>
On 30 September 2013 at 21:02, Martyn Plummer wrote:<br>
| I've already been asked to look into this by Brian.  I can't say I've rolled up my sleeves very far, but the problem appears to be in Eigen itself.<br>
<br>
</div>Cool!  Good to know.<br>
<div><br>
| I should clarify that I'm not volunteering to maintain RcppEigen, but I will provide patches as and when.<br>
<br>
</div>That too is appreciated!<br>
<br>
The Eigen code on CRAN has not been updated in a while. The version number<br>
(0.3.1.2.1) reflect both upstream Eigen 3.1.2 (two patch releases behind<br>
Eigen in the 3.1.* series which is now at 3.1.4) and the mandatory change for<br>
JSS once the paper was accepted.<br>
<br>
Changing 3.1.2 to 3.1.3/3.1.4 may help. Or it may not.  Maybe 3.2.0 upstream<br>
will help.<br>
<br></blockquote></div><div>Changing to 3.1.4 or, better, 3.2, would be a good idea.  In particular, the 3.2 version provides Cholesky, LU and QR sparse factorizations that do not depend on SuiteSparse (the package of C libraries of which CHOLMOD is part).</div>

<div><br></div><div>Unfortunately, the current maintainer played a little fast and loose by adding a couple of wrinkles to the Eigen code itself and using these wrinkles in lme4.</div><div><br></div><div>I had been hoping to put off the switch to Eigen 3.2 for a while longer as I am facing some deadline pressures on another project.  But perhaps it would be best to create an experimental branch using Eigen 3.2.  Best case, the fixes for lme4 are easy.  Worst case ...</div>

<div><br></div><div>Ben's message came through while I was writing this.  If you are interested look at the last link he gives</div><div><br></div><div><a href="http://permalink.gmane.org/gmane.comp.lib.eigen/991" target="_blank">http://permalink.gmane.org/gmane.comp.lib.eigen/991</a></div>

<div><br></div><div>and you will see that the sort of error messages showing up in the RcppEigen log were present in 2010 and the reaction of the Eigen developers is that the compiler is complaining about code that follows the standard and is compilable on other systems.</div>
</div></div></div></blockquote><div></div></div><br></div><div class="gmail_extra">It seems that a good place for someone with access to a Solaris system to start is to download versions 3.1.4 and 3.2 of Eigen from <a href="http://eigen.tuxfamily.org">http://eigen.tuxfamily.org</a> and see if the simple program </div>
<div class="gmail_extra"><pre style="line-height:19px;text-align:justify">$ cat test.cc
#include "eigen/Eigen/Eigen"
using namespace Eigen;
int main() {
 return 0;
}
$
$ CC test.cc</pre><pre style="line-height:19px;text-align:justify"><span style="font-family:arial;line-height:normal;text-align:start;white-space:normal">compiles with one of those versions.  If that doesn't compile the way forward will be difficult.</span></pre>
</div></div>