1 // Copyright (C) 2007-2020 CEA/DEN, EDF R&D
3 // This library is free software; you can redistribute it and/or
4 // modify it under the terms of the GNU Lesser General Public
5 // License as published by the Free Software Foundation; either
6 // version 2.1 of the License, or (at your option) any later version.
8 // This library is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 // Lesser General Public License for more details.
13 // You should have received a copy of the GNU Lesser General Public
14 // License along with this library; if not, write to the Free Software
15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
19 // Author : Anthony Geay (CEA/DEN)
21 #include "MEDCouplingPointSet.hxx"
23 #include "MEDCoupling1GTUMesh.hxx"
24 #include "MEDCouplingUMesh.hxx"
25 #include "MEDCouplingMemArray.hxx"
26 #include "PlanarIntersector.txx"
27 #include "InterpKernelGeo2DQuadraticPolygon.hxx"
28 #include "InterpKernelGeo2DNode.hxx"
29 #include "DirectedBoundingBox.hxx"
30 #include "InterpKernelAutoPtr.hxx"
37 using namespace MEDCoupling;
39 MEDCouplingPointSet::MEDCouplingPointSet():_coords(0)
43 MEDCouplingPointSet::MEDCouplingPointSet(const MEDCouplingPointSet& other, bool deepCpy):MEDCouplingMesh(other),_coords(0)
46 _coords=other._coords->performCopyOrIncrRef(deepCpy);
49 MEDCouplingPointSet::~MEDCouplingPointSet()
55 mcIdType MEDCouplingPointSet::getNumberOfNodes() const
58 return _coords->getNumberOfTuples();
60 throw INTERP_KERNEL::Exception("Unable to get number of nodes because no coordinates specified !");
63 int MEDCouplingPointSet::getSpaceDimension() const
66 return (int)_coords->getNumberOfComponents();
68 throw INTERP_KERNEL::Exception("Unable to get space dimension because no coordinates specified !");
71 void MEDCouplingPointSet::updateTime() const
75 updateTimeWith(*_coords);
79 std::size_t MEDCouplingPointSet::getHeapMemorySizeWithoutChildren() const
81 return MEDCouplingMesh::getHeapMemorySizeWithoutChildren();
84 std::vector<const BigMemoryObject *> MEDCouplingPointSet::getDirectChildrenWithNull() const
86 std::vector<const BigMemoryObject *> ret;
87 ret.push_back(_coords);
91 void MEDCouplingPointSet::setCoords(const DataArrayDouble *coords)
93 if( coords != _coords )
97 _coords=const_cast<DataArrayDouble *>(coords);
105 * Returns a pointer to the array of point coordinates held by \a this.
106 * \return DataArrayDouble * - the pointer to the array of point coordinates. The
107 * caller is to delete this array using decrRef() as it is no more needed.
109 DataArrayDouble *MEDCouplingPointSet::getCoordinatesAndOwner() const
117 * Copies string attributes from an \a other mesh. The copied strings are
121 * - textual data of the coordinates array (name and components info)
123 * \param [in] other - the mesh to copy string attributes from.
125 void MEDCouplingPointSet::copyTinyStringsFrom(const MEDCouplingMesh *other)
127 MEDCouplingMesh::copyTinyStringsFrom(other);
128 const MEDCouplingPointSet *otherC=dynamic_cast<const MEDCouplingPointSet *>(other);
130 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::copyTinyStringsFrom : meshes have not same type !");
131 if(_coords && otherC->_coords)
132 _coords->copyStringInfoFrom(*otherC->_coords);
135 bool MEDCouplingPointSet::isEqualIfNotWhy(const MEDCouplingMesh *other, double prec, std::string& reason) const
138 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::isEqualIfNotWhy : null mesh instance in input !");
139 const MEDCouplingPointSet *otherC=dynamic_cast<const MEDCouplingPointSet *>(other);
142 reason="mesh given in input is not castable in MEDCouplingPointSet !";
145 if(!MEDCouplingMesh::isEqualIfNotWhy(other,prec,reason))
147 if(!areCoordsEqualIfNotWhy(*otherC,prec,reason))
153 * Checks equality of point coordinates with coordinates of an \a other mesh.
154 * None textual data is considered.
155 * \param [in] other - the mesh to compare coordinates with \a this one.
156 * \param [in] prec - precision value to compare coordinates.
157 * \return bool - \a true if coordinates of points are equal, \a false else.
159 bool MEDCouplingPointSet::isEqualWithoutConsideringStr(const MEDCouplingMesh *other, double prec) const
161 const MEDCouplingPointSet *otherC=dynamic_cast<const MEDCouplingPointSet *>(other);
164 if(!areCoordsEqualWithoutConsideringStr(*otherC,prec))
169 bool MEDCouplingPointSet::areCoordsEqualIfNotWhy(const MEDCouplingPointSet& other, double prec, std::string& reason) const
171 if(_coords==0 && other._coords==0)
173 if(_coords==0 || other._coords==0)
175 reason="Only one PointSet between the two this and other has coordinate defined !";
178 if(_coords==other._coords)
180 bool ret=_coords->isEqualIfNotWhy(*other._coords,prec,reason);
182 reason.insert(0,"Coordinates DataArray do not match : ");
187 * Checks equality of point coordinates with \a other point coordinates.
188 * Textual data (name and components info) \b is compared as well.
189 * \param [in] other - the point coordinates to compare with \a this one.
190 * \param [in] prec - precision value to compare coordinates.
191 * \return bool - \a true if coordinates of points are equal, \a false else.
193 bool MEDCouplingPointSet::areCoordsEqual(const MEDCouplingPointSet& other, double prec) const
196 return areCoordsEqualIfNotWhy(other,prec,tmp);
200 * Checks equality of point coordinates with \a other point coordinates.
201 * None textual data is considered.
202 * \param [in] other - the point coordinates to compare with \a this one.
203 * \param [in] prec - precision value to compare coordinates.
204 * \return bool - \a true if coordinates of points are equal, \a false else.
206 bool MEDCouplingPointSet::areCoordsEqualWithoutConsideringStr(const MEDCouplingPointSet& other, double prec) const
208 if(_coords==0 && other._coords==0)
210 if(_coords==0 || other._coords==0)
212 if(_coords==other._coords)
214 return _coords->isEqualWithoutConsideringStr(*other._coords,prec);
218 * Returns coordinates of \a nodeId-th node.
219 * \param [in] nodeId - the ID of the node of interest.
220 * \param [in, out] coo - the array filled with coordinates of the \a nodeId-th
221 * node. This array is not cleared before filling in, the coordinates are
222 * appended to its end.
223 * \throw If the coordinates array is not set.
224 * \throw If \a nodeId is not a valid index for the coordinates array.
226 * \if ENABLE_EXAMPLES
227 * \ref cpp_mcpointset_getcoordinatesofnode "Here is a C++ example".<br>
228 * \ref py_mcpointset_getcoordinatesofnode "Here is a Python example".
231 void MEDCouplingPointSet::getCoordinatesOfNode(mcIdType nodeId, std::vector<double>& coo) const
234 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::getCoordinatesOfNode : no coordinates array set !");
235 mcIdType nbNodes=getNumberOfNodes();
236 if(nodeId>=0 && nodeId<nbNodes)
238 const double *cooPtr=_coords->getConstPointer();
239 std::size_t spaceDim=getSpaceDimension();
240 coo.insert(coo.end(),cooPtr+spaceDim*nodeId,cooPtr+spaceDim*(nodeId+1));
244 std::ostringstream oss; oss << "MEDCouplingPointSet::getCoordinatesOfNode : request of nodeId \"" << nodeId << "\" but it should be in [0,"<< nbNodes << ") !";
245 throw INTERP_KERNEL::Exception(oss.str().c_str());
250 * Finds nodes equal within \a precision and returns an array describing the
251 * permutation to remove duplicated nodes.
252 * \param [in] precision - minimal absolute distance between two nodes at which they are
253 * considered not coincident.
254 * \param [in] limitNodeId - limit node id. If all nodes within a group of coincident
255 * nodes have id strictly lower than \a limitTupleId then they are not
256 * returned. Put -1 to this parameter to have all nodes returned.
257 * \param [out] areNodesMerged - is set to \a true if any coincident nodes found.
258 * \param [out] newNbOfNodes - returns number of unique nodes.
259 * \return DataArrayIdType * - the permutation array in "Old to New" mode. For more
260 * info on "Old to New" mode see \ref numbering. The caller
261 * is to delete this array using decrRef() as it is no more needed.
262 * \throw If the coordinates array is not set.
264 DataArrayIdType *MEDCouplingPointSet::buildPermArrayForMergeNode(double precision, mcIdType limitNodeId, bool& areNodesMerged, mcIdType& newNbOfNodes) const
266 DataArrayIdType *comm,*commI;
267 findCommonNodes(precision,limitNodeId,comm,commI);
268 mcIdType oldNbOfNodes=getNumberOfNodes();
269 MCAuto<DataArrayIdType> ret=buildNewNumberingFromCommonNodesFormat(comm,commI,newNbOfNodes);
270 areNodesMerged=(oldNbOfNodes!=newNbOfNodes);
277 * Finds nodes coincident within \a prec tolerance.
278 * Ids of coincident nodes are stored in output arrays in the \ref numbering-indirect format.
279 * \param [in] prec - minimal absolute distance (using infinite norm) between two nodes at which they are
280 * considered not coincident.
281 * \param [in] limitNodeId - limit node id. If all nodes within a group of coincident
282 * nodes have id strictly lower than \a limitTupleId then they are not
283 * returned. Put -1 to this parameter to have all nodes treated.
284 * \param [out] comm - the array holding ids of coincident nodes.
285 * \a comm->getNumberOfComponents() == 1.
286 * \a comm->getNumberOfTuples() == \a commIndex->back(). The caller
287 * is to delete this array using decrRef() as it is no more needed.
288 * \param [out] commIndex - the array dividing all ids stored in \a comm into
289 * groups of (ids of) coincident nodes (\ref numbering-indirect). Its every value is a tuple
290 * index where a next group of nodes begins. For example the second
291 * group of nodes in \a comm is described by following range of indices:
292 * [ \a commIndex[1], \a commIndex[2] ). \a commIndex->getNumberOfTuples()-1
293 * gives the number of groups of coincident nodes. The caller
294 * is to delete this array using decrRef() as it is no more needed.
295 * \throw If the coordinates array is not set.
297 * \if ENABLE_EXAMPLES
298 * \ref cpp_mcpointset_findcommonnodes "Here is a C++ example".<br>
299 * \ref py_mcpointset_findcommonnodes "Here is a Python example".
302 void MEDCouplingPointSet::findCommonNodes(double prec, mcIdType limitNodeId, DataArrayIdType *&comm, DataArrayIdType *&commIndex) const
305 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::findCommonNodes : no coords specified !");
306 _coords->findCommonTuples(prec,limitNodeId,comm,commIndex);
310 * Finds nodes located at distances lower that \a eps from a given point.
311 * \param [in] pos - pointer to coordinates of the point. This array is expected to
312 * be of length \a this->getSpaceDimension() at least, else the
313 * behavior is not warranted.
314 * \param [in] eps - the lowest distance between a point and a node (using infinite norm) at which the node is
315 * not returned by this method.
316 * \return DataArrayIdType * - a new instance of DataArrayIdType holding ids of nodes
317 * close to the point. The caller is to delete this
318 * array using decrRef() as it is no more needed.
319 * \throw If the coordinates array is not set.
321 * \if ENABLE_EXAMPLES
322 * \ref cpp_mcpointset_getnodeidsnearpoint "Here is a C++ example".<br>
323 * \ref py_mcpointset_getnodeidsnearpoint "Here is a Python example".
326 DataArrayIdType *MEDCouplingPointSet::getNodeIdsNearPoint(const double *pos, double eps) const
328 DataArrayIdType *c=0,*cI=0;
329 getNodeIdsNearPoints(pos,1,eps,c,cI);
330 MCAuto<DataArrayIdType> cITmp(cI);
335 * Finds nodes located at distances lower that \a eps from given points.
336 * \param [in] pos - pointer to coordinates of the points. This array is expected to
337 * be of length \a nbOfPoints * \a this->getSpaceDimension() at least, else the
338 * behavior is not warranted.
339 * \param [in] nbOfPoints - number of points whose coordinates are given by \a pos
341 * \param [in] eps - the lowest distance between (using infinite norm) a point and a node at which the node is
342 * not returned by this method.
343 * \param [out] c - array (\ref numbering-indirect) returning ids of nodes located closer than \a eps to the
344 * given points. The caller
345 * is to delete this array using decrRef() as it is no more needed.
346 * \param [out] cI - for each i-th given point, the array specifies tuples of \a c
347 * holding ids of nodes close to the i-th point (\ref numbering-indirect). <br>The i-th value of \a cI is an
348 * index of tuple of \a c holding id of a first (if any) node close to the
349 * i-th given point. Difference between the i-th and (i+1)-th value of \a cI
350 * (i.e. \a cI[ i+1 ] - \a cI[ i ]) defines number of nodes close to the i-th
351 * point (that can be zero!). For example, the group of nodes close to the
352 * second point is described by following range of indices [ \a cI[1], \a cI[2] ).
353 * The caller is to delete this array using decrRef() as it is no more needed.
354 * \throw If the coordinates array is not set.
356 * \if ENABLE_EXAMPLES
357 * \ref cpp_mcpointset_getnodeidsnearpoints "Here is a C++ example".<br>
358 * \ref py_mcpointset_getnodeidsnearpoints "Here is a Python example".
361 void MEDCouplingPointSet::getNodeIdsNearPoints(const double *pos, mcIdType nbOfPoints, double eps, DataArrayIdType *& c, DataArrayIdType *& cI) const
364 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::getNodeIdsNearPoint : no coordiantes set !");
365 std::size_t spaceDim=getSpaceDimension();
366 MCAuto<DataArrayDouble> points=DataArrayDouble::New();
367 points->useArray(pos,false,DeallocType::CPP_DEALLOC,nbOfPoints,spaceDim);
368 _coords->computeTupleIdsNearTuples(points,eps,c,cI);
372 * @param comm in param in the same format than one returned by findCommonNodes method (\ref numbering-indirect).
373 * @param commI in param in the same format than one returned by findCommonNodes method (\ref numbering-indirect).
374 * @return the old to new correspondence array.
376 DataArrayIdType *MEDCouplingPointSet::buildNewNumberingFromCommonNodesFormat(const DataArrayIdType *comm, const DataArrayIdType *commIndex,
377 mcIdType& newNbOfNodes) const
380 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::buildNewNumberingFromCommonNodesFormat : no coords specified !");
381 return DataArrayIdType::ConvertIndexArrayToO2N(getNumberOfNodes(),comm->begin(),commIndex->begin(),commIndex->end(),newNbOfNodes);
385 * Permutes and possibly removes nodes as specified by \a newNodeNumbers array.
386 * If \a newNodeNumbers[ i ] < 0 then the i-th node is removed,
387 * else \a newNodeNumbers[ i ] is a new id of the i-th node. The nodal connectivity
388 * array is modified accordingly.
389 * \param [in] newNodeNumbers - a permutation array, of length \a
390 * this->getNumberOfNodes(), in "Old to New" mode.
391 * See \ref numbering for more info on renumbering modes.
392 * \param [in] newNbOfNodes - number of nodes remaining after renumbering.
393 * \throw If the coordinates array is not set.
394 * \throw If the nodal connectivity of cells is not defined.
396 * \if ENABLE_EXAMPLES
397 * \ref cpp_mcumesh_renumberNodes "Here is a C++ example".<br>
398 * \ref py_mcumesh_renumberNodes "Here is a Python example".
401 void MEDCouplingPointSet::renumberNodes(const mcIdType *newNodeNumbers, mcIdType newNbOfNodes)
404 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::renumberNodes : no coords specified !");
405 MCAuto<DataArrayDouble> newCoords=_coords->renumberAndReduce(newNodeNumbers,newNbOfNodes);
406 renumberNodesInConn(newNodeNumbers);
407 setCoords(newCoords);//let it here not before renumberNodesInConn because old number of nodes is sometimes used...
411 * Permutes and possibly removes nodes as specified by \a newNodeNumbers array.
412 * If \a newNodeNumbers[ i ] < 0 then the i-th node is removed,
413 * else \a newNodeNumbers[ i ] is a new id of the i-th node. The nodal connectivity
414 * array is modified accordingly. In contrast to renumberNodes(), location
415 * of merged nodes (whose new ids coincide) is changed to be at their barycenter.
416 * \param [in] newNodeNumbers - a permutation array, of length \a
417 * this->getNumberOfNodes(), in "Old to New" mode.
418 * See \ref numbering for more info on renumbering modes.
419 * \param [in] newNbOfNodes - number of nodes remaining after renumbering, which is
420 * actually one more than the maximal id in \a newNodeNumbers.
421 * \throw If the coordinates array is not set.
422 * \throw If the nodal connectivity of cells is not defined.
424 * \if ENABLE_EXAMPLES
425 * \ref cpp_mcumesh_renumberNodes "Here is a C++ example".<br>
426 * \ref py_mcumesh_renumberNodes "Here is a Python example".
429 void MEDCouplingPointSet::renumberNodesCenter(const mcIdType *newNodeNumbers, mcIdType newNbOfNodes)
431 DataArrayDouble *newCoords=DataArrayDouble::New();
432 std::vector<mcIdType> div(newNbOfNodes);
433 std::size_t spaceDim=getSpaceDimension();
434 newCoords->alloc(newNbOfNodes,spaceDim);
435 newCoords->copyStringInfoFrom(*_coords);
436 newCoords->fillWithZero();
437 mcIdType oldNbOfNodes=getNumberOfNodes();
438 double *ptToFill=newCoords->getPointer();
439 const double *oldCoordsPtr=_coords->getConstPointer();
440 for(mcIdType i=0;i<oldNbOfNodes;i++)
442 std::transform(oldCoordsPtr+i*spaceDim,oldCoordsPtr+(i+1)*spaceDim,ptToFill+newNodeNumbers[i]*spaceDim,
443 ptToFill+newNodeNumbers[i]*spaceDim,std::plus<double>());
444 div[newNodeNumbers[i]]++;
446 for(mcIdType i=0;i<newNbOfNodes;i++)
447 ptToFill=std::transform(ptToFill,ptToFill+spaceDim,ptToFill,std::bind2nd(std::multiplies<double>(),1./(double)div[i]));
448 setCoords(newCoords);
449 newCoords->decrRef();
450 renumberNodesInConn(newNodeNumbers);
454 * Computes the minimum box bounding all nodes. The edges of the box are parallel to
455 * the Cartesian coordinate axes. The bounding box is described by coordinates of its
456 * two extremum points with minimal and maximal coordinates.
457 * \param [out] bbox - array filled with coordinates of extremum points in "no
458 * interlace" mode, i.e. xMin, xMax, yMin, yMax, zMin, zMax (if in 3D). This
459 * array, of length 2 * \a this->getSpaceDimension() at least, is to be
460 * pre-allocated by the caller.
461 * \throw If the coordinates array is not set.
463 * \if ENABLE_EXAMPLES
464 * \ref cpp_mcpointset_getBoundingBox "Here is a C++ example".<br>
465 * \ref py_mcpointset_getBoundingBox "Here is a Python example".
468 void MEDCouplingPointSet::getBoundingBox(double *bbox) const
471 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::getBoundingBox : Coordinates not set !");
472 _coords->getMinMaxPerComponent(bbox);
476 * Removes "free" nodes, i.e. nodes not used to define any element.
477 * \throw If the coordinates array is not set.
478 * \throw If the elements are not defined.
480 void MEDCouplingPointSet::zipCoords()
483 DataArrayIdType *traducer=zipCoordsTraducer();
487 /*! \cond HIDDEN_ITEMS */
488 struct MEDCouplingCompAbs
490 bool operator()(double x, double y) { return std::abs(x)<std::abs(y);}
495 * Returns the carateristic dimension of \a this point set, that is a maximal
496 * absolute values of node coordinates.
497 * \throw If the coordinates array is not set.
499 double MEDCouplingPointSet::getCaracteristicDimension() const
502 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::getCaracteristicDimension : Coordinates not set !");
503 const double *coords=_coords->getConstPointer();
504 std::size_t nbOfValues=_coords->getNbOfElems();
505 return std::abs(*std::max_element(coords,coords+nbOfValues,MEDCouplingCompAbs()));
509 * This method recenter coordinates of nodes in \b this in order to be centered at the origin to benefit about the advantages of the precision to be around the box
510 * around origin of 'radius' 1.
512 * \warning this method is non const and alterates coordinates in \b this without modifying.
513 * \param [in] eps absolute epsilon. under that value of delta between max and min no scale is performed.
516 void MEDCouplingPointSet::recenterForMaxPrecision(double eps)
519 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::recenterForMaxPrecision : Coordinates not set !");
520 _coords->recenterForMaxPrecision(eps);
525 * Rotates \a this set of nodes by \a angle around either an axis (in 3D) or a point
527 * \param [in] center - coordinates either of an origin of rotation axis (in 3D) or
528 * of center of rotation (in 2D). This array is to be of size \a
529 * this->getSpaceDimension() at least.
530 * \param [in] vector - 3 components of a vector defining direction of the rotation
531 * axis in 3D. In 2D this parameter is not used.
532 * \param [in] angle - the rotation angle in radians.
533 * \throw If the coordinates array is not set.
534 * \throw If \a this->getSpaceDimension() != 2 && \a this->getSpaceDimension() != 3.
535 * \throw If \a center == NULL
536 * \throw If \a vector == NULL && \a this->getSpaceDimension() == 3.
537 * \throw If Magnitude of \a vector is zero.
539 * \if ENABLE_EXAMPLES
540 * \ref cpp_mcpointset_rotate "Here is a C++ example".<br>
541 * \ref py_mcpointset_rotate "Here is a Python example".
544 void MEDCouplingPointSet::rotate(const double *center, const double *vector, double angle)
546 std::size_t spaceDim=getSpaceDimension();
548 rotate3D(center,vector,angle);
550 rotate2D(center,angle);
552 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::rotate : invalid space dim for rotation must be 2 or 3");
553 _coords->declareAsNew();
558 * Translates \a this set of nodes.
559 * \param [in] vector - components of a translation vector. This array is to be of
560 * size \a this->getSpaceDimension() at least.
561 * \throw If the coordinates array is not set.
562 * \throw If \a vector == NULL.
564 * \if ENABLE_EXAMPLES
565 * \ref cpp_mcpointset_translate "Here is a C++ example".<br>
566 * \ref py_mcpointset_translate "Here is a Python example".
569 void MEDCouplingPointSet::translate(const double *vector)
572 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::translate : NULL input vector !");
574 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::translate : no coordinates set !");
575 double *coords=_coords->getPointer();
576 mcIdType nbNodes=getNumberOfNodes();
577 std::size_t dim=getSpaceDimension();
578 for(mcIdType i=0; i<nbNodes; i++)
579 for(std::size_t idim=0; idim<dim;idim++)
580 coords[i*dim+idim]+=vector[idim];
581 _coords->declareAsNew();
587 * Applies scaling transformation to \a this set of nodes.
588 * \param [in] point - coordinates of a scaling center. This array is to be of
589 * size \a this->getSpaceDimension() at least.
590 * \param [in] factor - a scale factor.
591 * \throw If the coordinates array is not set.
592 * \throw If \a point == NULL.
594 * \if ENABLE_EXAMPLES
595 * \ref cpp_mcpointset_scale "Here is a C++ example".<br>
596 * \ref py_mcpointset_scale "Here is a Python example".
599 void MEDCouplingPointSet::scale(const double *point, double factor)
602 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::scale : NULL input point !");
604 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::scale : no coordinates set !");
605 double *coords=_coords->getPointer();
606 mcIdType nbNodes=getNumberOfNodes();
607 std::size_t dim=getSpaceDimension();
608 for(mcIdType i=0;i<nbNodes;i++)
610 std::transform(coords+i*dim,coords+(i+1)*dim,point,coords+i*dim,std::minus<double>());
611 std::transform(coords+i*dim,coords+(i+1)*dim,coords+i*dim,std::bind2nd(std::multiplies<double>(),factor));
612 std::transform(coords+i*dim,coords+(i+1)*dim,point,coords+i*dim,std::plus<double>());
614 _coords->declareAsNew();
619 * Converts \a this set of points to an other dimension by changing number of
620 * components of point coordinates. If the dimension increases, added components
621 * are filled with \a dftValue. If the dimension decreases, last components are lost.
622 * If the new dimension is same as \a this->getSpaceDimension(), nothing is done.
623 * \param [in] newSpaceDim - the new space dimension.
624 * \param [in] dftValue - the value to assign to added components of point coordinates
625 * (if the dimension increases).
626 * \throw If the coordinates array is not set.
627 * \throw If \a newSpaceDim < 1.
629 void MEDCouplingPointSet::changeSpaceDimension(int newSpaceDim, double dftValue)
632 throw INTERP_KERNEL::Exception("changeSpaceDimension must be called on an MEDCouplingPointSet instance with coordinates set !");
634 throw INTERP_KERNEL::Exception("changeSpaceDimension must be called a newSpaceDim >=1 !");
635 int oldSpaceDim=getSpaceDimension();
636 if(newSpaceDim==oldSpaceDim)
638 DataArrayDouble *newCoords=getCoords()->changeNbOfComponents(newSpaceDim,dftValue);
639 setCoords(newCoords);
640 newCoords->decrRef();
645 * Substitutes \a this->_coords with \a other._coords provided that coordinates of
646 * the two point sets match with a specified precision, else an exception is thrown.
647 * \param [in] other - the other point set whose coordinates array will be used by
648 * \a this point set in case of their equality.
649 * \param [in] epsilon - the precision used to compare coordinates.
650 * \throw If the coordinates array of \a this is not set.
651 * \throw If the coordinates array of \a other is not set.
652 * \throw If the coordinates of \a this and \a other do not match.
654 void MEDCouplingPointSet::tryToShareSameCoords(const MEDCouplingPointSet& other, double epsilon)
656 if(_coords==other._coords)
659 throw INTERP_KERNEL::Exception("Current instance has no coords whereas other has !");
661 throw INTERP_KERNEL::Exception("Other instance has no coords whereas current has !");
662 if(!_coords->isEqualWithoutConsideringStr(*other._coords,epsilon))
663 throw INTERP_KERNEL::Exception("Coords are not the same !");
664 setCoords(other._coords);
668 * This method duplicates the nodes whose ids are in [\b nodeIdsToDuplicateBg, \b nodeIdsToDuplicateEnd) and put the result of their duplication at the end
669 * of existing node ids.
671 * \param [in] nodeIdsToDuplicateBg begin of node ids (included) to be duplicated in connectivity only
672 * \param [in] nodeIdsToDuplicateEnd end of node ids (excluded) to be duplicated in connectivity only
674 void MEDCouplingPointSet::duplicateNodesInCoords(const mcIdType *nodeIdsToDuplicateBg, const mcIdType *nodeIdsToDuplicateEnd)
677 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::duplicateNodesInCoords : no coords set !");
678 MCAuto<DataArrayDouble> newCoords=_coords->selectByTupleIdSafe(nodeIdsToDuplicateBg,nodeIdsToDuplicateEnd);
679 MCAuto<DataArrayDouble> newCoords2=DataArrayDouble::Aggregate(_coords,newCoords);
680 setCoords(newCoords2);
684 * Finds nodes located at distance lower that \a eps from a specified plane.
685 * \param [in] pt - 3 components of a point defining location of the plane.
686 * \param [in] vec - 3 components of a normal vector to the plane. Vector magnitude
687 * must be greater than 10*\a eps.
688 * \param [in] eps - maximal distance of a node from the plane at which the node is
689 * considered to lie on the plane.
690 * \param [in,out] nodes - a vector returning ids of found nodes. This vector is not
691 * cleared before filling in.
692 * \throw If the coordinates array is not set.
693 * \throw If \a pt == NULL.
694 * \throw If \a vec == NULL.
695 * \throw If the magnitude of \a vec is zero.
696 * \throw If \a this->getSpaceDimension() != 3.
698 void MEDCouplingPointSet::findNodesOnPlane(const double *pt, const double *vec, double eps, std::vector<mcIdType>& nodes) const
700 if(getSpaceDimension()!=3)
701 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::findNodesOnPlane : Invalid spacedim to be applied on this ! Must be equal to 3 !");
703 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::findNodesOnPlane : NULL point pointer specified !");
705 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::findNodesOnPlane : NULL vector pointer specified !");
706 mcIdType nbOfNodes=getNumberOfNodes();
707 double a=vec[0],b=vec[1],c=vec[2],d=-pt[0]*vec[0]-pt[1]*vec[1]-pt[2]*vec[2];
708 double deno=sqrt(a*a+b*b+c*c);
709 if(deno<std::numeric_limits<double>::min())
710 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::findNodesOnPlane : vector pointer specified has norm equal to 0. !");
711 const double *work=_coords->getConstPointer();
712 for(mcIdType i=0;i<nbOfNodes;i++)
714 if(std::abs(a*work[0]+b*work[1]+c*work[2]+d)/deno<eps)
721 * Finds nodes located at distance lower that \a eps from a specified line in 2D and 3D.
722 * \param [in] pt - components of coordinates of an initial point of the line. This
723 * array is to be of size \a this->getSpaceDimension() at least.
724 * \param [in] vec - components of a vector defining the line direction. This array
725 * is to be of size \a this->getSpaceDimension() at least. Vector magnitude
726 * must be greater than 10*\a eps.
727 * \param [in] eps - maximal distance of a node from the line at which the node is
728 * considered to lie on the line.
729 * \param [in,out] nodes - a vector returning ids of found nodes. This vector is not
730 * cleared before filling in.
731 * \throw If the coordinates array is not set.
732 * \throw If \a pt == NULL.
733 * \throw If \a vec == NULL.
734 * \throw If the magnitude of \a vec is zero.
735 * \throw If ( \a this->getSpaceDimension() != 3 && \a this->getSpaceDimension() != 2 ).
737 void MEDCouplingPointSet::findNodesOnLine(const double *pt, const double *vec, double eps, std::vector<mcIdType>& nodes) const
739 std::size_t spaceDim=getSpaceDimension();
740 if(spaceDim!=2 && spaceDim!=3)
741 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::findNodesOnLine : Invalid spacedim to be applied on this ! Must be equal to 2 or 3 !");
743 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::findNodesOnLine : NULL point pointer specified !");
745 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::findNodesOnLine : NULL vector pointer specified !");
746 mcIdType nbOfNodes=getNumberOfNodes();
748 for(std::size_t i=0;i<spaceDim;i++)
750 double deno=sqrt(den);
752 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::findNodesOnLine : Invalid given direction vector ! Norm is too small !");
753 INTERP_KERNEL::AutoPtr<double> vecn=new double[spaceDim];
754 for(std::size_t i=0;i<spaceDim;i++)
756 const double *work=_coords->getConstPointer();
759 for(mcIdType i=0;i<nbOfNodes;i++)
761 if(std::abs(vecn[0]*(work[1]-pt[1])-vecn[1]*(work[0]-pt[0]))<eps)
768 for(mcIdType i=0;i<nbOfNodes;i++)
770 double a=vecn[0]*(work[1]-pt[1])-vecn[1]*(work[0]-pt[0]);
771 double b=vecn[1]*(work[2]-pt[2])-vecn[2]*(work[1]-pt[1]);
772 double c=vecn[2]*(work[0]-pt[0])-vecn[0]*(work[2]-pt[2]);
773 if(std::sqrt(a*a+b*b+c*c)<eps)
781 * Returns a new array of node coordinates by concatenating node coordinates of two
782 * given point sets, so that (1) the number of nodes in the result array is a sum of the
783 * number of nodes of given point sets and (2) the number of component in the result array
784 * is same as that of each of given point sets. Info on components is copied from the first
785 * of the given point set. Space dimension of the given point sets must be the same.
786 * \param [in] m1 - a point set whose coordinates will be included in the result array.
787 * \param [in] m2 - another point set whose coordinates will be included in the
789 * \return DataArrayDouble * - the new instance of DataArrayDouble.
790 * The caller is to delete this result array using decrRef() as it is no more
792 * \throw If both \a m1 and \a m2 are NULL.
793 * \throw If \a m1->getSpaceDimension() != \a m2->getSpaceDimension().
795 DataArrayDouble *MEDCouplingPointSet::MergeNodesArray(const MEDCouplingPointSet *m1, const MEDCouplingPointSet *m2)
797 int spaceDim=m1->getSpaceDimension();
798 if(spaceDim!=m2->getSpaceDimension())
799 throw INTERP_KERNEL::Exception("Mismatch in SpaceDim during call of MergeNodesArray !");
800 return DataArrayDouble::Aggregate(m1->getCoords(),m2->getCoords());
803 DataArrayDouble *MEDCouplingPointSet::MergeNodesArray(const std::vector<const MEDCouplingPointSet *>& ms)
806 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::MergeNodesArray : input array must be NON EMPTY !");
807 std::vector<const MEDCouplingPointSet *>::const_iterator it=ms.begin();
808 std::vector<const DataArrayDouble *> coo(ms.size());
809 int spaceDim=(*it)->getSpaceDimension();
810 coo[0]=(*it++)->getCoords();
811 if(!coo[0]->isAllocated())
812 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::MergeNodesArray : first element in coordinates is not allocated !");
813 for(int i=1;it!=ms.end();it++,i++)
815 const DataArrayDouble *tmp=(*it)->getCoords();
818 if(tmp->isAllocated())
820 if((*it)->getSpaceDimension()==spaceDim)
823 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::MergeNodesArray : Mismatch in SpaceDim !");
826 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::MergeNodesArray : Presence of a non allocated array !");
829 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::MergeNodesArray : Empty coords detected !");
831 return DataArrayDouble::Aggregate(coo);
835 * Factory to build new instance of instanciable subclasses of MEDCouplingPointSet.
836 * This method is used during unserialization process.
838 MEDCouplingPointSet *MEDCouplingPointSet::BuildInstanceFromMeshType(MEDCouplingMeshType type)
843 return MEDCouplingUMesh::New();
844 case SINGLE_STATIC_GEO_TYPE_UNSTRUCTURED:
845 return MEDCoupling1SGTUMesh::New();
846 case SINGLE_DYNAMIC_GEO_TYPE_UNSTRUCTURED:
847 return MEDCoupling1DGTUMesh::New();
849 throw INTERP_KERNEL::Exception("Invalid type of mesh specified");
854 * First step of serialization process. Used by ParaMEDMEM and MEDCouplingCorba to transfert data between process.
856 void MEDCouplingPointSet::getTinySerializationInformation(std::vector<double>& tinyInfoD, std::vector<mcIdType>& tinyInfo, std::vector<std::string>& littleStrings) const
859 double time=getTime(it,order);
862 int spaceDim=getSpaceDimension();
863 littleStrings.resize(spaceDim+4);
864 littleStrings[0]=getName();
865 littleStrings[1]=getDescription();
866 littleStrings[2]=_coords->getName();
867 littleStrings[3]=getTimeUnit();
868 for(int i=0;i<spaceDim;i++)
869 littleStrings[i+4]=getCoords()->getInfoOnComponent(i);
871 tinyInfo.push_back(getType());
872 tinyInfo.push_back(spaceDim);
873 tinyInfo.push_back(getNumberOfNodes());
874 tinyInfo.push_back(it);
875 tinyInfo.push_back(order);
876 tinyInfoD.push_back(time);
880 littleStrings.resize(3);
881 littleStrings[0]=getName();
882 littleStrings[1]=getDescription();
883 littleStrings[2]=getTimeUnit();
885 tinyInfo.push_back(getType());
886 tinyInfo.push_back(-1);
887 tinyInfo.push_back(-1);
888 tinyInfo.push_back(it);
889 tinyInfo.push_back(order);
890 tinyInfoD.push_back(time);
895 * Third and final step of serialization process.
897 void MEDCouplingPointSet::serialize(DataArrayIdType *&a1, DataArrayDouble *&a2) const
901 a2=const_cast<DataArrayDouble *>(getCoords());
909 * Second step of serialization process.
910 * @param tinyInfo must be equal to the result given by getTinySerializationInformation method.
912 void MEDCouplingPointSet::resizeForUnserialization(const std::vector<mcIdType>& tinyInfo, DataArrayIdType *a1, DataArrayDouble *a2, std::vector<std::string>& littleStrings) const
914 if(tinyInfo[2]>=0 && tinyInfo[1]>=1)
916 a2->alloc(tinyInfo[2],tinyInfo[1]);
917 littleStrings.resize(tinyInfo[1]+4);
921 littleStrings.resize(3);
926 * Second and final unserialization process.
927 * @param tinyInfo must be equal to the result given by getTinySerializationInformation method.
929 void MEDCouplingPointSet::unserialization(const std::vector<double>& tinyInfoD, const std::vector<mcIdType>& tinyInfo, const DataArrayIdType *a1, DataArrayDouble *a2, const std::vector<std::string>& littleStrings)
931 if(tinyInfo[2]>=0 && tinyInfo[1]>=1)
934 setName(littleStrings[0]);
935 setDescription(littleStrings[1]);
936 a2->setName(littleStrings[2]);
937 setTimeUnit(littleStrings[3]);
938 for(mcIdType i=0;i<tinyInfo[1];i++)
939 getCoords()->setInfoOnComponent(i,littleStrings[i+4]);
940 setTime(tinyInfoD[0],FromIdType<int>(tinyInfo[3]),FromIdType<int>(tinyInfo[4]));
944 setName(littleStrings[0]);
945 setDescription(littleStrings[1]);
946 setTimeUnit(littleStrings[2]);
947 setTime(tinyInfoD[0],FromIdType<int>(tinyInfo[3]),FromIdType<int>(tinyInfo[4]));
951 void MEDCouplingPointSet::checkConsistencyLight() const
954 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::checkConsistencyLight : no coordinates set !");
958 * Intersect Bounding Box given 2 Bounding Boxes.
960 bool MEDCouplingPointSet::intersectsBoundingBox(const double* bb1, const double* bb2, int dim, double eps)
962 double* bbtemp = new double[2*dim];
965 for (int i=0; i< dim; i++)
967 double delta = bb1[2*i+1]-bb1[2*i];
968 if ( delta > deltamax )
973 for (int i=0; i<dim; i++)
975 bbtemp[i*2]=bb1[i*2]-deltamax*eps;
976 bbtemp[i*2+1]=bb1[i*2+1]+deltamax*eps;
979 for (int idim=0; idim < dim; idim++)
981 bool intersects = (bbtemp[idim*2]<bb2[idim*2+1])
982 && (bb2[idim*2]<bbtemp[idim*2+1]) ;
994 * Intersect 2 given Bounding Boxes.
996 bool MEDCouplingPointSet::intersectsBoundingBox(const INTERP_KERNEL::DirectedBoundingBox& bb1, const double* bb2, int dim, double eps)
998 double* bbtemp(new double[2*dim]);
999 double deltamax(0.0);
1001 for (int i=0; i< dim; i++)
1003 double delta = bb2[2*i+1]-bb2[2*i];
1004 if ( delta > deltamax )
1009 for (int i=0; i<dim; i++)
1011 bbtemp[i*2]=bb2[i*2]-deltamax*eps;
1012 bbtemp[i*2+1]=bb2[i*2+1]+deltamax*eps;
1015 bool intersects(!bb1.isDisjointWith(bbtemp));
1021 * 'This' is expected to be of spaceDim==3. Idem for 'center' and 'vect'
1023 void MEDCouplingPointSet::rotate3D(const double *center, const double *vect, double angle)
1025 double *coords(_coords->getPointer());
1026 mcIdType nbNodes(getNumberOfNodes());
1027 DataArrayDouble::Rotate3DAlg(center,vect,angle,nbNodes,coords,coords);
1031 * This method allows to give for each cell in \a trgMesh, how much it interacts with cells of \a srcMesh.
1032 * The returned array can be seen as a weighted array on the target cells of \a trgMesh input parameter.
1034 * \param [in] srcMesh - source mesh
1035 * \param [in] trgMesh - target mesh
1036 * \param [in] eps - precision of the detection
1037 * \return DataArrayIdType * - An array that gives for each cell of \a trgMesh, how many cells in \a srcMesh (regarding the precision of detection \a eps) can interacts.
1039 * \throw If \a srcMesh and \a trgMesh have not the same space dimension.
1041 DataArrayIdType *MEDCouplingPointSet::ComputeNbOfInteractionsWithSrcCells(const MEDCouplingPointSet *srcMesh, const MEDCouplingPointSet *trgMesh, double eps)
1043 if(!srcMesh || !trgMesh)
1044 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::ComputeNbOfInteractionsWithSrcCells : the input meshes must be not NULL !");
1045 MCAuto<DataArrayDouble> sbbox(srcMesh->getBoundingBoxForBBTree(eps)),tbbox(trgMesh->getBoundingBoxForBBTree(eps));
1046 return tbbox->computeNbOfInteractionsWith(sbbox,eps);
1050 * Creates a new MEDCouplingMesh containing a part of cells of \a this mesh. The new
1051 * mesh shares a coordinates array with \a this one. The cells to include to the
1052 * result mesh are specified by an array of cell ids.
1053 * \param [in] start - an array of cell ids to include to the result mesh.
1054 * \param [in] end - specifies the end of the array \a start, so that
1055 * the last value of \a start is \a end[ -1 ].
1056 * \return MEDCouplingMesh * - a new instance of MEDCouplingMesh. The caller is to
1057 * delete this mesh using decrRef() as it is no more needed.
1059 MEDCouplingMesh *MEDCouplingPointSet::buildPart(const mcIdType *start, const mcIdType *end) const
1061 return buildPartOfMySelf(start,end,true);
1065 * Creates a new MEDCouplingMesh containing a part of cells of \a this mesh. The
1066 * cells to include to the result mesh are specified by an array of cell ids.
1067 * <br> This method additionally returns a renumbering map in "Old to New" mode
1068 * which allows the caller to know the mapping between nodes in \a this and the result mesh.
1069 * \param [in] start - an array of cell ids to include to the result mesh.
1070 * \param [in] end - specifies the end of the array \a start, so that
1071 * the last value of \a start is \a end[ -1 ].
1072 * \param [out] arr - a new DataArrayIdType that is the "Old to New" renumbering
1073 * map. The caller is to delete this array using decrRef() as it is no more needed.
1074 * \return MEDCouplingMesh * - a new instance of MEDCouplingMesh. The caller is to
1075 * delete this mesh using decrRef() as it is no more needed.
1077 MEDCouplingMesh *MEDCouplingPointSet::buildPartAndReduceNodes(const mcIdType *start, const mcIdType *end, DataArrayIdType*& arr) const
1079 MCAuto<MEDCouplingPointSet> ret=buildPartOfMySelf(start,end,true);
1080 arr=ret->zipCoordsTraducer();
1085 * This method specialized the MEDCouplingMesh::buildPartRange.
1086 * This method is equivalent to MEDCouplingMesh::buildPart method except that here the cell ids are specified using slice
1087 * \a beginCellIds \a endCellIds and \a stepCellIds.
1088 * \b WARNING , there is a big difference compared to MEDCouplingMesh::buildPart method.
1089 * If the input range is equal all cells in \a this, \a this is returned !
1091 * \return a new ref to be managed by the caller. Warning this ref can be equal to \a this if input slice is exactly equal to the whole cells in the same order.
1093 * \sa MEDCouplingUMesh::buildPartOfMySelfSlice
1095 MEDCouplingMesh *MEDCouplingPointSet::buildPartRange(mcIdType beginCellIds, mcIdType endCellIds, mcIdType stepCellIds) const
1097 if(beginCellIds==0 && endCellIds==ToIdType(getNumberOfCells()) && stepCellIds==1)
1099 MEDCouplingMesh *ret(const_cast<MEDCouplingPointSet *>(this));
1105 return buildPartOfMySelfSlice(beginCellIds,endCellIds,stepCellIds,true);
1110 * This method specialized the MEDCouplingMesh::buildPartRangeAndReduceNodes
1112 * \param [out] beginOut valid only if \a arr not NULL !
1113 * \param [out] endOut valid only if \a arr not NULL !
1114 * \param [out] stepOut valid only if \a arr not NULL !
1115 * \param [out] arr correspondence old to new in node ids.
1117 * \sa MEDCouplingUMesh::buildPartOfMySelfSlice
1119 MEDCouplingMesh *MEDCouplingPointSet::buildPartRangeAndReduceNodes(mcIdType beginCellIds, mcIdType endCellIds, mcIdType stepCellIds, mcIdType& beginOut, mcIdType& endOut, mcIdType& stepOut, DataArrayIdType*& arr) const
1121 MCAuto<MEDCouplingPointSet> ret(buildPartOfMySelfSlice(beginCellIds,endCellIds,stepCellIds,true));
1122 arr=ret->zipCoordsTraducer();
1127 * 'This' is expected to be of spaceDim==2. Idem for 'center' and 'vect'
1129 void MEDCouplingPointSet::rotate2D(const double *center, double angle)
1131 double *coords(_coords->getPointer());
1132 mcIdType nbNodes(getNumberOfNodes());
1133 DataArrayDouble::Rotate2DAlg(center,angle,nbNodes,coords,coords);
1141 static const int MY_SPACEDIM=3;
1142 static const int MY_MESHDIM=2;
1143 typedef mcIdType MyConnType;
1144 static const INTERP_KERNEL::NumberingPolicy My_numPol=INTERP_KERNEL::ALL_C_MODE;
1150 * res should be an empty vector before calling this method.
1151 * This method returns all the node coordinates included in _coords which ids are in [startConn;endConn) and put it into 'res' vector.
1152 * If spaceDim==3 a projection will be done for each nodes on the middle plane containing these all nodes in [startConn;endConn).
1153 * And after each projected nodes are moved to Oxy plane in order to consider these nodes as 2D nodes.
1155 void MEDCouplingPointSet::project2DCellOnXY(const mcIdType *startConn, const mcIdType *endConn, std::vector<double>& res) const
1157 const double *coords(_coords->getConstPointer());
1158 std::size_t spaceDim(getSpaceDimension());
1159 for(const mcIdType *it=startConn;it!=endConn;it++)
1160 res.insert(res.end(),coords+spaceDim*(*it),coords+spaceDim*(*it+1));
1165 std::vector<double> cpy(res);
1166 mcIdType nbNodes=ToIdType(std::distance(startConn,endConn));
1167 INTERP_KERNEL::PlanarIntersector<DummyClsMCPS,mcIdType>::Projection(&res[0],&cpy[0],nbNodes,nbNodes,1.e-12,0./*max distance*/,-1./*min dot*/,0.,true);
1168 res.resize(2*nbNodes);
1169 for(mcIdType i=0;i<nbNodes;i++)
1172 res[2*i+1]=cpy[3*i+1];
1176 throw INTERP_KERNEL::Exception("Invalid spacedim for project2DCellOnXY !");
1180 * low level method that checks that the 2D cell is not a butterfly cell.
1182 bool MEDCouplingPointSet::isButterfly2DCell(const std::vector<double>& res, bool isQuad, double eps)
1184 INTERP_KERNEL::QuadraticPlanarPrecision prec(eps);
1186 std::size_t nbOfNodes(res.size()/2);
1187 std::vector<INTERP_KERNEL::Node *> nodes(nbOfNodes);
1188 for(std::size_t i=0;i<nbOfNodes;i++)
1190 INTERP_KERNEL::Node *tmp=new INTERP_KERNEL::Node(res[2*i],res[2*i+1]);
1193 INTERP_KERNEL::QuadraticPolygon *pol=0;
1195 pol=INTERP_KERNEL::QuadraticPolygon::BuildArcCirclePolygon(nodes);
1197 pol=INTERP_KERNEL::QuadraticPolygon::BuildLinearPolygon(nodes);
1198 bool ret(pol->isButterflyAbs());
1204 * This method compares 2 cells coming from two unstructured meshes : \a this and \a other.
1205 * This method compares 2 cells having the same id 'cellId' in \a this and \a other.
1207 bool MEDCouplingPointSet::areCellsFrom2MeshEqual(const MEDCouplingPointSet *other, mcIdType cellId, double prec) const
1209 if(getTypeOfCell(cellId)!=other->getTypeOfCell(cellId))
1211 std::vector<mcIdType> c1,c2;
1212 getNodeIdsOfCell(cellId,c1);
1213 other->getNodeIdsOfCell(cellId,c2);
1214 std::size_t sz(c1.size());
1217 for(std::size_t i=0;i<sz;i++)
1219 std::vector<double> n1,n2;
1220 getCoordinatesOfNode(c1[0],n1);
1221 other->getCoordinatesOfNode(c2[0],n2);
1222 std::transform(n1.begin(),n1.end(),n2.begin(),n1.begin(),std::minus<double>());
1223 std::transform(n1.begin(),n1.end(),n1.begin(),std::ptr_fun<double,double>(fabs));
1224 if(*std::max_element(n1.begin(),n1.end())>prec)
1231 * Substitutes node coordinates array of \a this mesh with that of \a other mesh
1232 * (i.e. \a this->_coords with \a other._coords) provided that coordinates of the two
1233 * meshes match with a specified precision, else an exception is thrown and \a this
1234 * remains unchanged. In case of success the nodal connectivity of \a this mesh
1235 * is permuted according to new order of nodes.
1236 * Contrary to tryToShareSameCoords() this method makes a deeper analysis of
1237 * coordinates (and so more expensive) than simple equality.
1238 * \param [in] other - the other mesh whose node coordinates array will be used by
1239 * \a this mesh in case of their equality.
1240 * \param [in] epsilon - the precision used to compare coordinates (using infinite norm).
1241 * \throw If the coordinates array of \a this is not set.
1242 * \throw If the coordinates array of \a other is not set.
1243 * \throw If the coordinates of \a this and \a other do not match.
1245 void MEDCouplingPointSet::tryToShareSameCoordsPermute(const MEDCouplingPointSet& other, double epsilon)
1247 const DataArrayDouble *coords=other.getCoords();
1249 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::tryToShareSameCoordsPermute : No coords specified in other !");
1251 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::tryToShareSameCoordsPermute : No coords specified in this whereas there is any in other !");
1252 mcIdType otherNbOfNodes=other.getNumberOfNodes();
1253 MCAuto<DataArrayDouble> newCoords=MergeNodesArray(&other,this);
1255 MCAuto<DataArrayDouble> oldCoords=_coords;
1256 setCoords(newCoords);
1257 bool areNodesMerged;
1258 mcIdType newNbOfNodes;
1259 MCAuto<DataArrayIdType> da=buildPermArrayForMergeNode(epsilon,otherNbOfNodes,areNodesMerged,newNbOfNodes);
1262 setCoords(oldCoords);
1263 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::tryToShareSameCoordsPermute fails : no nodes are mergeable with specified given epsilon !");
1265 mcIdType maxId=*std::max_element(da->getConstPointer(),da->getConstPointer()+otherNbOfNodes);
1266 const mcIdType *pt=std::find_if(da->getConstPointer()+otherNbOfNodes,da->getConstPointer()+da->getNbOfElems(),std::bind2nd(std::greater<mcIdType>(),maxId));
1267 if(pt!=da->getConstPointer()+da->getNbOfElems())
1269 setCoords(oldCoords);
1270 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::tryToShareSameCoordsPermute fails : some nodes in this are not in other !");
1272 setCoords(oldCoords);
1273 renumberNodesInConn(da->getConstPointer()+otherNbOfNodes);
1277 MEDCouplingPointSet *MEDCouplingPointSet::buildPartOfMySelf(const mcIdType *begin, const mcIdType *end, bool keepCoords) const
1279 MCAuto<MEDCouplingPointSet> ret=buildPartOfMySelfKeepCoords(begin,end);
1285 MEDCouplingPointSet *MEDCouplingPointSet::buildPartOfMySelfSlice(mcIdType start, mcIdType end, mcIdType step, bool keepCoords) const
1287 MCAuto<MEDCouplingPointSet> ret=buildPartOfMySelfKeepCoordsSlice(start,end,step);
1294 Creates a new MEDCouplingUMesh containing some cells of \a this mesh. The cells to
1295 copy are selected basing on specified node ids and the value of \a fullyIn
1296 parameter. If \a fullyIn ==\c true, a cell is copied if its all nodes are in the
1297 array \a begin of node ids. If \a fullyIn ==\c false, a cell is copied if any its
1298 node is in the array of node ids. The created mesh shares the node coordinates array
1300 * \param [in] begin - the array of node ids.
1301 * \param [in] end - a pointer to the (last+1)-th element of \a begin.
1302 * \param [in] fullyIn - if \c true, then cells whose all nodes are in the
1303 * array \a begin are copied, else cells whose any node is in the
1304 * array \a begin are copied.
1305 * \return MEDCouplingPointSet * - new instance of MEDCouplingUMesh. The caller is
1306 * to delete this mesh using decrRef() as it is no more needed.
1307 * \throw If the coordinates array is not set.
1308 * \throw If the nodal connectivity of cells is not defined.
1309 * \throw If any node id in \a begin is not valid.
1311 * \if ENABLE_EXAMPLES
1312 * \ref cpp_mcumesh_buildPartOfMySelfNode "Here is a C++ example".<br>
1313 * \ref py_mcumesh_buildPartOfMySelfNode "Here is a Python example".
1316 MEDCouplingPointSet *MEDCouplingPointSet::buildPartOfMySelfNode(const mcIdType *begin, const mcIdType *end, bool fullyIn) const
1318 DataArrayIdType *cellIdsKept=0;
1319 fillCellIdsToKeepFromNodeIds(begin,end,fullyIn,cellIdsKept);
1320 MCAuto<DataArrayIdType> cellIdsKept2(cellIdsKept);
1321 return buildPartOfMySelf(cellIdsKept->begin(),cellIdsKept->end(),true);
1325 * Removes duplicates of cells from \a this mesh and returns an array mapping between
1326 * new and old cell ids in "Old to New" mode. Nothing is changed in \a this mesh if no
1327 * equal cells found.
1328 * \warning Cells of the result mesh are \b not sorted by geometric type, hence,
1329 * to write this mesh to the MED file, its cells must be sorted using
1330 * sortCellsInMEDFileFrmt().
1331 * \param [in] compType - specifies a cell comparison technique. Meaning of its
1332 * valid values [0,1,2] is as follows.
1333 * - 0 : "exact". Two cells are considered equal \c iff they have exactly same nodal
1334 * connectivity and type. This is the strongest policy.
1335 * - 1 : "permuted same orientation". Two cells are considered equal \c iff they
1336 * are based on same nodes and have the same type and orientation.
1337 * - 2 : "nodal". Two cells are considered equal \c iff they
1338 * are based on same nodes and have the same type. This is the weakest
1339 * policy, it can be used by users not sensitive to cell orientation.
1340 * \param [in] startCellId - specifies the cell id at which search for equal cells
1341 * starts. By default it is 0, which means that all cells in \a this will be
1343 * \return DataArrayIdType - a new instance of DataArrayIdType, of length \a
1344 * this->getNumberOfCells() before call of this method. The caller is to
1345 * delete this array using decrRef() as it is no more needed.
1346 * \throw If the coordinates array is not set.
1347 * \throw If the nodal connectivity of cells is not defined.
1348 * \throw If the nodal connectivity includes an invalid id.
1350 * \if ENABLE_EXAMPLES
1351 * \ref cpp_mcumesh_zipConnectivityTraducer "Here is a C++ example".<br>
1352 * \ref py_mcumesh_zipConnectivityTraducer "Here is a Python example".
1355 DataArrayIdType *MEDCouplingPointSet::zipConnectivityTraducer(int compType, mcIdType startCellId)
1357 DataArrayIdType *commonCells(nullptr),*commonCellsI(nullptr);
1358 findCommonCells(compType,startCellId,commonCells,commonCellsI);
1359 MCAuto<DataArrayIdType> commonCellsTmp(commonCells),commonCellsITmp(commonCellsI);
1360 mcIdType newNbOfCells=-1;
1361 MCAuto<DataArrayIdType> ret=DataArrayIdType::ConvertIndexArrayToO2N(ToIdType(getNumberOfCells()),commonCells->begin(),commonCellsI->begin(),commonCellsI->end(),newNbOfCells);
1362 MCAuto<DataArrayIdType> ret2=ret->invertArrayO2N2N2O(newNbOfCells);
1363 MCAuto<MEDCouplingPointSet> self=buildPartOfMySelf(ret2->begin(),ret2->end(),true);
1364 shallowCopyConnectivityFrom(self);
1369 * This const method states if the nodal connectivity of this fetches all nodes in \a this.
1370 * In other words, this method looks is there are no orphan nodes in \a this.
1371 * \sa zipCoordsTraducer, getNodeIdsInUse, computeFetchedNodeIds.
1373 bool MEDCouplingPointSet::areAllNodesFetched() const
1375 checkFullyDefined();
1376 mcIdType nbNodes(getNumberOfNodes());
1377 std::vector<bool> fetchedNodes(nbNodes,false);
1378 computeNodeIdsAlg(fetchedNodes);
1379 return std::find(fetchedNodes.begin(),fetchedNodes.end(),false)==fetchedNodes.end();
1383 * Checks if \a this and \a other meshes are geometrically equivalent, else an
1384 * exception is thrown. The meshes are
1385 * considered equivalent if (1) \a this mesh contains the same nodes as the \a other
1386 * mesh (with a specified precision) and (2) \a this mesh contains the same cells as
1387 * the \a other mesh (with use of a specified cell comparison technique). The mapping
1388 * from \a other to \a this for nodes and cells is returned via out parameters.
1390 * If \a cellCor is null (or Py_None) it means that for all #i cell in \a other is equal to cell # i in \a this.
1392 * If \a nodeCor is null (or Py_None) it means that for all #i node in \a other is equal to node # i in \a this.
1394 * So null (or Py_None) returned in \a cellCor and/or \a nodeCor means identity array. This is for optimization reason to avoid building
1395 * useless arrays for some \a levOfCheck (for example 0).
1397 * \param [in] other - the mesh to compare with.
1398 * \param [in] cellCompPol - id [0-2] of cell comparison method. See meaning of
1399 * each method in description of MEDCouplingPointSet::zipConnectivityTraducer().
1400 * \param [in] prec - the precision used to compare nodes of the two meshes.
1401 * \param [out] cellCor - a cell permutation array in "Old to New" mode. The caller is
1402 * to delete this array using decrRef() as it is no more needed.
1403 * \param [out] nodeCor - a node permutation array in "Old to New" mode. The caller is
1404 * to delete this array using decrRef() as it is no more needed.
1405 * \throw If the two meshes do not match.
1407 * \if ENABLE_EXAMPLES
1408 * \ref cpp_mcumesh_checkDeepEquivalWith "Here is a C++ example".<br>
1409 * \ref py_mcumesh_checkDeepEquivalWith "Here is a Python example".
1412 void MEDCouplingPointSet::checkDeepEquivalWith(const MEDCouplingMesh *other, int cellCompPol, double prec,
1413 DataArrayIdType *&cellCor, DataArrayIdType *&nodeCor) const
1416 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::checkDeepEquivalWith : input is null !");
1417 const MEDCouplingPointSet *otherC=dynamic_cast<const MEDCouplingPointSet *>(other);
1419 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::checkDeepEquivalWith : other is not a PointSet mesh !");
1420 MCAuto<MEDCouplingPointSet> m=dynamic_cast<MEDCouplingPointSet *>(mergeMyselfWith(otherC));
1421 bool areNodesMerged;
1422 mcIdType newNbOfNodes;
1423 mcIdType oldNbOfNodes=getNumberOfNodes();
1424 MCAuto<DataArrayIdType> da=m->buildPermArrayForMergeNode(prec,oldNbOfNodes,areNodesMerged,newNbOfNodes);
1426 if(!areNodesMerged && oldNbOfNodes != 0)
1427 throw INTERP_KERNEL::Exception("checkDeepEquivalWith : Nodes are incompatible ! ");
1428 const mcIdType *pt=std::find_if(da->getConstPointer()+oldNbOfNodes,da->getConstPointer()+da->getNbOfElems(),std::bind2nd(std::greater<mcIdType>(),oldNbOfNodes-1));
1429 if(pt!=da->getConstPointer()+da->getNbOfElems())
1430 throw INTERP_KERNEL::Exception("checkDeepEquivalWith : some nodes in other are not in this !");
1431 m->renumberNodes(da->getConstPointer(),newNbOfNodes);
1433 MCAuto<DataArrayIdType> nodeCor2=da->subArray(oldNbOfNodes);
1434 da=m->mergeNodes(prec,areNodesMerged,newNbOfNodes);
1436 da=m->zipConnectivityTraducer(cellCompPol);
1437 mcIdType nbCells=ToIdType(getNumberOfCells());
1438 if (nbCells != ToIdType(other->getNumberOfCells()))
1439 throw INTERP_KERNEL::Exception("checkDeepEquivalWith : some cells in other are not in this !");
1440 mcIdType dan(da->getNumberOfTuples());
1443 MCAuto<DataArrayIdType> da1(DataArrayIdType::New()),da2(DataArrayIdType::New());
1444 da1->alloc(dan/2,1); da2->alloc(dan/2,1);
1445 std::copy(da->getConstPointer(), da->getConstPointer()+dan/2, da1->getPointer());
1446 std::copy(da->getConstPointer()+dan/2, da->getConstPointer()+dan, da2->getPointer());
1447 da1->sort(); da2->sort();
1448 if (!da1->isEqualWithoutConsideringStr(*da2))
1449 throw INTERP_KERNEL::Exception("checkDeepEquivalWith : some cells in other are not in this !");
1451 MCAuto<DataArrayIdType> cellCor2=da->selectByTupleIdSafeSlice(nbCells,ToIdType(da->getNbOfElems()),1);
1452 nodeCor=nodeCor2->isIota(nodeCor2->getNumberOfTuples())?0:nodeCor2.retn();
1453 cellCor=cellCor2->isIota(cellCor2->getNumberOfTuples())?0:cellCor2.retn();
1457 * Checks if \a this and \a other meshes are geometrically equivalent, else an
1458 * exception is thrown. The meshes are considered equivalent if (1) they share one
1459 * node coordinates array and (2) they contain the same cells (with use of a specified
1460 * cell comparison technique). The mapping from cells of the \a other to ones of \a this
1461 * is returned via an out parameter.
1463 * If \a cellCor is null (or Py_None) it means that for all #i cell in \a other is equal to cell # i in \a this.
1465 * \param [in] other - the mesh to compare with.
1466 * \param [in] cellCompPol - id [0-2] of cell comparison method. See the meaning of
1467 * each method in description of MEDCouplingPointSet::zipConnectivityTraducer().
1468 * \param [in] prec - a not used parameter.
1469 * \param [out] cellCor - the permutation array in "Old to New" mode. The caller is
1470 * to delete this array using decrRef() as it is no more needed.
1471 * \throw If the two meshes do not match.
1473 * \if ENABLE_EXAMPLES
1474 * \ref cpp_mcumesh_checkDeepEquivalWith "Here is a C++ example".<br>
1475 * \ref py_mcumesh_checkDeepEquivalWith "Here is a Python example".
1478 void MEDCouplingPointSet::checkDeepEquivalOnSameNodesWith(const MEDCouplingMesh *other, int cellCompPol, double prec,
1479 DataArrayIdType *&cellCor) const
1482 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::checkDeepEquivalOnSameNodesWith : input is null !");
1483 const MEDCouplingPointSet *otherC=dynamic_cast<const MEDCouplingPointSet *>(other);
1485 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::checkDeepEquivalOnSameNodesWith : other is not a PointSet mesh !");
1486 if(_coords!=otherC->_coords)
1487 throw INTERP_KERNEL::Exception("checkDeepEquivalOnSameNodesWith : meshes do not share the same coordinates ! Use tryToShareSameCoordinates or call checkDeepEquivalWith !");
1488 MCAuto<MEDCouplingPointSet> m=mergeMyselfWithOnSameCoords(otherC);
1489 MCAuto<DataArrayIdType> da=m->zipConnectivityTraducer(cellCompPol);
1490 mcIdType maxId=*std::max_element(da->getConstPointer(),da->getConstPointer()+getNumberOfCells());
1491 const mcIdType *pt=std::find_if(da->getConstPointer()+getNumberOfCells(),da->getConstPointer()+da->getNbOfElems(),std::bind2nd(std::greater<mcIdType>(),maxId));
1492 if(pt!=da->getConstPointer()+da->getNbOfElems())
1494 throw INTERP_KERNEL::Exception("checkDeepEquivalOnSameNodesWith : some cells in other are not in this !");
1496 MCAuto<DataArrayIdType> cellCor2=da->selectByTupleIdSafeSlice(ToIdType(getNumberOfCells()),ToIdType(da->getNbOfElems()),1);
1497 cellCor=cellCor2->isIota(cellCor2->getNumberOfTuples())?0:cellCor2.retn();
1500 void MEDCouplingPointSet::checkFastEquivalWith(const MEDCouplingMesh *other, double prec) const
1502 MEDCouplingMesh::checkFastEquivalWith(other,prec);
1503 //other not null checked by the line before
1504 const MEDCouplingPointSet *otherC=dynamic_cast<const MEDCouplingPointSet *>(other);
1506 throw INTERP_KERNEL::Exception("MEDCouplingPointSet::checkFastEquivalWith : fails because other is not a pointset mesh !");
1507 mcIdType nbOfCells=ToIdType(getNumberOfCells());
1511 status&=areCellsFrom2MeshEqual(otherC,0,prec);
1512 status&=areCellsFrom2MeshEqual(otherC,nbOfCells/2,prec);
1513 status&=areCellsFrom2MeshEqual(otherC,nbOfCells-1,prec);
1515 throw INTERP_KERNEL::Exception("checkFastEquivalWith : Two meshes are not equal because on 3 test cells some difference have been detected !");
1519 * Finds cells whose all or some nodes are in a given array of node ids.
1520 * \param [in] begin - the array of node ids.
1521 * \param [in] end - a pointer to the (last+1)-th element of \a begin.
1522 * \param [in] fullyIn - if \c true, then cells whose all nodes are in the
1523 * array \a begin are returned only, else cells whose any node is in the
1524 * array \a begin are returned.
1525 * \return DataArrayIdType * - a new instance of DataArrayIdType holding ids of found
1526 * cells. The caller is to delete this array using decrRef() as it is no more
1528 * \throw If the coordinates array is not set.
1529 * \throw If the nodal connectivity of cells is not defined.
1530 * \throw If any cell id in \a begin is not valid.
1532 * \sa MEDCouplingPointSet::getCellIdsFullyIncludedInNodeIds
1534 * \if ENABLE_EXAMPLES
1535 * \ref cpp_mcumesh_getCellIdsLyingOnNodes "Here is a C++ example".<br>
1536 * \ref py_mcumesh_getCellIdsLyingOnNodes "Here is a Python example".
1539 DataArrayIdType *MEDCouplingPointSet::getCellIdsLyingOnNodes(const mcIdType *begin, const mcIdType *end, bool fullyIn) const
1541 DataArrayIdType *cellIdsKept=0;
1542 fillCellIdsToKeepFromNodeIds(begin,end,fullyIn,cellIdsKept);
1543 cellIdsKept->setName(getName());
1548 * Finds cells whose all nodes are in a given array of node ids.
1549 * This method is a specialization of MEDCouplingPointSet::getCellIdsLyingOnNodes (true
1550 * as last input argument).
1551 * \param [in] partBg - the array of node ids.
1552 * \param [in] partEnd - a pointer to a (last+1)-th element of \a partBg.
1553 * \return DataArrayIdType * - a new instance of DataArrayIdType holding ids of found
1554 * cells. The caller is to delete this array using decrRef() as it is no
1556 * \throw If the coordinates array is not set.
1557 * \throw If the nodal connectivity of cells is not defined.
1558 * \throw If any cell id in \a partBg is not valid.
1560 * \sa MEDCouplingPointSet::getCellIdsLyingOnNodes
1562 * \if ENABLE_EXAMPLES
1563 * \ref cpp_mcumesh_getCellIdsFullyIncludedInNodeIds "Here is a C++ example".<br>
1564 * \ref py_mcumesh_getCellIdsFullyIncludedInNodeIds "Here is a Python example".
1567 DataArrayIdType *MEDCouplingPointSet::getCellIdsFullyIncludedInNodeIds(const mcIdType *partBg, const mcIdType *partEnd) const
1569 return getCellIdsLyingOnNodes(partBg,partEnd,true);
1573 * Removes unused nodes (the node coordinates array is shorten) and returns an array
1574 * mapping between new and old node ids in "Old to New" mode. -1 values in the returned
1575 * array mean that the corresponding old node is no more used.
1576 * \return DataArrayIdType * - a new instance of DataArrayIdType of length \a
1577 * this->getNumberOfNodes() before call of this method. The caller is to
1578 * delete this array using decrRef() as it is no more needed.
1579 * \throw If the coordinates array is not set.
1580 * \throw If the nodal connectivity of cells is not defined.
1581 * \throw If the nodal connectivity includes an invalid id.
1583 * \if ENABLE_EXAMPLES
1584 * \ref cpp_mcumesh_zipCoordsTraducer "Here is a C++ example".<br>
1585 * \ref py_mcumesh_zipCoordsTraducer "Here is a Python example".
1588 DataArrayIdType *MEDCouplingPointSet::zipCoordsTraducer()
1590 mcIdType newNbOfNodes=-1;
1591 MCAuto<DataArrayIdType> traducer=getNodeIdsInUse(newNbOfNodes);
1592 renumberNodes(traducer->getConstPointer(),newNbOfNodes);
1593 return traducer.retn();
1597 * Merges nodes equal within \a precision and returns an array describing the
1598 * permutation used to remove duplicate nodes.
1599 * \param [in] precision - minimal absolute distance between two nodes at which they are
1600 * considered not coincident.
1601 * \param [out] areNodesMerged - is set to \c true if any coincident nodes removed.
1602 * \param [out] newNbOfNodes - number of nodes remaining after the removal.
1603 * \return DataArrayIdType * - the permutation array in "Old to New" mode. For more
1604 * info on "Old to New" mode see \ref numbering. The caller
1605 * is to delete this array using decrRef() as it is no more needed.
1606 * \throw If the coordinates array is not set.
1607 * \throw If the nodal connectivity of cells is not defined.
1609 * \if ENABLE_EXAMPLES
1610 * \ref cpp_mcumesh_mergeNodes "Here is a C++ example".<br>
1611 * \ref py_mcumesh_mergeNodes "Here is a Python example".
1614 DataArrayIdType *MEDCouplingPointSet::mergeNodes(double precision, bool& areNodesMerged, mcIdType& newNbOfNodes)
1616 MCAuto<DataArrayIdType> ret=buildPermArrayForMergeNode(precision,-1,areNodesMerged,newNbOfNodes);
1618 renumberNodes(ret->begin(),newNbOfNodes);
1623 * Merges nodes equal within \a precision and returns an array describing the
1624 * permutation used to remove duplicate nodes. In contrast to mergeNodes(), location
1625 * of merged nodes is changed to be at their barycenter.
1626 * \param [in] precision - minimal absolute distance between two nodes at which they are
1627 * considered not coincident.
1628 * \param [out] areNodesMerged - is set to \c true if any coincident nodes removed.
1629 * \param [out] newNbOfNodes - number of nodes remaining after the removal.
1630 * \return DataArrayIdType * - the permutation array in "Old to New" mode. For more
1631 * info on "Old to New" mode see \ref numbering. The caller
1632 * is to delete this array using decrRef() as it is no more needed.
1633 * \throw If the coordinates array is not set.
1634 * \throw If the nodal connectivity of cells is not defined.
1636 * \if ENABLE_EXAMPLES
1637 * \ref cpp_mcumesh_mergeNodes "Here is a C++ example".<br>
1638 * \ref py_mcumesh_mergeNodes "Here is a Python example".
1641 DataArrayIdType *MEDCouplingPointSet::mergeNodesCenter(double precision, bool& areNodesMerged, mcIdType& newNbOfNodes)
1643 DataArrayIdType *ret=buildPermArrayForMergeNode(precision,-1,areNodesMerged,newNbOfNodes);
1645 renumberNodesCenter(ret->getConstPointer(),newNbOfNodes);