def test_figure3b_a3_equals_point_five(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0]) * 20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0, A3=4.53514739e-06, H=20, rs=0.175, ks=2e-8 / 1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[0.81168093, 0.13560054], [0.8355792, 0.18756625], [0.86669814, 0.27523356], [0.90708096, 0.42885195], [0.92011769, 0.48939465], [0.92219713, 0.49958604]]))
def test_figure3a_A2_equals_point_nine_nine(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0]) * 20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0.99, A3=0, H=20, rs=0.175, ks=2e-8 / 1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[0.81168066, 0.12412292], [0.83668852, 0.16812724], [0.87096244, 0.25118534], [0.92073203, 0.43785721], [0.94257216, 0.55353625], [0.95050775, 0.60194473]]))
def test_figure3b_a3_equals_point_one(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0]) * 20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0, A3=9.07029478e-07, H=20, rs=0.175, ks=2e-8 / 1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[0.81110629, 0.12524912], [0.83472394, 0.16807639], [0.86563054, 0.24288244], [0.90599734, 0.38274085], [0.91909172, 0.44098945], [0.92118321, 0.45096744]]))
def test_figure3c_A2_equals_point_nin_nine_a3_equals_one(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0]) * 20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0.99, A3=9.07029478e-06, H=20, rs=0.175, ks=2e-8 / 1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[0.81315577, 0.15423227], [0.83889148, 0.22420948], [0.87368554, 0.34204751], [0.92325897, 0.55703101], [0.9446685, 0.66714553], [0.95239683, 0.70945769]]))
def test_figure3a_A2_equals_point_one(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0]) * 20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0.1, A3=0, H=20, rs=0.175, ks=2e-8 / 1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[0.81103472, 0.12313867], [0.83472101, 0.16421523], [0.86588137, 0.23691031], [0.90695675, 0.3765874], [0.92048265, 0.43667275], [0.92267961, 0.44720765]]))
def test_figure3b_a3_equals_point_five(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0])*20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0, A3=4.53514739e-06, H=20, rs=0.175, ks=2e-8/1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[ 0.81168093, 0.13560054], [ 0.8355792 , 0.18756625], [ 0.86669814, 0.27523356], [ 0.90708096, 0.42885195], [ 0.92011769, 0.48939465], [ 0.92219713, 0.49958604]]))
def test_figure3b_a3_equals_one(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0]) * 20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0, A3=9.07029478e-06, H=20, rs=0.175, ks=2e-8 / 1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[0.81241413, 0.15242832], [0.83666321, 0.21757866], [0.86803975, 0.32108789], [0.90842794, 0.4868681], [0.92138857, 0.54790671], [0.92345241, 0.55797917]]))
def test_figure3b_a3_equals_point_one(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0])*20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0, A3=9.07029478e-07, H=20, rs=0.175, ks=2e-8/1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[ 0.81110629, 0.12524912], [ 0.83472394, 0.16807639], [ 0.86563054, 0.24288244], [ 0.90599734, 0.38274085], [ 0.91909172, 0.44098945], [ 0.92118321, 0.45096744]]))
def test_figure3a_A2_equals_point_nine_nine(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0])*20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0.99, A3=0, H=20, rs=0.175, ks=2e-8/1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[ 0.81168066, 0.12412292], [ 0.83668852, 0.16812724], [ 0.87096244, 0.25118534], [ 0.92073203, 0.43785721], [ 0.94257216, 0.55353625], [ 0.95050775, 0.60194473]]))
def test_figure3a_A2_equals_point_one(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0])*20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0.1, A3=0, H=20, rs=0.175, ks=2e-8/1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[ 0.81103472, 0.12313867], [ 0.83472101, 0.16421523], [ 0.86588137, 0.23691031], [ 0.90695675, 0.3765874 ], [ 0.92048265, 0.43667275], [ 0.92267961, 0.44720765]]))
def test_figure3c_A2_equals_point_nin_nine_a3_equals_one(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0])*20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0.99, A3=9.07029478e-06, H=20, rs=0.175, ks=2e-8/1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[ 0.81315577, 0.15423227], [ 0.83889148, 0.22420948], [ 0.87368554, 0.34204751], [ 0.92325897, 0.55703101], [ 0.9446685 , 0.66714553], [ 0.95239683, 0.70945769]]))
def test_figure3b_a3_equals_one(self): assert_allclose( dengetal2013(z=np.array([0.05, 0.1, 0.2, 0.5, 0.8, 1.0])*20, t=[11025., 110250.], rw=0.035, re=0.525, A1=1, A2=0, A3=9.07029478e-06, H=20, rs=0.175, ks=2e-8/1.8, kw0=1e-3, kh=2e-8, mv=0.2e-3, gamw=10, ui=1), np.array([[ 0.81241413, 0.15242832], [ 0.83666321, 0.21757866], [ 0.86803975, 0.32108789], [ 0.90842794, 0.4868681 ], [ 0.92138857, 0.54790671], [ 0.92345241, 0.55797917]]))