The R Project SVN R

Rev

Rev 76131 | Rev 81517 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 76131 Rev 76206
Line 380... Line 380...
380 
	}
 380 
	}
381 
    }
 381 
    }
382 
    for (i = 0; i < n; i++) q[i] += i;
 382 
    for (i = 0; i < n; i++) q[i] += i;
383 
 
 383 
 
384 
    /* generate sample */
 384 
    /* generate sample */
- 
 
 385 
    Sampletype Sample_kind = R_sample_kind();
385 
    for (i = 0; i < nans; i++) {
 386 
    for (i = 0; i < nans; i++) {
- 
 
 387 
	if (Sample_kind == ROUNDING) {
386 
	rU = R_unif_index(n);
 388 
	    rU = unif_rand() * n;
387 
	k = (int) rU;
 389 
	    k = (int) rU;
- 
 
 390 
	}
- 
 
 391 
	else {
- 
 
 392 
	    k = (int) R_unif_index(n);
- 
 
 393 
	    rU = k + unif_rand();
- 
 
 394 
	}
388 
	ans[i] = (rU < q[k]) ? k+1 : a[k]+1;
 395 
	ans[i] = (rU < q[k]) ? k+1 : a[k]+1;
389 
    }
 396 
    }
390 
    if(n > SMALL) {
 397 
    if(n > SMALL) {
391 
	Free(HL);
 398 
	Free(HL);
392 
	Free(q);
 399 
	Free(q);