MOAB: Mesh Oriented datABase  (version 5.4.1)
Public Member Functions | Private Member Functions | Private Attributes
IdealTargetTest Class Reference
+ Inheritance diagram for IdealTargetTest:
+ Collaboration diagram for IdealTargetTest:

Public Member Functions

void test_tri_corner ()
void test_tri_edge ()
void test_tri_center ()
void test_hex_corner ()
void test_hex_edge ()
void test_hex_face ()
void test_hex_center ()

Private Member Functions

 CPPUNIT_TEST_SUITE (IdealTargetTest)
 CPPUNIT_TEST (test_tri_corner)
 CPPUNIT_TEST (test_tri_edge)
 CPPUNIT_TEST (test_tri_center)
 CPPUNIT_TEST (test_hex_corner)
 CPPUNIT_TEST (test_hex_edge)
 CPPUNIT_TEST (test_hex_face)
 CPPUNIT_TEST (test_hex_center)
 CPPUNIT_TEST_SUITE_END ()
void get_calc_target (EntityTopology type, Sample sample, MsqMatrix< 3, 3 > &, MsqMatrix< 2, 2 > &)
void get_ideal_target (EntityTopology type, Sample sample, MsqMatrix< 3, 3 > &, MsqMatrix< 2, 2 > &)
void do_test (EntityTopology type, Sample location)

Private Attributes

Settings settings

Detailed Description

Definition at line 47 of file IdealTargetTest.cpp.

Member Function Documentation

IdealTargetTest::CPPUNIT_TEST ( test_tri_corner  ) [private]
IdealTargetTest::CPPUNIT_TEST ( test_tri_edge  ) [private]
IdealTargetTest::CPPUNIT_TEST ( test_tri_center  ) [private]
IdealTargetTest::CPPUNIT_TEST ( test_hex_corner  ) [private]
IdealTargetTest::CPPUNIT_TEST ( test_hex_edge  ) [private]
IdealTargetTest::CPPUNIT_TEST ( test_hex_face  ) [private]
IdealTargetTest::CPPUNIT_TEST ( test_hex_center  ) [private]
IdealTargetTest::CPPUNIT_TEST_SUITE ( IdealTargetTest  ) [private]
IdealTargetTest::CPPUNIT_TEST_SUITE_END ( ) [private]
void IdealTargetTest::do_test ( EntityTopology  type,
Sample  location 
) [private]

Definition at line 205 of file IdealTargetTest.cpp.

References ASSERT_MATRICES_EQUAL, CPPUNIT_ASSERT_DOUBLES_EQUAL, MBMesquite::det(), MBMesquite::inverse(), moab::R, and MBMesquite::transpose().

{
    MsqMatrix< 3, 3 > w3_calc, w3_exp;
    MsqMatrix< 2, 2 > w2_calc, w2_exp;
    get_calc_target( type, location, w3_calc, w2_calc );
    get_ideal_target( type, location, w3_exp, w2_exp );
    if( TopologyInfo::dimension( type ) == 2 )
    {
        CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0, det( w2_calc ), 1e-6 );
        CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0, det( w2_exp ), 1e-6 );
        // a rotation of the expected matrix is acceptable.
        MsqMatrix< 2, 2 > R = inverse( w2_calc ) * w2_exp;
        ASSERT_MATRICES_EQUAL( transpose( R ), inverse( R ), 1e-6 );
    }
    else
    {
        MsqMatrix< 3, 3 > R = inverse( w3_calc ) * w3_exp;
        ASSERT_MATRICES_EQUAL( transpose( R ), inverse( R ), 1e-6 );
    }
}
void IdealTargetTest::get_calc_target ( EntityTopology  type,
Sample  sample,
MsqMatrix< 3, 3 > &  w3,
MsqMatrix< 2, 2 > &  w2 
) [private]

Definition at line 131 of file IdealTargetTest.cpp.

References MBMesquite::arrptr(), MBMesquite::PatchData::attach_settings(), conn, CPPUNIT_ASSERT, MBMesquite::PatchData::fill(), MBMesquite::IdealShapeTarget::get_2D_target(), MBMesquite::IdealShapeTarget::get_3D_target(), MSQ_CHKERR, and settings.

{
    MsqPrintError err( std::cout );
    const int elem_dim = TopologyInfo::dimension( type );

    // create a patch -- actual coords and such don't really matter
    std::vector< double > coords( 24, 0.0 );
    const size_t conn[] = { 0, 1, 2, 3, 4, 5, 6, 7 };
    PatchData pd;
    pd.fill( 8, arrptr( coords ), 1, type, conn, 0, err );
    CPPUNIT_ASSERT( !MSQ_CHKERR( err ) );
    pd.attach_settings( &settings );

    IdealShapeTarget tc;
    if( elem_dim == 2 )
        tc.get_2D_target( pd, 0, location, w2, err );
    else
        tc.get_3D_target( pd, 0, location, w3, err );
    CPPUNIT_ASSERT( !MSQ_CHKERR( err ) );
}
void IdealTargetTest::get_ideal_target ( EntityTopology  type,
Sample  sample,
MsqMatrix< 3, 3 > &  w3,
MsqMatrix< 2, 2 > &  w2 
) [private]

