<div>Hi,<br></div><div><br></div><div>Thanks for sharing your experiments.</div><div><br></div><div>It's possible that Eigen detect error on "Optimization".</div><div><br></div><div>I want say that they are not true error ...</div> <div><br></div><div>One other possibility is to report the valgrind error the the Eigen developement team in order to isolate potential problem that they could fix ;)</div><div><br></div><div>Regards,</div><div>Pierre</div> <br><div class="gmail_quote">2011/5/20 José Luis Lisani <span dir="ltr"><<a href="mailto:joseluis.lisani@uib.es">joseluis.lisani@uib.es</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"> Hi,<br> <br> this is a follow up of the discussion "SVD for IPOL: which library, and how?"<br> <br> I just wanted to share my experience with the various options that were proposed.<br> <br> I wanted to find a substitute for a PCA computation function from Numerical Recipes<br> <br> <br> - CImg:<br> I tried what David proposed and it worked fine<br> (except for a bug in the compilation line, it should be:<br> g++ -o toto toto.cpp -Dcimg_display_type=0)<br> <br> The only drawback is that I don't need a whole image processing library,<br> just a linear algebra library<br> <br> - GSL:<br> It also worked fine, but then I tried to follow Enric's suggestion<br> (ie. 'The desired code from GSL can be easily copy-pasted into your program').<br> I spent two days trying to isolate the SVD functions that I needed but<br> I was unable to do it (too many dependencies between the source codes<br> and on global variables that were set during the GSL installation process).<br> So I gave up this option.<br> <br> - Eigen:<br> it works fine and it is very easy to use.<br> I added the 'Eigen' folder to my source code (it adds 'only' 2.16MB)<br> and now it can be compiled without problem with any C++ compiler<br> (it even works fine with Visual C++ Express).<br> The code is now slightly slower than my original code (the one using<br> Numerical Recipes) but it is OK.<br> The main drawback is that valgrind detects several errors in Eigen,<br> which raises doubts about the quality of the code.<br> <br> If you want, you can take a look at the code, available from the web page:<br> <a href="http://www.ipol.im/pub/algo/blm_color_dimensional_filtering/" target="_blank">http://www.ipol.im/pub/algo/blm_color_dimensional_filtering/</a><br> <br> <br> Best,<br> José Luis<br> <br> <br> <br> <br> <br> <br> _______________________________________________<br> discuss mailing list<br> <a href="mailto:discuss@list.ipol.im" target="_blank">discuss@list.ipol.im</a><br> <a href="http://tools.ipol.im/mailman/listinfo/discuss" target="_blank">http://tools.ipol.im/mailman/listinfo/discuss</a><br> </blockquote></div><br>