MOAB: Mesh Oriented datABase  (version 5.2.1)
MappingFunctionTest Class Reference
+ Inheritance diagram for MappingFunctionTest:
+ Collaboration diagram for MappingFunctionTest:

Public Member Functions

void test_jacobian_2d ()
void test_jacobian_3d ()
void test_ideal_2d ()
void test_ideal_3d ()
MsqMatrix< 3, 2 > calculate_jacobian_2d (Sample s, const Vector3D *coords)
MsqMatrix< 3, 3 > calculate_jacobian_3d (Sample s, const Vector3D *coords)

Private Member Functions

 CPPUNIT_TEST_SUITE (MappingFunctionTest)
 CPPUNIT_TEST (test_jacobian_2d)
 CPPUNIT_TEST (test_jacobian_3d)
 CPPUNIT_TEST (test_ideal_2d)
 CPPUNIT_TEST (test_ideal_3d)
 CPPUNIT_TEST_SUITE_END ()

Private Attributes

LinearTriangle tri
LinearTetrahedron tet

Detailed Description

Definition at line 42 of file MappingFunctionTest.cpp.


Member Function Documentation

Definition at line 185 of file MappingFunctionTest.cpp.

References ASSERT_MATRICES_EQUAL, ASSERT_NO_ERROR, CPPUNIT_ASSERT, CPPUNIT_ASSERT_DOUBLES_EQUAL, MBMesquite::MsqMatrix< R, C >::data(), MBMesquite::det(), MBMesquite::inverse(), moab::R, MBMesquite::transpose(), MBMesquite::TRIANGLE, and MBMesquite::unit_edge_element().

{
    MsqError err;
    MsqMatrix< 3, 2 > W_prime;
    tri.MappingFunction2D::ideal( Sample( 2, 0 ), W_prime, err );
    ASSERT_NO_ERROR( err );

    // for this test that everything is in the xy-plane
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.0, W_prime( 2, 0 ), 1e-12 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.0, W_prime( 2, 1 ), 1e-12 );
    MsqMatrix< 2, 2 > W( W_prime.data() );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0, det( W ), 1e-6 );

    const Vector3D* verts = unit_edge_element( TRIANGLE );
    CPPUNIT_ASSERT( verts );

    size_t indices[3], num;
    MsqVector< 2 > coeff[3];
    tri.derivatives( Sample( 2, 0 ), NodeSet(), indices, coeff, num, err );
    ASSERT_NO_ERROR( err );

    MsqMatrix< 3, 2 > J_exp( 0.0 );  // zero
    for( size_t j = 0; j < num; ++j )
        J_exp += MsqVector< 3 >( verts[j].to_array() ) * transpose( coeff[j] );

    // for this test that everything is in the xy-plane
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.0, W_prime( 2, 0 ), 1e-12 );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.0, W_prime( 2, 1 ), 1e-12 );
    MsqMatrix< 2, 2 > W_exp( J_exp.data() );
    W_exp /= sqrt( det( W_exp ) );

    // Matrices should be a rotation of each other.
    // First, calculate tentative rotation matrix
    MsqMatrix< 2, 2 > R = inverse( W_exp ) * W;
    // next check that it is a rotation
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0, det( R ), 1e-6 );          // no scaling
    ASSERT_MATRICES_EQUAL( transpose( R ), inverse( R ), 1e-6 );  // orthogonal
}

Definition at line 224 of file MappingFunctionTest.cpp.

References ASSERT_MATRICES_EQUAL, ASSERT_NO_ERROR, MBMesquite::cbrt(), CPPUNIT_ASSERT, CPPUNIT_ASSERT_DOUBLES_EQUAL, MBMesquite::det(), MBMesquite::inverse(), moab::R, MBMesquite::TETRAHEDRON, MBMesquite::transpose(), and MBMesquite::unit_edge_element().

