<p>Thanks. I guess the error in<br>
double Radical::DoRadical2D(const mat&amp; matX)<br>
{<br>
 .....<br>
    matJacobi(0, 0) = cosTheta;<br>
 ...<br>
 }<br>
is the same. The workaround for this is to use<br>
matJacobi.at(0,0) = cosTheta</p>

<p style="font-size:small;-webkit-text-size-adjust:none;color:#666;">&mdash;<br>Reply to this email directly or <a href="https://github.com/mlpack/mlpack/issues/424#issuecomment-84268273">view it on GitHub</a>.<img alt="" height="1" src="https://github.com/notifications/beacon/AJ4bFLNQ66qXRCH5oVO8QwbSi02igZNaks5n3Q0bgaJpZM4DyYOG.gif" width="1" /></p>
<div itemscope itemtype="http://schema.org/EmailMessage">
  <div itemprop="action" itemscope itemtype="http://schema.org/ViewAction">
    <link itemprop="url" href="https://github.com/mlpack/mlpack/issues/424#issuecomment-84268273"></link>
    <meta itemprop="name" content="View Issue"></meta>
  </div>
  <meta itemprop="description" content="View this Issue on GitHub"></meta>
</div>