complex_overloads.h

Go to the documentation of this file.

// -----------------------------------------------------------------------------
//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception OR LGPL-2.1-or-later
// Copyright (C) 2015 - 2025 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Detailed license information governing the source code and contributions
// can be found in LICENSE.md and CONTRIBUTING.md at the top level directory.
//
// -----------------------------------------------------------------------------
 
#ifndef dealii_complex_overloads_h
#define dealii_complex_overloads_h
 
#include <deal.II/base/config.h>
 
#include <complex>
#include <type_traits>
 
 
DEAL_II_NAMESPACE_OPEN
 
#ifndef DOXYGEN
 
// Forward declarations
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> &&
    !std::is_same_v<T, 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> &&
    !std::is_same_v<T, 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> && !std::is_same_v<T, 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> && !std::is_same_v<T, 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> && !std::is_same_v<T, 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> && !std::is_same_v<T, 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