Den 29.03.2011 14:56, skrev Sturla Molden:
> import numpy as np
> def randombits(n, p):
> b = (np.random.rand(n*8).reshape((n,8))< p).astype(int)
> return (b<< range(8)).sum(axis=1).astype(np.uint8)
n is the number of bytes. If you prefer to count in bits:
def randombits(n, p):
assert(n%8 == 0)
b = (np.random.rand(n).reshape((n>>3,8))< p).astype(int)
return (b<< range(8)).sum(axis=1).astype(np.uint8)
Sturla
_______________________________________________
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion
