|
MOAB: Mesh Oriented datABase
(version 5.4.1)
|
|
MOAB: Mesh Oriented datABase
(version 5.4.1)
|
#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) |
| MeshDomain * | get_itaps_domain () |
| MeshDomain * | get_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 = "unittest/mesquite/2D/vtk/quads/untangled/quads_on_sphere_529.vtk" |
| #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().
| 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;
}
| MeshDomain * get_itaps_domain | ( | ) |
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().
{
return new SphericalDomain( SPHERE_CENTER, SPHERE_RADIUS );
}
| 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::unique_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 );
}
| std::string default_file_name = "unittest/mesquite/2D/vtk/quads/untangled/quads_on_sphere_529.vtk" |
Definition at line 60 of file test/mesquite/igeom/main.cpp.
| const Vector3D SPHERE_CENTER(2.0, 2.0, 0.0) |
Referenced by check_results(), get_itaps_domain(), and get_mesquite_domain().
| const double SPHERE_RADIUS = 3.0 |
Definition at line 46 of file test/mesquite/igeom/main.cpp.
Referenced by check_results(), get_itaps_domain(), and get_mesquite_domain().
1.7.6.1