expm1 & log1p

[David Bateman]

David Bateman wrote:
> Jaroslav Hajek wrote:
>   
>>> Also the matlab
>>>  functions accept complex arguments whereas your version seems to be
>>>  limited to real arguments
>>>
>>>     
>>>       
>> I missed that. Apparently there are only real versions in GNU C
>> library. I'll try to make
>> a complex version using the real version.
>>
>>   
>>
>>     
>
> Complex expm1 (const Complex& x)
> {
>     double xr = x.real ();
>     Complex xc = Complex (0., x.imag ());
>     return exp (xc) * expm1 (xr);
> }
>
> D.
>    
>
>   
Opps

Complex expm1 (const Complex& x)
{
  return ((exp (Complex (0., x.imag ())) + 1.0) * expm1 (x.real ()));
}

D.


-- 
David Bateman                                David.Bateman at motorola.com
Motorola Labs - Paris                        +33 1 69 35 48 04 (Ph) 
Parc Les Algorithmes, Commune de St Aubin    +33 6 72 01 06 33 (Mob) 
91193 Gif-Sur-Yvette FRANCE                  +33 1 69 35 77 01 (Fax) 

The information contained in this communication has been classified as: 

[x] General Business Information 
[ ] Motorola Internal Use Only 
[ ] Motorola Confidential Proprietary