def test_noinfo2(self) : """ Run ALR eqn 1 under noinfo with tol 10e-7 """ ss = Chou2006() ss.exptype = "noinfo" ss.equations = [1] ar = ALRSolver(ss) ar.solve(l1penalty=0.0,maxiter=10000,tol=10e-7) assert len(ar.all_exp_art) == 1 assert len(ar.all_exp_art[0]['eqns']) == 1 art = ar.all_exp_art[0]['eqns'][0] params = art.params[-1] assert params is None
def test_noinfo2(self):
""" Run ALR eqn 1 under noinfo with tol 10e-7 """
ss = Chou2006()
ss.exptype = "noinfo"
ss.equations = [1]
ar = ALRSolver(ss)
ar.solve(l1penalty=0.0, maxiter=10000, tol=10e-7)
assert len(ar.all_exp_art) == 1
assert len(ar.all_exp_art[0]['eqns']) == 1
art = ar.all_exp_art[0]['eqns'][0]
params = art.params[-1]
assert params is None
def test_noinfo1(self) :
""" Run ALR eqn 1 under no info with tol 10e-6 """
ss = Chou2006()
ss.exptype = "noinfo"
ss.equations = [1]
ar = ALRSolver(ss)
ar.solve(l1penalty=0.0,maxiter=10000,tol=10e-6)
assert len(ar.all_exp_art) == 1
assert len(ar.all_exp_art[0]['eqns']) == 1
art = ar.all_exp_art[0]['eqns'][0]
sol = {'alpha': 12,
'beta': 10,
'g': np.array([ 0,0,-0.8,0]),
'h': np.array([ 0.5,0,0,0])}
params = art.params[-1]
assert np.allclose(sol['alpha'],params['alpha'],atol=1.0)
assert np.allclose(sol['beta'],params['beta'],atol=1.0)
assert np.allclose(sol['g'],params['g'],atol=0.5)
assert np.allclose(sol['h'],params['h'],atol=0.5)
def test_fullinfo1(self) :
""" RUN ALR (slow) eqn 1 under full info and alpha=0.0"""
ss = Chou2006()
ss.exptype = "fullinfo"
ss.equations = [1]
ar = ALRSolver(ss)
ar.solve(l1penalty=0.0,maxiter=10000,tol=10e-6)
assert len(ar.all_exp_art) == 1
assert len(ar.all_exp_art[0]['eqns']) == 1
art = ar.all_exp_art[0]['eqns'][0]
sol = {'alpha': 11.177018643003306,
'beta': 8.7027543289527056,
'g': np.array([-1.0236953]),
'h': np.array([ 0.62268492])}
params = art.params[-1]
assert np.allclose(sol['alpha'],params['alpha'],atol=1)
assert np.allclose(sol['beta'],params['beta'],atol=1.5)
assert np.allclose(sol['g'],params['g'],atol=0.5)
assert np.allclose(sol['h'],params['h'],atol=0.5)
def test_noinfo1(self):
""" Run ALR eqn 1 under no info with tol 10e-6 """
ss = Chou2006()
ss.exptype = "noinfo"
ss.equations = [1]
ar = ALRSolver(ss)
ar.solve(l1penalty=0.0, maxiter=10000, tol=10e-6)
assert len(ar.all_exp_art) == 1
assert len(ar.all_exp_art[0]['eqns']) == 1
art = ar.all_exp_art[0]['eqns'][0]
sol = {
'alpha': 12,
'beta': 10,
'g': np.array([0, 0, -0.8, 0]),
'h': np.array([0.5, 0, 0, 0])
}
params = art.params[-1]
assert np.allclose(sol['alpha'], params['alpha'], atol=1.0)
assert np.allclose(sol['beta'], params['beta'], atol=1.0)
assert np.allclose(sol['g'], params['g'], atol=0.5)
assert np.allclose(sol['h'], params['h'], atol=0.5)
def test_fullinfo1(self):
""" RUN ALR (slow) eqn 1 under full info and alpha=0.0"""
ss = Chou2006()
ss.exptype = "fullinfo"
ss.equations = [1]
ar = ALRSolver(ss)
ar.solve(l1penalty=0.0, maxiter=10000, tol=10e-6)
assert len(ar.all_exp_art) == 1
assert len(ar.all_exp_art[0]['eqns']) == 1
art = ar.all_exp_art[0]['eqns'][0]
sol = {
'alpha': 11.177018643003306,
'beta': 8.7027543289527056,
'g': np.array([-1.0236953]),
'h': np.array([0.62268492])
}
params = art.params[-1]
assert np.allclose(sol['alpha'], params['alpha'], atol=1)
assert np.allclose(sol['beta'], params['beta'], atol=1.5)
assert np.allclose(sol['g'], params['g'], atol=0.5)
assert np.allclose(sol['h'], params['h'], atol=0.5)