MiPointProbeIjk.h

/*=======================================================================
 *** THE CONTENT OF THIS WORK IS PROPRIETARY TO FEI S.A.S, (FEI S.A.S.),            ***
 ***              AND IS DISTRIBUTED UNDER A LICENSE AGREEMENT.                     ***
 ***                                                                                ***
 ***  REPRODUCTION, DISCLOSURE,  OR USE,  IN WHOLE OR IN PART,  OTHER THAN AS       ***
 ***  SPECIFIED  IN THE LICENSE ARE  NOT TO BE  UNDERTAKEN  EXCEPT WITH PRIOR       ***
 ***  WRITTEN AUTHORIZATION OF FEI S.A.S.                                           ***
 ***                                                                                ***
 ***                        RESTRICTED RIGHTS LEGEND                                ***
 ***  USE, DUPLICATION, OR DISCLOSURE BY THE GOVERNMENT OF THE CONTENT OF THIS      ***
 ***  WORK OR RELATED DOCUMENTATION IS SUBJECT TO RESTRICTIONS AS SET FORTH IN      ***
 ***  SUBPARAGRAPH (C)(1) OF THE COMMERCIAL COMPUTER SOFTWARE RESTRICTED RIGHT      ***
 ***  CLAUSE  AT FAR 52.227-19  OR SUBPARAGRAPH  (C)(1)(II)  OF  THE RIGHTS IN      ***
 ***  TECHNICAL DATA AND COMPUTER SOFTWARE CLAUSE AT DFARS 52.227-7013.             ***
 ***                                                                                ***
 ***                   COPYRIGHT (C) 1996-2019 BY FEI S.A.S,                        ***
 ***                        BORDEAUX, FRANCE                                        ***
 ***                      ALL RIGHTS RESERVED                                       ***
**=======================================================================*/
/*=======================================================================
** Author      : VSG (MMM YYYY)
**=======================================================================*/
#ifndef  _MiPointProbeIjk_h
#define  _MiPointProbeIjk_h
 
#include <MeshVizXLM/extractors/MiPointProbeUnstructured.h>
 
#include <MeshVizXLM/mesh/data/MiDataSetIjk.h>
#include <MeshVizXLM/mesh/cell/MiCell.h>
 
class MiCellFilterIjk;
class MiVolumeMeshRegular;
class MiVolumeMeshRectilinear;
class MiVolumeMeshCurvilinear;
 
class MESHVIZXLM_EXTR_API MiPointProbeIjk 
{
public:
  virtual ~MiPointProbeIjk() {}
 
 
  static MiPointProbeIjk* getNewInstance(const MiVolumeMeshRegular& mesh, bool parallel = true);
  static MiPointProbeIjk* getNewInstance(const MiVolumeMeshRectilinear& mesh, bool parallel = true);
  static MiPointProbeIjk* getNewInstance(const MiVolumeMeshCurvilinear& mesh, bool parallel = true);
 
  virtual bool setLocation(const MbVec3d& point, const MiCellFilterIjk* cellFilter=NULL) = 0;
 
  virtual bool moveLocation(const MbVec3d& point, size_t cellIdI, size_t cellIdJ, size_t cellIdK, 
                            const MiCellFilterIjk* cellFilter=NULL) = 0;
 
  virtual bool isFound() const = 0;
 
  template <typename _T>
  _T getValue(const MiDataSetIjk<_T>& dataset) const;
 
  virtual void getCellId(size_t &i, size_t &j, size_t& k) const = 0;
 
private: protected:
  virtual void getWeight(std::vector<double>& weight) const = 0;
};
 
//-----------------------------------------------------------------------------
template <typename _T>
inline _T MiPointProbeIjk::getValue(const MiDataSetIjk<_T>& dataset) const
{
  _T val(0);
  if (isFound())
  {
    size_t icell,jcell,kcell;
    getCellId(icell,jcell,kcell);
    if(dataset.getBinding()==MiDataSetIjk<_T>::PER_NODE)
    {
      std::vector<double> weight(8);
      getWeight(weight);
      val += weight[0]*dataset.get(icell,jcell,kcell);
      val += weight[1]*dataset.get(icell+1,jcell,kcell);
      val += weight[2]*dataset.get(icell+1,jcell+1,kcell);
      val += weight[3]*dataset.get(icell,jcell+1,kcell);
      val += weight[4]*dataset.get(icell,jcell,kcell+1);
      val += weight[5]*dataset.get(icell+1,jcell,kcell+1);
      val += weight[6]*dataset.get(icell+1,jcell+1,kcell+1);
      val += weight[7]*dataset.get(icell,jcell+1,kcell+1);
    }
    else 
      val = dataset.get(icell,jcell,kcell);
  }
  return val;
}
 
 
 
#endif