Definition at line 155 of file IdealTargetTest.cpp.

References CPPUNIT_ASSERT, MBMesquite::MappingFunction2D::derivatives(), MBMesquite::MappingFunction3D::derivatives(), MSQ_CHKERR, settings, MBMesquite::TargetCalculator::skew(), and MBMesquite::unit_element().

{
    MsqPrintError err( std::cout );
    const unsigned elem_dim = TopologyInfo::dimension( type );

    // get the target matrix for an ideal element
    size_t indices[100];
    size_t num_vtx;
    const Vector3D* coords = unit_element( type );
    Vector3D c[3];
    if( elem_dim == 2 )
    {
        MsqVector< 2 > derivs[100];
        const MappingFunction2D* func = settings.get_mapping_function_2D( type );
        func->derivatives( location, NodeSet(), indices, derivs, num_vtx, err );
        CPPUNIT_ASSERT( !MSQ_CHKERR( err ) );

        MsqMatrix< 3, 2 > J;
        for( size_t i = 0; i < num_vtx; ++i )
            for( unsigned j = 0; j < 2; ++j )
                c[j] += derivs[i][j] * coords[indices[i]];

        for( unsigned i = 0; i < 3; ++i )
            for( unsigned j = 0; j < 2; ++j )
                J( i, j ) = c[j][i];

        w2 = TargetCalculator::skew( J );
    }
    else
    {
        MsqVector< 3 > derivs[100];
        const MappingFunction3D* func = settings.get_mapping_function_3D( type );
        func->derivatives( location, NodeSet(), indices, derivs, num_vtx, err );
        CPPUNIT_ASSERT( !MSQ_CHKERR( err ) );

        for( size_t i = 0; i < num_vtx; ++i )
            for( unsigned j = 0; j < 3; ++j )
                c[j] += derivs[i][j] * coords[indices[i]];

        for( unsigned i = 0; i < 3; ++i )
            for( unsigned j = 0; j < 3; ++j )
                w3( i, j ) = c[j][i];

        w3 = TargetCalculator::skew( w3 );
    }
}
void IdealTargetTest::test_hex_center ( )

Definition at line 126 of file IdealTargetTest.cpp.

References do_test(), and MBMesquite::HEXAHEDRON.

{
    do_test( HEXAHEDRON, Sample( 3, 0 ) );
}
void IdealTargetTest::test_hex_corner ( )

Definition at line 101 of file IdealTargetTest.cpp.

References do_test(), and MBMesquite::HEXAHEDRON.

{
    do_test( HEXAHEDRON, Sample( 0, 0 ) );
    do_test( HEXAHEDRON, Sample( 0, 1 ) );
    do_test( HEXAHEDRON, Sample( 0, 6 ) );
    do_test( HEXAHEDRON, Sample( 0, 7 ) );
}
void IdealTargetTest::test_hex_edge ( )

Definition at line 109 of file IdealTargetTest.cpp.

References do_test(), and MBMesquite::HEXAHEDRON.

{
    do_test( HEXAHEDRON, Sample( 1, 1 ) );
    do_test( HEXAHEDRON, Sample( 1, 4 ) );
    do_test( HEXAHEDRON, Sample( 1, 11 ) );
}
void IdealTargetTest::test_hex_face ( )

Definition at line 116 of file IdealTargetTest.cpp.

References do_test(), and MBMesquite::HEXAHEDRON.

{
    do_test( HEXAHEDRON, Sample( 2, 0 ) );
    do_test( HEXAHEDRON, Sample( 2, 1 ) );
    do_test( HEXAHEDRON, Sample( 2, 2 ) );
    do_test( HEXAHEDRON, Sample( 2, 3 ) );
    do_test( HEXAHEDRON, Sample( 2, 4 ) );
    do_test( HEXAHEDRON, Sample( 2, 5 ) );
}
void IdealTargetTest::test_tri_center ( )

Definition at line 96 of file IdealTargetTest.cpp.

References do_test(), and MBMesquite::TRIANGLE.

{
    do_test( TRIANGLE, Sample( 2, 0 ) );
}
void IdealTargetTest::test_tri_corner ( )

Definition at line 82 of file IdealTargetTest.cpp.

References do_test(), and MBMesquite::TRIANGLE.

{
    do_test( TRIANGLE, Sample( 0, 0 ) );
    do_test( TRIANGLE, Sample( 0, 1 ) );
    do_test( TRIANGLE, Sample( 0, 2 ) );
}
void IdealTargetTest::test_tri_edge ( )

Definition at line 89 of file IdealTargetTest.cpp.

References do_test(), and MBMesquite::TRIANGLE.

{
    do_test( TRIANGLE, Sample( 1, 0 ) );
    do_test( TRIANGLE, Sample( 1, 1 ) );
    do_test( TRIANGLE, Sample( 1, 2 ) );
}

Member Data Documentation

Settings IdealTargetTest::settings [private]

Definition at line 76 of file IdealTargetTest.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