Leonid Keselman leonidk
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| import numpy as np | |
| class QuasiRandom(): | |
| def __init__(self,dim=1,seed=None): | |
| self.dim = dim | |
| self.x = np.random.rand(dim) if seed is None else seed | |
| root_sys = [1] +[0 for i in range(dim-1)] + [-1,-1] | |
| self.const = sorted(np.roots(root_sys))[-1].real | |
| self.phi = np.array([1/(self.const)**(i+1) for i in range(dim)]) | |
| def generate(self,n_points=1): |
OlderNewer