return pfl->getName();
}
-int MEDMeshMultiLev::getNumberOfCells() const throw(INTERP_KERNEL::Exception)
+/*!
+ * Returns the number of cells having geometric type \a t.
+ * The profiles are **NOT** taken into account here.
+ */
+int MEDMeshMultiLev::getNumberOfCells(INTERP_KERNEL::NormalizedCellType t) const throw(INTERP_KERNEL::Exception)
{
- std::size_t sz(_nb_entities.size()/3);
- if(_nb_entities.size()%3!=0)
- throw INTERP_KERNEL::Exception("MEDMeshMultiLev::getNumberOfCells : invalid code !");
- int ret(0);
+ std::size_t sz(_nb_entities.size());
for(std::size_t i=0;i<sz;i++)
- {
- if(_nb_entities[3*i+2]==-1)
- ret+=_nb_entities[3*i+1];
- else
- {
- const DataArrayInt *pfl(_pfls[_nb_entities[3*i+2]]);
- if(!pfl)
- throw INTERP_KERNEL::Exception("MEDMeshMultiLev::getNumberOfCells : internal error !");
- ret+=pfl->getNumberOfTuples();
- }
- }
- return ret;
+ if(_geo_types[i]==t)
+ return _nb_entities[i];
+ throw INTERP_KERNEL::Exception("MEDMeshMultiLev::getNumberOfCells : not existing geometric type in this !");
}
int MEDMeshMultiLev::getNumberOfNodes() const throw(INTERP_KERNEL::Exception)
else
{
std::size_t sz(fst.getNumberOfItems());
- std::vector< MEDCouplingAutoRefCountObjectPtr<DataArray> > arrSafe(sz);
- std::vector< const DataArray *> arr(sz);
- for(std::size_t i=0;i<sz;i++)
+ std::set<INTERP_KERNEL::NormalizedCellType> s(_geo_types.begin(),_geo_types.end());
+ if(s.size()!=_geo_types.size())
+ throw INTERP_KERNEL::Exception("MEDMeshMultiLev::constructDataArray : unexpected situation for cells 2 !");
+ std::vector< const DataArray *> arr(s.size());
+ std::vector< MEDCouplingAutoRefCountObjectPtr<DataArray> > arrSafe(s.size());
+ int iii(0);
+ int nc(vals->getNumberOfComponents());
+ std::vector<std::string> compInfo(vals->getInfoOnComponents());
+ for(std::vector< INTERP_KERNEL::NormalizedCellType >::const_iterator it=_geo_types.begin();it!=_geo_types.end();it++,iii++)
{
- const MEDFileField1TSStructItem2& p(fst[i]);
- const std::pair<int,int>& strtStop(p.getStartStop());
- std::vector< INTERP_KERNEL::NormalizedCellType >::const_iterator it(std::find(_geo_types.begin(),_geo_types.end(),p.getGeo()));
- if(it==_geo_types.end())
- throw INTERP_KERNEL::Exception("MEDMeshMultiLev::constructDataArray : unexpected situation for cells 1 !");
- if(std::find(it+1,_geo_types.end(),p.getGeo())!=_geo_types.end())
- throw INTERP_KERNEL::Exception("MEDMeshMultiLev::constructDataArray : unexpected situation for cells 2 !");
- std::size_t pos(std::distance(_geo_types.begin(),it));
- const DataArrayInt *thisP(_pfls[pos]),*otherP(p.getPfl(globs));
- MEDCouplingAutoRefCountObjectPtr<DataArray> ret(vals->selectByTupleId2(strtStop.first,strtStop.second,1));
- if(!thisP && !otherP)
+ const DataArrayInt *thisP(_pfls[iii]);
+ std::vector<const MEDFileField1TSStructItem2 *> ps;
+ for(std::size_t i=0;i<sz;i++)
{
- arrSafe[i]=ret; arr[i]=ret;
- continue;
- }
- int nbi(p.getNbOfIntegrationPts(globs)),nc(ret->getNumberOfComponents());
- if(!thisP && otherP)
- {
- MEDCouplingAutoRefCountObjectPtr<DataArrayInt> p1(otherP->deepCpy());
- p1->sort(true);
- if(!p1->isIdentity() || p1->getNumberOfTuples()!=getNumberOfCells())
- throw INTERP_KERNEL::Exception("MEDMeshMultiLev::constructDataArray : unexpected situation for cells 3 !");
- ret->rearrange(nbi*nc); ret->renumberInPlace(otherP->begin()); ret->rearrange(nc);
- arrSafe[i]=ret; arr[i]=ret;
- continue;
+ const MEDFileField1TSStructItem2& p(fst[i]);
+ if(p.getGeo()==*it)
+ ps.push_back(&p);
}
- if(thisP && otherP)
+ if(ps.empty())
+ throw INTERP_KERNEL::Exception("MEDMeshMultiLev::constructDataArray : unexpected situation for cells 1 !");
+ if(ps.size()==1)
{
- MEDCouplingAutoRefCountObjectPtr<DataArrayInt> p1(otherP->deepCpy());
- MEDCouplingAutoRefCountObjectPtr<DataArrayInt> p2(thisP->deepCpy());
- p1->sort(true); p2->sort(true);
- if(!p1->isEqualWithoutConsideringStr(*p2))
- throw INTERP_KERNEL::Exception("MEDMeshMultiLev::constructDataArray : unexpected situation for cells 4 !");
- p1=DataArrayInt::FindPermutationFromFirstToSecond(otherP,thisP);
- ret->rearrange(nbi*nc); ret->renumberInPlace(p1->begin()); ret->rearrange(nc);
- continue;
+ int nbi(ps[0]->getNbOfIntegrationPts(globs));
+ const DataArrayInt *otherP(ps[0]->getPfl(globs));
+ const std::pair<int,int>& strtStop(ps[0]->getStartStop());
+ MEDCouplingAutoRefCountObjectPtr<DataArray> ret(vals->selectByTupleId2(strtStop.first,strtStop.second,1));
+ if(!thisP && !otherP)
+ {
+ arrSafe[iii]=ret; arr[iii]=ret;
+ continue;
+ }
+ if(thisP && otherP)
+ {
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> p1(otherP->invertArrayN2O2O2N(getNumberOfCells(ps[0]->getGeo())));
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> p2(thisP->deepCpy());
+ p2->transformWithIndArr(p1->begin(),p1->end());
+ //p1=p2->getIdsNotEqual(-1);
+ //p1=p2->selectByTupleIdSafe(p1->begin(),p1->end());
+ ret->rearrange(nbi*nc); ret=ret->selectByTupleIdSafe(p2->begin(),p2->end()); ret->rearrange(nc); ret->setInfoOnComponents(compInfo);
+ arrSafe[iii]=ret; arr[iii]=ret;
+ continue;
+ }
+ if(!thisP && otherP)
+ {
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> p1(otherP->deepCpy());
+ p1->sort(true);
+ p1->checkAllIdsInRange(0,getNumberOfCells(ps[0]->getGeo()));
+ p1=DataArrayInt::FindPermutationFromFirstToSecond(otherP,p1);
+ ret->rearrange(nbi*nc); ret->renumberInPlace(p1->begin()); ret->rearrange(nc); ret->setInfoOnComponents(compInfo);
+ arrSafe[iii]=ret; arr[iii]=ret;
+ continue;
+ }
+ throw INTERP_KERNEL::Exception("MEDMeshMultiLev::constructDataArray : unexpected situation for cells 3 !");
}
- if(thisP && !otherP)
+ else
{
- MEDCouplingAutoRefCountObjectPtr<DataArrayInt> p1(thisP->deepCpy());
- p1->sort(true);
- if(!p1->isIdentity() || p1->getNumberOfTuples()!=p.getNbEntity())
- throw INTERP_KERNEL::Exception("MEDMeshMultiLev::constructDataArray : unexpected situation for cells 3 !");
- ret->rearrange(nbi*nc); ret->renumberInPlaceR(otherP->begin()); ret->rearrange(nc);
- arrSafe[i]=ret; arr[i]=ret;
+ std::vector< const DataArrayInt * >otherPS(ps.size());
+ std::vector< const DataArray * > arr2(ps.size());
+ std::vector< MEDCouplingAutoRefCountObjectPtr<DataArray> > arr2Safe(ps.size());
+ std::vector< const DataArrayInt * > nbis(ps.size());
+ std::vector< MEDCouplingAutoRefCountObjectPtr<DataArrayInt> > nbisSafe(ps.size());
+ int jj(0);
+ for(std::vector<const MEDFileField1TSStructItem2 *>::const_iterator it2=ps.begin();it2!=ps.end();it2++,jj++)
+ {
+ int nbi((*it2)->getNbOfIntegrationPts(globs));
+ const DataArrayInt *otherPfl((*it2)->getPfl(globs));
+ const std::pair<int,int>& strtStop((*it2)->getStartStop());
+ MEDCouplingAutoRefCountObjectPtr<DataArray> ret2(vals->selectByTupleId2(strtStop.first,strtStop.second,1));
+ if(!otherPfl)
+ throw INTERP_KERNEL::Exception("MEDMeshMultiLev::constructDataArray : unexpected situation for cells 4 !");
+ arr2[jj]=ret2; arr2Safe[jj]=ret2; otherPS[jj]=otherPfl;
+ nbisSafe[jj]=DataArrayInt::New(); nbisSafe[jj]->alloc(otherPfl->getNumberOfTuples(),1); nbisSafe[jj]->fillWithValue(nbi);
+ nbis[jj]=nbisSafe[jj];
+ }
+ MEDCouplingAutoRefCountObjectPtr<DataArray> arr3(DataArray::Aggregate(arr2));
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> otherP(DataArrayInt::Aggregate(otherPS));
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> zenbis(DataArrayInt::Aggregate(nbis));
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> otherPN(otherP->invertArrayN2O2O2N(getNumberOfCells(*it)));
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> p1;
+ if(thisP)
+ p1=DataArrayInt::FindPermutationFromFirstToSecond(otherP,thisP);
+ else
+ p1=otherP->deepCpy();
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> zenbisN(zenbis->renumber(p1->begin()));
+ zenbisN->computeOffsets2();
+ jj=0;
+ for(std::vector<const MEDFileField1TSStructItem2 *>::const_iterator it2=ps.begin();it2!=ps.end();it2++,jj++)
+ {
+ int nbi((*it2)->getNbOfIntegrationPts(globs));
+ const DataArrayInt *otherPfl((*it2)->getPfl(globs));
+ const std::pair<int,int>& strtStop((*it2)->getStartStop());
+ MEDCouplingAutoRefCountObjectPtr<DataArray> ret2(vals->selectByTupleId2(strtStop.first,strtStop.second,1));
+ //
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> p2(otherPfl->deepCpy());
+ p2->transformWithIndArr(otherPN->begin(),otherPN->end());
+ p2->transformWithIndArr(p1->begin(),p1->end());
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> idsN(p2->buildExplicitArrByRanges(zenbisN));
+ arr3->setPartOfValuesBase3(ret2,idsN->begin(),idsN->end(),0,nc,1);
+ }
+ arrSafe[iii]=arr3; arr[iii]=arr3;
continue;
}
- throw INTERP_KERNEL::Exception("MEDMeshMultiLev::constructDataArray : unexpected situation for cells 6 !");
}
return DataArray::Aggregate(arr);
}
MEDCouplingAutoRefCountObjectPtr<DataArrayInt> cellIdsSafe(cellIds);
MEDCouplingAutoRefCountObjectPtr<MEDCouplingPointSet> m2(m->buildPartOfMySelfKeepCoords(cellIds->begin(),cellIds->end()));
int tmp=-1;
- a[i]=m2->getNodeIdsInUse(tmp); aa[i]=a[i];
+ MEDCouplingAutoRefCountObjectPtr<DataArrayInt> o2n(m2->getNodeIdsInUse(tmp));
+ a[i]=o2n->invertArrayO2N2N2O(tmp); aa[i]=a[i];
if(pfl)
_pfls[i]=pfl->selectByTupleIdSafe(cellIds->begin(),cellIds->end());
else
{
for(int i=0;i<curNbCells;i++,connIPtr++)
{
- std::set<int> s(connPtr+connIPtr[0],connPtr+connIPtr[1]);
+ std::set<int> s(connPtr+connIPtr[0],connPtr+connIPtr[1]); s.erase(-1);
*dPtr++=(int)s.size();
dPtr=std::copy(s.begin(),s.end(),dPtr);
*cPtr=k+(int)s.size(); k=*cPtr++;
fPtr=std::copy(work,work2,fPtr);
work=work2+1;
}
- *ePtr=kk; kk+=connIPtr[1]-connIPtr[0]+2;
+ *ePtr++=kk; kk+=connIPtr[1]-connIPtr[0]+2;
}
}
}