[ptx] Help on SVD decompoistion code.

Xianyong Fang FANGXIANYONG at CAD.ZJU.EDU.CN
Fri Apr 9 04:51:57 BST 2004


Hi Alexandre,

  I appreciate your kindly help.
  Actually I have integrated one svd decomposition into my code. In addition, you have suggested me to check the homography computation in the http://www.cs.unc.edu/~blloyd/comp290-089/fmatrix/. I find one strange phenomena, that is, my current svd code generate a different V(eigen vector matrix). Sometimes(not always!) the V of mine is negative to the V generated in Matlab(that is they have the same absolute value, but all the elements are negative to the other ones. Of course at these times the U is different from that geneated in Matlab). So the final F(the fundamental matrix) of my program is different from that generated in Matlab. I think it is very strange. 

  Of course I don't check the principle of SVD decomposition.

  Have you ever encountered such phenomena? 
	

Xianyong

========================

>BTW you can also use my old highly optimized c++ matrix class :
>http://www.le-geo.com/computer/matrix/index.html
>...
>Unfortunately, when I wrote it, I relied on the Intel MKL which was free.
>Anyway. It's a great matrice implementation (with SVD, QR, etc).
>
>Alexandre
>
>> -----Message d'origine-----
>> De : ptx-bounces at email-lists.org 
>> [mailto:ptx-bounces at email-lists.org] De la part de Xianyong Fang
>> Envoy?: jeudi 8 avril 2004 15:01
>> ?: ptx at email-lists.org
>> Objet : Re: Re: [ptx] Help on SVD decompoistion code.
>> 
>> 
>> Hi Rafal,
>> 
>>   Wonderful. It is really very useful. Thanks very much, :)
>> 
>> 	
>> 
>> Xianyong
>> 
>> ========================
>> 
>> >On Thu, 8 Apr 2004, alexandre jenny wrote:
>> >
>> >> For any numerical algorithm, there's only one thing to know : 
>> >> Numerical Recipies : http://www.nr.com/ You can get the full book 
>> >> online and the svd decomposition is inside "in c or fortran code".
>> >> (the c++ version isn't free but doesn't worth it).
>> >
>> >Unfortunately, the code is not free. You can't distribute 
>> it. I'm not 
>> >sure
>> >if you can use it for your own if you don't buy the book, either. 
>> >
>> >On the other hand www.netlib.org is full of free (as beer and as in
>> >speech) programs.  SVD is a part of LAPACK which is the 
>> standard set of 
>> >routines for linear algebra applications.  Individual routine with 
>> >dependencies in Fortran is here
>> >
>> >http://www.netlib.org/cgi-bin/netlibfiles.pl?filename=/lapack
>/double/dg
>>esvd.f
>>
>>There is a C interface to lapack (http://www.netlib.org/clapack/). 
>>There is also a C++ interface, http://math.nist.gov/lapack++/ which 
>>according to the page is superseded by the Template Numerical Toolkit 
>>(TNT) so you may also look here http://math.nist.gov/tnt/
>>
>>LAPACK is used in many scientific programs and it's well optimized and
>>documented. 
>>
>>Regards,
>>
>>Rafal
>>
>>
>>
>
>= = = = = = = = = = = = = = = = = = = =

= = = = = = = = = = = = = = = = = = = =





More information about the ptX mailing list