feat(unitproc/bkt-get_path_fork_level):Add Bash function
[BK-2020-03.git] / unitproc / python / bktemp-subsetExpandReduce.py
1 #!/usr/bin/env python3
2
3 # Desc: Expand or reduce list of int/floats to specified length.
4
5 import numpy as np
6
7 def medianSubset(listIn: list = [], listOutLen: int = 0) -> list:
8 # Input: int: listOutLen: quantity of elements in output list
9 # list: listIn: input list consisting of integers or floats
10 # Output: list: ints/floats of specified size
11 # Ref/Attrib: PEP 3107 typing https://stackoverflow.com/a/21384492
12 # Version: 0.1.0
13 #print('DEBUG:listOutLen:' + str(listOutLen))
14 #print('DEBUG:listIn:' + str(listIn))
15
16 # Exit for invalid input
17 if not isinstance(listOutLen, int):
18 raise ValueError('ERROR:Not a valid int:' + str(listOutLen))
19 else:
20 if not listOutLen > 0:
21 raise ValueError('ERROR:Invalid value:' + str(listOutLen))
22 if not isinstance(listIn, list):
23 raise ValueError('ERROR:Not a valid list:' + str(listOutLen))
24 if not all([( (isinstance(x,int)) or (isinstance(x,float)) ) for x in listIn]):
25 raise ValueError('ERROR:Input list contains something besides integers or floating point numbers.')
26
27 # Initialize listOut
28 listOut = [None] * listOutLen
29 #print('DEBUG:listOut:' + str(listOut))
30
31 # Calc listIn length
32 listInLen = len(listIn)
33 #print('DEBUG:listInLen:' + str(listInLen))
34
35 # Calc subset length float
36 subsetLenFloat = ( (max([listInLen,listOutLen]) - 1) /min([listInLen,listOutLen]))
37 subsetIndRatio = ( (listInLen)/(listOutLen) )
38 #print('DEBUG:subsetLenFloat: %.5f' % subsetLenFloat)
39 #print('DEBUG:subsetLenFloat2: %.5f' % subsetIndRatio)
40
41 # Iterate for each element in listOut
42 for i_out in range(listOutLen):
43 #print('DEBUG:i_out:' + str(i_out))
44 ## Decide to expand or reduce listIn to produce listOut
45 if listInLen > listOutLen:
46 ### reduce listIn to listOut
47 #print('DEBUG:listOutLen:' + str(listOutLen))
48 #print('DEBUG:listInLen:' + str(listInLen))
49 if i_out == 0:
50 #### Initialize subsetIndLo in first loop
51 subsetIndLo = int(0)
52 #print('DEBUG:subsetIndLo:' + str(subsetIndLo))
53 #print('DEBUG:i_out:' + str(i_out))
54 #### Calc indices of i_out'th subset of listIn
55 subsetIndHi = (listInLen - 1) * (i_out + 1) // listOutLen
56 subsetLen = subsetIndHi - subsetIndLo + 1
57 #print('DEBUG:subsetIndLo:' + str(subsetIndLo))
58 #print('DEBUG:subsetIndHi:' + str(subsetIndHi))
59 #print('DEBUG:subsetLen:' + str(subsetLen))
60 #### Extract subset from listIn using indices inclusively
61 subset = listIn[ int(subsetIndLo) : int(subsetIndHi)+1 ]
62 #print('DEBUG:subset:' + str(subset))
63 #### Calculate median for subset
64 subsetMedian = np.median(subset)
65 #print('DEBUG:subset median:' + str(subsetMedian))
66 #### Set listOut element
67 listOut[i_out] = subsetMedian
68 #### Housekeeping
69 ##### Update subsetIndLo for next loop
70 subsetIndLo = subsetIndHi + 1
71 #print('DEBUG:Updated subsetIndLo:' + str(subsetIndLo))
72 elif listOutLen > listInLen:
73 ### Expand listIn to listOut
74 #print('DEBUG:listOutLen:' + str(listOutLen))
75 #print('DEBUG:listInLen:' + str(listInLen))
76 #### Identify index list of lists mapping listIn to ListOut
77 expandIndex = int(i_out / subsetLenFloat)
78 expandIndex = min([expandIndex,(listInLen - 1)])
79 #print('DEBUG:expandIndex:' + str(expandIndex))
80 listOut[i_out] = listIn[expandIndex]
81 #print('DEBUG:listOut[i_out]:' + str(listOut[i_out]))
82 elif listOutLen == listInLen:
83 listOut = listIn
84 #print('DEBUG:end for loop===========')
85 return listOut
86
87 example_list1 = [738,253,104,543,501,352,770,945,992,920,977,435,783,592,261,176,888]
88 example_list2 = [641,679,371]
89 example_list3 = [738,253,104,543,501,352,770,945,992,920,977,435,783,592,261,176,888,767,169,495,382,291,116,480,8,52,761,455,743,592,54,490,596,8,816,134,641,705,2,675,414,511,220,60,116,563,19,316,605,735,196,210,107,848,630,392,700,94,824,329,307,693,374,511,17,488,626,730,606,273,875,466,790,40,681,954,947,981,775,498,49,327,268,492,799,273,403,839,895,200,329,373,819,302,456,23,944,543,997,605,338,919,385,783,339,825,457,761,994,530,862,504,690,720,10,102,662,579,835,747,337,471,409,904,321,928,908,573,984,853,244,499,940,468,480,668,128,570,584,325,940,230,738,624,672,601,384,242,800,172,125,150,124,958,242,560,829,22,614,707,677,413,800,536,269,748,209,821,212,939,196,2,830,815,235,847,456,181,920,413,153,935,626,215,821,813,847,875,550,656,237,760,718,179,549,63,839,204,631,917,95,733,341,881,539,561,227,492,549,423,565,856,384,376,809,920,116,609,205,644,791,187,15,55,776,376,882,905,937,304,626,312,263,16,163,651,422,656,116,495,158,156,981,167,786,536,987,697,88,354,335,281,718,562,169,783,483,16,314,344,760,189,965,494,68,314,806,26,866,741,283,302,761,518,713,264,570,239,410,629,47,547,270,94,228,303,811,795,658,102,96,438,157,62,826,27,382,847,341,953,671,147,746,408,201,536,809,235,778,41,266,682,686,620,279,167,183,533,471,173,145,676,774,191,885,583,694,1000,484,115,159,31,704,248,852,441,901,705,829,496,784,163,88,1,503,986,929,631,839,124,55,136,907,342,863,328,417,674,881,567,876,30,812,515,606,469,725,478,100,149,881,12,221,839,138,354,712,550,698,931,914,610,704,513,477,527,388,341,373,77,673,970,612,384,613,321,411,382,273,100,55,760,246,161,268,348,22,533,851,515,807,76,396,254,554,503,322,292,821,345,147,303,225,389,356,119,530,502,407,474,824,900,561,624,292,893,597,60,393,856,675,27,257,535,5,879,809,944,904,327,209,626,440,102,528,970,839,77,0,987,422,414,337,499,266,640,194,454,414,517,221,83,4,473,84,39,831,640,557,685,45,91,410,327,734,470,928,22,837,583,570,671,773,94,317,935,365,319,965,465,723,471,77,420,683,581,604,17,373,767,722,15,742,533,497,847,8,84,842,28,526,591,466,396,807,676,564,695,967,47,886,266,638,767,58,763,775,778,25,699,895,64,616,11,369,521,856,277,726,165,281,692,736,256,545,912,954,325,620,220,407,304,191,363,884,818,138,783,559,164,516,880,791,884,246,963,522,547,356,475,302,582,687,399,862,109,524,570,867,961,780,611,491,677,977,879,803,651,41,400,634,415,750,129,89,861,265,669,833,833,333,221,102,505,353,264,31,782,595,838,365,141,475,697,46,503,71,205,164,811,46,57,278,278,12,297,33,326,106,916,352,528,616,21,129,534,485,753,410,952,758,26,152,123,578,42,33,881,145,71,662,902,782,634,256,598,265,953,135,778,407,990,643,309,448,381,620,559,208,648,830,332,59,420,456,558,439,503,182,728,486,568,986,172,652,176,380,542,92,551,335,602,854,485,543,752,860,163,120,607,207,710,298,711,32,908,228,41,481,313,581,874,443,273,204,265,567,792,326,173,29,762,398,923,18,312,144,965,764,342,884,395,624,853,709,7,424,11,836,677,210,717,950,525,782,39,770,972,993,517,495,323,772,133,941,300,809,517,574,243,664,772,721,232,304,429,864,135,667,871,371,894,485,86,42,445,426,556,588,205,533,392,900,962,112,837,961,617,302,779,793,969,199,199,591,747,321,453,861,145,451,149,887,726,414,696,421,260,92,159,56,528,879,844,408,568,603,182,492,796,656,876,566,468,237,400,629,366,868,145,150,151,561,738,241,891,813,235,400,476,224,241,169,157,969,540,901,925,723,536,421,369,24,560,882,817,677,234,16,241,464,923,849,300,144,261,916,529,586,905,268,950,673,697,664,840,725,926,603,96,845,410,321,9,996,844,952,380,272,125,140,719,402,383,499,9,48,955,429,458,965,584,381,739,984,10,282,165,486,33,353,616,304,215,293,362,148,542,815,659,874,890,502,369,499,533,144,876,169,375,951,784,320,782,763,473,38,871,790,933,522,88,702,387,705,544,575,719,757,867,470,807,628,914,719,764,536,748,766,419,42,604,678,742,321,280,147,975,918,982,440,842,723,910,893,111,560,922,837,116,858,374,262,878,354,397,638,338,853]
90 example_list4 = [641,679,371,3.141592]
91 example_list5 = [641,679,371,3.141592, 'jan']
92 #print(medianSubset(example_list,3))
93 #print(medianSubset(example_list2,17))
94 #print(medianSubset(example_list3,3))
95 #print(medianSubset(example_list3,1))
96 #print(medianSubset(example_list4,1))
97 print(medianSubset(example_list3,47))