#include <Tqdcfr.hpp>

Collaboration diagram for moab::Tqdcfr::BlockHeader:

Public Member Functions
void	print ()

	BlockHeader ()

Static Public Member Functions
static ErrorCode	read_info_header (const double data_version, const unsigned int model_offset, const FEModelHeader::ArrayInfo &info, Tqdcfr instance, BlockHeader &block_headers)

Public Attributes
unsigned int	blockID

unsigned int	blockElemType

unsigned int	memCt

unsigned int	memOffset

unsigned int	memTypeCt

unsigned int	attribOrder

unsigned int	blockCol

unsigned int	blockMixElemType

unsigned int	blockPyrType

unsigned int	blockMat

unsigned int	blockLength

unsigned int	blockDim

EntityHandle	setHandle

EntityType	blockEntityType

int	hasMidNodes [4]

Detailed Description

Definition at line 154 of file Tqdcfr.hpp.

Constructor & Destructor Documentation

◆ BlockHeader()

moab::Tqdcfr::BlockHeader::BlockHeader ( )

Definition at line 2939 of file Tqdcfr.cpp.

     : blockID( 0 ), blockElemType( 0 ), memCt( 0 ), memOffset( 0 ), memTypeCt( 0 ), attribOrder( 0 ), blockCol( 0 ),
       blockMixElemType( 0 ), blockPyrType( 0 ), blockMat( 0 ), blockLength( 0 ), blockDim( 0 ), setHandle( 0 ),
       blockEntityType( MBMAXTYPE )
 {
 }

Member Function Documentation

◆ print()

void moab::Tqdcfr::BlockHeader::print ( )

Definition at line 2946 of file Tqdcfr.cpp.

 {
     std::cout << "blockID = " << blockID << std::endl;
     std::cout << "blockElemType = " << blockElemType << std::endl;
     std::cout << "memCt = " << memCt << std::endl;
     std::cout << "memOffset = " << memOffset << std::endl;
     std::cout << "memTypeCt = " << memTypeCt << std::endl;
     std::cout << "attribOrder = " << attribOrder << std::endl;
     std::cout << "blockCol = " << blockCol << std::endl;
     std::cout << "blockMixElemType = " << blockMixElemType << std::endl;
     std::cout << "blockPyrType = " << blockPyrType << std::endl;
     std::cout << "blockMat = " << blockMat << std::endl;
     std::cout << "blockLength = " << blockLength << std::endl;
     std::cout << "blockDim = " << blockDim << std::endl;
     std::cout << "setHandle = " << setHandle << std::endl;
     std::cout << "blockEntityType = " << blockEntityType << std::endl;
 }

◆ read_info_header()

ErrorCode moab::Tqdcfr::BlockHeader::read_info_header	(	const double	data_version,
		const unsigned int	model_offset,
		const FEModelHeader::ArrayInfo &	info,
		Tqdcfr *	instance,
		Tqdcfr::BlockHeader *&	block_headers
	)

static

