Intrepid2_PointTools.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//                           Intrepid2 Package
//
// Copyright 2007 NTESS and the Intrepid2 contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER

#ifndef __INTREPID2_POINTTOOLS_HPP__
#define __INTREPID2_POINTTOOLS_HPP__

#include "Intrepid2_ConfigDefs.hpp"

#include "Intrepid2_Types.hpp"
#include "Intrepid2_Utils.hpp"

#include "Shards_CellTopology.hpp"

#include "Intrepid2_Polylib.hpp"
#include "Intrepid2_CellTools.hpp"

#include "Kokkos_Core.hpp"

namespace Intrepid2 {

class PointTools {
public:
  inline
  static ordinal_type
  getLatticeSize( const shards::CellTopology cellType,
      const ordinal_type order,
      const ordinal_type offset = 0 );

  template<typename pointValueType, class ...pointProperties>
  static void
  getLattice(     Kokkos::DynRankView<pointValueType,pointProperties...> points,
      const shards::CellTopology cellType,
      const ordinal_type order,
      const ordinal_type offset = 0 ,
      const EPointType pointType = POINTTYPE_EQUISPACED );

  template<typename pointValueType, class ...pointProperties>
  static void
  getLatticeLine(     Kokkos::DynRankView<pointValueType,pointProperties...> points,
      const ordinal_type order,
      const ordinal_type offset = 0 ,
      const EPointType pointType = POINTTYPE_EQUISPACED );

  template<typename pointValueType, class ...pointProperties>
  static void
  getLatticeTriangle(     Kokkos::DynRankView<pointValueType,pointProperties...> points,
      const ordinal_type order,
      const ordinal_type offset = 0 ,
      const EPointType pointType = POINTTYPE_EQUISPACED );

  template<typename pointValueType, class ...pointProperties>
  static void
  getLatticeTetrahedron(     Kokkos::DynRankView<pointValueType,pointProperties...> points,
      const ordinal_type order,
      const ordinal_type offset = 0 ,
      const EPointType pointType = POINTTYPE_EQUISPACED );

  template<typename pointValueType, class ...pointProperties>
  static void
  getLatticePyramid(     Kokkos::DynRankView<pointValueType,pointProperties...> points,
      const ordinal_type order,
      const ordinal_type offset = 0 ,
      const EPointType pointType = POINTTYPE_EQUISPACED );

  template<typename pointValueType, class ...pointProperties>
  static void getGaussPoints(      Kokkos::DynRankView<pointValueType,pointProperties...> points,
      const ordinal_type order );

private:

  template<typename pointValueType, class ...pointProperties>
  static void
  getEquispacedLatticeLine(     Kokkos::DynRankView<pointValueType,pointProperties...> points,
      const ordinal_type order,
      const ordinal_type offset = 0 );


  template<typename pointValueType, class ...pointProperties>
  static void
  getWarpBlendLatticeLine(       Kokkos::DynRankView<pointValueType,pointProperties...> points,
      const ordinal_type order ,
      const ordinal_type offset = 0 );


  // /* \brief Converts Cartesian coordinates to barycentric coordinates
  //     on a batch of triangles.
  //     The input array cartValues is (C,P,2)
  //     The output array baryValues is (C,P,3).
  //     The input array vertices is (C,3,2), where
  //     \code
  //     C - num. integration domains
  //     P - number of points per cell
  //     \endcode

  //     \param  baryValues      [out] - Output array of barycentric coords
  //     \param  cartValues      [in]  - Input array of Cartesian coords
  //     \param  vertices        [out] - Vertices of each cell.

  // */
  // template<class Scalar, class ArrayTypeOut, class ArrayTypeIn1, class ArrayTypeIn2>
  // static void cartToBaryTriangle( ArrayTypeOut & baryValues ,
  //                                 const ArrayTypeIn1 & cartValues ,
  //                                 const ArrayTypeIn2 & vertices );


  // /* \brief Converts barycentric coordinates to Cartesian coordinates
  //     on a batch of triangles.
  //     The input array baryValues is (C,P,3)
  //     The output array cartValues is (C,P,2).
  //     The input array vertices is (C,3,2), where
  //     \code
  //     C - num. integration domains
  //     P - number of points per cell
  //     D - is the spatial dimension
  //     \endcode

  //     \param  baryValues      [out] - Output array of barycentric coords
  //     \param  cartValues      [in]  - Input array of Cartesian coords
  //     \param  vertices        [out] - Vertices of each cell.

  // */
  // template<class Scalar, class ArrayTypeOut, class ArrayTypeIn1, class ArrayTypeIn2>
  // static void baryToCartTriangle( ArrayTypeOut & cartValues ,
  //                                 const ArrayTypeIn1 & baryValues ,
  //                                 const ArrayTypeIn2 & vertices );


