The test to validate the multi discr per geo type management.

[modules/med.git] / src / MEDLoader / MEDFileField.cxx

diff --git a/src/MEDLoader/MEDFileField.cxx b/src/MEDLoader/MEDFileField.cxx
index 4e461f329d6ca2232588b78f1d9c5db0030e1a28..aa007f4052600bb038776e6edd95381d2a2de548 100644 (file)
--- a/src/MEDLoader/MEDFileField.cxx
+++ b/src/MEDLoader/MEDFileField.cxx
@@ -1537,6 +1537,7 @@ bool MEDFileFieldPerMeshPerType::keepOnlyGaussDiscretization(std::size_t idOfDis
   elt->setNewStart(globalNum);
   globalNum=elt->getEnd();
   its.push_back(bgEnd);
+  _field_pm_pt_pd=newPmPtPd;
   return true;
 }
 
@@ -4248,7 +4249,7 @@ std::vector< MEDCouplingAutoRefCountObjectPtr<MEDFileAnyTypeField1TSWithoutSDA>
       return ret0;
     }
   std::vector< MEDCouplingAutoRefCountObjectPtr<MEDFileAnyTypeField1TSWithoutSDA> > ret(nbOfMDPGT);
-  for(std::size_t i=0;i<allEnt.size();i++)
+  for(std::size_t i=0;i<nbOfMDPGT;i++)
     {
       std::vector< std::pair<int,int> > its;
       ret[i]=shallowCpy();
@@ -7714,7 +7715,7 @@ std::vector< MEDCouplingAutoRefCountObjectPtr<MEDFileAnyTypeFieldMultiTSWithoutS
     {
       MEDCouplingAutoRefCountObjectPtr<MEDFileAnyTypeFieldMultiTSWithoutSDA> elt(createNew());
       for(std::size_t j=0;j<sz;j++)
-        elt->pushBackTimeStep(items[i][j]);
+        elt->pushBackTimeStep(items[j][i]);
       ret[i]=elt;
       elt->MEDFileFieldNameScope::operator=(*this);
     }