Open Inventor Release 2024.2.1
 
Loading...
Searching...
No Matches
MxTriangle6CellExtract.h
1/*=======================================================================
2 *** THE CONTENT OF THIS WORK IS PROPRIETARY TO FEI S.A.S, (FEI S.A.S.), ***
3 *** AND IS DISTRIBUTED UNDER A LICENSE AGREEMENT. ***
4 *** ***
5 *** REPRODUCTION, DISCLOSURE, OR USE, IN WHOLE OR IN PART, OTHER THAN AS ***
6 *** SPECIFIED IN THE LICENSE ARE NOT TO BE UNDERTAKEN EXCEPT WITH PRIOR ***
7 *** WRITTEN AUTHORIZATION OF FEI S.A.S. ***
8 *** ***
9 *** RESTRICTED RIGHTS LEGEND ***
10 *** USE, DUPLICATION, OR DISCLOSURE BY THE GOVERNMENT OF THE CONTENT OF THIS ***
11 *** WORK OR RELATED DOCUMENTATION IS SUBJECT TO RESTRICTIONS AS SET FORTH IN ***
12 *** SUBPARAGRAPH (C)(1) OF THE COMMERCIAL COMPUTER SOFTWARE RESTRICTED RIGHT ***
13 *** CLAUSE AT FAR 52.227-19 OR SUBPARAGRAPH (C)(1)(II) OF THE RIGHTS IN ***
14 *** TECHNICAL DATA AND COMPUTER SOFTWARE CLAUSE AT DFARS 52.227-7013. ***
15 *** ***
16 *** COPYRIGHT (C) 1996-2017 BY FEI S.A.S, ***
17 *** BORDEAUX, FRANCE ***
18 *** ALL RIGHTS RESERVED ***
19**=======================================================================*/
20/*=======================================================================
21** Author : VSG (MMM YYYY)
22**=======================================================================*/
23#ifndef _MxTriangle6CellExtract_h
24#define _MxTriangle6CellExtract_h
25
26#include <MeshVizXLM/mesh/geometry/MiGeometryI.h>
27
28#include <Inventor/STL/vector>
29#include <Inventor/STL/utility>
30
31class MiSurfaceCell;
32
57class MESHVIZXLM_EXTR_API MxTriangle6CellExtract
58{
59public:
60
74 static MbVec3d getIsoParametricCoord(const MiGeometryI& meshGeometry, const MiSurfaceCell* triangleCell, const MbVec3d &point);
75
81 static MbVec3d getIsoParametricCoord(size_t nodeIndex)
82 {
83 return MbVec3d(s_nodesIsoParametricCoords[nodeIndex]);
84 }
85
105 static void getWeight(const MiGeometryI& meshGeometry, const MiSurfaceCell* triangleCell, const MbVec3d &point, std::vector<double>& weights)
106 {
107 getWeight(getIsoParametricCoord(meshGeometry,triangleCell,point),weights);
108 }
109
122 static void getWeight(const MbVec3d &ipcoord, std::vector<double>& weights);
123
135 static void getDerivs(const MbVec3d &ipcoord, std::vector<double>& derivs);
136
150 static bool isPointInsideCell(const MiGeometryI& meshGeometry, const MiSurfaceCell* triangleCell, const MbVec3d &point, std::vector<double>& weights)
151 {
152 MbVec3d pcoord = getIsoParametricCoord(meshGeometry,triangleCell,point);
153 bool inside = ( -1.E-5 < pcoord[0] && pcoord[0] < 1 + 1.E-5 &&
154 -1.E-5 < pcoord[1] && pcoord[1] < 1 + 1.E-5 &&
155 -1.E-5 < (pcoord[0]+pcoord[1]) && (pcoord[0]+pcoord[1]) < 1 + 1.E-5);
156 if (inside)
157 getWeight(pcoord,weights);
158 return inside;
159 }
160
168 static size_t getSubTriangleNodesIndex (std::vector<size_t>& triangleNodeIds)
169 {
170 triangleNodeIds.assign(s_subTriangleNodeIds,s_subTriangleNodeIds+(4*3));
171 return 4;
172 }
173
174 friend std::ostream& operator << (std::ostream& s, const MxTriangle6CellExtract& cell);
175
176private:
177 static double s_nodesIsoParametricCoords[6][3];
178 static size_t s_subTriangleNodeIds[4*3];
179};
180
181//-----------------------------------------------------------------------------
182inline void
183MxTriangle6CellExtract::getWeight(const MbVec3d &ipcoord, std::vector<double>& weight)
184{
185 double r = ipcoord[0];
186 double s = ipcoord[1];
187 weight[0] = (1-r-s) * ( 2*(1-r-s)-1 );
188 weight[1] = r * (2*r-1);
189 weight[2] = s * (2*s-1);
190 weight[3] = 4 * r * (1-r-s);
191 weight[4] = 4 * r * s;
192 weight[5] = 4 * s * (1-r-s);
193}
194
195
196#endif
197
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Generic geometry interfa...
Definition MiGeometryI.h:39
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Defines an abstract cell...
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Utility class for quadra...
static void getDerivs(const MbVec3d &ipcoord, std::vector< double > &derivs)
Gets the value of the derivatives of the shape functions (aka weights) at the point given by its iso ...
static size_t getSubTriangleNodesIndex(std::vector< size_t > &triangleNodeIds)
Gets a triangulation of a quadratic triangle of 6 nodes.
static MbVec3d getIsoParametricCoord(size_t nodeIndex)
Returns the iso parametric coordinates of one of the 6 nodes of a triangle cell.
static MbVec3d getIsoParametricCoord(const MiGeometryI &meshGeometry, const MiSurfaceCell *triangleCell, const MbVec3d &point)
Computes the iso parametric coordinates of the given point in the given cell with the given geometry.
static bool isPointInsideCell(const MiGeometryI &meshGeometry, const MiSurfaceCell *triangleCell, const MbVec3d &point, std::vector< double > &weights)
Checks if a point is inside or outside a triangle cell.
static void getWeight(const MiGeometryI &meshGeometry, const MiSurfaceCell *triangleCell, const MbVec3d &point, std::vector< double > &weights)
Gets the weights in the given cell of the given point.
MbVec3< double > MbVec3d
Vector of 3 double coordinates.
Definition MbVec3.h:532