Amesos2_ShyLUBasker_FunctionMap.hpp

// @HEADER
// *****************************************************************************
//           Amesos2: Templated Direct Sparse Solver Package
//
// Copyright 2011 NTESS and the Amesos2 contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef AMESOS2_SHYLUBASKER_FUNCTIONMAP_HPP
#define AMESOS2_SHYLUBASKER_FUNCTIONMAP_HPP
 
#ifdef HAVE_TEUCHOS_COMPLEX
#include <complex>
#endif
 
#include "Amesos2_FunctionMap.hpp"
#include "Amesos2_ShyLUBasker_TypeMap.hpp"
 
//#include "shylubasker_decl.hpp"
//#include "shylubasker_def.hpp"
//#include "shylubasker_trilinos_decl.hpp"
 
namespace Amesos2 {
 
  /* ==================== Specializations ====================
   *
   * \cond ShyLUBasker_function_specializations
   */
 
#ifdef HAVE_TEUCHOS_INST_COMPLEX_DOUBLE
  template <>
  struct FunctionMap<ShyLUBasker,Kokkos::complex<double>>
  {
    static std::complex<double> * convert_scalar(Kokkos::complex<double> * pData) {
      return reinterpret_cast<std::complex<double> *>(pData);
    }
  };
#endif // HAVE_TEUCHOS_COMPLEX_DOUBLE
 
#ifdef HAVE_TEUCHOS_INST_COMPLEX_FLOAT
  template <>
  struct FunctionMap<ShyLUBasker,Kokkos::complex<float>>
  {
    static std::complex<float> * convert_scalar(Kokkos::complex<float> * pData) {
      return reinterpret_cast<std::complex<float> *>(pData);
    }
  };
#endif // HAVE_TEUCHOS_INST_COMPLEX_FLOAT
 
  // if not specialized, then assume generic conversion is fine
  template <typename scalar_t>
  struct FunctionMap<ShyLUBasker,scalar_t>
  {
    static scalar_t * convert_scalar(scalar_t * pData) {
      return pData; // no conversion necessary
    }
  };
 
 
  /* \endcond ShyLUBasker_function_specializations */
 
 
} // end namespace Amesos2
 
#endif  // AMESOS2_SHYLUBASKER_FUNCTIONMAP_HPP