MOAB: Mesh Oriented datABase  (version 5.2.1)
main.cpp File Reference
#include "Mesquite.hpp"
#include "MsqIGeom.hpp"
#include "FacetModifyEngine.hpp"
#include "MeshImpl.hpp"
#include "MsqError.hpp"
#include "ShapeImprovementWrapper.hpp"
#include "MsqVertex.hpp"
#include "QualityAssessor.hpp"
#include "SphericalDomain.hpp"
#include "TestUtil.hpp"
#include <memory>
#include <iostream>
#include "MsqIBase.hpp"
+ Include dependency graph for test/mesquite/igeom/main.cpp:

Go to the source code of this file.

Defines

#define CHKIGEOM   if( chk_igeom_error( ierr, __FILE__, __LINE__ ) ) return 0

Functions

bool chk_igeom_error (int ierr, const char *file, int line)
void usage ()
void run_smoother (Mesh &mesh, MeshDomain *domain, MsqError &err)
bool check_results (Mesh &mesh, MeshDomain *domain, MsqError &err)
MeshDomainget_itaps_domain ()
MeshDomainget_mesquite_domain ()
int main (int argc, char *argv[])

Variables

const double SPHERE_RADIUS = 3.0
const Vector3D SPHERE_CENTER (2.0, 2.0, 0.0)
std::string default_file_name = "/2D/vtk/quads/untangled/quads_on_sphere_529.vtk"

Define Documentation

#define CHKIGEOM   if( chk_igeom_error( ierr, __FILE__, __LINE__ ) ) return 0

Definition at line 49 of file test/mesquite/igeom/main.cpp.

Referenced by get_itaps_domain().


Function Documentation

bool check_results ( Mesh mesh,
MeshDomain domain,
MsqError err 
)

Definition at line 167 of file test/mesquite/igeom/main.cpp.

References MBMesquite::arrptr(), dof(), MBMesquite::MeshDomain::domain_DoF(), EPS, MBMesquite::Mesh::get_all_vertices(), MBMesquite::Vector3D::length(), MBMesquite::length(), MSQ_ERRZERO, SPHERE_CENTER, SPHERE_RADIUS, MBMesquite::MeshDomain::vertex_normal_at(), and MBMesquite::Mesh::vertices_get_coordinates().

Referenced by main().

{
    double EPS = 1e-2;

    std::vector< Mesh::VertexHandle > handles;
    mesh.get_all_vertices( handles, err );
    MSQ_ERRZERO( err );
    std::vector< MsqVertex > coords( handles.size() );
    mesh.vertices_get_coordinates( arrptr( handles ), arrptr( coords ), handles.size(), err );
    MSQ_ERRZERO( err );

    bool valid = true;
    size_t c   = 0;
    for( size_t i = 0; i < coords.size(); ++i )
    {
        double d = ( coords[i] - SPHERE_CENTER ).length();
        if( fabs( d - SPHERE_RADIUS ) > EPS ) ++c;
    }
    if( c )
    {
        std::cerr << c << " vertices not on domain" << std::endl;
        valid = false;
    }

    std::vector< unsigned short > dof( handles.size() ), exp_dof( handles.size(), 2 );
    dom->domain_DoF( arrptr( handles ), arrptr( dof ), handles.size(), err );
    MSQ_ERRZERO( err );
    if( dof != exp_dof )
    {
        std::cerr << "Invalid domain dimension for one or more vertices" << std::endl;
        valid = false;
    }

    c = 0;
    std::vector< Vector3D > normals( coords.begin(), coords.end() );
    dom->vertex_normal_at( arrptr( handles ), arrptr( normals ), handles.size(), err );
    MSQ_ERRZERO( err );
    for( size_t i = 0; i < handles.size(); ++i )
    {
        Vector3D exp_norm = ~( coords[i] - SPHERE_CENTER );
        Vector3D act_norm = ~normals[i];
        Vector3D diff     = exp_norm - act_norm;
        if( diff.length() > EPS ) ++c;
    }
    if( c )
    {
        std::cerr << c << " invalid vertex normals" << std::endl;
        valid = false;
    }

    return valid;
}
bool chk_igeom_error ( int  ierr,
const char *  file,
int  line 
)

Definition at line 52 of file test/mesquite/igeom/main.cpp.

References iBase_SUCCESS, and MBMesquite::process_itaps_error().

{
    if( iBase_SUCCESS == ierr ) return false;
    std::cerr << "iGeom call failed at " << file << ":" << line << std::endl;
    std::cerr << process_itaps_error( ierr ) << std::endl;
    return true;
}

Definition at line 220 of file test/mesquite/igeom/main.cpp.

References center(), CHKIGEOM, EPS, iBase_FACE, ierr, size, SPHERE_CENTER, and SPHERE_RADIUS.

Referenced by main().

{
    const double EPS = 1e-2;
    int ierr;
    iGeom_Instance igeom;
    iGeom_newGeom( "", &igeom, &ierr, 0 );
    CHKIGEOM;

#ifdef MOAB_HAVE_CGM_FACET
    FacetModifyEngine::set_modify_enabled( CUBIT_TRUE );
#endif

    iBase_EntityHandle sphere_vol;
    iGeom_createSphere( igeom, SPHERE_RADIUS, &sphere_vol, &ierr );
    CHKIGEOM;
    iGeom_moveEnt( igeom, sphere_vol, SPHERE_CENTER[0], SPHERE_CENTER[1], SPHERE_CENTER[2], &ierr );
    CHKIGEOM;

    iBase_EntityHandle sphere_surf;
    iBase_EntityHandle* ptr = &sphere_surf;
    int size = 0, alloc = 1;
    iGeom_getEntAdj( igeom, sphere_vol, iBase_FACE, &ptr, &alloc, &size, &ierr );
    CHKIGEOM;
    if( size != 1 )
    {
        std::cerr << "Failed to get spherical surface from iGeom instance" << std::endl;
        return 0;
    }

    Vector3D bmin, bmax;
    iGeom_getEntBoundBox( igeom, sphere_surf, &bmin[0], &bmin[1], &bmin[2], &bmax[0], &bmax[1], &bmax[2], &ierr );
    CHKIGEOM;
    Vector3D center = 0.5 * ( bmin + bmax );
    Vector3D rad    = 0.5 * ( bmax - bmin );
    if( ( center - SPHERE_CENTER ).length() > EPS || fabs( rad[0] - SPHERE_RADIUS ) > EPS ||
        fabs( rad[1] - SPHERE_RADIUS ) > EPS || fabs( rad[2] - SPHERE_RADIUS ) > EPS )
    {
        std::cerr << "iGeom implementation returned invalid box for sphere" << std::endl
                  << "  Expected Sphere Center: " << SPHERE_CENTER << std::endl
                  << "  Expected Sphere Radius: " << SPHERE_RADIUS << std::endl
                  << "  Actual Box Minimum    : " << bmin << std::endl
                  << "  Actual Box Maximum    : " << bmax << std::endl;
        return 0;
    }

    return new MsqIGeom( igeom, sphere_surf );
}

Definition at line 268 of file test/mesquite/igeom/main.cpp.

References SPHERE_CENTER, and SPHERE_RADIUS.

Referenced by main().

int main ( int  argc,
char *  argv[] 
)

Definition at line 75 of file test/mesquite/igeom/main.cpp.

References check_results(), default_file_name, get_itaps_domain(), get_mesquite_domain(), MBMesquite::MsqError::INTERNAL_ERROR, mesh, MSQ_CHKERR, MSQ_SETERR, MBMesquite::MeshImpl::read_vtk(), run_smoother(), usage, and MBMesquite::MeshImpl::write_vtk().

{
    // command line arguments
    const char* file_name   = 0;
    const char* output_file = 0;
    bool use_native = false, opts_done = false;
    for( int arg = 1; arg < argc; ++arg )
    {
        if( !opts_done && argv[arg][0] == '-' )
        {
            if( !strcmp( argv[arg], "-N" ) )
                use_native = true;
            else if( !strcmp( argv[arg], "-o" ) )
                output_file = argv[arg++];
            else if( !strcmp( argv[arg], "--" ) )
                opts_done = true;
            else
                usage();
        }
        else if( !file_name )
            file_name = argv[arg];
        else
            usage();
    }
    if( !file_name )
    {
        file_name = default_file_name.c_str();
        std::cout << "No file specified: using default: " << default_file_name << std::endl;
    }

    MsqError err;
    MeshImpl mesh;
    mesh.read_vtk( file_name, err );
    if( MSQ_CHKERR( err ) )
    {
        std::cerr << err << std::endl;
        std::cerr << "Failed to read input file: " << file_name << std::endl;
        return 1;
    }

    std::auto_ptr< MeshDomain > dom( use_native ? get_mesquite_domain() : get_itaps_domain() );

    if( !dom.get() )
    {
        MSQ_SETERR( err )( "Domain creation failed", MsqError::INTERNAL_ERROR );
        std::cerr << err << std::endl;
        return 1;
    }

    run_smoother( mesh, dom.get(), err );
    if( MSQ_CHKERR( err ) )
    {
        std::cerr << err << std::endl;
        std::cerr << "Smoother failed." << std::endl;
        return 1;
    }

    if( output_file )
    {
        mesh.write_vtk( output_file, err );
        if( MSQ_CHKERR( err ) )
        {
            std::cerr << err << std::endl;
            std::cerr << "Failed to write file: " << output_file << std::endl;
            return 1;
        }
    }

    bool valid = check_results( mesh, dom.get(), err );
    if( MSQ_CHKERR( err ) )
    {
        std::cerr << err << std::endl;
        std::cerr << "Error while validating results." << std::endl;
        return 1;
    }

    return !valid;
}
void run_smoother ( Mesh mesh,
MeshDomain domain,
MsqError err 
)

Definition at line 154 of file test/mesquite/igeom/main.cpp.

References MBMesquite::QualityAssessor::invalid_elements(), MBMesquite::MsqError::INVALID_MESH, MSQ_CHKERR, MSQ_SETERR, MBMesquite::Wrapper::quality_assessor(), and MBMesquite::IQInterface::run_instructions().

{
    ShapeImprovementWrapper smoother;
    MeshDomainAssoc mesh_and_domain = MeshDomainAssoc( &mesh, dom );
    smoother.run_instructions( &mesh_and_domain, err );MSQ_CHKERR( err );

    if( smoother.quality_assessor().invalid_elements() )
    {
        MSQ_SETERR( err )( "Smoothing produced invalid elements", MsqError::INVALID_MESH );
        return;
    }
}
void usage ( )

Definition at line 62 of file test/mesquite/igeom/main.cpp.

References default_file_name.

{
    std::cout << "main [-N] [filename] [-o <output_file>]" << std::endl;
    std::cout << "  -N : Use native representation instead of iGeom implementation\n" << std::endl;
    std::cout << "  If no file name is specified, will use \"" << default_file_name << '"' << std::endl;
    exit( 1 );
}

Variable Documentation

std::string default_file_name = "/2D/vtk/quads/untangled/quads_on_sphere_529.vtk"

Definition at line 60 of file test/mesquite/igeom/main.cpp.

Referenced by main(), and usage().

const double SPHERE_RADIUS = 3.0
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines