[Date Index] [Thread Index] [Author Index]

Re: Clip[] doesn't work as expected

• To: mathgroup at smc.vnet.net
• Subject: [mg115383] Re: Clip[] doesn't work as expected
• From: "Alexey Popkov" <lehin.p at gmail.com>
• Date: Mon, 10 Jan 2011 02:33:52 -0500 (EST)
• References: <igbnil$hi0$1@smc.vnet.net>

"Leonid Shifrin" <lshifr at gmail.com> news:igbnil$hi0$1 at smc.vnet.net...
> This behavior is correct, and there is nothing to warn about. You use
> machine-precision arithmetic,
> and your result is within machine epsilon of 1:
>
> In[1]:=Chop[Clip[1.00000000000000036,{0,1}]-1]
>
> Out[1] = 0
>
> If you want exact result, you need fixed higher  precision, and exact
bounds
> in Clip, such as:
>
> In[2]:= Clip[N[1+36*10^(-17),20],{0,1}]
>
> Out[2]= 1
>
>
> Regards,
> Leonid


I think that there is no excuse for this:
In[6]:= Clip[1.00000000000000003`18, {0, 1}]
Out[6]= 1.00000000000000003