/home/jgoppert/Projects/uvsim/uvsim/src/test/areTest.cc

Go to the documentation of this file.

#include <iostream>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/matrix_proxy.hpp>
#include "uvsim/utilities/utilities.h"
#include <boost/numeric/ublas/io.hpp>
#include <boost/numeric/bindings/lapack/lapack.h>


using namespace boost::numeric::ublas;


int main ()
{
//matrix<double>A,matrix<double>B,matrix<double>R,matrix<double>Q
matrix<double> A;
matrix<double> B;
matrix<double> R;
matrix<double> Q;

A(0,0)=1;
A(0,1)=2;
A(1,0)=2;
A(1,1)=1;

B(0,0)=1;
B(0,1)=1;
B(1,0)=1;
B(1,1)=1;

R(0,0)=1;
R(0,1)=0;
R(1,0)=0;
R(1,1)=1;

Q(0,0)=1;
Q(0,1)=0;
Q(1,0)=0;
Q(1,1)=1;

//std::cout<<dare(A,B,R,Q);


//std::cout<<A<<std::endl;

//matrix<double> B(3,3)
//matrix<double> R(3,3)

  //  matrix<double> m(3, 3);
        //    for (unsigned i = 0; i < m.size1(); ++ i)
        //                      for (unsigned j = 0; j < m.size2(); ++ j)
        //                                                  m (i, j) = 3 * i + j;               
        //identity_matrix<double> B(3);
        //std::cout<<prod(m,B);
        //      std::cout<<m<<std::endl;
        //      std::cout<<std::endl;
        
        return 0;

}

---