<font><font face="verdana,sans-serif">I think it is useful to put this finding on the mailing list for future reference.<br></font></font><br><div class="gmail_quote">On Mon, Feb 6, 2012 at 18:28, Maarten Beek <span dir="ltr">&lt;<a href="mailto:beekmaarten@yahoo.com">beekmaarten@yahoo.com</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-size:12pt;font-family:times new roman,new york,times,serif"><div><span>Hate to admit, but I found a bug in my code as well. Same bug as Dzenan. Both vtkPoints objects were created in similar for-loops. So one of the arrays was not totally filled. Gives very unpredictable behavior (sometimes not even crashing).</span></div>

<div><span>Switching the point sets however always caused crashes and should have been the clue I shouldn&#39;t have missed.</span></div><div><br></div><div>Remain the facts that vtkMath::vtkJacobiN doesn&#39;t really handle this error nicely and that the algorithm used is potentially unstable.<span class="HOEnZb"><font color="#888888"><br>

</font></span></div><span class="HOEnZb"><font color="#888888"><div><br><span></span></div><div><span>Maarten<br></span></div></font></span><div class="hm HOEnZb"><br></div><div style="font-family:times new roman,new york,times,serif;font-size:12pt">

<div class="hm HOEnZb"> </div><div style="font-family:times new roman,new york,times,serif;font-size:12pt"><div class="hm HOEnZb"> </div><div dir="ltr"><div class="hm HOEnZb"> </div><font face="Arial"><div class="hm HOEnZb">

 <hr size="1">  <b><span style="font-weight:bold">From:</span></b> Dženan Zukić &lt;<a href="mailto:dzenanz@gmail.com" target="_blank">dzenanz@gmail.com</a>&gt;<br> <b><span style="font-weight:bold">To:</span></b> David Gobbi &lt;<a href="mailto:david.gobbi@gmail.com" target="_blank">david.gobbi@gmail.com</a>&gt; <br>

<b><span style="font-weight:bold">Cc:</span></b> Maarten Beek &lt;<a href="mailto:beekmaarten@yahoo.com" target="_blank">beekmaarten@yahoo.com</a>&gt; <br> <b><span style="font-weight:bold">Sent:</span></b> Friday, February 3, 2012 10:55:58 AM</div>

<div><div class="h5"><br> <b><span style="font-weight:bold">Subject:</span></b> Re: [vtkusers] ICP Error: vtkMath::Jacobi: Error extracting eigenfunctions<br> </div></div></font> </div><div><div class="h5"> <br><div><font><font face="verdana,sans-serif">In trying to create a self-contained test program, I discovered a bug in my own code. I was supplying some junk points to the vtkLandmarkTransform (one points array not completely filled with values).<br>



<br>Now the execution in release mode doesn&#39;t give any warnings either. Maybe something similar happened to you Maarten?<br></font></font><br><div>2012/2/1 David Gobbi <span dir="ltr">&lt;<a rel="nofollow" href="mailto:david.gobbi@gmail.com" target="_blank">david.gobbi@gmail.com</a>&gt;</span><br>



<blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Would it be possible for you to send me a short, simple program that<br>
includes vtkLandmarkTransform that triggers the error?  Or can you try<br>
to modify vtkMath::JacobiN in such a way that the error goes away,<br>
using the Numerical Recipes discussion on the following link as a<br>
guideline? <a rel="nofollow" href="http://www.nr.com/forum/showthread.php?p=4913" target="_blank">http://www.nr.com/forum/showthread.php?p=4913</a><br>
<br>
There isn&#39;t much that I can do on my side unless I have some way of<br>
reproducing the error...<br>
<br>
2012/2/1 Dženan Zukić &lt;<a rel="nofollow" href="mailto:dzenanz@gmail.com" target="_blank">dzenanz@gmail.com</a>&gt;:<br>
<div><div>&gt; Same: problem is not encountered.<br>
&gt;<br>
&gt; 2012/2/1 David Gobbi &lt;<a rel="nofollow" href="mailto:david.gobbi@gmail.com" target="_blank">david.gobbi@gmail.com</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; Can you modify my test program so that main() uses vtkMath::JacobiN<br>
&gt;&gt; instead of using its own JacobiN, and then try it again?<br>
&gt;&gt;<br>
&gt;&gt; 2012/2/1 Dženan Zukić &lt;<a rel="nofollow" href="mailto:dzenanz@gmail.com" target="_blank">dzenanz@gmail.com</a>&gt;:<br>
&gt;&gt; &gt; This program does not crash in either debug or release mode. The<br>
&gt;&gt; &gt; optimization options are the same (as CMake sets them).<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; I have modified it to read input for the files I sent you, but it goes<br>
&gt;&gt; &gt; through them without errors. I guess that whatever is happening, is not<br>
&gt;&gt; &gt; triggered in this simplified example.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; 2012/1/31 David Gobbi &lt;<a rel="nofollow" href="mailto:david.gobbi@gmail.com" target="_blank">david.gobbi@gmail.com</a>&gt;<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; I haven&#39;t been able to reproduce the problem, but at work all I have is<br>
&gt;&gt; &gt;&gt; a Mac and the optimizations it does are probably different from MSVC.<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; Can you try different matrices with the attached program and see if<br>
&gt;&gt; &gt;&gt; you can get it to fail?<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;  - David<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; 2012/1/31 Dženan Zukić &lt;<a rel="nofollow" href="mailto:dzenanz@gmail.com" target="_blank">dzenanz@gmail.com</a>&gt;:<br>
&gt;&gt; &gt;&gt; &gt; Here is an output from non-release version, same input as before.<br>
&gt;&gt; &gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt; &gt; 2012/1/31 Dženan Zukić &lt;<a rel="nofollow" href="mailto:dzenanz@gmail.com" target="_blank">dzenanz@gmail.com</a>&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt; Here are 3 different runs of my program, identical input.<br>
&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt; 2012/1/31 David Gobbi &lt;<a rel="nofollow" href="mailto:david.gobbi@gmail.com" target="_blank">david.gobbi@gmail.com</a>&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; Thanks for the info.  Can you post your compiler optimization flags<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; to<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; the list?  There is definitely a problem with JacobiN that needs to<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; be<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; fixed.  The fix should be very easy to do if we can get a matrix<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; that<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; reliably causes JacobiN to fail when certain optimization flags are<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; used, so if either of you can print out such a matrix (with all<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; numbers printed with a precision of 20 digits), it would be a big<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; help.<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt;  - David<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; 2012/1/31 Dženan Zukić &lt;<a rel="nofollow" href="mailto:dzenanz@gmail.com" target="_blank">dzenanz@gmail.com</a>&gt;:<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; I also saw this warning recently, from the block that &quot;never gets<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; called&quot;.<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; And it only happened in release version of my program.<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; On Tue, Jan 31, 2012 at 00:48, Maarten Beek<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; &lt;<a rel="nofollow" href="mailto:beekmaarten@yahoo.com" target="_blank">beekmaarten@yahoo.com</a>&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt; wrote:<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Hi David,<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Ai, seems these crashes are hard to get rid off....<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I currently use VTK within Amira modules.<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I changed the type of vtkPoints from the default VTK_FLOAT to<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; VTK_DOUBLE<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; when I convert an Amira object to a VTK object. This reduces the<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; occurrences<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; of the errors, but doesn&#39;t solve the issue.<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; (Ps: why is the default in vtkPoints still float, while<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; everywhere<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; else<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; floats were changed into doubles a long while ago?)<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I don&#39;t really know much about different optimizations, so I<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; just<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; copy-paste what I found in MSVS:<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I build my VTK with:<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /O2 /Ob2 /FD /EHsc /MD /Fo&quot;vtkCommon.dir\Release\\&quot;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /Fd&quot;G:\vtk-5.6.0\build_vs90_noQt\bin\Release/vtkCommon.pdb&quot; /W4<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /nologo /c<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /TP /errorReport:prompt/Zm 1000 and a couple of &#39;/I&#39;-s and<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &#39;/D&#39;-s<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I build my Amira modules with:<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /O2 /Ob1 /Oi /Oy /GF /FD /EHsc /MD /GS- /Gy /fp:fast /openmp<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /Fo&quot;./../../obj/arch-Win64VC9-Optimize/GraphPackage/&quot;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /Fd&quot;./../.././bin/arch-Win64VC9-Optimize/GraphPackage.pdb&quot; /W3<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /nologo<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /c<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /Zi /wd4068 /wd4305 /wd4018 /wd4244 /wd4800 /wd4275 /wd4251<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; /errorReport:prompt and a couple of &#39;/I&#39;-s and &#39;/D&#39;-s<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I am working on getting vtkMath::JacobiN to print the matrix in<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; the<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; if-loop that &#39;gets never called&#39; using vtkGenericWarningMacro to<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; avoid<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; having to include Amira stuff. Having a hard time getting that<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; to<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; work. (<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; char* test = &quot;hello&quot;;vtkGenericWarningMacro( test ); won&#39;t do<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; the<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; trick).<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Any suggestions?<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Maarten<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; ________________________________<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; From: David Gobbi &lt;<a rel="nofollow" href="mailto:david.gobbi@gmail.com" target="_blank">david.gobbi@gmail.com</a>&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; To: Maarten Beek &lt;<a rel="nofollow" href="mailto:beekmaarten@yahoo.com" target="_blank">beekmaarten@yahoo.com</a>&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Cc: VTK list &lt;<a rel="nofollow" href="mailto:vtkusers@vtk.org" target="_blank">vtkusers@vtk.org</a>&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Sent: Monday, January 30, 2012 2:31:38 PM<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Subject: Re: [vtkusers] ICP Error: vtkMath::Jacobi: Error<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; extracting<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; eigenfunctions<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Hi Maarten,<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; I think there are two issues here, one is that<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; vtkLandmarkTransform<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; should have a check as you have described, but another<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; potentially<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; serious issue is that, as far as I understand it, JacobiN should<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; never<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; fail.  I did a google search and found a message from 2002 by<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; the<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Numerical Recipes authors stating that the NR Jacobi code can<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; become<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; numerically unstable due to certain compiler optimizations.  I&#39;m<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; guessing that certain compiler optimizations could similarly<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; make<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; the<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; VTK Jacobi code become unstable.<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; <a rel="nofollow" href="http://www.nr.com/forum/showthread.php?p=4913" target="_blank">http://www.nr.com/forum/showthread.php?p=4913</a><br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; Can you provide info about your compiler/flags etc?  Also, is it<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; possible to get your code to print a matrix that causes JacobiN<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; to<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; fail?<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; - David<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; On Mon, Jan 30, 2012 at 9:44 AM, Maarten Beek<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &lt;<a rel="nofollow" href="mailto:beekmaarten@yahoo.com" target="_blank">beekmaarten@yahoo.com</a>&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; wrote:<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; Hi all,<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; I get the following warning when using ICP:<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; vtkMath::Jacobi: Error extracting eigenfunctions<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; It comes from vtkMath::vtkJacobiN that calculates the<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; eigenvalues,<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; -vectors<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; of a matrix.<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; The function returns 0, when this error occurs (otherwise 1),<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; but<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; vtkLandmarkTransform doesn&#39;t check for the return value and<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; thus<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; uses<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; undefined values for the eigenvalues, -vectors.<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; This means my point set becomes invalid.<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; I believe the way to avoid this error is to &#39;giggle&#39; the<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; coordinates<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; a<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; little when it happens, but how can I intervene when<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; vtkLandmarkTransform<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; stubbornly uses whatever values it gets from<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; vtkMath::vtkJacobiN?<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; Is there a way to intercept the error message?<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; I also noticed that if I apply ICP on A and B1 resulting in<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; the<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; error<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; message, applying ICP on A and B2 (loading the same data file<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; that<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; gives<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; me<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; B1) works fine. I don&#39;t see vtkMath::vtkJacobiN using a random<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; number or<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; doing something similar, which makes me believe the error<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; message<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; I<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; get<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; is<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; due to floating point accuracy? Even more a reason to think<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; that<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; &#39;giggling&#39;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; might work (if I were able to intervene before the invalid<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; transform<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; is<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; applied to my data...)<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;&gt; &gt;&gt; &gt; Thanks - Maarten<br>
&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;&gt;<br>
&gt;&gt; &gt;&gt; &gt;<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt;<br>
&gt;<br>
&gt;<br>
</div></div></blockquote></div><br>
</div><br><br> </div></div></div> </div>  </div></div></blockquote></div><br>