{
    MsqError err;
    MsqMatrix< 3, 3 > J;
    tet.MappingFunction3D::ideal( Sample( 3, 0 ), J, err );
    ASSERT_NO_ERROR( err );
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0, det( J ), 1e-6 );

    const Vector3D* verts = unit_edge_element( TETRAHEDRON );
    CPPUNIT_ASSERT( verts );

    size_t indices[4], num;
    MsqVector< 3 > coeff[4];
    tet.derivatives( Sample( 2, 0 ), NodeSet(), indices, coeff, num, err );
    ASSERT_NO_ERROR( err );

    MsqMatrix< 3, 3 > J_exp( 0.0 );  // zero
    for( size_t j = 0; j < num; ++j )
        J_exp += MsqVector< 3 >( verts[j].to_array() ) * transpose( coeff[j] );
    J_exp /= MBMesquite::cbrt( det( J_exp ) );

    // Matrices should be a rotation of each other.
    // First, calculate tentative rotation matrix
    MsqMatrix< 3, 3 > R = inverse( J_exp ) * J;
    // next check that it is a rotation
    CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0, det( R ), 1e-6 );          // no scaling
    ASSERT_MATRICES_EQUAL( transpose( R ), inverse( R ), 1e-6 );  // orthogonal
}

Definition at line 68 of file MappingFunctionTest.cpp.

References ASSERT_MATRICES_EQUAL, ASSERT_NO_ERROR, conn, CPPUNIT_ASSERT, CPPUNIT_ASSERT_EQUAL, MBMesquite::PatchData::element_by_index(), MBMesquite::PatchData::fill(), fixed, MBMesquite::MsqMeshEntity::get_vertex_index_array(), n, MBMesquite::Vector3D::to_array(), MBMesquite::transpose(), MBMesquite::TRIANGLE, and MBMesquite::PatchData::vertex_by_index().

{
    MsqError err;
    PatchData pd;
    const size_t conn[3]      = { 0, 1, 2 };
    const bool fixed[3]       = { false, false, false };
    const double coords[3][3] = { { 0.5, 0.5, 0.0 }, { 2.0, 0.2, 0.0 }, { 0.7, 0.8, 0.0 } };

    size_t indices_exp[3], indices_act[3], num_exp, num_act;
    MsqVector< 2 > coeff_exp[3], coeff_act[3];
    MsqMatrix< 3, 2 > J_exp, J_act;
    for( int r = 0; r < 3; ++r )
    {
        const double verts[] = { coords[r][0],
                                 coords[r][1],
                                 coords[r][2],
                                 coords[( r + 1 ) % 3][0],
                                 coords[( r + 1 ) % 3][1],
                                 coords[( r + 1 ) % 3][2],
                                 coords[( r + 2 ) % 3][0],
                                 coords[( r + 2 ) % 3][1],
                                 coords[( r + 2 ) % 3][2] };
        pd.fill( 3, verts, 1, TRIANGLE, conn, fixed, err );
        ASSERT_NO_ERROR( err );

        for( int d = 0; d < 3; ++d )
        {
            int n = ( d == 2 ) ? 1 : 3;
            for( int s = 0; s < n; ++s )
            {
                tri.MappingFunction2D::jacobian( pd, 0, NodeSet(), Sample( d, s ), indices_act, coeff_act, num_act,
                                                 J_act, err );
                ASSERT_NO_ERROR( err );
                CPPUNIT_ASSERT( num_act <= 3 );
                CPPUNIT_ASSERT( num_act > 0 );

                tri.derivatives( Sample( d, s ), NodeSet(), indices_exp, coeff_exp, num_exp, err );
                ASSERT_NO_ERROR( err );
                CPPUNIT_ASSERT_EQUAL( num_exp, num_act );

                J_exp = MsqMatrix< 3, 2 >( 0.0 );  // zero
                for( size_t j = 0; j < num_exp; ++j )
                {
                    size_t idx = pd.element_by_index( 0 ).get_vertex_index_array()[j];
                    size_t k   = std::find( indices_act, indices_act + num_act, idx ) - indices_act;
                    CPPUNIT_ASSERT( k < num_act );  // found
                    ASSERT_MATRICES_EQUAL( coeff_exp[j], coeff_act[k], 1e-10 );

                    J_exp += MsqVector< 3 >( pd.vertex_by_index( idx ).to_array() ) * transpose( coeff_exp[j] );
                }

                ASSERT_MATRICES_EQUAL( J_exp, J_act, 1e-6 );
            }
        }
    }
}