Definition at line 2006 of file Tqdcfr.cpp.

 {
     block_headers = new BlockHeader[info.numEntities];
     instance->FSEEK( model_offset + info.tableOffset );
     ErrorCode result;
  
     if( 0 == instance->categoryTag )
     {
         static const char val[CATEGORY_TAG_SIZE] = { 0 };
         result = instance->mdbImpl->tag_get_handle( CATEGORY_TAG_NAME, CATEGORY_TAG_SIZE, MB_TYPE_OPAQUE,
                                                     instance->categoryTag, MB_TAG_SPARSE | MB_TAG_CREAT, val );
         if( MB_SUCCESS != result && MB_ALREADY_ALLOCATED != result ) return result;
     }
  
     for( unsigned int i = 0; i < info.numEntities; i++ )
     {
         // Create an entity set for this entity
         result = instance->create_set( block_headers[i].setHandle );
         if( MB_SUCCESS != result ) return result;
         static const char material_category[CATEGORY_TAG_SIZE] = "Material Set\0";
  
         instance->FREADI( 12 );
         block_headers[i].blockID          = instance->uint_buf[0];
         block_headers[i].blockElemType    = instance->uint_buf[1];
         block_headers[i].memCt            = instance->uint_buf[2];
         block_headers[i].memOffset        = instance->uint_buf[3];
         block_headers[i].memTypeCt        = instance->uint_buf[4];
         block_headers[i].attribOrder      = instance->uint_buf[5];  // Attrib order
         block_headers[i].blockCol         = instance->uint_buf[6];
         block_headers[i].blockMixElemType = instance->uint_buf[7];  // Mixed elem type
         block_headers[i].blockPyrType     = instance->uint_buf[8];
         block_headers[i].blockMat         = instance->uint_buf[9];
         block_headers[i].blockLength      = instance->uint_buf[10];
         block_headers[i].blockDim         = instance->uint_buf[11];
  
         Tag bhTag_header;
         {
             std::vector< int > def_uint_zero( 3, 0 );
             result = instance->mdbImpl->tag_get_handle( BLOCK_HEADER, 3 * sizeof( unsigned int ), MB_TYPE_INTEGER,
                                                         bhTag_header, MB_TAG_CREAT | MB_TAG_SPARSE | MB_TAG_BYTES,
                                                         &def_uint_zero[0] );
             if( MB_SUCCESS != result ) return result;
             int block_header_data[] = { static_cast< int >( block_headers[i].blockCol ),
                                         static_cast< int >( block_headers[i].blockMat ),
                                         static_cast< int >( block_headers[i].blockDim ) };
             result =
                 instance->mdbImpl->tag_set_data( bhTag_header, &( block_headers[i].setHandle ), 1, block_header_data );
         }
  
         if( MB_SUCCESS != result ) return result;
  
         // Adjust element type for data version; older element types didn't include
         // 4 new trishell element types
         if( data_version <= 1.0 && block_headers[i].blockElemType >= 15 ) block_headers[i].blockElemType += 4;
  
         if( block_headers[i].blockElemType >= (unsigned)cub_elem_num_verts_len )
         {
             // Block element type unassigned, will have to infer from verts/element; make sure it's
             // the expected value of 52
             if( ( 14 == major && 2 < minor ) || 15 <= major )
             {
                 if( 55 != block_headers[i].blockElemType )
                     MB_SET_ERR( MB_FAILURE, "Invalid block element type: " << block_headers[i].blockElemType );
             }
             else
             {
                 if( 52 != block_headers[i].blockElemType )
                     MB_SET_ERR( MB_FAILURE, "Invalid block element type: " << block_headers[i].blockElemType );
             }
         }
  
         // Set the material set tag and id tag both to id
         result = instance->mdbImpl->tag_set_data( instance->blockTag, &( block_headers[i].setHandle ), 1,
                                                   &( block_headers[i].blockID ) );
         if( MB_SUCCESS != result ) return result;
         result = instance->mdbImpl->tag_set_data( instance->globalIdTag, &( block_headers[i].setHandle ), 1,
                                                   &( block_headers[i].blockID ) );
         if( MB_SUCCESS != result ) return result;
         result = instance->mdbImpl->tag_set_data( instance->categoryTag, &( block_headers[i].setHandle ), 1,
                                                   material_category );
         if( MB_SUCCESS != result ) return result;
  
         // If this block is empty, continue
         if( !block_headers[i].memCt ) continue;
  
         // Check the number of vertices in the element type, and set the has mid nodes tag
         // accordingly; if element type wasn't set, they're unlikely to have mid nodes
         // 52 is for CUBIT versions below 14.1, 55 for CUBIT version 14.9 and above
         if( 52 != block_headers[i].blockElemType && 55 != block_headers[i].blockElemType )
         {
             int num_verts                    = cub_elem_num_verts[block_headers[i].blockElemType];
             block_headers[i].blockEntityType = block_type_to_mb_type[block_headers[i].blockElemType];
             if( ( block_headers[i].blockEntityType < MBMAXTYPE ) &&
                 ( num_verts != CN::VerticesPerEntity( block_headers[i].blockEntityType ) ) )
             {
                 // Not a linear element; try to find hasMidNodes values
                 for( int j = 0; j < 4; j++ )
                     block_headers[i].hasMidNodes[j] = 0;
                 if( 0 == instance->hasMidNodesTag )
                 {
                     result = instance->mdbImpl->tag_get_handle( HAS_MID_NODES_TAG_NAME, 4, MB_TYPE_INTEGER,
                                                                 instance->hasMidNodesTag, MB_TAG_SPARSE | MB_TAG_CREAT,
                                                                 block_headers[i].hasMidNodes );
                     if( MB_SUCCESS != result ) return result;
                 }
  
                 CN::HasMidNodes( block_headers[i].blockEntityType, num_verts, block_headers[i].hasMidNodes );
  
                 // Now set the tag on this set
                 result = instance->mdbImpl->tag_set_data( instance->hasMidNodesTag, &block_headers[i].setHandle, 1,
                                                           block_headers[i].hasMidNodes );
                 if( MB_SUCCESS != result ) return result;
             }
         }
     }
  
     return MB_SUCCESS;
 }