 |
Mesh Oriented datABase
(version 5.5.1)
An array-based unstructured mesh library
|
|
Mesh Oriented datABase
(version 5.5.1)
An array-based unstructured mesh library
|
#include <ReadDamsel.hpp>
Inheritance diagram for moab::ReadDamsel: Collaboration diagram for moab::ReadDamsel:Classes | |
| class | subrange |
Public Member Functions | |
| ErrorCode | load_file (const char *file_name, const EntityHandle *file_set, const FileOptions &opts, const SubsetList *subset_list=0, const Tag *file_id_tag=0) |
| Load an NC file. More... | |
| ReadDamsel (Interface *impl=NULL) | |
| Constructor. More... | |
| virtual | ~ReadDamsel () |
| Destructor. More... | |
| virtual ErrorCode | read_tag_values (const char *file_name, const char *tag_name, const FileOptions &opts, std::vector< int > &tag_values_out, const SubsetList *subset_list=0) |
| Read tag values from a file. More... | |
| Public Member Functions inherited from moab::ReaderIface | |
| virtual | ~ReaderIface () |
Static Public Member Functions | |
| static ReaderIface * | factory (Interface *) |
Private Member Functions | |
| ErrorCode | get_contents (damsel_model m, damsel_container c, Range &ents) |
| Get contents of the container (containing file-side handles) and translate to moab-side handles. More... | |
| ErrorCode | get_contents (damsel_model m, damsel_container c, EntityHandle *ents) |
| Get contents of the container (containing file-side handles) and translate to moab-side handles ents argument should point to already-allocated space. More... | |
| ErrorCode | init () |
| ErrorCode | parse_options (const FileOptions &opts, bool ¶llel) |
| ErrorCode | process_tags (std::vector< damsel_tag_buf_type > &tag_infos) |
| ErrorCode | process_ent_info (const damsel_entity_buf_type &einfo) |
| ErrorCode | process_entity_tags (int count, damsel_container tag_container, damsel_container app_cont, Range &these_ents) |
| ErrorCode | process_coll_infos (std::vector< damsel_collection_buf_type > &coll_infos) |
| ErrorCode | container_to_handle_pairs (damsel_container &cont, std::vector< damsel_handle > &handle_pairs) |
| Convert handles in a container into handle pairs, one pair per contiguous sequence of handles in the container. More... | |
| ErrorCode | insert_into_map (std::vector< damsel_handle > &handle_pairs, EntityHandle start_handle) |
| Store MOAB handles starting from start_handle, corresponding to store handles in handle_pairs, into the entity map. More... | |
Private Attributes | |
| Interface * | mbImpl |
| Interface instance. More... | |
| ReadUtilIface * | readMeshIface |
| UtilIface. More... | |
| std::string | fileName |
| File name. More... | |
| bool | nativeParallel |
| Whether this reader natively supports parallel semantics. More... | |
| ParallelComm * | myPcomm |
| Parallel info. More... | |
| Tag | mGlobalIdTag |
| Used to track file handles. More... | |
| RangeMap< damsel_handle, EntityHandle, 0 > | dmHandleRMap |
| map from damsel to moab handles More... | |
| DamselUtil | dU |
| Keep various damsel data. More... | |
Definition at line 36 of file ReadDamsel.hpp.
| moab::ReadDamsel::ReadDamsel | ( | Interface * | impl = NULL | ) |
Constructor.
Definition at line 30 of file ReadDamsel.cpp.
31 : mbImpl( impl ), readMeshIface( NULL ), nativeParallel( false ), myPcomm( NULL ), mGlobalIdTag( 0 ), dU()
32 {
33 init();
34 }
References init().
Referenced by factory().
|
virtual |
Destructor.
Definition at line 36 of file ReadDamsel.cpp.
37 { 38 if( readMeshIface ) mbImpl->release_interface( readMeshIface ); 39 40 DMSLlib_finalize( dU.dmslLib ); 41 }
References moab::DamselUtil::dmslLib, dU, mbImpl, readMeshIface, and moab::Interface::release_interface().
|
private |
Convert handles in a container into handle pairs, one pair per contiguous sequence of handles in the container.
|
static |
Definition at line 25 of file ReadDamsel.cpp.
26 {
27 return new ReadDamsel( iface );
28 }
References iface, and ReadDamsel().
Referenced by moab::ReaderWriterSet::ReaderWriterSet().
|
private |
Get contents of the container (containing file-side handles) and translate to moab-side handles ents argument should point to already-allocated space.
Definition at line 489 of file ReadDamsel.cpp.
490 {
491 EntityHandle eh;
492 int ind = 0;
493
494 if( DMSLcontainer_get_type( c ) == DAMSEL_HANDLE_CONTAINER_TYPE_SEQUENCE )
495 {
496 damsel_handle start;
497 size_t count, stride;
498 damsel_err_t err = DMSLcontainer_sequence_get_contents( m, c, &start, &count, &stride );
499 CHK_DMSL_ERR( err, "Problem calling DMSLcontainer_sequence_get_contents" );
500 if( stride == 1 )
501 {
502 while( count )
503 {
504 // Get start in rangemap
505 RangeMap< damsel_handle, EntityHandle, 0 >::iterator beg = dmHandleRMap.lower_bound( start );
506 if( beg == dmHandleRMap.end() ) return MB_SUCCESS;
507 unsigned int diff = std::max( ( *beg ).begin - start, (damsel_handle)0 );
508 unsigned int num = std::min( count - diff, (size_t)( *beg ).count );
509 for( EntityHandle hdl = ( *beg ).begin + diff; hdl <= (int)( *beg ).begin + diff + num - 1; hdl++ )
510 ents[ind++] = hdl;
511 count -= ( diff + num );
512 ++beg;
513 }
514 }
515 else
516 {
517 for( int i = count - 1; i >= 0; i-- )
518 {
519 if( dmHandleRMap.find( start + i, eh ) ) ents[i] = eh;
520 }
521 }
522 }
523 else if( DMSLcontainer_get_type( c ) == DAMSEL_HANDLE_CONTAINER_TYPE_VECTOR )
524 {
525 damsel_handle_ptr handle_ptr;
526 size_t count;
527 damsel_err_t err = DMSLcontainer_vector_get_contents( m, c, &handle_ptr, &count );
528 CHK_DMSL_ERR( err, "Failed to get vector contents" );
529 for( int i = count - 1; i >= 0; i-- )
530 {
531 if( dmHandleRMap.find( handle_ptr[i], eh ) ) ents[i] = eh;
532 }
533 }
534 else if( DMSLcontainer_get_type( c ) == DAMSEL_HANDLE_CONTAINER_TYPE_TREE )
535 {
536 damsel_handle_ptr node_ptr = NULL;
537 damsel_container cont = NULL;
538 while( DMSLcontainer_tree_get_contents( m, c, &node_ptr, &cont ) == DMSL_OK && cont )
539 {
540 ErrorCode rval = get_contents( m, cont, ents );
541 if( MB_SUCCESS != rval ) return rval;
542 unsigned int num = DMSLcontainer_count( cont );
543 ents += num;
544 }
545 }
546
547 return MB_SUCCESS;
548 }
References CHK_DMSL_ERR, dmHandleRMap, moab::RangeMap< KeyType, ValType, NullVal >::end(), ErrorCode, moab::RangeMap< KeyType, ValType, NullVal >::find(), get_contents(), moab::RangeMap< KeyType, ValType, NullVal >::lower_bound(), and MB_SUCCESS.
|
private |
Get contents of the container (containing file-side handles) and translate to moab-side handles.
Definition at line 433 of file ReadDamsel.cpp.
434 {
435 EntityHandle eh;
436 if( DMSLcontainer_get_type( c ) == DAMSEL_HANDLE_CONTAINER_TYPE_SEQUENCE )
437 {
438 damsel_handle start;
439 size_t count, stride;
440 damsel_err_t err = DMSLcontainer_sequence_get_contents( m, c, &start, &count, &stride );
441 CHK_DMSL_ERR( err, "DMSLcontainer_sequence_get_contents" );
442 if( stride == 1 )
443 {
444 while( count )
445 {
446 // Get start in rangemap
447 RangeMap< damsel_handle, EntityHandle, 0 >::iterator beg = dmHandleRMap.lower_bound( start );
448 if( beg == dmHandleRMap.end() ) return MB_SUCCESS;
449 unsigned long diff = std::max( ( *beg ).begin - start, (damsel_handle)0 );
450 unsigned long num = std::min( count - diff, (size_t)( *beg ).count );
451 ents.insert( ( *beg ).begin + diff, ( *beg ).begin + diff + num - 1 );
452 count -= ( diff + num );
453 ++beg;
454 }
455 }
456 else
457 {
458 for( int i = count - 1; i >= 0; i-- )
459 {
460 if( dmHandleRMap.find( start + i, eh ) ) ents.insert( eh );
461 }
462 }
463 }
464 else if( DMSLcontainer_get_type( c ) == DAMSEL_HANDLE_CONTAINER_TYPE_VECTOR )
465 {
466 damsel_handle* handle_ptr;
467 size_t count;
468 damsel_err_t err = DMSLcontainer_vector_get_contents( m, c, &handle_ptr, &count );
469 CHK_DMSL_ERR( err, "Trouble getting vector contents" );
470 for( int i = count - 1; i >= 0; i-- )
471 {
472 if( dmHandleRMap.find( handle_ptr[i], eh ) ) ents.insert( eh );
473 }
474 }
475 else if( DMSLcontainer_get_type( c ) == DAMSEL_HANDLE_CONTAINER_TYPE_TREE )
476 {
477 damsel_handle_ptr node_ptr = NULL;
478 damsel_container cont = NULL;
479 while( DMSLcontainer_tree_get_contents( m, c, &node_ptr, &cont ) == DMSL_OK && cont )
480 {
481 ErrorCode rval = get_contents( m, c, ents );
482 if( MB_SUCCESS != rval ) return rval;
483 }
484 }
485
486 return MB_SUCCESS;
487 }
References CHK_DMSL_ERR, dmHandleRMap, moab::RangeMap< KeyType, ValType, NullVal >::end(), ErrorCode, moab::RangeMap< KeyType, ValType, NullVal >::find(), moab::Range::insert(), moab::RangeMap< KeyType, ValType, NullVal >::lower_bound(), and MB_SUCCESS.
Referenced by get_contents(), and process_ent_info().
|
private |
Definition at line 43 of file ReadDamsel.cpp.
44 {
45 mbImpl->query_interface( readMeshIface );
46 assert( readMeshIface );
47
48 return MB_SUCCESS;
49 }
References MB_SUCCESS, mbImpl, moab::Interface::query_interface(), and readMeshIface.
Referenced by ReadDamsel().
|
private |
Store MOAB handles starting from start_handle, corresponding to store handles in handle_pairs, into the entity map.
|
virtual |
Load an NC file.
Implements moab::ReaderIface.
Definition at line 70 of file ReadDamsel.cpp.
75 {
76 ErrorCode rval;
77
78 rval = parse_options( opts, nativeParallel );
79 if( MB_SUCCESS != rval ) return rval;
80
81 // Initialize damsel
82 dU.dmslLib = DMSLlib_init();
83
84 // Create a damsel model
85 dU.dmslModel =
86 DMSLmodel_create( sizeof( EntityHandle ) == 8 ? DAMSEL_HANDLE_TYPE_HANDLE64 : DAMSEL_HANDLE_TYPE_HANDLE32 );
87
88 // Model attach - need model id from make model, filename
89 #ifdef MOAB_HAVE_MPI
90 MPI_Comm comm = MPI_COMM_WORLD;
91 if( nativeParallel )
92 {
93 comm = myPcomm->proc_config().proc_comm();
94 }
95 #endif
96
97 damsel_err_t err;
98 err = DMSLmodel_attach( dU.dmslModel, filename, comm, NULL );
99 CHK_DMSL_ERR( err, "DMSLmodel_attach failed" );
100 err = DMSLmodel_populate( dU.dmslModel );
101 CHK_DMSL_ERR( err, "DMSLmodel_populate failed" );
102
103 // STEP 0: GET COLLECTION, TAG, ENTITY INFOS FOR GLOBAL MODEL
104 int num_containers = 0, num_tag_infos = 0, num_ent_infos = 0;
105 DMSLmodel_get_tuple_count( dU.dmslModel, &num_containers, &num_tag_infos );
106 num_ent_infos = DMSLmodel_get_entity_count( dU.dmslModel );
107 int num_coll_infos = DMSLmodel_get_collection_count( dU.dmslModel );
108 CHK_DMSL_ERR( err, "DMSLmodel_get_collection_count failed" );
109 if( -1 == num_containers || -1 == num_tag_infos || -1 == num_ent_infos )
110 MB_SET_ERR( MB_FAILURE, "Bad count for containers/tags/ents" );
111
112 std::vector< damsel_entity_buf_type > ent_infos( num_ent_infos );
113 std::vector< damsel_collection_buf_type > coll_infos( num_coll_infos );
114 std::vector< damsel_tag_buf_type > tag_infos( num_tag_infos );
115 std::vector< damsel_container_buf_type > cont_infos( num_containers );
116 err = DMSLmodel_get_entity_infos( dU.dmslModel, &ent_infos[0] );
117 CHK_DMSL_ERR( err, "Failure getting entity infos" );
118 err = DMSLmodel_get_collection_infos( dU.dmslModel, &coll_infos[0] );
119 CHK_DMSL_ERR( err, "Failure getting collection infos" );
120 err = DMSLmodel_get_tag_infos( dU.dmslModel, &tag_infos[0] );
121 CHK_DMSL_ERR( err, "Failure getting tag infos" );
122 err = DMSLmodel_get_container_infos( dU.dmslModel, &cont_infos[0] );
123 CHK_DMSL_ERR( err, "Failure getting container infos" );
124
125 // Create MOAB-side tags for all damsel tags except pre-defined ones
126 rval = process_tags( tag_infos );MB_CHK_SET_ERR( rval, "Error processing tags" );
127
128 /*
129 if (nativeParallel) {
130 // STEP 1: GET COLLECTION(S) REPRESENTING PARTITION:
131 // Input: tag name, optionally value;
132 // Output: container with file-side handles of collections satisfying those criteria
133 // - Get all handles/values for tag
134 // - Select handles matching criteria for tag value (will be collection handles)
135 std::string partn_tag_name("PARALLEL_PARTITION");
136 damsel_handle partn_tag = DMSLselect_tag_by_name(dU.dmslModel, partn_tag_name.c_str());
137 // Get all the parts with that tag regardless of value
138 damsel_container part_handles = DMSLselect_handles_with_values(dU.dmslModel, partn_tag);
139
140 // STEP 2: GET HANDLES FOR TAGS WE NEED FOR THIS READER:
141 // - "SET_CHARACTERISTIC"
142 damsel_handle setchar_tag = DMSLselect_tag_by_name(dU.dmslModel, "SET_CHARACTERISTIC");
143 // - "PARENT_LIST"
144 //damsel_handle plist_tag = DMSLselect_tag_by_name(dU.dmslModel, "PARENT_LIST");
145 // - "CHILD_LIST"
146 //damsel_handle clist_tag = DMSLselect_tag_by_name(dU.dmslModel, "CHILD_LIST");
147
148 // STEP 3: GET VALUES FOR "SET_CHARACTERISTIC" TAG ON PARTITION COLLECTIONS,
149 // GET VECTOR- OR SET-TYPE FLAGS FOR PARTITION COLLECTIONS
150 // (gives tracking flag for moab)
151 int num_handles = DMSLcontainer_count(part_handles);
152 std::vector<unsigned> char_tagvals(num_handles);
153 // Map the set chars
154 err = DMSLmodel_map_tag(&char_tagvals[0], part_handles, &setchar_tag);CHK_DMSL_ERR(err,
155 "Problem calling DMSLmodel_map_tag");
156
157 // Execute the transfer
158 err = DMSLmodel_transfer_sync(dU.dmslModel, DAMSEL_TRANSFER_TYPE_READ);CHK_DMSL_ERR(err,
159 "Problem calling DMSLmodel_transfer_sync");
160
161 // STEP 4: READ/PROCESS PARTITION COLLECTION(S)
162 // Decide the handles I am responsible using round-robin for now
163 // - GET TYPE, CONTENTS OF COLLECTION CONTENTS CONTAINER
164 // - Allocate moab-side container (using count from container)
165 // - MAP storage TO CONTAINER
166 // - EXECUTE
167 // ==> have list of all handles (entities + collections) represented on this proc
168
169 int tmp_num = num_handles / proc_size, extra = num_handles % proc_size;
170 if (extra) tmp_num++;
171 int my_num_handles = tmp_num;
172 if (proc_rank >= extra) my_num_handles--;
173 int first_ind = std::min(proc_rank, extra) * tmp_num +
174 std::max(proc_rank - extra, 0) * (tmp_num - 1);
175
176 // - Create moab entity sets for partition collection(s)
177 EntityHandle start_handle;
178 rval = readMeshIface->create_entity_sets(my_num_handles, &char_tagvals[first_ind], 0,
179 start_handle);MB_CHK_SET_ERR(rval, "Problem creating entity sets");
180 }
181 else {
182 */
183 // Initialize just by entity; each call to process_ent_info will:
184 // a. Create moab-side representation to read into
185 // b. Map those handles to damsel handles
186 // c. Map coords / connectivity storage to damsel equivalent
187 // d. For each tag, map moab storage to damsel storage
188 std::vector< damsel_entity_buf_type >::iterator eiit;
189
190 // Process verts info first
191 for( eiit = ent_infos.begin(); eiit != ent_infos.end(); ++eiit )
192 {
193 if( ( *eiit ).entity_type == DAMSEL_ENTITY_TYPE_VERTEX )
194 {
195 rval = process_ent_info( *eiit );MB_CHK_ERR( rval );
196 }
197 }
198
199 for( eiit = ent_infos.begin(); eiit != ent_infos.end(); ++eiit )
200 {
201 if( ( *eiit ).entity_type != DAMSEL_ENTITY_TYPE_VERTEX )
202 {
203 rval = process_ent_info( *eiit );MB_CHK_ERR( rval );
204 }
205 }
206
207 /*
208 }
209
210 // Process collections
211 rval = process_coll_infos(coll_infos);MB_CHK_ERR(rval);
212
213 // STEP 5: Process into list of local info structs, each represents file-side struct and
214 // portion of that struct
215 // ASSUMPTION: Each local info struct represents single entity type & # vertices or collection
216
217 // STEP 6: For each local info struct:
218
219 // STEP 6b: READ CONTAINER INTO LOCAL BUFFER
220 // STEP 6c: Create app representation of entities/vertices/collection, and damsel container
221 for them,
222 // and MAP APP HANDLE CONTAINER TO DAMSEL CONTAINER
223 // STEP 6d: Process vertices/entities/collection
224 // 6d1: If vertices, continue
225 // 6d2: If entities:
226 // - MAP LOCAL CONNECTIVITY REP'N TO DAMSEL (might be tag, don't know yet)
227 // 6d3: If collection:
228 // - (everything in STEP 4 for this collection except for EXECUTE)
229 // - might need to filter out things not represented on this rank
230 // 6d4: If sparse tag:
231 // - Create app-side representation of sparse tag
232 // - READ CONTAINER OF MODEL HANDLES INTO LOCAL BUFFER
233 // - Allocate app-side storage for tag values
234 // - MAP APP STORAGE TO MODEL TAG + (implicit?) CONTAINER
235
236 // STEP 6e: Process dense tags for the local info struct; for each dense tag:
237 // - Get app tag handle for model tag handle
238 // - Get app storage for app tag handle + container
239 // - MAP APP STORAGE TO MODEL TAG + CONTAINER
240
241 // STEP 7: EXECUTE
242 // - Assign all mapped data
243 // - Translate all damsel handles to app handles
244 // uninit damsel
245
246 */
247
248 return MB_SUCCESS;
249 }
References CHK_DMSL_ERR, moab::DamselUtil::dmslLib, moab::DamselUtil::dmslModel, dU, ErrorCode, MB_CHK_ERR, MB_CHK_SET_ERR, MB_SET_ERR, MB_SUCCESS, myPcomm, nativeParallel, parse_options(), moab::ProcConfig::proc_comm(), moab::ParallelComm::proc_config(), process_ent_info(), and process_tags().
|
private |
Definition at line 51 of file ReadDamsel.cpp.
52 {
53 // Handle parallel options
54 bool use_mpio = ( MB_SUCCESS == opts.get_null_option( "USE_MPIO" ) );
55 ErrorCode rval = opts.match_option( "PARALLEL", "READ_PART" );
56 parallel = ( rval != MB_ENTITY_NOT_FOUND );
57 nativeParallel = ( rval == MB_SUCCESS );
58 if( use_mpio && !parallel )
59 {
60 MB_SET_ERR( MB_NOT_IMPLEMENTED, "'USE_MPIO' option specified w/out 'PARALLEL' option" );
61 }
62
63 return MB_SUCCESS;
64 }
References ErrorCode, moab::FileOptions::get_null_option(), moab::FileOptions::match_option(), MB_ENTITY_NOT_FOUND, MB_NOT_IMPLEMENTED, MB_SET_ERR, MB_SUCCESS, and nativeParallel.
Referenced by load_file().
|
private |
|
private |
Definition at line 317 of file ReadDamsel.cpp.
318 {
319 // Create this chunk of entities
320 EntityHandle *connect, start_handle;
321 ErrorCode rval;
322 damsel_err_t err;
323 damsel_container app_cont;
324 Range these_ents;
325
326 // Check that there's only one contiguous run of file-side handles, fail if there isn't
327 #ifndef NDEBUG
328 Range fside_handles;
329 rval = DamselUtil::container_to_range( dU.dmslModel, const_cast< damsel_container& >( einfo.entity_container ),
330 fside_handles );
331 if( MB_SUCCESS != rval || fside_handles.size() != einfo.count || fside_handles.psize() != 1 ) return MB_FAILURE;
332 #endif
333
334 if( einfo.entity_type != DAMSEL_ENTITY_TYPE_VERTEX )
335 {
336 // Create the moab entities
337 rval = readMeshIface->get_element_connect( einfo.count, einfo.vertices_per_entity,
338 DamselUtil::dtom_entity_type[einfo.entity_type], 0, start_handle,
339 connect );MB_CHK_ERR( rval );
340 these_ents.insert( start_handle, start_handle + einfo.count - 1 );
341
342 // Create an app-side sequence and map to file-side container
343 app_cont = DMSLcontainer_create_sequence( dU.dmslModel, einfo.count, start_handle, 1 );
344 err = DMSLmodel_map_handles( app_cont, einfo.entity_container );
345 CHK_DMSL_ERR( err, "Error returned mapping entity handles" );
346
347 // Map connectivity
348 assert( DMSLcontainer_count( einfo.vertex_container ) == (int)( einfo.vertices_per_entity * einfo.count ) );
349 rval = get_contents( dU.dmslModel, einfo.vertex_container, connect );MB_CHK_SET_ERR( rval, "Error returned mapping connectivity" );
350 }
351 else
352 {
353 // Get the number of coordinate arrays
354 int num_ctags = 0;
355 damsel_handle xcoord_dtag = DMSLselect_tag_by_name( dU.dmslModel, "mbdmsl_XCOORDS" );
356 if( xcoord_dtag ) num_ctags++;
357 damsel_handle ycoord_dtag = DMSLselect_tag_by_name( dU.dmslModel, "mbdmsl_YCOORDS" );
358 if( ycoord_dtag ) num_ctags++;
359 damsel_handle zcoord_dtag = DMSLselect_tag_by_name( dU.dmslModel, "mbdmsl_ZCOORDS" );
360 if( zcoord_dtag ) num_ctags++;
361
362 // Should have one vertex per entity
363 assert( einfo.vertices_per_entity == 1 );
364 std::vector< double* > coord_arrays;
365 rval = readMeshIface->get_node_coords( num_ctags, einfo.count, 0, start_handle, coord_arrays );MB_CHK_ERR( rval );
366
367 these_ents.insert( start_handle, start_handle + einfo.count - 1 );
368
369 // Create an app-side sequence and map to file-side container
370 app_cont = DMSLcontainer_create_sequence( dU.dmslModel, start_handle, einfo.count, 1 );
371 err = DMSLmodel_map_handles( app_cont, einfo.entity_container );
372 CHK_DMSL_ERR( err, "Trouble mapping entity handles" );
373
374 // Map the coords storage
375 if( xcoord_dtag != 0 )
376 {
377 err = DMSLmodel_map_tag( coord_arrays[0], app_cont, (damsel_handle_ptr)&dU.xcoordsTag.mTagh );
378 CHK_DMSL_ERR( err, "Trouble mapping x coordinate tag" );
379 }
380 if( ycoord_dtag != 0 )
381 {
382 err = DMSLmodel_map_tag( coord_arrays[1], app_cont, (damsel_handle_ptr)&dU.ycoordsTag.mTagh );
383 CHK_DMSL_ERR( err, "Trouble mapping y coordinate tag" );
384 }
385 if( zcoord_dtag != 0 )
386 {
387 err = DMSLmodel_map_tag( coord_arrays[2], app_cont, (damsel_handle_ptr)&dU.zcoordsTag.mTagh );
388 CHK_DMSL_ERR( err, "Trouble mapping z coordinate tag" );
389 }
390 }
391
392 // Save mapping from moab entity to einfo
393 dmHandleRMap.insert( DMSLcontainer_handle_at_position( einfo.entity_container, 0 ), start_handle, einfo.count );
394
395 rval = process_entity_tags( einfo.tag_count, einfo.tag_handle_container, app_cont, these_ents );
396
397 return rval;
398 }
References CHK_DMSL_ERR, moab::DamselUtil::container_to_range(), dmHandleRMap, moab::DamselUtil::dmslModel, moab::DamselUtil::dtom_entity_type, dU, ErrorCode, get_contents(), moab::ReadUtilIface::get_element_connect(), moab::ReadUtilIface::get_node_coords(), moab::Range::insert(), moab::RangeMap< KeyType, ValType, NullVal >::insert(), MB_CHK_ERR, MB_CHK_SET_ERR, MB_SUCCESS, moab::DamselUtil::tinfo::mTagh, process_entity_tags(), moab::Range::psize(), readMeshIface, moab::Range::size(), moab::DamselUtil::xcoordsTag, moab::DamselUtil::ycoordsTag, and moab::DamselUtil::zcoordsTag.
Referenced by load_file().
|
private |
Definition at line 400 of file ReadDamsel.cpp.
404 {
405 // Process tags on these entities
406 ErrorCode rval = MB_SUCCESS;
407 for( int i = 0; i < count; i++ )
408 {
409 damsel_handle dtagh = DMSLcontainer_handle_at_position( tag_container, i );
410
411 // Don't do conventional tags here
412 std::vector< DamselUtil::tinfo >::iterator vit =
413 std::find_if( dU.tagMap.begin(), dU.tagMap.end(), DamselUtil::DtagP< DamselUtil::tinfo >( dtagh ) );
414
415 if( ( *vit ).tagType == MB_TAG_ANY )
416 continue;
417 else if( vit == dU.tagMap.end() )
418 MB_SET_ERR( MB_FAILURE, "Failed to find tag" );
419
420 Tag tagh = ( *vit ).mTagh;
421 assert( tagh );
422 void* tag_data;
423 int ecount = these_ents.size();
424 rval = mbImpl->tag_iterate( tagh, these_ents.begin(), these_ents.end(), ecount, tag_data );MB_CHK_SET_ERR( rval, "Problem getting tag iterator" );
425 assert( ecount == (int)these_ents.size() );
426 damsel_err_t err = DMSLmodel_map_tag( tag_data, app_cont, (damsel_handle_ptr)&tagh );
427 CHK_DMSL_ERR( err, "Problem calling DMSLmodel_map_tag" );
428 }
429
430 return rval;
431 }
References moab::Range::begin(), CHK_DMSL_ERR, dU, moab::Range::end(), ErrorCode, MB_CHK_SET_ERR, MB_SET_ERR, MB_SUCCESS, MB_TAG_ANY, mbImpl, moab::Range::size(), moab::Interface::tag_iterate(), and moab::DamselUtil::tagMap.
Referenced by process_ent_info().
|
private |
Definition at line 260 of file ReadDamsel.cpp.
261 {
262 Tag tagh;
263 ErrorCode rval = MB_SUCCESS, tmp_rval;
264 for( std::vector< damsel_tag_buf_type >::iterator tit = tag_infos.begin(); tit != tag_infos.end(); ++tit )
265 {
266 if( DamselUtil::dtom_data_type[( *tit ).tag_datatype] == MB_TYPE_OPAQUE )
267 {
268 std::cout << "Damsel reader encountered opaque tag." << std::endl;
269 continue;
270 }
271
272 tmp_rval = mbImpl->tag_get_handle( ( *tit ).name, 1, DamselUtil::dtom_data_type[( *tit ).tag_datatype], tagh,
273 MB_TAG_CREAT | MB_TAG_DENSE );
274 if( MB_SUCCESS != tmp_rval )
275 rval = tmp_rval;
276 else
277 {
278 dU.tagMap.push_back( DamselUtil::tinfo( tagh, 0, MB_TAG_DENSE ) );
279 // Also store predefined tags specially...
280 if( !strncmp( ( *tit ).name, "mbdmsl_", 7 ) )
281 {
282 // Predefined tag name, store the handle
283 if( !strcmp( ( *tit ).name, "mbdmsl_XCOORDS" ) )
284 dU.xcoordsTag = dU.tagMap.back();
285 else if( !strcmp( ( *tit ).name, "mbdmsl_YCOORDS" ) )
286 {
287 dU.ycoordsTag = dU.tagMap.back();
288 }
289 else if( !strcmp( ( *tit ).name, "mbdmsl_ZCOORDS" ) )
290 {
291 dU.zcoordsTag = dU.tagMap.back();
292 }
293 else if( !strcmp( ( *tit ).name, "mbdmsl_COLL_FLAGS" ) )
294 {
295 dU.collFlagsTag = dU.tagMap.back();
296 }
297 else if( !strcmp( ( *tit ).name, "mbdmsl_PARENTS" ) )
298 {
299 dU.parentsTag = dU.tagMap.back();
300 }
301 else if( !strcmp( ( *tit ).name, "mbdmsl_CHILDREN" ) )
302 {
303 dU.childrenTag = dU.tagMap.back();
304 }
305 else
306 {
307 rval = MB_FAILURE;
308 continue;
309 }
310 }
311 }
312 }
313
314 return rval;
315 }
References moab::DamselUtil::childrenTag, moab::DamselUtil::collFlagsTag, moab::DamselUtil::dtom_data_type, dU, ErrorCode, MB_SUCCESS, MB_TAG_CREAT, MB_TAG_DENSE, MB_TYPE_OPAQUE, mbImpl, moab::DamselUtil::parentsTag, moab::Interface::tag_get_handle(), moab::DamselUtil::tagMap, moab::DamselUtil::xcoordsTag, moab::DamselUtil::ycoordsTag, and moab::DamselUtil::zcoordsTag.
Referenced by load_file().
|
virtual |
Read tag values from a file.
Read the list if all integer tag values from the file for a tag that is a single integer value per entity.
| file_name | The file to read. |
| tag_name | The tag for which to read values |
| tag_values_out | Output: The list of tag values. |
| subset_list | An array of tag name and value sets specifying the subset of the file to read. If multiple tags are specified, the sets that match all tags (intersection) should be read. |
| subset_list_length | The length of the 'subset_list' array. |
Implements moab::ReaderIface.
Definition at line 251 of file ReadDamsel.cpp.
256 {
257 return MB_FAILURE;
258 }
|
private |
map from damsel to moab handles
Definition at line 122 of file ReadDamsel.hpp.
Referenced by get_contents(), and process_ent_info().
|
private |
Keep various damsel data.
Definition at line 125 of file ReadDamsel.hpp.
Referenced by load_file(), process_ent_info(), process_entity_tags(), process_tags(), and ~ReadDamsel().
|
private |
File name.
Definition at line 110 of file ReadDamsel.hpp.
|
private |
Interface instance.
Definition at line 104 of file ReadDamsel.hpp.
Referenced by init(), process_entity_tags(), process_tags(), and ~ReadDamsel().
|
private |
Used to track file handles.
Definition at line 119 of file ReadDamsel.hpp.
|
private |
|
private |
Whether this reader natively supports parallel semantics.
Definition at line 113 of file ReadDamsel.hpp.
Referenced by load_file(), and parse_options().
|
private |
UtilIface.
Definition at line 107 of file ReadDamsel.hpp.
Referenced by init(), process_ent_info(), and ~ReadDamsel().
1.9.1.