-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathSectors.py
executable file
·172 lines (154 loc) · 6.98 KB
/
Sectors.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
import numpy as np
from Configuration import Configuration
## Las matrices (TMat, EMat) almacenan los incrementos relativos al primer año
class Sector:
STypes = { 'Interpolated': "InterpolatedSector",
'Predetermined': "PredeterminedSector",
'Dependent': "DependentSector" }
def __init__(self, sectorID, sectorData):
self.name = sectorID
self.indV = sectorData["indV"]
self.constrained = [False for i in range(Configuration.nYears)]
self.e0 = sectorData['e0']
self.t0 = sectorData['t0']
if self.t0==0:
raise Exception (sectorID + ": Invalid t0 = 0.")
self.eTConstant = self.e0/self.t0
if self.indV=='t':
self.tf = sectorData['mf']
elif self.indV=='e':
self.tf = sectorData['mf'] / self.eTConstant
else:
raise Exception ("Unrecognized indV in sector " + sectorID)
self.currentStep = -1
def __getitem__(self, year):
return self.getRelYear(year)
def getAbsYear(self, year):
absE = self.e0 + self.e0*self.m[0];
absT = self.t0 + self.t0*self.m[0];
for i in range(year-Configuration.firstYear):
absE = absE + absE*self.m[i+1]
absT = absT + absT*self.m[i+1]
return [absE, absT]
def getRelYear(self, year):
i=year-Configuration.firstYear
return [self.m[i], self.m[i]]
def getData(self, ini, fin):
if (ini<Configuration.firstYear) or (ini>=Configuration.firstYear+Configuration.nYears):
raise Exception("ini("+str(ini)+") must be between "+
str(Configuration.firstYear) + " and " + str(Configuration.firstYear+Configuration.nYears-1) +".")
if (fin<Configuration.firstYear) or (fin>=Configuration.firstYear+Configuration.nYears):
raise Exception("fin("+str(fin)+") must be between "+
str(Configuration.firstYear) + " and " + str(Configuration.firstYear+Configuration.nYears-1) +".")
if (ini>fin):
raise Exception("fin("+fin+") must be greater than ini("+
str(ini) + ".")
data = np.zeros((2, fin-ini+1))
for i in range(ini, fin+1):
dat = self.getAbsYear(i)
data[0][i-Configuration.firstYear]=dat[0]
data[1][i-Configuration.firstYear]=dat[1]
return data
def getRelData(self, ini, fin):
if (ini<Configuration.firstYear) or (ini>=Configuration.firstYear+Configuration.nYears):
raise Exception("ini("+str(ini)+") must be between "+
str(Configuration.firstYear) + " and " + str(Configuration.firstYear+Configuration.nYears-1) +".")
if (fin<Configuration.firstYear) or (fin>=Configuration.firstYear+Configuration.nYears):
raise Exception("fin("+str(fin)+") must be between "+
str(Configuration.firstYear) + " and " + str(Configuration.firstYear+Configuration.nYears-1) +".")
if (ini>fin):
raise Exception("fin("+fin+") must be greater than ini("+
str(ini) + ".")
data = np.zeros(fin-ini+1)
for i in range(ini, fin+1):
dat = self.getRelYear(i)
data[i-Configuration.firstYear]=dat[0]
return data
def _detectConstraints(self, i, dat):
if hasattr(self,"maxDecr") and dat < self.maxDecr:
self.constrained[i] = True
return self.maxDecr
if hasattr(self,"maxIncr") and dat > self.maxIncr:
self.constrained[i] = True
return self.maxIncr
return dat
def step(self):
#@todo: Detect constraint violation!
self.currentStep+=1
self.m[self.currentStep] = self._detectConstraints(self.currentStep, self.m[self.currentStep])
#return self.m[self.currentStep]
class PredeterminedSector(Sector):
## data in percentages: 5(%) => +0.05
def __init__(self, sectorID, sectorData):
super().__init__(sectorID, sectorData)
IncrS = [0] + [sectorData['data'][i] for i in range(2020,2030)]
self.m = np.array(IncrS)/100
class InterpolatedSector(Sector):
def __init__(self, sectorID, sectorData):
super().__init__(sectorID, sectorData)
self.maxIncr = sectorData["maxIncr"]/100
self.maxDecr = sectorData["maxDecr"]/100
self.m = np.zeros(Configuration.nYears)
ini = Configuration.firstYear-1
data = sectorData["data"]
if ini in data:
iniData = eval("self."+sectorData["indV"]+'0')
if data[ini]!=eval("self."+sectorData["indV"]+'0'):
raise Exception(sectorID + ": Non conformant "+
sectorData["indV"]+
"0 from data ("+ str(data[ini])
+") and sector data ("+ str(iniData) +")")
if sectorData["indV"]!='e' and sectorData["indV"]!='t':
raise Exception(sectorID + ": indV is not well defined (" + sectorData["indV"] +")")
data[ini] = self.t0
keys = list(data)
if ini in keys:
keys.remove(ini)
if Configuration.lastYear not in keys:
data[Configuration.lastYear] = self.tf
keys.append(Configuration.lastYear)
i = 0
for dataEntry in sorted(keys):
fin=dataEntry
if ini>fin:
raise Exception(sectorID + ": ini > fin keys in Interpolation.")
interp = self.__linearInterp(data[ini], data[fin], int(fin)-int(ini))
for dat in interp:
year = Configuration.firstYear + i
if (year > Configuration.firstYear-1):
self.m[i] += dat
i += 1
ini=fin
def __linearInterp(self, ini, fin, steps):
interp = np.zeros(steps)
difStepAbs=(fin-ini)/steps
val = ini
if val==0:
print("__linearInterp: Divide by zero avoided (val=0.0001; ini="+ str(ini)+"; fin=" + str(fin) + ")")
val = 0.0001
for i in range(steps):
interp[i] += difStepAbs/val
val += val*interp[i]
return interp
class DependentSector(Sector):
## AssSector must be the ID of the sector
def __init__(self, sectorID, sectorData):
super().__init__(sectorID, sectorData)
cfg = Configuration()
data = sectorData["data"]
self.maxIncr = sectorData["maxIncr"]/100
self.maxDecr = sectorData["maxDecr"]/100
self.assSectorIDs = data["assSector"]
self.assSectorObjs = []
self.g = data["g"]
self.m = np.zeros(Configuration.nYears);
def step(self):
if (self.assSectorObjs == None) or (self.assSectorObjs==[]):
raise Exception(self.sectorID + ": associated sector not defined.")
self.currentStep+=1
year = Configuration.firstYear + self.currentStep
incr = 0
for assSectorObj in self.assSectorObjs:
incr += assSectorObj[year][0]
self.m[self.currentStep]=self._detectConstraints(self.currentStep, self.g*incr)
#return self.m[self.currentStep]