sl@0: /*
sl@0:  * © Portions copyright (c) 2006-2007 Nokia Corporation.  All rights reserved.
sl@0:  *
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: # include "stlport_prefix.h"
sl@0: 
sl@0: 
sl@0: // Trigonometric and hyperbolic functions for complex<float>, 
sl@0: // complex<double>, and complex<long double>
sl@0: 
sl@0: 
sl@0: #include "complex_impl.h"
sl@0: 
sl@0: #include <cfloat>
sl@0: #include <cmath>
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: //----------------------------------------------------------------------
sl@0: // helpers
sl@0: 
sl@0: #ifdef __sgi
sl@0:   static const union { unsigned int i; float f; } float_ulimit = { 0x42b2d4fc };
sl@0:   static const float float_limit = float_ulimit.f;
sl@0:   static union {
sl@0:     struct { unsigned int h; unsigned int l; } w;
sl@0:     double d;
sl@0:   } double_ulimit = { 0x408633ce, 0x8fb9f87d };
sl@0:   static const double double_limit = double_ulimit.d;
sl@0:   static union {
sl@0:     struct { unsigned int h[2]; unsigned int l[2]; } w;
sl@0:     long double ld;
sl@0:   } ldouble_ulimit = {0x408633ce, 0x8fb9f87e, 0xbd23b659, 0x4e9bd8b1};
sl@0: # ifndef _STLP_NO_LONG_DOUBLE
sl@0:   static const long double ldouble_limit = ldouble_ulimit.ld;
sl@0: # endif
sl@0: #else
sl@0: #if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
sl@0: void complex_trig_limit_init()
sl@0: {
sl@0: 	get_complex_trig_float_limit() = _STLP_LOGF(FLT_MAX);
sl@0: 	get_complex_trig_double_limit() = _STLP_DO_LOG(double)(DBL_MAX);    
sl@0: } 
sl@0: # else
sl@0:   static const float float_limit = _STLP_LOGF(FLT_MAX);
sl@0:   static const double double_limit = _STLP_DO_LOG(double)(DBL_MAX);  
sl@0: # endif //__LIBSTD_CPP_SYMBIAN32_WSD__
sl@0: # ifndef _STLP_NO_LONG_DOUBLE
sl@0:   static const long double ldouble_limit = _STLP_LOGL(LDBL_MAX);
sl@0: # endif
sl@0: #endif
sl@0: 
sl@0: 
sl@0: //----------------------------------------------------------------------
sl@0: // sin
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float>  _STLP_CALL sin(const complex<float>& z) {
sl@0:   return complex<float>(_STLP_SINF(z._M_re) * _STLP_COSHF(z._M_im),
sl@0:                         _STLP_COSF(z._M_re) * _STLP_SINHF(z._M_im));
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL sin(const complex<double>& z) {
sl@0:   return complex<double>(_STLP_SIN(z._M_re) * _STLP_COSH(z._M_im),
sl@0:                          _STLP_COS(z._M_re) * _STLP_SINH(z._M_im));
sl@0: }
sl@0: 
sl@0: #ifndef _STLP_NO_LONG_DOUBLE
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL sin(const complex<long double>& z) {
sl@0:   return complex<long double>(_STLP_SINL(z._M_re) * _STLP_COSHL(z._M_im),
sl@0:                               _STLP_COSL(z._M_re) * _STLP_SINHL(z._M_im));
sl@0: }
sl@0: #endif
sl@0: 
sl@0: //----------------------------------------------------------------------
sl@0: // cos
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL cos(const complex<float>& z) {
sl@0:   return complex<float>(_STLP_COSF(z._M_re) * _STLP_COSHF(z._M_im),
sl@0:                         -_STLP_SINF(z._M_re) * _STLP_SINHF(z._M_im));
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL cos(const complex<double>& z) {
sl@0:   return complex<double>(_STLP_COS(z._M_re) * _STLP_COSH(z._M_im),
sl@0:                         -_STLP_SIN(z._M_re) * _STLP_SINH(z._M_im));
sl@0: }
sl@0: 
sl@0: #ifndef _STLP_NO_LONG_DOUBLE
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL cos(const complex<long double>& z) {
sl@0:   return complex<long double>(_STLP_COSL(z._M_re) * _STLP_COSHL(z._M_im),
sl@0:                               -_STLP_SINL(z._M_re) * _STLP_SINHL(z._M_im));
sl@0: }
sl@0: # endif
sl@0: 
sl@0: //----------------------------------------------------------------------
sl@0: // tan
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL tan(const complex<float>& z) {
sl@0:   float re2 = 2.f * z._M_re;
sl@0:   float im2 = 2.f * z._M_im;
sl@0: 
sl@0: #if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
sl@0:   if (_STLP_ABSF(im2) > get_complex_trig_float_limit())  
sl@0: # else
sl@0:   if (_STLP_ABSF(im2) > float_limit)	
sl@0: # endif  
sl@0:     return complex<float>(0.f, (im2 > 0 ? 1.f : -1.f));
sl@0:   else {
sl@0:     float den = _STLP_COSF(re2) + _STLP_COSHF(im2);
sl@0:     return complex<float>(_STLP_SINF(re2) / den, _STLP_SINHF(im2) / den);
sl@0:   }
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL tan(const complex<double>& z) {
sl@0:   double re2 = 2. * z._M_re;
sl@0:   double im2 = 2. * z._M_im;
sl@0: 
sl@0: #if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
sl@0:   if (fabs(float(im2)) > get_complex_trig_double_limit())
sl@0: # else
sl@0:   if (fabs(float(im2)) > double_limit)
sl@0: # endif  //__LIBSTD_CPP_SYMBIAN32_WSD__
sl@0:     return complex<double>(0., (im2 > 0 ? 1. : -1.));
sl@0:   else {
sl@0:     double den = _STLP_COS(re2) + _STLP_COSH(im2);
sl@0:     return complex<double>(_STLP_SIN(re2) / den, _STLP_SINH(im2) / den);
sl@0:   }
sl@0: }
sl@0: 
sl@0: #ifndef _STLP_NO_LONG_DOUBLE
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL tan(const complex<long double>& z) {
sl@0:   long double re2 = 2.l * z._M_re;
sl@0:   long double im2 = 2.l * z._M_im;
sl@0:   if (_STLP_ABSL(im2) > ldouble_limit)
sl@0:     return complex<long double>(0.l, (im2 > 0 ? 1.l : -1.l));
sl@0:   else {
sl@0:     long double den = _STLP_COSL(re2) + _STLP_COSHL(im2);
sl@0:     return complex<long double>(_STLP_SINL(re2) / den, _STLP_SINHL(im2) / den);
sl@0:   }
sl@0: }
sl@0: 
sl@0: # endif
sl@0: 
sl@0: //----------------------------------------------------------------------
sl@0: // sinh
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL sinh(const complex<float>& z) {
sl@0:   return complex<float>(_STLP_SINHF(z._M_re) * _STLP_COSF(z._M_im),
sl@0:                         _STLP_COSHF(z._M_re) * _STLP_SINF(z._M_im));
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL sinh(const complex<double>& z) {
sl@0:   return complex<double>(_STLP_SINH(z._M_re) * _STLP_COS(z._M_im),
sl@0:                          _STLP_COSH(z._M_re) * _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 sinh(const complex<long double>& z) {
sl@0:   return complex<long double>(_STLP_SINHL(z._M_re) * _STLP_COSL(z._M_im),
sl@0:                               _STLP_COSHL(z._M_re) * _STLP_SINL(z._M_im));
sl@0: }
sl@0: #endif
sl@0: 
sl@0: //----------------------------------------------------------------------
sl@0: // cosh
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL cosh(const complex<float>& z) {
sl@0:   return complex<float>(_STLP_COSHF(z._M_re) * _STLP_COSF(z._M_im),
sl@0:                         _STLP_SINHF(z._M_re) * _STLP_SINF(z._M_im));
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL cosh(const complex<double>& z) {
sl@0:   return complex<double>(_STLP_COSH(z._M_re) * _STLP_COS(z._M_im),
sl@0:                          _STLP_SINH(z._M_re) * _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 cosh(const complex<long double>& z) {
sl@0:   return complex<long double>(_STLP_COSHL(z._M_re) * _STLP_COSL(z._M_im),
sl@0:                               _STLP_SINHL(z._M_re) * _STLP_SINL(z._M_im));
sl@0: }
sl@0: #endif
sl@0: 
sl@0: //----------------------------------------------------------------------
sl@0: // tanh
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<float> _STLP_CALL tanh(const complex<float>& z) {
sl@0:   float re2 = 2.f * z._M_re;
sl@0:   float im2 = 2.f * z._M_im;
sl@0: #if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
sl@0:    if (_STLP_ABSF(re2) > get_complex_trig_float_limit())
sl@0: # else
sl@0:    if (_STLP_ABSF(re2) > float_limit)
sl@0: # endif  //__LIBSTD_CPP_SYMBIAN32_WSD__
sl@0:     return complex<float>((re2 > 0 ? 1.f : -1.f), 0.f);
sl@0:   else {
sl@0:     float den = _STLP_COSHF(re2) + _STLP_COSF(im2);
sl@0:     return complex<float>(_STLP_SINHF(re2) / den, _STLP_SINF(im2) / den);
sl@0:   }
sl@0: }
sl@0: 
sl@0: _STLP_EXP_DECLSPEC complex<double> _STLP_CALL tanh(const complex<double>& z) {
sl@0:   double re2 = 2. * z._M_re;
sl@0:   double im2 = 2. * z._M_im;  
sl@0:  
sl@0: #if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
sl@0:   if (fabs(float(re2)) > get_complex_trig_double_limit())
sl@0: # else
sl@0:   if (fabs(float(re2)) > double_limit)
sl@0: # endif //__LIBSTD_CPP_SYMBIAN32_WSD__  
sl@0:     return complex<double>((re2 > 0 ? 1. : -1.), 0.);
sl@0:   else {
sl@0:     double den = _STLP_COSH(re2) + _STLP_COS(im2);
sl@0:     return complex<double>(_STLP_SINH(re2) / den, _STLP_SIN(im2) / den);
sl@0:   }
sl@0: }
sl@0: 
sl@0: #ifndef _STLP_NO_LONG_DOUBLE
sl@0: _STLP_EXP_DECLSPEC complex<long double> _STLP_CALL tanh(const complex<long double>& z) {
sl@0:   long double re2 = 2.l * z._M_re;
sl@0:   long double im2 = 2.l * z._M_im;
sl@0:   if (_STLP_ABSL(re2) > ldouble_limit)
sl@0:     return complex<long double>((re2 > 0 ? 1.l : -1.l), 0.l);
sl@0:   else {
sl@0:     long double den = _STLP_COSHL(re2) + _STLP_COSL(im2);
sl@0:     return complex<long double>(_STLP_SINHL(re2) / den, _STLP_SINL(im2) / den);
sl@0:   }
sl@0: }
sl@0: #endif
sl@0: _STLP_END_NAMESPACE