fork download
copy 
  1. def Halton_sequence(i0):
  2. 	xbase = 2
  3. 	ybase = 3
  4.  
  5. 	invxbase = 1.0 / xbase
  6. 	facx = 1.0 / xbase
  7.  
  8. 	invybase = 1.0 / ybase
  9. 	facy = 1.0 / ybase
  10.  
  11. 	inp = i0
  12. 	x0 = 0.0
  13. 	while inp > 0:
  14. 		x0 = x0 + (inp % xbase) * invxbase
  15. 		inp = inp / xbase
  16. 		invxbase = invxbase * facx
  17.  
  18. 	inp = i0
  19. 	y0 = 0.0
  20. 	while inp > 0:
  21. 		y0 = y0 + (inp % ybase) * invybase
  22. 		inp = inp / ybase
  23. 		invybase = invybase * facy
  24.  
  25. 	return x0, y0
  26.  
  27. for i0 in range(10):
  28. 	x0,y0 = Halton_sequence(i0)
  29. 	print x0,y0
  30.  
  31.  
  32.  
Success #stdin #stdout 0.01s 9024KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
0.0 0.0
0.5 0.333333333333
0.25 0.666666666667
0.75 0.111111111111
0.125 0.444444444444
0.625 0.777777777778
0.375 0.222222222222
0.875 0.555555555556
0.0625 0.888888888889
0.5625 0.037037037037
https://ideone.com/9LC7WO
Python > Halton_sequence
language:
Python (cpython 2.7.16)
created:
7 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!