// File : DriverCGNS_Read.cxx
// Created : Thu Jun 30 10:33:31 2011
// Author : Edward AGAPOV (eap)
+#define _DEBUG_
+#include <utilities.h>
#include "DriverCGNS_Read.hxx"
#include <map>
+
#if CGNS_VERSION < 3100
# define cgsize_t int
#endif
Driver_Mesh::Status DriverCGNS_Read::Perform()
{
+ MESSAGE("DriverCGNS_Read::Perform");
myErrorMessages.clear();
Status aResult;
// read nb of meshes (CGNSBase_t)
if ( myMeshId < 0 || myMeshId >= GetNbMeshes(aResult))
return addMessage( SMESH_Comment("Invalid mesh index :") << myMeshId );
+ MESSAGE("NbMeshes: " << GetNbMeshes(aResult));
// read a name and a dimension of the mesh
const int cgnsBase = myMeshId + 1;
<< " in mesh '" << meshName << "'");
myMeshName = meshName;
+ MESSAGE("myMeshName: " << myMeshName);
+
// read nb of domains (Zone_t) in the mesh
int nbZones = 0;
if ( nbZones < 1 )
return addMessage( SMESH_Comment("Empty mesh: '") << meshName << "'");
+ MESSAGE("nbZones: " << nbZones);
// read the domains (zones)
// ------------------------
zone._nodeIdShift = meshInfo.NbNodes();
zone._elemIdShift = meshInfo.NbElements();
zone.SetSizeAndDim( sizes, meshDim );
+ MESSAGE(" zone name: " << name);
// mesh type of the zone
if ( cg_zone_type ( _fn, cgnsBase, iZone, &zone._type) != CG_OK) {
switch ( zone._type )
{
case CGNS_ENUMV( Unstructured ):
+ MESSAGE(" zone type: unstructured");
+ break;
case CGNS_ENUMV( Structured ):
+ MESSAGE(" zone type: structured");
break;
case CGNS_ENUMV( ZoneTypeNull ):
addMessage( "Meshes with ZoneTypeNull are not supported");
// -----------
// Read nodes
// -----------
-
+ MESSAGE(" Read nodes");
if ( cg_ncoords( _fn, cgnsBase, iZone, &spaceDim) != CG_OK ) {
addMessage( cg_get_error() );
continue;
}
// read coordinates
+ MESSAGE(" Read coordinates");
cgsize_t rmin[3] = {1,1,1}; // range of nodes to read
cgsize_t rmax[3] = {1,1,1};
int nbNodes = rmax[0] = zone._sizes[0];
coords[ c-1 ].resize( nbNodes, 0.0 );
// create nodes
+ MESSAGE(" create nodes");
try {
for ( int i = 0; i < nbNodes; ++i )
myMesh->AddNodeWithID( coords[0][i], coords[1][i], coords[2][i], i+1+zone._nodeIdShift );
// --------------
// Read elements
// --------------
+ MESSAGE(" read elements");
if ( zone.IsStructured())
{
int nbI = zone._sizeX - 1, nbJ = zone._sizeY - 1, nbK = zone._sizeZ - 1;
// read element data
+ MESSAGE(" read element data");
CGNS_ENUMT( ElementType_t ) elemType;
cgsize_t start, end; // range of ids of elements of a zone
cgsize_t eDataSize = 0;
int nbBnd, parent_flag;
for ( int iSec = 1; iSec <= nbSections; ++iSec )
{
+ MESSAGE(" section " << iSec << " of " << nbSections);
if ( cg_section_read( _fn, cgnsBase, iZone, iSec, name, &elemType,
&start, &end, &nbBnd, &parent_flag) != CG_OK ||
cg_ElementDataSize( _fn, cgnsBase, iZone, iSec, &eDataSize ) != CG_OK )
}
// store elements
+ MESSAGE(" store elements");
int pos = 0, cgnsNbNodes = 0, elemID = start + zone._elemIdShift;
cg_npe( elemType, &cgnsNbNodes ); // get nb nodes by element type
curAddElemFun = addElemFuns[ elemType ];
// -------------------------------------------
// Read Boundary Conditions into SMESH groups
// -------------------------------------------
+
+ MESSAGE(" read Boundary Conditions");
int nbBC = 0;
if ( cg_nbocos( _fn, cgnsBase, iZone, &nbBC) == CG_OK )
{
CGNS_ENUMT( DataType_t ) normDataType;
cgsize_t nbPnt, normFlag;
int normIndex[3], nbDS;
+ MESSAGE(" nbBC: " << nbBC);
for ( int iBC = 1; iBC <= nbBC; ++iBC )
{
+ MESSAGE(" iBC: " << iBC);
if ( cg_boco_info( _fn, cgnsBase, iZone, iBC, name, &bcType, &psType,
&nbPnt, normIndex, &normFlag, &normDataType, &nbDS ) != CG_OK )
{
addMessage( cg_get_error() );
continue;
}
+ MESSAGE(" iBC info OK: " << iBC);
vector< cgsize_t > ids( nbPnt * zone.IndexSize() );
CGNS_ENUMT( GridLocation_t ) location;
if ( cg_boco_read( _fn, cgnsBase, iZone, iBC, &ids[0], NULL ) != CG_OK ||
} // loop on BCs of the zone
}
- else
+ else addMessage( cg_get_error() );
+
+
+ MESSAGE(" read flow solutions");
+ int nsols = 0;
+ if ( cg_nsols( _fn, cgnsBase, iZone, &nsols) == CG_OK )
{
- addMessage( cg_get_error() );
+ MESSAGE(" nb flow solutions: " << nsols);
+ }
+ else addMessage( cg_get_error() );
+
+ MESSAGE(" read discrete data");
+ int nbdiscrete = 0;
+ if ( cg_ndiscrete( _fn, cgnsBase, iZone, &nbdiscrete) == CG_OK )
+ {
+ MESSAGE(" nb discrete data: " << nbdiscrete);
+ char nameDiscrete[CGNS_NAME_SIZE];
+ for (int idisc = 1; idisc <= nbdiscrete; idisc++)
+ {
+ if ( cg_discrete_read( _fn, cgnsBase, iZone, idisc, nameDiscrete) == CG_OK )
+ {
+ MESSAGE(" discrete data #"<< idisc << " name: " << nameDiscrete);
+ PointSetType_t ptset_type;
+ cgsize_t npnts;
+ if ( cg_discrete_ptset_info( _fn, cgnsBase, iZone, idisc, &ptset_type, &npnts) == CG_OK )
+ {
+ MESSAGE(" discrete data #"<< idisc << " npnts: " << npnts);
+ }
+ else addMessage( cg_get_error() );
+ }
+ else addMessage( cg_get_error() );
+ }
+ }
+ else addMessage( cg_get_error() );
+
+
+ MESSAGE(" read subregions");
+ int nbSubrg = 0;
+ if ( cg_nsubregs( _fn, cgnsBase, iZone, &nbSubrg) == CG_OK )
+ {
+ MESSAGE(" nb subregions: " << nbSubrg);
}
+ else addMessage( cg_get_error() );
+
+ MESSAGE(" end zone");
} // loop on the zones of a mesh
+ MESSAGE("read families");
+ int nbFam = 0;
+ if ( cg_nfamilies( _fn, cgnsBase, &nbFam) == CG_OK )
+ {
+ MESSAGE("nb families: " << nbFam);
+ }
+ else addMessage( cg_get_error() );
+
+
// ------------------------------------------------------------------------
// Make groups for multiple zones and remove free nodes at zone interfaces
map< string, TZoneData >::iterator nameZoneIt = zonesByName.begin();
for ( ; nameZoneIt != zonesByName.end(); ++nameZoneIt )
{
+ MESSAGE("nameZone: " << nameZoneIt->first);
TZoneData& zone = nameZoneIt->second;
if ( zone._nbElems == 0 ) continue;
if ( zone._nbElems == meshInfo.NbElements() ) break; // there is only one non-empty zone
myMesh->Modified();
myMesh->CompactMesh();
-
+ MESSAGE("end perform");
return aResult;
}