exp

kyosu::exp = {}

inlineconstexpr

Computes the exponential of the argument.

Defined in Header

#include <kyosu/functions.hpp>

Callable Signatures

namespace kyosu
{
   template<eve::floating_ordered_value T>     constexpr T exp(T z) noexcept; //1
   template<kyosu::concepts::ccomplex T>       constexpr T exp(T z) noexcept; //2
   template<kyosu::concepts::cayley_dickson T> constexpr T exp(T z) noexcept; //3
}

Parameters

• z: Value to process.

Return value

1. Returns the exponential of the argument, calling eve::exp.
2. Returns the exponential of the complex input following IEEE standards:
   • for every z: kyosu::exp(eve::conj(z)) == kyosu::conj(std::exp(z))
   • If z is \(\pm0\), the result is \(1\)
   • If z is \(x+i \infty\) (for any finite x), the result is \(NaN+i NaN\).
   • If z is \(x+i NaN\) (for any finite x), the result is \(NaN+i NaN\).
   • If z is \(+\infty+i 0\), the result is \(+\infty\)
   • If z is \(-\infty+i y\) (for any finite y), the result is \(+0 \mathrm{cis}(y)\).
   • If z is \(+\infty+i y\) (for any finite nonzero y), the result is \(+\infty \mathrm{cis}(y)\).
   • If z is \(-\infty+i \infty\), the result is \(\pm 0+i \pm 0\) (signs are unspecified)
   • If z is \(+\infty+i \pm\infty\), the result is \(\pm \infty+i NaN\) (the sign of the real part is unspecified).
   • If z is \(-\infty+i NaN\), the result is \(\pm 0+i \pm 0\) (signs are unspecified).
   • If z is \(\pm\infty+i NaN\), the result is \(\pm \infty+i NaN\) (the sign of the real part is unspecified).
   • If z is \(NaN\), the result is \(NaN\).
   • If z is \(NaN+i y\) (for any nonzero y), the result is \(NaN+i NaN\).
   • If z is \(NaN+i NaN\), the result is \(NaN+i NaN\).
3. Returns \(e^{z_0}(\cos|\underline{z}|+\underline{z}\; \mathop{sinc}|\underline{z}|)\)

Example

#include <kyosu/kyosu.hpp>
#include <eve/wide.hpp>
#include <iostream>
 
int main()
{
  using kyosu::exp;
  using kyosu::complex_t;
  using kyosu::quaternion_t;
  using e_t = float;
  using c_t = kyosu::complex_t<float>;
  using q_t = kyosu::quaternion_t<float>;
  using we_t = eve::wide<e_t, eve::fixed<2>>;
  using wc_t = eve::wide<c_t, eve::fixed<2>>;
  using wq_t = eve::wide<q_t, eve::fixed<2>>;
 
  std::cout << "Real:        \n";
  e_t e(2.9f);
  we_t we = we_t(e);
  std::cout << e << " -> " << exp(e) << "\n";
  std::cout << we << " -> " << exp(we) << "\n";
  std::cout               << exp(c_t(e))<< "\n";
  std::cout               << exp(q_t(e))<< "\n";
  std::cout               << exp(wc_t(e))<< "\n";
  std::cout               << exp(wq_t(e))<< "\n";
 
  std::cout << "Complex:     \n";
  c_t c(3.5f,-2.9f);
  wc_t wc = wc_t(c);
  std::cout << c << " -> " << exp(c) << "\n";
  std::cout << wc << " -> " << exp(wc) << "\n";
  std::cout               << exp(q_t(c))<< "\n";
  std::cout               << exp(wq_t(c))<< "\n";
 
  std::cout << "Quaternion:  \n";
  q_t q(3.5f,-2.9f, 2.1f, 3.2f);
  wq_t wq = wq_t(q);
  std::cout << q << " -> " << exp(q) << "\n";
  std::cout << wq << " -> " << exp(wq) << "\n";
 
  return 0;
}