sl@0: /*
sl@0:  * © Portions copyright (c) 2006-2007 Nokia Corporation.  All rights reserved.
sl@0:  * Copyright (c) 1999
sl@0:  * Silicon Graphics Computer Systems, Inc.
sl@0:  *
sl@0:  * Copyright (c) 1999 
sl@0:  * Boris Fomitchev
sl@0:  *
sl@0:  * This material is provided "as is", with absolutely no warranty expressed
sl@0:  * or implied. Any use is at your own risk.
sl@0:  *
sl@0:  * Permission to use or copy this software for any purpose is hereby granted 
sl@0:  * without fee, provided the above notices are retained on all copies.
sl@0:  * Permission to modify the code and to distribute modified code is granted,
sl@0:  * provided the above notices are retained, and a notice that the code was
sl@0:  * modified is included with the above copyright notice.
sl@0:  *
sl@0:  */ 
sl@0: 
sl@0: # include "stlport_prefix.h"
sl@0: // exp, log, pow for complex<float>, complex<double>, and complex<long double>
sl@0: 
sl@0: #include <numeric>
sl@0: #include "complex_impl.h"
sl@0: 
sl@0: #if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
sl@0: #include "libstdcppwsd.h"
sl@0: # endif
sl@0: 
sl@0: _STLP_BEGIN_NAMESPACE
sl@0: 
sl@0: //----------------------------------------------------------------------
sl@0: // exp
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float>  _STLP_CALL
sl@0: exp(const complex<float>& z)
sl@0: {
sl@0:   float expx = _STLP_EXPF(z._M_re);
sl@0:   return complex<float>(expx * _STLP_COSF(z._M_im),
sl@0:                         expx * _STLP_SINF(z._M_im));
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL exp(const complex<double>& z)
sl@0: {
sl@0:   double expx = _STLP_EXP(z._M_re);
sl@0:   return complex<double>(expx * _STLP_COS(z._M_im),
sl@0:                          expx * _STLP_SIN(z._M_im));
sl@0: }
sl@0: 
sl@0: # ifndef _STLP_NO_LONG_DOUBLE
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL exp(const complex<long double>& z)
sl@0: {
sl@0:   long double expx = _STLP_EXPL(z._M_re);
sl@0:   return complex<long double>(expx * _STLP_COSL(z._M_im),
sl@0:                               expx * _STLP_SINL(z._M_im));
sl@0: }
sl@0: # endif
sl@0: 
sl@0: //----------------------------------------------------------------------
sl@0: // log10
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL log10(const complex<float>& z)
sl@0: {
sl@0:   complex<float> r;
sl@0:   
sl@0: #if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
sl@0:   get_complex_exp_float_ln10_inv() = 1.f / _STLP_LOGF(10.f);  
sl@0:   r._M_im = _STLP_ATAN2F(z._M_im, z._M_re) * get_complex_exp_float_ln10_inv();  	
sl@0: # else
sl@0:   static float ln10_inv = 1.f / _STLP_LOGF(10.f);
sl@0:   r._M_im = _STLP_ATAN2F(z._M_im, z._M_re) * ln10_inv;  
sl@0: # endif //__LIBSTD_CPP_SYMBIAN32_WSD__
sl@0:   
sl@0:   r._M_re = _STLP_LOG10F(_STLP_HYPOTF(z._M_re, z._M_im));
sl@0:   return r;
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL log10(const complex<double>& z)
sl@0: {
sl@0:   complex<double> r;
sl@0: 
sl@0: #if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
sl@0: 	get_complex_exp_double_ln10_inv() = 1. / _STLP_LOG(10.);
sl@0: 	r._M_im = _STLP_ATAN2(z._M_im, z._M_re) * get_complex_exp_double_ln10_inv();
sl@0: # else
sl@0: 	static double ln10_inv = 1. / _STLP_LOG(10.);
sl@0: 	r._M_im = _STLP_ATAN2(z._M_im, z._M_re) * ln10_inv;
sl@0: # endif //__LIBSTD_CPP_SYMBIAN32_WSD__
sl@0:   
sl@0:   r._M_re = _STLP_LOG10(_STLP_HYPOT(z._M_re, z._M_im));
sl@0:   return r;
sl@0: }
sl@0: 
sl@0: #ifndef _STLP_NO_LONG_DOUBLE
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL log10(const complex<long double>& z)
sl@0: {
sl@0:   complex<long double> result;
sl@0: #if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
sl@0:   get_complex_exp_long_double_ln10_inv() = 1.l / _STLP_LOGL(10.l);
sl@0:   result._M_im = _STLP_ATAN2L(z._M_im, z._M_re) * get_complex_exp_long_double_ln10_inv();
sl@0: # else
sl@0:   static long double ln10_inv = 1.l / _STLP_LOGL(10.l);
sl@0:   result._M_im = _STLP_ATAN2L(z._M_im, z._M_re) * ln10_inv;
sl@0: # endif //__LIBSTD_CPP_SYMBIAN32_WSD__
sl@0:  
sl@0:     result._M_re = _STLP_LOG10L(_STLP_HYPOTL(z._M_re, z._M_im));
sl@0:   return result;
sl@0: }
sl@0: # endif
sl@0: 
sl@0: //----------------------------------------------------------------------
sl@0: // log
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL log(const complex<float>& z)
sl@0: {
sl@0:   complex<float> r;
sl@0: 
sl@0:   r._M_im = _STLP_ATAN2F(z._M_im, z._M_re);
sl@0:   r._M_re = _STLP_LOGF(_STLP_HYPOTF(z._M_re, z._M_im));
sl@0:   return r;
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL log(const complex<double>& z)
sl@0: {
sl@0:   complex<double> r;
sl@0: 
sl@0:   r._M_im = _STLP_ATAN2(z._M_im, z._M_re);
sl@0:   r._M_re = _STLP_LOG(_STLP_HYPOT(z._M_re, z._M_im));
sl@0:   return r;
sl@0: }
sl@0: 
sl@0: #ifndef _STLP_NO_LONG_DOUBLE
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL log(const complex<long double>& z)
sl@0: {
sl@0:   complex<long double> result;
sl@0: 
sl@0:   result._M_im = _STLP_ATAN2L(z._M_im, z._M_re);
sl@0:   result._M_re = _STLP_LOGL(_STLP_HYPOTL(z._M_re, z._M_im));
sl@0:   return result;
sl@0: }
sl@0: # endif
sl@0: 
sl@0: //----------------------------------------------------------------------
sl@0: // pow
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL pow(const float& a, const complex<float>& b) {
sl@0:   float logr = _STLP_LOGF(a);
sl@0:   float x = _STLP_EXPF(logr*b._M_re);
sl@0:   float y = logr*b._M_im;
sl@0: 
sl@0:   return complex<float>(x * _STLP_COSF(y), x * _STLP_SINF(y));
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL pow(const complex<float>& z_in, int n) {
sl@0:   complex<float> z = z_in;
sl@0:   z = __power(z, (n < 0 ? -n : n), multiplies< complex<float> >());
sl@0:   if (n < 0)
sl@0:     return 1.f / z;
sl@0:   else
sl@0:     return z;
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL pow(const complex<float>& a, const float& b) {
sl@0:   float logr = _STLP_LOGF(_STLP_HYPOTF(a._M_re,a._M_im));
sl@0:   float logi = _STLP_ATAN2F(a._M_im, a._M_re);
sl@0:   float x = _STLP_EXPF(logr * b);
sl@0:   float y = logi * b;
sl@0: 
sl@0:   return complex<float>(x * _STLP_COSF(y), x * _STLP_SINF(y));
sl@0: }  
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL pow(const complex<float>& a, const complex<float>& b) {
sl@0:   float logr = _STLP_LOGF(_STLP_HYPOTF(a._M_re,a._M_im));
sl@0:   float logi = _STLP_ATAN2F(a._M_im, a._M_re);
sl@0:   float x = _STLP_EXPF(logr*b._M_re - logi*b._M_im);
sl@0:   float y = logr*b._M_im + logi*b._M_re;
sl@0: 
sl@0:   return complex<float>(x * _STLP_COSF(y), x * _STLP_SINF(y));
sl@0: }
sl@0: 
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL pow(const double& a, const complex<double>& b) {
sl@0:   double logr = _STLP_LOG(a);
sl@0:   double x = _STLP_EXP(logr*b._M_re);
sl@0:   double y = logr*b._M_im;
sl@0: 
sl@0:   return complex<double>(x * _STLP_COS(y), x * _STLP_SIN(y));
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL pow(const complex<double>& z_in, int n) {
sl@0:   complex<double> z = z_in;
sl@0:   z = __power(z, (n < 0 ? -n : n), multiplies< complex<double> >());
sl@0:   if (n < 0)
sl@0: #if !defined(__SC__)			//*TY 04/15/2000 - 
sl@0:     return 1. / z;
sl@0: #else							//*TY 04/15/2000 - added workaround for SCpp compiler
sl@0: 	return double(1.0) / z;		//*TY 04/15/2000 - it incorrectly assign long double attribute to floating point literals
sl@0: #endif							//*TY 04/15/2000 - 
sl@0:   else
sl@0:     return z;
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL pow(const complex<double>& a, const double& b) {
sl@0:   double logr = _STLP_LOG(_STLP_HYPOT(a._M_re,a._M_im));
sl@0:   double logi = _STLP_ATAN2(a._M_im, a._M_re);
sl@0:   double x = _STLP_EXP(logr * b);
sl@0:   double y = logi * b;
sl@0: 
sl@0:   return complex<double>(x * _STLP_COS(y), x * _STLP_SIN(y));
sl@0: }  
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL pow(const complex<double>& a, const complex<double>& b) {
sl@0:   double logr = _STLP_LOG(_STLP_HYPOT(a._M_re,a._M_im));
sl@0:   double logi = _STLP_ATAN2(a._M_im, a._M_re);
sl@0:   double x = _STLP_EXP(logr*b._M_re - logi*b._M_im);
sl@0:   double y = logr*b._M_im + logi*b._M_re;
sl@0: 
sl@0:   return complex<double>(x * _STLP_COS(y), x * _STLP_SIN(y));
sl@0: }
sl@0: 
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL pow(const long double& a,
sl@0:                                                    const complex<long double>& b) {
sl@0:   long double logr = _STLP_LOGL(a);
sl@0:   long double x = _STLP_EXPL(logr*b._M_re);
sl@0:   long double y = logr*b._M_im;
sl@0: 
sl@0:   return complex<long double>(x * _STLP_COSL(y), x * _STLP_SINL(y));
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL pow(const complex<long double>& z_in, int n) {
sl@0:   complex<long double> z = z_in;
sl@0:   z = __power(z, (n < 0 ? -n : n), multiplies< complex<long double> >());
sl@0:   if (n < 0)
sl@0:     return 1.l / z;
sl@0:   else
sl@0:     return z;
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL pow(const complex<long double>& a,
sl@0:                          const long double& b) {
sl@0:   long double logr = _STLP_LOGL(_STLP_HYPOTL(a._M_re,a._M_im));
sl@0:   long double logi = _STLP_ATAN2L(a._M_im, a._M_re);
sl@0:   long double x = _STLP_EXPL(logr * b);
sl@0:   long double y = logi * b;
sl@0: 
sl@0:   return complex<long double>(x * _STLP_COSL(y), x * _STLP_SINL(y));
sl@0: }  
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL pow(const complex<long double>& a,
sl@0:                          const complex<long double>& b) {
sl@0:   long double logr = _STLP_LOGL(_STLP_HYPOTL(a._M_re,a._M_im));
sl@0:   long double logi = _STLP_ATAN2L(a._M_im, a._M_re);
sl@0:   long double x = _STLP_EXPL(logr*b._M_re - logi*b._M_im);
sl@0:   long double y = logr*b._M_im + logi*b._M_re;
sl@0: 
sl@0:   return complex<long double>(x * _STLP_COSL(y), x * _STLP_SINL(y));
sl@0: }
sl@0: 
sl@0: _STLP_END_NAMESPACE
sl@0: