complex_overloads.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2015 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Part of the source code is dual licensed under Apache-2.0 WITH
// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
// governing the source code and code contributions can be found in
// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
//
// ------------------------------------------------------------------------
 
#ifndef dealii_complex_overloads_h
#define dealii_complex_overloads_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/template_constraints.h>
 
DEAL_II_NAMESPACE_OPEN
 
// Forward declarations
#ifndef DOXYGEN
template <typename T, typename U>
struct ProductType;
#endif
 
#ifndef DEAL_II_HAVE_COMPLEX_OPERATOR_OVERLOADS
template <typename T, typename U>
inline std::enable_if_t<
  std::is_floating_point_v<T> && std::is_floating_point_v<U>,
  typename ProductType<std::complex<T>, std::complex<U>>::type>
operator*(const std::complex<T> &left, const std::complex<U> &right)
{
  using result_type =
    typename ProductType<std::complex<T>, std::complex<U>>::type;
  return static_cast<result_type>(left) * static_cast<result_type>(right);
}
 
template <typename T, typename U>
inline std::enable_if_t<
  std::is_floating_point_v<T> && std::is_floating_point_v<U>,
  typename ProductType<std::complex<T>, std::complex<U>>::type>
operator/(const std::complex<T> &left, const std::complex<U> &right)
{
  using result_type =
    typename ProductType<std::complex<T>, std::complex<U>>::type;
  return static_cast<result_type>(left) / static_cast<result_type>(right);
}
 
template <typename T, typename U>
inline std::enable_if_t<std::is_floating_point_v<T> &&
                          std::is_floating_point_v<U>,
                        typename ProductType<std::complex<T>, U>::type>
operator*(const std::complex<T> &left, const U &right)
{
  using result_type = typename ProductType<std::complex<T>, U>::type;
  return static_cast<result_type>(left) * static_cast<result_type>(right);
}
 
template <typename T, typename U>
inline std::enable_if_t<std::is_floating_point_v<T> &&
                          std::is_floating_point_v<U>,
                        typename ProductType<std::complex<T>, U>::type>
operator/(const std::complex<T> &left, const U &right)
{
  using result_type = typename ProductType<std::complex<T>, U>::type;
  return static_cast<result_type>(left) / static_cast<result_type>(right);
}
 
template <typename T, typename U>
inline std::enable_if_t<std::is_floating_point_v<T> &&
                          std::is_floating_point_v<U>,
                        typename ProductType<T, std::complex<U>>::type>
operator*(const T &left, const std::complex<U> &right)
{
  using result_type = typename ProductType<T, std::complex<U>>::type;
  return static_cast<result_type>(left) * static_cast<result_type>(right);
}
 
template <typename T, typename U>
inline std::enable_if_t<std::is_floating_point_v<T> &&
                          std::is_floating_point_v<U>,
                        typename ProductType<T, std::complex<U>>::type>
operator/(const T &left, const std::complex<U> &right)
{
  using result_type = typename ProductType<T, std::complex<U>>::type;
  return static_cast<result_type>(left) / static_cast<result_type>(right);
}
#endif /* DEAL_II_HAVE_COMPLEX_OPERATOR_OVERLOADS */
 
DEAL_II_NAMESPACE_CLOSE
 
#endif