[ublas] LU Matrix Inversion

From: Mei, Longyu (LMei_at_[hidden])
Date: 2007-09-17 10:09:29

I am trying to use the LU Matrix Inversion
<http://www.crystalclearsoftware.com/cgi-bin/boost_wiki/wiki.pl?search=L
U+Matrix+Inversion>  to do matrix inversion. I have latest boost
(boost_1_34_1). I am use this InvertMatrix method:
http://www.crystalclearsoftware.com/cgi-bin/boost_wiki/wiki.pl?LU_Matrix
_Inversion
 
 
When I test it with the following piece of code, I cannot get the right
answer, and the result is "inverseResult: Failed". Anybody can help me
out?
 
Do I need to update "lu.hpp" header includes from boost-CVS? if yest,
anybody can email me one or tell me how to get from there?
 
 
matrix<float>inverseTest(3,3);
    for(int i = 0; i < inverseTest.size1(); ++i)
    {
        for(int j = 0; j < inverseTest.size2(); ++j)
        {
            inverseTest(i, j) = j + 1;
        }
    }
    cout << "inverseTest: " << inverseTest << endl;
    
    matrix<float>inverseResult(3,3);
    bool status = InvertMatrix(inverseTest, inverseResult);
    if(status)
    {
        cout << "inverseResult: " << inverseResult << endl;
    }
    else
    {
        cout << "inverseResult: Failed" << endl;
    } 
 
thanks,
 
James 

****************************************************************
The information transmitted in this message is confidential and may be privileged.  Any review, retransmission, dissemination, or other use of this information by persons or entities other than the intended recipient is prohibited.  If you are not the intended recipient, please notify Analogic Corporation immediately - by replying to this message or by sending an email to DeliveryErrors_at_[hidden] - and destroy all copies of this information, including any attachments, without reading or disclosing them.

Thank you.