Re: precision of floor(a/b)

[Martin Helm]

Am Dienstag, 6. April 2010 12:57:19 schrieb CdeMills:
> To avoid implementation-dependent problems, I suggest to compute your loop
> as
> alpha_min = 0;
> alpha_max = 6;
> alpha_step = 0.1;
> for alpha = alpha_min : alpha_max,
>         # <do some stuff with alpha*alpha_step>
> endfor
> This way, the latest processed value will be 6*0.1, and not either .5 or
>  .6, depending on the way the equality test is performed. It's less
>  elegant, requires one more multiplication or variable, but it's
>  predictable.
> 
> Regards
> 
> Pascal
> 
why not simply

n = number of intervals

for alpha = linspace(alpha_min, alpha_max, n+1),
.. work with alpha
endfor

- mh