BarycentricGeomShapeComputer.h

// -*- tab-width: 2; indent-tabs-mode: nil; coding: utf-8-with-signature -*-
//-----------------------------------------------------------------------------
// Copyright 2000-2022 CEA (www.cea.fr) IFPEN (www.ifpenergiesnouvelles.com)
// See the top-level COPYRIGHT file for details.
// SPDX-License-Identifier: Apache-2.0
//-----------------------------------------------------------------------------
/*---------------------------------------------------------------------------*/
/* BarycentricGeomShapeComputer.h                              (C) 2000-2016 */
/*                                                                           */
/* Calcul des GeomShape en utilisant les barycentres.                        */
/*---------------------------------------------------------------------------*/
#ifndef ARCANE_GEOMETRIC_BARYCENTRICGEOMSHAPECOMPUTER_H
#define ARCANE_GEOMETRIC_BARYCENTRICGEOMSHAPECOMPUTER_H
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
#include "arcane/VariableTypes.h"
 
#include "arcane/geometric/GeomShapeView.h"
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
ARCANE_BEGIN_NAMESPACE
GEOMETRIC_BEGIN_NAMESPACE
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
class ARCANE_CEA_GEOMETRIC_EXPORT BarycentricGeomShapeComputer
{
 public:

  static void computeAll(GeomShapeMutableView elem,const VariableNodeReal3& coords,Cell cell);

  static void computeAll(GeomShapeMng& shape_mng,VariableNodeReal3& coords,const CellGroup& cells);

  static void computeTriangle3(GeomShapeMutableView elem);
  static void computeQuad4(GeomShapeMutableView elem);
  static void computeTetraedron4(GeomShapeMutableView elem);
  static void computePyramid5(GeomShapeMutableView elem);
  static void computePentaedron6(GeomShapeMutableView elem);
  static void computeHexaedron8(GeomShapeMutableView elem);
  static void computeHeptaedron10(GeomShapeMutableView elem);
  static void computeOctaedron12(GeomShapeMutableView elem);

  template<GeomType ItemType> static
  void compute(GeomShapeMutableView elem);

  static void setNodes(GeomShapeMutableView elem,const VariableNodeReal3& node_coord,Cell cell)
  {
    Integer nb_node = cell.nbNode();
    for( Integer node_id=0; node_id<nb_node; ++node_id){
      elem.setNode(node_id,node_coord[cell.node(node_id)]);
    }
  }
 
 private:
 
  inline static void
  _setFace3D(Integer fid,GeomShapeMutableView& elem,Integer id1,Integer id2,Integer id3,Integer id4)
  {
    elem.setFace(fid, 0.25 * ( elem.node(id1) + elem.node(id2) + elem.node(id3) + elem.node(id4) ));
  }
 
  inline static void
  _setFace3D(Integer fid,GeomShapeMutableView& elem,Integer id1,Integer id2,Integer id3)
  {
    elem.setFace(fid, (1.0/3.0) * ( elem.node(id1) + elem.node(id2) + elem.node(id3) ));
  }
 
  inline static void
  _setFace2D(Integer fid,GeomShapeMutableView& elem,Integer id1,Integer id2)
  {
    elem.setFace(fid,
                 Real3( 0.5 * ( elem.node(id1).x + elem.node(id2).x ),
                        0.5 * ( elem.node(id1).y + elem.node(id2).y ),
                        0.0));
  }
};
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
GEOMETRIC_END_NAMESPACE
ARCANE_END_NAMESPACE
 
/*---------------------------------------------------------------------------*/
/*---------------------------------------------------------------------------*/
 
#endif