[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
RE: [Help-gsl] Exponential calculation
From: |
Ashish Negi |
Subject: |
RE: [Help-gsl] Exponential calculation |
Date: |
Mon, 16 Feb 2009 11:51:29 +0530 |
Hi Kevin,
Thanks.
You are right exp(-100) returns 3.7e-44.
Actually, during mathematical programming. We are encountering such
scenarios while calculating probability density function.
As structure of probability density function is as follows.
Pdf = (expression_1) * (exponential (expression_2).
Now during our calculation we are encountering the situation where
expression_2 is getting evaluated to such values like -800.0.
In this case exponential (-800.0) is zero.
So how do we deal with such situations what will be pdf in this case.
Thanks
Ashish
-----Original Message-----
From: address@hidden [mailto:address@hidden On Behalf
Of Kevin Channon
Sent: Wednesday, February 11, 2009 6:49 PM
To: Ashish Negi
Cc: address@hidden
Subject: Re: [Help-gsl] Exponential calculation
Hi Ashish,
> Is there any known mathematical method to calculate : exponential( -
> val), when val is too small number like -100
You don't need gsl for this, just use the standard c library math.h (or
cmath in c++)
printf("%g\n", exp(-100));
returns the correct answer of 3.7e-44
There might be some penalty in terms of speed and accuracy, I don't
know, but it does return a value that should be at least in the ballpark
>
> If there is any such instance during probability density function
> calculation when exp( - x2/ 2) is calculated, How to handle such
events?
You might need to clarify this a bit more, I'm not sure that I
understand
what you're asking
Cheers,
Kevin
--
Kevin Channon
School of Chemistry
University of Bristol
UK