def test_find_intersection_array(self):
# SA == SAcap at index 1.5, therefore SD=2.5
SD = array([[1, 2, 3.0], [1, 2, 3.0], [1, 2, 3.0]])
SA = array([[2, 4, 6.0], [2, 4, 6.0], [2, 4, 6.0]])
SAcap = array([[0, 3, 7.0], [0, 3, 4.0], [0, 6, 7.0]])
answer = [2.5, 3, 1.5]
assert allclose(find_intersection(SD, SA, SAcap), answer)
def test_find_intersection_array(self):
# SA == SAcap at index 1.5, therefore SD=2.5
SD = array([[1, 2, 3.0], [1, 2, 3.0], [1, 2, 3.0]])
SA = array([[2, 4, 6.0], [2, 4, 6.0], [2, 4, 6.0]])
SAcap = array([[0, 3, 7.0], [0, 3, 4.0], [0, 6, 7.0]])
answer = [2.5, 3, 1.5]
assert allclose(find_intersection(SD, SA, SAcap), answer)
def test_find_intersection(self):
# SA == SAcap at index 1.5, therefore SD=2.5
SD = array([1, 2, 3.0]) # x-axis
SA = array([2, 4, 6.0]) # y-axis
SAcap = array([0, 3, 7.0]) # y-axis
assert (find_intersection(SD, SA, SAcap) == 2.5)
# SA exceeds SAcap the whole way, so SD = max(SD)
SD = array([1, 2, 3.0])
SA = array([2, 4, 6.0])
SAcap = array([0, 3, 4.0])
assert (find_intersection(SD, SA, SAcap) == 3.0)
# SA is exceeded by SAcap in the first bin
SD = array([1, 2, 3.0])
SA = array([2, 4, 6.0])
SAcap = array([0, 6, 7.0])
assert (find_intersection(SD, SA, SAcap) == 1.5)
def test_find_intersection(self):
# SA == SAcap at index 1.5, therefore SD=2.5
SD = array([1, 2, 3.0]) # x-axis
SA = array([2, 4, 6.0]) # y-axis
SAcap = array([0, 3, 7.0]) # y-axis
assert (find_intersection(SD, SA, SAcap) == 2.5)
# SA exceeds SAcap the whole way, so SD = max(SD)
SD = array([1, 2, 3.0])
SA = array([2, 4, 6.0])
SAcap = array([0, 3, 4.0])
assert (find_intersection(SD, SA, SAcap) == 3.0)
# SA is exceeded by SAcap in the first bin
SD = array([1, 2, 3.0])
SA = array([2, 4, 6.0])
SAcap = array([0, 6, 7.0])
assert (find_intersection(SD, SA, SAcap) == 1.5)
def test_find_intersection_again(self):
'''
[SDcrNew, SAcrNew] = find_intersection( SDnew, ...
SAnew, ...
SAcapNew, ...
BCAP_PARAMS_T.('DyV'), ...
BCAP_PARAMS_T.('AyV'), ...
BCAP_PARAMS_T.('DuV'), ...
BCAP_PARAMS_T.('AuV'), ...
BCAP_PARAMS_T.('aaV'), ...
BCAP_PARAMS_T.('bbV'), ...
BCAP_PARAMS_T.('ccV') );'''
SAnew = array([0.289001251823406, 0.645044680932497,
0.552493097099059, 0.448383584870071,
0.378214709055794, 0.334886136749447,
0.303739197972188, 0.283211327986589,
0.267078826474521, 0.243891573578266,
0.21891385215188, 0.140684880819053,
0.0970201359460657, 0.0706874801683211,
0.0485714053533709, 0.0350210210488472,
0.0256594254042202, 0.0188746569449423,
0.0138838913645612])
SDnew = array([0, 1.6012, 5.486, 10.017, 15.022, 20.783, 27.144, 34.449,
42.431, 49.04, 54.342, 78.577, 96.336, 109.67, 108.51,
106.5, 101.91, 94.879, 86.162])
SAcap = array([0, 0.0716716312147219, 0.209842038805638,
0.25543506124516, 0.265904181516395,
0.26797788578358, 0.26829075866169,
0.268329611544255, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333])
SDcr = find_intersection(SDnew, SAnew, SAcap, axis=-1)
SDcr_m = 41.810644943005
assert allclose(SDcr_m, SDcr)
def test_find_intersection_again(self):
'''
[SDcrNew, SAcrNew] = find_intersection( SDnew, ...
SAnew, ...
SAcapNew, ...
BCAP_PARAMS_T.('DyV'), ...
BCAP_PARAMS_T.('AyV'), ...
BCAP_PARAMS_T.('DuV'), ...
BCAP_PARAMS_T.('AuV'), ...
BCAP_PARAMS_T.('aaV'), ...
BCAP_PARAMS_T.('bbV'), ...
BCAP_PARAMS_T.('ccV') );'''
SAnew = array([0.289001251823406, 0.645044680932497,
0.552493097099059, 0.448383584870071,
0.378214709055794, 0.334886136749447,
0.303739197972188, 0.283211327986589,
0.267078826474521, 0.243891573578266,
0.21891385215188, 0.140684880819053,
0.0970201359460657, 0.0706874801683211,
0.0485714053533709, 0.0350210210488472,
0.0256594254042202, 0.0188746569449423,
0.0138838913645612])
SDnew = array([0, 1.6012, 5.486, 10.017, 15.022, 20.783, 27.144, 34.449,
42.431, 49.04, 54.342, 78.577, 96.336, 109.67, 108.51,
106.5, 101.91, 94.879, 86.162])
SAcap = array([0, 0.0716716312147219, 0.209842038805638,
0.25543506124516, 0.265904181516395,
0.26797788578358, 0.26829075866169,
0.268329611544255, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333])
SDcr = find_intersection(SDnew, SAnew, SAcap, axis=-1)
SDcr_m = 41.810644943005
assert allclose(SDcr_m, SDcr)
def test_find_intersection(self):
SA = array([0.289739794645192, 0.646693092888598, 0.553904993443363,
0.449529429311655, 0.379120948142093, 0.335578413347301,
0.304315540473696, 0.283748719032318, 0.267585606238139,
0.24435435573003, 0.219329239293237, 0.140951829163807,
0.0972042308150227, 0.0708216090506566, 0.0486635691749981,
0.0350874731334362, 0.0257081139420298, 0.018910471443273,
0.0139102359283618])
SD = array([0,1.60533088580753, 5.49998633698612, 10.0430812761127,
15.05788941807, 20.8257362526643, 27.1952446916894,
34.5140821612566, 42.5117576683308, 49.1327918769187,
54.4456104247827, 78.726251497955, 96.5187076679668,
109.878354239657, 108.720841961029, 106.697568070816,
102.107240125687, 95.0591053452788, 86.3258461016045])
SAcap = array([0, 0.0718547880987331, 0.210114969058595,
0.255544749303803, 0.26593317666683,
0.267982944023487, 0.268291483967979,
0.268329691827397, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333])
'''[SDcrNew, SAcrNew] = find_intersection( SDnew, ...
SAnew, ...
SAcapNew, ...
BCAP_PARAMS_T.('DyV'), ...
BCAP_PARAMS_T.('AyV'), ...
BCAP_PARAMS_T.('DuV'), ...
BCAP_PARAMS_T.('AuV'), ...
BCAP_PARAMS_T.('aaV'), ...
BCAP_PARAMS_T.('bbV'), ...
BCAP_PARAMS_T.('ccV') );'''
SDcr = find_intersection(SD, SA, SAcap, axis=-1)
SDcr_m = 42.1418578996947
assert allclose(SDcr_m, SDcr)
def test_find_intersection(self):
SA = array([0.289739794645192, 0.646693092888598, 0.553904993443363,
0.449529429311655, 0.379120948142093, 0.335578413347301,
0.304315540473696, 0.283748719032318, 0.267585606238139,
0.24435435573003, 0.219329239293237, 0.140951829163807,
0.0972042308150227, 0.0708216090506566, 0.0486635691749981,
0.0350874731334362, 0.0257081139420298, 0.018910471443273,
0.0139102359283618])
SD = array([0,1.60533088580753, 5.49998633698612, 10.0430812761127,
15.05788941807, 20.8257362526643, 27.1952446916894,
34.5140821612566, 42.5117576683308, 49.1327918769187,
54.4456104247827, 78.726251497955, 96.5187076679668,
109.878354239657, 108.720841961029, 106.697568070816,
102.107240125687, 95.0591053452788, 86.3258461016045])
SAcap = array([0, 0.0718547880987331, 0.210114969058595,
0.255544749303803, 0.26593317666683,
0.267982944023487, 0.268291483967979,
0.268329691827397, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333,
0.268333333333333, 0.268333333333333])
'''[SDcrNew, SAcrNew] = find_intersection( SDnew, ...
SAnew, ...
SAcapNew, ...
BCAP_PARAMS_T.('DyV'), ...
BCAP_PARAMS_T.('AyV'), ...
BCAP_PARAMS_T.('DuV'), ...
BCAP_PARAMS_T.('AuV'), ...
BCAP_PARAMS_T.('aaV'), ...
BCAP_PARAMS_T.('bbV'), ...
BCAP_PARAMS_T.('ccV') );'''
SDcr = find_intersection(SD, SA, SAcap, axis=-1)
SDcr_m = 42.1418578996947
assert allclose(SDcr_m, SDcr)
def test_find_intersectionII(self):
# If there is 2 intersections, are they found? No
SD = array([1, 2, 3.0])
SA = array([2, 4, 6.0])
SAcap = array([0, 6, 4.0])
assert (find_intersection(SD, SA, SAcap) == 1.5)
def test_find_intersectionII(self):
# If there is 2 intersections, are they found? No
SD = array([1, 2, 3.0])
SA = array([2, 4, 6.0])
SAcap = array([0, 6, 4.0])
assert (find_intersection(SD, SA, SAcap) == 1.5)
