It probably would be useful to have a linspace-ish interface for 
constructing FloatRanges: I wrote a quick one for KernelDensity.jl:
https://github.com/JuliaStats/KernelDensity.jl/blob/master/src/univariate.jl#L60
but there's probably a better way to do this.


On Thursday, 24 April 2014 19:05:51 UTC+1, Hans W Borchers wrote:
>
> There is a strange difference in results between Matlab and Julia that 
> could be the reason why  [0:pi/100:pi]  ends up differently:
>
>     matlab> x = ones(100, 1) * (pi/100);
>     matlab> sum(x) - pi
>     ans =
>       -4.4409e-15
>
>
>     julia> x = ones(100)*(pi/100);
>     julia> sum(x) - pi
>     1.3322676295501878e-15
>
>
> So it appears that summing pi/100 100-times is smaller than pi in Matlab 
> and larger than pi in Julia.
>
>
> On Thursday, April 24, 2014 5:47:34 PM UTC+2, Stefan Karpinski wrote:
>>
>> On Thu, Apr 24, 2014 at 11:43 AM, andrew cooke <[email protected]> wrote:
>>
>>>
>>> sorry, rational is a stupid suggestion for irrational numbers...
>>>
>>> what i am trying to say is that there is no perfect solution to this 
>>> problem (that i know of).  if matlab works better then it is either luck or 
>>> because they are actually fudging things.  fudging things is appealing but 
>>> usually bites you somewhere down the line.
>>>
>>
>> This is exactly how I feel about it. I still think it may be possible to 
>> improve our FloatRange behavior even further, but I'd like to avoid 
>> "fudging it".
>>
>

Reply via email to