[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: expm1 & log1p
From: |
Jaroslav Hajek |
Subject: |
Re: expm1 & log1p |
Date: |
Tue, 18 Mar 2008 15:35:51 +0100 |
On Tue, Mar 18, 2008 at 3:21 PM, David Bateman
<address@hidden> wrote:
> Jaroslav Hajek wrote:
> > another pair of mappers.
> > These are in GNU C library, so usually they will be in math.h; that's
> > why I supplied "quick hack" replacements from top of my head
> > (actually, pencil, paper & computer) instead of stealing the
> > complicated formulas from glibc or other software. The minor advantage
> > is that there is no need to sort out any copyright notices.
> >
> >
> <quote>
>
> +#if defined (HAVE_EXPM1)
> +double
> +expm1 (double x)
> +
>
> </quote>
>
> Shouldn't that read "#if ! defined (HAVE_EXPM1)" ?
Indeed. I'll fix it.
> 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.
> 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
>
>
--
RNDr. Jaroslav Hajek
computing expert
Aeronautical Research and Test Institute (VZLU)
Prague, Czech Republic
url: www.highegg.matfyz.cz