Re: [ublas] Lapack-solve linear equations

From: nisha kannookadan (nishak44_at_[hidden])
Date: 2008-01-12 13:01:10

Seems like half of the msg was dissappeared before, I give a second try with an easier example.
I downloaded Lapack now, but I believe I link it the wrong way.

Thanks

#include 
#include 
#include 
#include 

namespace ublas = boost::numeric::ublas;
namespace lapack = boost::numeric::bindings::lapack;

int main() {
        
        ublas::matrix A(3,3);
        A(0,0) = 1;
        A(0,1) = 1;
        A(0,2) = 1;
        A(1,0) = 2;
        A(1,1) = 3;
        A(1,2) = 1;
        A(2,0) = 1;
        A(2,1) = -1;
        A(2,2) = -1;
        
        std:: cout << A << std::endl;
        
        ublas::matrix B(3,1);
        B(0,0) = 4;
        B(0,1) = 9;
        B(0,2) = -2;
        
        std::cout << B << std::endl;
        
        lapack::gesv(A,B); 
                                 
        std::cout << B << std::endl;
}
__________________________
> Express yourself instantly with MSN Messenger! Download today it's FREE!
> http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/
> _______________________________________________
> ublas mailing list
> ublas_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/ublas

_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE!
http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/