Open Inventor
Release 2023.2.3
Loading...
Searching...
No Matches
MxTriangle6CellExtract.h
Go to the documentation of this file.
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
31
class
MiSurfaceCell
;
32
57
class
MESHVIZXLM_EXTR_API
MxTriangle6CellExtract
58
{
59
public
:
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
176
private
:
177
static
double
s_nodesIsoParametricCoords[6][3];
178
static
size_t
s_subTriangleNodeIds[4*3];
179
};
180
181
//-----------------------------------------------------------------------------
182
inline
void
183
MxTriangle6CellExtract::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
MiGeometryI.h
getWeight
static void getWeight(const MiGeometryI &meshGeometry, const MiVolumeCell *hexahedronCell, const MbVec3d &point, std::vector< double > &weights)
Gets the weights of a point in a cell.
Definition
MxHexahedron21CellExtract.h:70
getIsoParametricCoord
static MbVec3d getIsoParametricCoord(size_t nodeIndex)
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Utility class for quadra...
MESHVIZXLM_EXTR_API
#define MESHVIZXLM_EXTR_API
Definition
MxMeshVizXLM.h:48
operator<<
std::ostream & operator<<(std::ostream &os, const SoPathList &pathList)
Writes the list to the specified output stream.
MbVec3< double >
MiGeometryI
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Generic geometry interfa...
Definition
MiGeometryI.h:39
MiSurfaceCell
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Defines an abstract cell...
Definition
MiSurfaceCell.h:70
MxTriangle6CellExtract
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Utility class for quadra...
Definition
MxTriangle6CellExtract.h:58
MxTriangle6CellExtract::getDerivs
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 ...
MxTriangle6CellExtract::getSubTriangleNodesIndex
static size_t getSubTriangleNodesIndex(std::vector< size_t > &triangleNodeIds)
Gets a triangulation of a quadratic triangle of 6 nodes.
Definition
MxTriangle6CellExtract.h:168
MxTriangle6CellExtract::getIsoParametricCoord
static MbVec3d getIsoParametricCoord(size_t nodeIndex)
Returns the iso parametric coordinates of one of the 6 nodes of a triangle cell.
Definition
MxTriangle6CellExtract.h:81
MxTriangle6CellExtract::getIsoParametricCoord
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.
MxTriangle6CellExtract::isPointInsideCell
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.
Definition
MxTriangle6CellExtract.h:150
MxTriangle6CellExtract::getWeight
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.
Definition
MxTriangle6CellExtract.h:105
MbVec3d
MbVec3< double > MbVec3d
Vector of 3 double coordinates.
Definition
MbVec3.h:532
MeshVizXLM
extractors
MxTriangle6CellExtract.h
Generated by
1.9.8