Noise not working in JWST experiment
I created a demo to run CDRl on JWST with noise. I got the following error:
[nix-shell:~/work/git/ProxPython/demos]$ python3 JWST_CDRl.py
Loading data.
Traceback (most recent call last):
File "JWST_CDRl.py", line 6, in JWST = JWST_Experiment(algorithm = 'CDRl', lambda_0 = 0.5, lambda_max = 0.5, noise = 'Poisson', MAXIT = 1000)
File "/home/rluke1/work/git/ProxPython/demos/../proxtoolbox/experiments/experiment.py", line 39, in call obj.initialize() # initialize experiment object
File "/home/rluke1/work/git/ProxPython/demos/../proxtoolbox/experiments/experiment.py", line 223, in initialize self.loadData()
File "/home/rluke1/work/git/ProxPython/demos/../proxtoolbox/experiments/phase/JWST_Experiment.py", line 150, in
loadData k_j[ii,jj]= np.random.poisson(k_j[ii,jj])*snr #use built in numpy possion instead of
utils.PoissonRan(k[ii,jj,i])*snr File "mtrand.pyx", line 3556, in numpy.random.mtrand.RandomState.poisson File "_common.pyx", line 815, in numpy.random._common.disc File "_common.pyx", line 379, in numpy.random._common.check_constraint
ValueError: lam value too large