rrcpplinspace

Rcpp code for "linspace" function: vector length incorrect when increment 1e-5


I want to write a function like "linspace" to create equal interval vectors in R with Rcpp. num_dis here is the number of intervals I want for that vector, i.e., a=0,b=10,num_dis = (10-0)/(1e-5)+1=1000001. When I put linspace(0,10,1000001), the vector length is correct as 1000001; however, when I write like: linspace(0,10,10/(1e-5)+1), it returns a vector of length 100000. For increments like 1e-3,1e-6, etc., it seems good. I am not quite sure what happens?

 NumericVector linspace(double a,double b,int num_dis) {
        NumericVector u(num_dis);
        for (int i = 0; i < num_dis; i++) {
            u[i] = a + i * ((b - a) / (num_dis-1));
        }
        return u;
    }

Solution

  • I think this explains things:

    > as.integer(10/(1e-5)+1)
    [1] 1000000
    > as.integer(10L*1e5L+1L)
    [1] 1000001
    >