operator_multiplication.hpp

Go to the documentation of this file.

#ifndef STAN_MATH_PRIM_CORE_OPERATOR_MULTIPLICATION_HPP
#define STAN_MATH_PRIM_CORE_OPERATOR_MULTIPLICATION_HPP
 
#include <stan/math/prim/meta/return_type.hpp>
#include <stan/math/prim/meta/is_complex.hpp>
#include <stan/math/prim/meta/is_stan_scalar.hpp>
#include <complex>
 
namespace stan {
namespace math {
 
namespace internal {
template <typename U, typename V, require_any_complex_t<U, V>* = nullptr>
inline complex_return_t<U, V> complex_multiply(const U& lhs, const V& rhs) {
  complex_return_t<U, V> y(lhs);
  y *= rhs;
  return y;
}
}  // namespace internal
 
template <typename U, typename V, require_all_stan_scalar_t<U, V>* = nullptr>
inline complex_return_t<U, V> operator*(const std::complex<U>& x,
                                        const std::complex<V>& y) {
  return internal::complex_multiply(x, y);
}
 
template <typename U, typename V, require_stan_scalar_t<V>* = nullptr>
inline complex_return_t<U, V> operator*(const std::complex<U>& x, const V& y) {
  return internal::complex_multiply(x, y);
}
 
template <typename U, typename V, require_stan_scalar_t<U>* = nullptr>
inline complex_return_t<U, V> operator*(const U& x, const std::complex<V>& y) {
  return internal::complex_multiply(x, y);
}
 
}  // namespace math
}  // namespace stan
 
#endif