<p>Thanks. I guess the error in<br>
double Radical::DoRadical2D(const mat& 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;">—<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>