  // /* \brief Converts Cartesian coordinates to barycentric coordinates
  //     on a batch of tetrahedra.
  //     The input array cartValues is (C,P,3)
  //     The output array baryValues is (C,P,4).
  //     The input array vertices is (C,4,3), where
  //     \code
  //     C - num. integration domains
  //     P - number of points per cell
  //     D - is the spatial dimension
  //     \endcode

  //     \param  baryValues      [out] - Output array of barycentric coords
  //     \param  cartValues      [in]  - Input array of Cartesian coords
  //     \param  vertices        [out] - Vertices of each cell.

  // */
  // template<class Scalar, class ArrayTypeOut, class ArrayTypeIn1, class ArrayTypeIn2>
  // static void cartToBaryTetrahedron( ArrayTypeOut & baryValues ,
  //                                    const ArrayTypeIn1 & cartValues ,
  //                                    const ArrayTypeIn2 & vertices );

  // /* \brief Converts barycentric coordinates to Cartesian coordinates
  //     on a batch of tetrahedra.
  //     The input array baryValues is (C,P,4)
  //     The output array cartValues is (C,P,3).
  //     The input array vertices is (C,4,3), where
  //     \code
  //     C - num. integration domains
  //     P - number of points per cell
  //     D - is the spatial dimension
  //     \endcode

  //     \param  baryValues      [out] - Output array of barycentric coords
  //     \param  cartValues      [in]  - Input array of Cartesian coords
  //     \param  vertices        [out] - Vertices of each cell.

  // */
  // template<class Scalar, class ArrayTypeOut, class ArrayTypeIn1, class ArrayTypeIn2>
  // static void baryToCartTetrahedron( ArrayTypeOut & cartValues ,
  //                                    const ArrayTypeIn1 & baryValues ,
  //                                    const ArrayTypeIn2 & vertices );




  template<typename pointValueType, class ...pointProperties>
  static void
  getEquispacedLatticeTriangle(       Kokkos::DynRankView<pointValueType,pointProperties...> points,
      const ordinal_type order,
      const ordinal_type offset = 0 );

  template<typename pointValueType, class ...pointProperties>
  static void
  getEquispacedLatticeTetrahedron( Kokkos::DynRankView<pointValueType,pointProperties...> points,
                                   const ordinal_type order ,
                                   const ordinal_type offset = 0 );


 template<typename pointValueType, class ...pointProperties>
 static void
 getEquispacedLatticePyramid(       Kokkos::DynRankView<pointValueType,pointProperties...> points,
     const ordinal_type order,
     const ordinal_type offset = 0 );


  template<typename pointValueType, class ...pointProperties>
  static void
  warpFactor( Kokkos::DynRankView<pointValueType,pointProperties...>  warp ,
              const ordinal_type order ,
              const Kokkos::DynRankView<pointValueType,pointProperties...> xnodes ,
              const Kokkos::DynRankView<pointValueType,pointProperties...> xout
              );

  template<typename pointValueType, class ...pointProperties>
  static void
  getWarpBlendLatticeTriangle( Kokkos::DynRankView<pointValueType,pointProperties...> points,
                               const ordinal_type order ,
                               const ordinal_type offset = 0 );

  template<typename pointValueType, class ...pointProperties>
  static void
  getWarpBlendLatticeTetrahedron( Kokkos::DynRankView<pointValueType,pointProperties...> points,
                                  const ordinal_type order ,
                                  const ordinal_type offset = 0 );


  template<typename pointValueType, class ...pointProperties>
  static void
  warpShiftFace3D( Kokkos::DynRankView<pointValueType,pointProperties...> dxy,
                   const ordinal_type order ,
                   const pointValueType pval ,
                   const Kokkos::DynRankView<pointValueType,pointProperties...> L1,
                   const Kokkos::DynRankView<pointValueType,pointProperties...> L2,
                   const Kokkos::DynRankView<pointValueType,pointProperties...> L3,
                   const Kokkos::DynRankView<pointValueType,pointProperties...> L4
                   );

  template<typename pointValueType, class ...pointProperties>
  static void
  evalshift( Kokkos::DynRankView<pointValueType,pointProperties...> dxy ,
             const ordinal_type order ,
             const pointValueType pval ,
             const Kokkos::DynRankView<pointValueType,pointProperties...> L1 ,
             const Kokkos::DynRankView<pointValueType,pointProperties...> L2 ,
             const Kokkos::DynRankView<pointValueType,pointProperties...> L3
             );

  template<typename pointValueType, class ...pointProperties>
  static void
  evalwarp( Kokkos::DynRankView<pointValueType,pointProperties...> warp ,
            const ordinal_type order ,
            const Kokkos::DynRankView<pointValueType,pointProperties...> xnodes ,
            const Kokkos::DynRankView<pointValueType,pointProperties...> xout );
};

} 

#include "Intrepid2_PointToolsDef.hpp"

#endif