Intrepid2_PointTools.hpp

Go to the documentation of this file.

// @HEADER
// ************************************************************************
//
//                           Intrepid2 Package
//                 Copyright (2007) Sandia Corporation
//
// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
// license for use of this work by or on behalf of the U.S. Government.
//
// Redistribution and use in source and binary forms, with or without
// modification, are permitted provided that the following conditions are
// met:
//
// 1. Redistributions of source code must retain the above copyright
// notice, this list of conditions and the following disclaimer.
//
// 2. Redistributions in binary form must reproduce the above copyright
// notice, this list of conditions and the following disclaimer in the
// documentation and/or other materials provided with the distribution.
//
// 3. Neither the name of the Corporation nor the names of the
// contributors may be used to endorse or promote products derived from
// this software without specific prior written permission.
//
// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
//
// Questions? Contact Kyungjoo Kim  (kyukim@sandia.gov), or
//                    Mauro Perego  (mperego@sandia.gov)
//
// ************************************************************************
// @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 getGaussPoints(      Kokkos::DynRankView<pointValueType,pointProperties...> points,
      const ordinal_type order );

private:

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


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


  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
  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