Definition at line 125 of file MappingFunctionTest.cpp.

References ASSERT_MATRICES_EQUAL, ASSERT_NO_ERROR, conn, CPPUNIT_ASSERT, CPPUNIT_ASSERT_EQUAL, MBMesquite::PatchData::element_by_index(), MBMesquite::PatchData::fill(), fixed, MBMesquite::MsqMeshEntity::get_vertex_index_array(), n, MBMesquite::TETRAHEDRON, MBMesquite::Vector3D::to_array(), MBMesquite::transpose(), and MBMesquite::PatchData::vertex_by_index().

{
    MsqError err;
    PatchData pd;
    const size_t conn[4]      = { 0, 1, 2, 3 };
    const bool fixed[4]       = { false, false, false, false };
    const double coords[4][3] = { { 0.5, 0.5, 0.0 }, { 2.0, 0.2, 0.0 }, { 0.7, 0.8, 0.0 }, { 0.6, 1.1, 13. } };

    size_t indices_exp[4], indices_act[4], num_exp, num_act;
    MsqVector< 3 > coeff_exp[4], coeff_act[4];
    MsqMatrix< 3, 3 > J_exp, J_act;
    for( int r = 0; r < 3; ++r )
    {
        const double verts[] = { coords[r][0],
                                 coords[r][1],
                                 coords[r][2],
                                 coords[( r + 1 ) % 3][0],
                                 coords[( r + 1 ) % 3][1],
                                 coords[( r + 1 ) % 3][2],
                                 coords[( r + 2 ) % 3][0],
                                 coords[( r + 2 ) % 3][1],
                                 coords[( r + 2 ) % 3][2],
                                 coords[3][0],
                                 coords[3][1],
                                 coords[3][2] };
        pd.fill( 4, verts, 1, TETRAHEDRON, conn, fixed, err );
        ASSERT_NO_ERROR( err );

        for( int d = 0; d < 3; ++d )
        {
            int n = ( d == 3 ) ? 1 : ( d == 1 ) ? 6 : 4;
            for( int s = 0; s < n; ++s )
            {
                tet.MappingFunction3D::jacobian( pd, 0, NodeSet(), Sample( d, s ), indices_act, coeff_act, num_act,
                                                 J_act, err );
                ASSERT_NO_ERROR( err );
                CPPUNIT_ASSERT( num_act <= 4 );
                CPPUNIT_ASSERT( num_act > 0 );

                tet.derivatives( Sample( d, s ), NodeSet(), indices_exp, coeff_exp, num_exp, err );
                ASSERT_NO_ERROR( err );
                CPPUNIT_ASSERT_EQUAL( num_exp, num_act );

                J_exp = MsqMatrix< 3, 3 >( 0.0 );  // zero
                for( size_t j = 0; j < num_exp; ++j )
                {
                    size_t idx = pd.element_by_index( 0 ).get_vertex_index_array()[j];
                    size_t k   = std::find( indices_act, indices_act + num_act, idx ) - indices_act;
                    CPPUNIT_ASSERT( k < num_act );  // found
                    ASSERT_MATRICES_EQUAL( coeff_exp[j], coeff_act[k], 1e-10 );

                    J_exp += MsqVector< 3 >( pd.vertex_by_index( idx ).to_array() ) * transpose( coeff_exp[j] );
                }

                ASSERT_MATRICES_EQUAL( J_exp, J_act, 1e-6 );
            }
        }
    }
}

Member Data Documentation

Definition at line 52 of file MappingFunctionTest.cpp.

List of all members.


The documentation for this class was generated from the following file:
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines