2013-12-14 07:43:55 +00:00
|
|
|
(* vdc.m -- van der Corput script *)
|
|
|
|
|
(* Copyright (C) 2013 SheetJS *)
|
|
|
|
|
(* vim: set ft=mma: *)
|
|
|
|
|
(* Usage: math -script vdc.m <dispdigits> <base> <nelts> *)
|
|
|
|
|
|
|
|
|
|
(* read arguments from command line *)
|
|
|
|
|
argdigs=ToExpression[$CommandLine[[-3]]];
|
|
|
|
|
argbase=ToExpression[$CommandLine[[-2]]];
|
|
|
|
|
argelts=ToExpression[$CommandLine[[-1]]];
|
|
|
|
|
d={argdigs,argdigs};
|
|
|
|
|
|
|
|
|
|
(*
|
|
|
|
|
van der Corput sequence generation in mathematica:
|
|
|
|
|
http://reference.wolfram.com/mathematica/CUDALink/tutorial/Applications.html
|
|
|
|
|
|
|
|
|
|
this function returns a table of the form
|
2016-01-16 18:47:42 +00:00
|
|
|
[iteration, fraction]
|
2013-12-14 07:43:55 +00:00
|
|
|
*)
|
|
|
|
|
VanDerCorput[base_][len_]:=Table[
|
2016-01-16 18:47:42 +00:00
|
|
|
{n, With[{digits = Reverse@IntegerDigits[n, base]},
|
2013-12-14 07:43:55 +00:00
|
|
|
N[Sum[argbase^(-ii)*digits[[ii]], {ii, Length[digits]}],d]]}, {n, len}];
|
|
|
|
|
|
|
|
|
|
(* print table form *)
|
|
|
|
|
data=VanDerCorput[argbase][argelts];
|
|
|
|
|
Print[Grid[data]//OutputForm]
|
