aGauss->myGaussInfo = aGaussInfo;
aGauss->myName = aGaussInfo->GetName();
aGauss->myNbPoints = aGaussInfo->GetNbGauss();
- }
+ } else if ( aTimeStampInfo.GetNbGauss(aMGeom) > 1 )
+ EXCEPTION( runtime_error, "aGaussIter == aGeom2Gauss.end()"
+ " && aTimeStampInfo.GetNbGauss(aMGeom) > 1 !!!" );
aGaussSubMesh->myNbCells = aNbCells*aGauss->myNbPoints;
aGaussSubMesh->myCellsSize = aGaussSubMesh->myNbCells*2;
}
- //---------------------------------------------------------------
- void
- InitGaussProfile(const MED::PWrapper& theMEDWrapper,
- const MED::PMeshInfo& theMeshInfo,
- MED::TTimeStampVal& theTimeStampVal,
- VISU::TMEDMeshOnEntity& theMeshOnEntity,
- MED::EEntiteMaillage theMEntity,
- const MED::TGeom2Size& theGeom2Size,
- VISU::TMEDValForTime& theValForTime)
- {
- TTimerLog aTimerLog(MYDEBUG,"InitGaussProfile");
- INITMSG(MYDEBUG,"InitGaussProfile"<<endl);
-
- // The order of the function calls is important
- InitProfile(theMEDWrapper,
- theMeshInfo,
- theTimeStampVal,
- theMeshOnEntity,
- theMEntity,
- theGeom2Size,
- theValForTime);
-
- InitGaussMesh(theTimeStampVal,
- theMeshOnEntity,
- theGeom2Size,
- theValForTime);
- }
-
-
//---------------------------------------------------------------
void
BuildMeshOnEntityMap(PMEDMesh theMesh,
theMed->GetPTimeStampVal(aTimeStampInfo,
aMKey2Profile,
aKey2Gauss);
-
- InitGaussProfile(theMed,
- aMeshInfo,
- aTimeStampVal,
- theMeshOnEntity,
- aMEntity,
- aGeom2Size,
- theValForTime);
+
+ InitProfile(theMed,
+ aMeshInfo,
+ aTimeStampVal,
+ theMeshOnEntity,
+ aMEntity,
+ aGeom2Size,
+ theValForTime);
LoadProfile(theMed,
theMesh,
theValForTime,
theMeshOnEntity);
- if(theIsGauss)
+ if ( theIsGauss ) {
+ InitGaussMesh(aTimeStampVal,
+ theMeshOnEntity,
+ aGeom2Size,
+ theValForTime);
+
LoadGaussMesh(theMed,
theMesh,
aTimeStampVal,
theValForTime,
theMeshOnEntity);
-
+ }
+
PMEDProfile aProfile = theValForTime->myProfile;
TGeom2SubProfile& aGeom2SubProfile = aProfile->myGeom2SubProfile;