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)