Open Inventor Release 2023.2.3
 
Loading...
Searching...
No Matches
MiPointProbeIjk.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-2019 BY FEI S.A.S, ***
17 *** BORDEAUX, FRANCE ***
18 *** ALL RIGHTS RESERVED ***
19**=======================================================================*/
20/*=======================================================================
21** Author : VSG (MMM YYYY)
22**=======================================================================*/
23#ifndef _MiPointProbeIjk_h
24#define _MiPointProbeIjk_h
25
27
30
31class MiCellFilterIjk;
35
48{
49public:
50 virtual ~MiPointProbeIjk() {}
51
53
59 static MiPointProbeIjk* getNewInstance(const MiVolumeMeshRegular& mesh, bool parallel = true);
60 static MiPointProbeIjk* getNewInstance(const MiVolumeMeshRectilinear& mesh, bool parallel = true);
61 static MiPointProbeIjk* getNewInstance(const MiVolumeMeshCurvilinear& mesh, bool parallel = true);
63
67 virtual bool setLocation(const MbVec3d& point, const MiCellFilterIjk* cellFilter=NULL) = 0;
68
93 virtual bool moveLocation(const MbVec3d& point, size_t cellIdI, size_t cellIdJ, size_t cellIdK,
94 const MiCellFilterIjk* cellFilter=NULL) = 0;
95
99 virtual bool isFound() const = 0;
100
104 template <typename _T>
105 _T getValue(const MiDataSetIjk<_T>& dataset) const;
106
114 virtual void getCellId(size_t &i, size_t &j, size_t& k) const = 0;
115
116private: protected:
121 virtual void getWeight(std::vector<double>& weight) const = 0;
122};
123
124//-----------------------------------------------------------------------------
125template <typename _T>
126inline _T MiPointProbeIjk::getValue(const MiDataSetIjk<_T>& dataset) const
127{
128 _T val(0);
129 if (isFound())
130 {
131 size_t icell,jcell,kcell;
132 getCellId(icell,jcell,kcell);
134 {
135 std::vector<double> weight(8);
136 getWeight(weight);
137 val += weight[0]*dataset.get(icell,jcell,kcell);
138 val += weight[1]*dataset.get(icell+1,jcell,kcell);
139 val += weight[2]*dataset.get(icell+1,jcell+1,kcell);
140 val += weight[3]*dataset.get(icell,jcell+1,kcell);
141 val += weight[4]*dataset.get(icell,jcell,kcell+1);
142 val += weight[5]*dataset.get(icell+1,jcell,kcell+1);
143 val += weight[6]*dataset.get(icell+1,jcell+1,kcell+1);
144 val += weight[7]*dataset.get(icell,jcell+1,kcell+1);
145 }
146 else
147 val = dataset.get(icell,jcell,kcell);
148 }
149 return val;
150}
151
152
153
154#endif
155
156
#define MESHVIZXLM_EXTR_API
valueRef getValue() const
Definition SoSubField.h:234
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Interface class defining...
virtual DataBinding getBinding() const =0
Returns the data binding of this data set.
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Generic data set interfa...
virtual _T get(size_t i, size_t j, size_t k) const =0
Returns the ijk-th value of this set.
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Probe at a given locatio...
virtual void getCellId(size_t &i, size_t &j, size_t &k) const =0
Gets the ids of the cell of the mesh containing the current location of the probe.
virtual void getWeight(std::vector< double > &weight) const =0
Gets the weights at the current probe location.
static MiPointProbeIjk * getNewInstance(const MiVolumeMeshCurvilinear &mesh, bool parallel=true)
virtual bool setLocation(const MbVec3d &point, const MiCellFilterIjk *cellFilter=NULL)=0
Sets the current probe location at the given point and looks for the cell containing this point.
static MiPointProbeIjk * getNewInstance(const MiVolumeMeshRectilinear &mesh, bool parallel=true)
virtual ~MiPointProbeIjk()
virtual bool moveLocation(const MbVec3d &point, size_t cellIdI, size_t cellIdJ, size_t cellIdK, const MiCellFilterIjk *cellFilter=NULL)=0
Moves the current probe location at the given point.
_T getValue(const MiDataSetIjk< _T > &dataset) const
Gets the value at the current probe location in the given dataset.
virtual bool isFound() const =0
static MiPointProbeIjk * getNewInstance(const MiVolumeMeshRegular &mesh, bool parallel=true)
Factory method returning a new instance of a class implementing this abstract class.
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Curvilinear volume mesh ...
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Rectilinear volume mesh ...
<a href="IconLegend.html"><img src="extMV.gif" alt="MeshViz" border="0"></a> Regular volume mesh abst...