qrTest.cpp File Reference

#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/io.hpp>
#include "pwiz/utility/misc/Std.hpp"
#include "pwiz/utility/misc/unit.hpp"
#include "Types.hpp"
#include "qr.hpp"

Go to the source code of this file.

Functions
template<class matrix_type >
bool	isUpperTriangular (const matrix_type &A, double eps)
void	testReflector ()
void	testQR ()
void	testRectangularQR ()
int	main (int argc, char **argv)

Variables
ostream *	os_ = 0
const double	epsilon = 1e-12

---

Function Documentation

template<class matrix_type >

bool isUpperTriangular	(	const matrix_type &	A,
		double	eps
	)

Definition at line 39 of file qrTest.cpp.

References A.

Referenced by testQR(), and testRectangularQR().

{
    typedef typename matrix_type::size_type size_type;
    bool upper = true;
    
    for (size_type i=1; i<A.size2() && upper; i++)
    {
        for (size_type j=0; j<i && upper; j++)
        {
            upper = fabs(A(i,j)) < eps;
        }
    }

    return upper;
}

void testReflector

(

)

Definition at line 56 of file qrTest.cpp.

References epsilon, F, os_, pwiz::math::Reflector(), unit_assert_equal, and x.

Referenced by main().

{
    if (os_) *os_ << "testReflector() begin" << endl;

    dmatrix F(3,3);
    dvector x(3);

    x(0) = 1;
    x(1) = 1;
    x(2) = 0;

    Reflector(x, F);

    dvector v = prod(F, x);

    unit_assert_equal(fabs(v(0)), norm_2(x), epsilon);
    unit_assert_equal(v(1), 0, epsilon);
    unit_assert_equal(v(2), 0, epsilon);
    
    x(0) = -1;
    x(1) = 1;
    x(2) = 0;

    if (os_) *os_ << "testReflector() end" << endl;
}

void testQR

(

)

Definition at line 82 of file qrTest.cpp.

References A, diff(), epsilon, isUpperTriangular(), os_, unit_assert, and unit_assert_equal.

Referenced by main().

{
    if (os_) *os_ << "testQR() begin" << endl;
    
    dmatrix A(3,3);
    dmatrix Q(3,3);
    dmatrix R(3,3);

    for (dmatrix::size_type i=0; i< A.size1(); i++)
    {
        for (dmatrix::size_type j=0; j<A.size2(); j++)
        {
            A(i,j) = ((i==j) || (j == 0));
        }
    }

    try {
        qr<dmatrix>(A, Q, R);

        unit_assert(isUpperTriangular(R, epsilon));

        dmatrix diff = (A - prod(Q, R));

        unit_assert_equal(norm_1(diff), 0, epsilon);

        identity_matrix<dmatrix::value_type> eye(Q.size1());
        diff = prod(Q, herm(Q)) - eye;
        
        unit_assert_equal(norm_1(diff), 0, epsilon);        
    }
    catch (boost::numeric::ublas::bad_argument ba)
    {
        if (os_) *os_ << "exception: " << ba.what() << endl;
    }
    
    if (os_) *os_ << "testQR() end" << endl;
}

void testRectangularQR

(

)

Definition at line 120 of file qrTest.cpp.

References A, diff(), epsilon, isUpperTriangular(), os_, unit_assert, and unit_assert_equal.

Referenced by main().

{
    if (os_) *os_ << "testRectangularQR() begin" << endl;
    
    dmatrix A(5,3);
    dmatrix Q;
    dmatrix R;

    for (dmatrix::size_type i=0; i< A.size1(); i++)
    {
        for (dmatrix::size_type j=0; j<A.size2(); j++)
        {
            A(i,j) = ((i==j) || (j == 0));
        }
    }

    try {
        qr<dmatrix>(A, Q, R);

        unit_assert(isUpperTriangular(R, epsilon));

        dmatrix diff = (A - prod(Q, R));

        unit_assert_equal(norm_1(diff), 0, epsilon);

        identity_matrix<dmatrix::value_type> eye(Q.size1());
        diff = prod(Q, herm(Q)) - eye;
        
        unit_assert_equal(norm_1(diff), 0, epsilon);        
    }
    catch (boost::numeric::ublas::bad_argument ba)
    {
        if (os_) *os_ << "exception: " << ba.what() << endl;
    }
    
    if (os_) *os_ << "testRectangularQR() end" << endl;
}

int main	(	int	argc,
		char **	argv
	)

Definition at line 158 of file qrTest.cpp.

References os_, testQR(), testRectangularQR(), and testReflector().

{
    if (argc>1 && !strcmp(argv[1],"-v")) os_ = &cout;
    if (os_) *os_ << "qrTest\n";
    testReflector();
    testQR();
    testRectangularQR();

    return 0;
}

---

Variable Documentation

ostream* os_ = 0

Definition at line 34 of file qrTest.cpp.

const double epsilon = 1e-12

Definition at line 35 of file qrTest.cpp.