From d592d340b70baa1b23cd76f4eaac72fa7785cb18 Mon Sep 17 00:00:00 2001
From: apo <apo@opencascade.com>
Date: Fri, 22 Jul 2005 06:53:44 +0000
Subject: [PATCH] Gauss Point reading. Next step

---
 src/CONVERTOR/VISU_MedConvertor.cxx | 55 +++++++++++++++++++++++++----
 1 file changed, 48 insertions(+), 7 deletions(-)

diff --git a/src/CONVERTOR/VISU_MedConvertor.cxx b/src/CONVERTOR/VISU_MedConvertor.cxx
index f9d9dbdb..7f9d742f 100644
--- a/src/CONVERTOR/VISU_MedConvertor.cxx
+++ b/src/CONVERTOR/VISU_MedConvertor.cxx
@@ -291,6 +291,8 @@ namespace{
 	if(aGaussIter != aGeom2Gauss.end()){
 	  PGaussInfo aGaussInfo = aGaussIter->second;
 
+	  aGauss->myGaussInfo = aGaussInfo;
+
 	  aName = aGaussInfo->GetName();
 	  aGauss->myName = aName;
 
@@ -364,6 +366,7 @@ namespace{
 		   const MED::TGeom2Size& theGeom2Size,
 		   VISU::TMEDValForTime& theValForTime)
   {
+    INITMSG(MYDEBUG,"InitGaussProfile"<<endl);
     // The order of the function calls is important
     InitProfile(theTimeStampVal,theMeshOnEntity,theGeom2Size,theValForTime);
     InitGaussMesh(theTimeStampVal,theMeshOnEntity,theGeom2Size,theValForTime);
@@ -1114,7 +1117,9 @@ LoadProfile(MED::TTimeStampVal& theTimeStampVal,
 
 //---------------------------------------------------------------
 void
-LoadGaussMesh(MED::TTimeStampVal& theTimeStampVal,
+LoadGaussMesh(const MED::PWrapper& theMed,
+	      VISU::PMEDMesh theMesh,
+	      MED::TTimeStampVal& theTimeStampVal,
 	      VISU::TMEDValForTime& theValForTime,
 	      VISU::TMEDMeshOnEntity& theMeshOnEntity)
 {
@@ -1124,6 +1129,12 @@ LoadGaussMesh(MED::TTimeStampVal& theTimeStampVal,
   if(aGaussMesh->myIsInitialized)
     return;
 
+  const PMeshInfo& aMeshInfo = theMesh->myMeshInfo;
+  PNodeInfo aNodeInfo = theMed->GetPNodeInfo(aMeshInfo);
+
+  TEntity aVEntity = theMeshOnEntity.myEntity;
+  EEntiteMaillage aMEntity = VTKEntityToMED(aVEntity);
+
   const TGeom2GaussSubMesh& aGeom2GaussSubMesh = aGaussMesh->myGeom2GaussSubMesh;
   const TTimeStampInfo& aTimeStampInfo = theTimeStampVal.GetTimeStampInfo();
   const TGeom2Gauss& aGeom2Gauss = aTimeStampInfo.GetGeom2Gauss();
@@ -1136,13 +1147,41 @@ LoadGaussMesh(MED::TTimeStampVal& theTimeStampVal,
     TGeom2GaussSubMesh::const_iterator anIter2 = aGeom2GaussSubMesh.find(aVGeom);
     if(anIter2 != aGeom2GaussSubMesh.end()){
       PMEDGaussSubMesh aGaussSubMesh = anIter2->second;
+
+      PMEDGauss aGauss = aGaussSubMesh->myGauss;
+      MED::PGaussInfo aGaussInfo = aGauss->myGaussInfo;
       PMEDSubProfile aSubProfile = aGaussSubMesh->mySubProfile;
+      
+      if(aGaussInfo){
+	MED::PCellInfo aCellInfo = theMed->GetPCellInfo(aMeshInfo,
+							aMEntity, 
+							aMGeom);
+	TElemNum anElemNum = aSubProfile->mySubMeshID;
+	
+	MED::TGaussCoord aGaussCoord;
+	MED::GetGaussCoord3D(aGaussInfo,
+			     aCellInfo,
+			     aNodeInfo,
+			     aGaussCoord,
+			     anElemNum);
+
+	TPoints& aPoints = aGaussSubMesh->myPoints;
+
+	INITMSG(MYDEBUG,
+		"- aVGeom = "<<aVGeom<<
+		"; aStatus = "<<aGaussSubMesh->myStatus<<
+		"; aNbCells = "<<aSubProfile->mySubMeshID.size()<<
+		"; aNbElem = "<<aGaussCoord.size()<<
+		endl);
+	
+      }else{
 
-      INITMSG(MYDEBUG,
-	      "- aVGeom = "<<aVGeom<<
-	      "; aStatus = "<<aGaussSubMesh->myStatus<<
-	      "; aNbCells = "<<aSubProfile->mySubMeshID.size()<<
-	      endl);
+	INITMSG(MYDEBUG,
+		"- aVGeom = "<<aVGeom<<
+		"; aStatus = "<<aGaussSubMesh->myStatus<<
+		"; aNbCells = "<<aSubProfile->mySubMeshID.size()<<
+		endl);
+      }
     }
   }
 }
@@ -1194,7 +1233,9 @@ VISU_MedConvertor
 	      theValForTime,
 	      theMeshOnEntity);
   
-  LoadGaussMesh(aTimeStampVal,
+  LoadGaussMesh(theMed,
+		theMesh,
+		aTimeStampVal,
 		theValForTime,
 		theMeshOnEntity);
   
-- 
2.39.2