From 5f5464d34bb79b13070b1507e6912f11313e1cb6 Mon Sep 17 00:00:00 2001
From: eap <eap@opencascade.com>
Date: Mon, 18 Jan 2010 10:16:31 +0000
Subject: [PATCH] 0020634: [CEA 381] MEDMEM:REMAPPER : intersection error +   
 unsigned getNumberOfNodesConstituentTheSon2(unsigned sonId, const int
 *nodalConn, int lgth) const;

---
 src/INTERP_KERNEL/CellModel.cxx | 29 +++++++++++++++++++++++++++++
 src/INTERP_KERNEL/CellModel.hxx |  1 +
 2 files changed, 30 insertions(+)

diff --git a/src/INTERP_KERNEL/CellModel.cxx b/src/INTERP_KERNEL/CellModel.cxx
index 418936050..77a648582 100644
--- a/src/INTERP_KERNEL/CellModel.cxx
+++ b/src/INTERP_KERNEL/CellModel.cxx
@@ -313,4 +313,33 @@ namespace INTERP_KERNEL
           }
       }
   }
+
+  //================================================================================
+  /*!
+   * \brief Return number of nodes in sonId-th son of a Dynamic() cell
+   */
+  //================================================================================
+
+  unsigned CellModel::getNumberOfNodesConstituentTheSon2(unsigned sonId, const int *nodalConn, int lgth) const
+  {
+    if(!isDynamic())
+      return getNumberOfNodesConstituentTheSon(sonId);
+
+    if(_dim==2)//polygon
+      {
+        return 2;
+      }
+    else
+      {//polyedron
+        const int *where=nodalConn;
+        for(int i=0;i<sonId;i++)
+          {
+            where=std::find(where,nodalConn+lgth,-1);
+            where++;
+          }
+        const int *where2=std::find(where,nodalConn+lgth,-1);
+        return where2-where;
+      }
+  }
+
 }
diff --git a/src/INTERP_KERNEL/CellModel.hxx b/src/INTERP_KERNEL/CellModel.hxx
index 5ca9d343b..2acdc6ef6 100644
--- a/src/INTERP_KERNEL/CellModel.hxx
+++ b/src/INTERP_KERNEL/CellModel.hxx
@@ -50,6 +50,7 @@ namespace INTERP_KERNEL
     unsigned getNumberOfSons() const { return _nb_of_sons; }
     unsigned getNumberOfSons2(const int *conn, int lgth) const;
     unsigned getNumberOfNodesConstituentTheSon(unsigned sonId) const { return _nb_of_sons_con[sonId]; }
+    unsigned getNumberOfNodesConstituentTheSon2(unsigned sonId, const int *nodalConn, int lgth) const;
     NormalizedCellType getSonType(unsigned sonId) const { return _sons_type[sonId]; }
     NormalizedCellType getSonType2(unsigned sonId) const;
     unsigned fillSonCellNodalConnectivity(int sonId, const int *nodalConn, int *sonNodalConn) const;
-- 
2.39.2