Inheritance diagram for MsqMatrixTest:

Collaboration diagram for MsqMatrixTest:

Public Member Functions
void	test_initialize ()
void	test_assign ()
void	test_matrix_multiply ()
void	test_scalar_multiply ()
void	test_matrix_add ()
void	test_scalar_add ()
void	test_matrix_subtract ()
void	test_scalar_subtract ()
void	test_get_set_row ()
void	test_get_set_column ()
void	test_inverse ()
void	test_qr ()
void	test_frobenius ()
void	test_vec_length ()
void	test_determinant ()
void	test_vec_outer_product ()
Private Member Functions
	CPPUNIT_TEST_SUITE (MsqMatrixTest)
	CPPUNIT_TEST (test_initialize)
	CPPUNIT_TEST (test_assign)
	CPPUNIT_TEST (test_matrix_multiply)
	CPPUNIT_TEST (test_scalar_multiply)
	CPPUNIT_TEST (test_matrix_add)
	CPPUNIT_TEST (test_scalar_add)
	CPPUNIT_TEST (test_matrix_subtract)
	CPPUNIT_TEST (test_scalar_subtract)
	CPPUNIT_TEST (test_get_set_row)
	CPPUNIT_TEST (test_get_set_column)
	CPPUNIT_TEST (test_inverse)
	CPPUNIT_TEST (test_qr)
	CPPUNIT_TEST (test_frobenius)
	CPPUNIT_TEST (test_vec_length)
	CPPUNIT_TEST (test_determinant)
	CPPUNIT_TEST (test_vec_outer_product)
	CPPUNIT_TEST_SUITE_END ()

Detailed Description

Definition at line 41 of file MsqMatrixTest.cpp.

Member Function Documentation

MsqMatrixTest::CPPUNIT_TEST ( test_initialize ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_assign ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_matrix_multiply ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_scalar_multiply ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_matrix_add ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_scalar_add ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_matrix_subtract ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_scalar_subtract ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_get_set_row ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_get_set_column ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_inverse ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_qr ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_frobenius ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_vec_length ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_determinant ) [private]

MsqMatrixTest::CPPUNIT_TEST ( test_vec_outer_product ) [private]

MsqMatrixTest::CPPUNIT_TEST_SUITE ( MsqMatrixTest ) [private]

MsqMatrixTest::CPPUNIT_TEST_SUITE_END ( ) [private]

void MsqMatrixTest::test_assign ( )

Definition at line 130 of file MsqMatrixTest.cpp.

References CPPUNIT_ASSERT_DOUBLES_EQUAL.

{
    MsqMatrix< 3, 2 > m;
    for( int i = 0; i < 6; ++i )
    {
        m( i / 2, i % 2 ) = (double)i;
        CPPUNIT_ASSERT_DOUBLES_EQUAL( (double)i, m( i / 2, i % 2 ), DBL_EPSILON );
    }
}

void MsqMatrixTest::test_determinant ( )

Definition at line 315 of file MsqMatrixTest.cpp.

References CPPUNIT_ASSERT_DOUBLES_EQUAL, and MBMesquite::det().

{
    MsqMatrix< 5, 5 > m5( 1.0 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0, det( m5 ), 1e-8 );

    double m4vals[] = { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16 };
    MsqMatrix< 4, 4 > m4( m4vals );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.0, det( m4 ), 1e-8 );

    double m4vals2[] = { 2, 5, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1 };
    MsqMatrix< 4, 4 > m4b( m4vals2 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( -3.0, det( m4b ), 1e-8 );

    double m3vals[] = { 1, 2, -1, -1, 0, 1, 1, 1, 3 };
    MsqMatrix< 3, 3 > m3( m3vals );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 8.0, det( m3 ), 1e-8 );

    double m2vals[] = { sqrt( 2.0 ), 1, 1, sqrt( 2.0 ) };
    MsqMatrix< 2, 2 > m2( m2vals );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0, det( m2 ), 1e-8 );
}

void MsqMatrixTest::test_frobenius ( )

Definition at line 294 of file MsqMatrixTest.cpp.

References CPPUNIT_ASSERT_DOUBLES_EQUAL, and MBMesquite::sqr_Frobenius().

