[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: expm1 & log1p
From: |
David Bateman |
Subject: |
Re: expm1 & log1p |
Date: |
Tue, 18 Mar 2008 16:26:01 +0100 |
User-agent: |
Thunderbird 2.0.0.12 (X11/20080306) |
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 address@hidden
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