{
    double v3[] = { 0.1, 0.5, 0.3, 1, 2, 3, 0.5, 0.1, 10 };
    MsqMatrix< 3, 3 > M3( v3 );
    double exp = 0.0;
    for( unsigned i = 0; i < 9; ++i )
        exp += v3[i] * v3[i];
    CPPUNIT_ASSERT_DOUBLES_EQUAL( exp, sqr_Frobenius( M3 ), 1e-10 );
}

void MsqMatrixTest::test_get_set_column ( )

Definition at line 230 of file MsqMatrixTest.cpp.

References MBMesquite::MsqMatrix< R, C >::column(), CPPUNIT_ASSERT_DOUBLES_EQUAL, and MBMesquite::MsqMatrix< R, C >::set_column().

{
    double values[] = { -1, 1, -2, 2, -3, 3, -4, 4, -5 };
    MsqMatrix< 3, 3 > A( values );
    MsqMatrix< 3, 1 > c0 = A.column( 0 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[0], c0( 0, 0 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[3], c0( 1, 0 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[6], c0( 2, 0 ), DBL_EPSILON );
    MsqMatrix< 3, 1 > c1 = A.column( 1 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[1], c1( 0, 0 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[4], c1( 1, 0 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[7], c1( 2, 0 ), DBL_EPSILON );
    MsqMatrix< 3, 1 > c2 = A.column( 2 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[2], c2( 0, 0 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[5], c2( 1, 0 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[8], c2( 2, 0 ), DBL_EPSILON );

    c0( 0, 0 ) = 20;
    c0( 1, 0 ) = 30;
    c0( 2, 0 ) = 10;
    A.set_column( 0, c0 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( A( 0, 0 ), 20, DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( A( 1, 0 ), 30, DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( A( 2, 0 ), 10, DBL_EPSILON );
}

void MsqMatrixTest::test_get_set_row ( )

Definition at line 204 of file MsqMatrixTest.cpp.

References CPPUNIT_ASSERT_DOUBLES_EQUAL, MBMesquite::MsqMatrix< R, C >::row(), and MBMesquite::MsqMatrix< R, C >::set_row().

{
    double values[] = { -1, 1, -2, 2, -3, 3, -4, 4, -5 };
    MsqMatrix< 3, 3 > A( values );
    MsqMatrix< 1, 3 > r0 = A.row( 0 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[0], r0( 0, 0 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[1], r0( 0, 1 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[2], r0( 0, 2 ), DBL_EPSILON );
    MsqMatrix< 1, 3 > r1 = A.row( 1 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[3], r1( 0, 0 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[4], r1( 0, 1 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[5], r1( 0, 2 ), DBL_EPSILON );
    MsqMatrix< 1, 3 > r2 = A.row( 2 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[6], r2( 0, 0 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[7], r2( 0, 1 ), DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( values[8], r2( 0, 2 ), DBL_EPSILON );

    r0( 0, 0 ) = 20;
    r0( 0, 1 ) = 30;
    r0( 0, 2 ) = 10;
    A.set_row( 0, r0 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( A( 0, 0 ), 20, DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( A( 0, 1 ), 30, DBL_EPSILON );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( A( 0, 2 ), 10, DBL_EPSILON );
}

void MsqMatrixTest::test_initialize ( )

Definition at line 87 of file MsqMatrixTest.cpp.

References ASSERT_IDENTITY_MATRIX, and CPPUNIT_ASSERT_DOUBLES_EQUAL.

{
    MsqMatrix< 4, 4 > I44( 1.0 );
    ASSERT_IDENTITY_MATRIX( I44 );

    double values[] = { 1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0 };
    MsqMatrix< 3, 3 > I33( values );
    ASSERT_IDENTITY_MATRIX( I33 );

    double c0v[] = { 1.0, 0.0 }, c1v[] = { 0.0, 1.0 };
    MsqMatrix< 2, 1 > cols[] = { c0v, c1v };
    MsqMatrix< 2, 2 > I22( cols );
    ASSERT_IDENTITY_MATRIX( I22 );

    double r0v[] = { 1.0, 2.0 }, r1v[] = { 3.0, 4.0 };
    MsqMatrix< 1, 2 > rows[] = { r0v, r1v };
    MsqMatrix< 2, 2 > I( rows );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( r0v[0], I( 0, 0 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( r0v[1], I( 0, 1 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( r1v[0], I( 1, 0 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( r1v[1], I( 1, 1 ), 1e-6 );

    std::string str = "-1 -2 1 2 3 4";
    MsqMatrix< 2, 3 > S1( str );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( -1, S1( 0, 0 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( -2, S1( 0, 1 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 1, S1( 0, 2 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 2, S1( 1, 0 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 3, S1( 1, 1 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 4, S1( 1, 2 ), 1e-6 );

    MsqMatrix< 2, 3 > S2( str.c_str() );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( -1, S2( 0, 0 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( -2, S2( 0, 1 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 1, S2( 0, 2 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 2, S2( 1, 0 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 3, S2( 1, 1 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 4, S2( 1, 2 ), 1e-6 );

    MsqMatrix< 2, 2 > II( I33, 1, 1 );
    ASSERT_IDENTITY_MATRIX( II );
}

void MsqMatrixTest::test_inverse ( )

Definition at line 256 of file MsqMatrixTest.cpp.

References ASSERT_IDENTITY_MATRIX, and MBMesquite::inverse().

{
    double v2[] = { 4, 5, 6, 7 };
    MsqMatrix< 2, 2 > M2( v2 );
    MsqMatrix< 2, 2 > M2I = inverse( M2 );
    MsqMatrix< 2, 2 > M2P = M2 * M2I;
    ASSERT_IDENTITY_MATRIX( M2P );

    double v3[] = { 0.1, 0.5, 0.3, 1, 2, 3, 0.5, 0.1, 10 };
    MsqMatrix< 3, 3 > M3( v3 );
    MsqMatrix< 3, 3 > M3I = inverse( M3 );
    MsqMatrix< 3, 3 > M3P = M3 * M3I;
    ASSERT_IDENTITY_MATRIX( M3P );
}

void MsqMatrixTest::test_matrix_add ( )

Definition at line 164 of file MsqMatrixTest.cpp.

References ASSERT_IDENTITY_MATRIX, and moab::R.

{
    double v1[] = { 3, 4, 5, 2, 1, 0, 7, 8, 9 };
    double v2[] = { -2, -4, -5, -2, 0, 0, -7, -8, -8 };
    MsqMatrix< 3, 3 > A( v1 ), B( v2 );
    MsqMatrix< 3, 3 > R = A + B;
    ASSERT_IDENTITY_MATRIX( R );
    A += B;
    ASSERT_IDENTITY_MATRIX( A );
}

void MsqMatrixTest::test_matrix_multiply ( )

Definition at line 140 of file MsqMatrixTest.cpp.

References CPPUNIT_ASSERT_DOUBLES_EQUAL, and moab::R.

{
    double v1[] = { 1, 2, 3, -3, -2, -1 };
    double v2[] = { 4, 3, 5, 10, 2, 0.5 };

    MsqMatrix< 2, 3 > A( v1 );
    MsqMatrix< 3, 2 > B( v2 );
    MsqMatrix< 2, 2 > R = A * B;

    CPPUNIT_ASSERT_DOUBLES_EQUAL( 20.0, R( 0, 0 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 24.5, R( 0, 1 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( -24.0, R( 1, 0 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( -29.5, R( 1, 1 ), 1e-6 );
}

void MsqMatrixTest::test_matrix_subtract ( )

Definition at line 184 of file MsqMatrixTest.cpp.

References ASSERT_IDENTITY_MATRIX, and moab::R.

{
    double v1[] = { 3, 4, 5, 2, 1, 0, 7, 8, 9 };
    double v2[] = { 2, 4, 5, 2, 0, 0, 7, 8, 8 };
    MsqMatrix< 3, 3 > A( v1 ), B( v2 );
    MsqMatrix< 3, 3 > R = A - B;
    ASSERT_IDENTITY_MATRIX( R );
    A -= B;
    ASSERT_IDENTITY_MATRIX( A );
}

void MsqMatrixTest::test_qr ( )

Definition at line 271 of file MsqMatrixTest.cpp.

References ASSERT_IDENTITY_MATRIX, ASSERT_MATRICES_EQUAL, CPPUNIT_ASSERT_DOUBLES_EQUAL, MBMesquite::QR(), and MBMesquite::transpose().

{
    double v2[] = { 4, 5, 6, 7 };
    MsqMatrix< 2, 2 > M2( v2 ), Q2, R2;
    QR( M2, Q2, R2 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.0, R2( 1, 0 ), 1e-6 );
    MsqMatrix< 2, 2 > P2 = Q2 * R2;
    ASSERT_MATRICES_EQUAL( M2, P2, 1e-6 );
    MsqMatrix< 2, 2 > I2 = Q2 * transpose( Q2 );
    ASSERT_IDENTITY_MATRIX( I2 );

    double v3[] = { 0.1, 0.5, 0.3, 1, 2, 3, 0.5, 0.1, 10 };
    MsqMatrix< 3, 3 > M3( v3 ), Q3, R3;
    QR( M3, Q3, R3 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.0, R3( 1, 0 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.0, R3( 2, 0 ), 1e-6 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.0, R3( 2, 1 ), 1e-6 );
    MsqMatrix< 3, 3 > P3 = Q3 * R3;
    ASSERT_MATRICES_EQUAL( M3, P3, 1e-6 );
    MsqMatrix< 3, 3 > I3 = Q3 * transpose( Q3 );
    ASSERT_IDENTITY_MATRIX( I3 );
}

void MsqMatrixTest::test_scalar_add ( )

Definition at line 175 of file MsqMatrixTest.cpp.

References ASSERT_IDENTITY_MATRIX.

{
    double v[] = { -6, -7, -7, -6 };
    MsqMatrix< 2, 2 > A( v ), B( A + 7 );
    A += 7;
    ASSERT_IDENTITY_MATRIX( A );
    ASSERT_IDENTITY_MATRIX( B );
}

void MsqMatrixTest::test_scalar_multiply ( )

Definition at line 155 of file MsqMatrixTest.cpp.

References ASSERT_IDENTITY_MATRIX.

{
    MsqMatrix< 4, 4 > A( 33.0 );
    MsqMatrix< 4, 4 > B( 1.0 / 33 * A );
    A *= 1.0 / 33.0;
    ASSERT_IDENTITY_MATRIX( A );
    ASSERT_IDENTITY_MATRIX( B );
}

void MsqMatrixTest::test_scalar_subtract ( )

Definition at line 195 of file MsqMatrixTest.cpp.

References ASSERT_IDENTITY_MATRIX.

{
    double v[] = { -6, -7, -7, -6 };
    MsqMatrix< 2, 2 > A( v ), B( A - -7.0 );
    A -= -7.0;
    ASSERT_IDENTITY_MATRIX( A );
    ASSERT_IDENTITY_MATRIX( B );
}

void MsqMatrixTest::test_vec_length ( )

Definition at line 304 of file MsqMatrixTest.cpp.

References CPPUNIT_ASSERT_DOUBLES_EQUAL, MBMesquite::length(), and MBMesquite::sqr_length().

{
    double values[] = { 2, 4, 4 };
    MsqMatrix< 3, 1 > c( values );
    MsqMatrix< 1, 3 > r( values );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 36.0, sqr_length( c ), 1e-12 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 36.0, sqr_length( r ), 1e-12 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 6.0, length( c ), 1e-12 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 6.0, length( r ), 1e-12 );
}

void MsqMatrixTest::test_vec_outer_product ( )

Definition at line 337 of file MsqMatrixTest.cpp.

References ASSERT_MATRICES_EQUAL, MBMesquite::outer(), and MBMesquite::transpose().

{
    double v1[] = { 1, 3, -1, 12 };
    double v2[] = { -10, 4, 0.5, -9, 2 };
    MsqMatrix< 4, 1 > V1( v1 );
    MsqMatrix< 5, 1 > V2( v2 );

    MsqMatrix< 4, 5 > M1 = V1 * transpose( V2 );
    MsqMatrix< 4, 5 > M2 = outer( V1, V2 );
    ASSERT_MATRICES_EQUAL( M1, M2, 1e-8 );
}

List of all members.

The documentation for this class was generated from the following file:

MsqMatrixTest.cpp

Public Member Functions

Private Member Functions

Detailed Description

Member Function Documentation