def get_coefficients(self):
log.log('get_coefficients function called')
self.coefficients = {}
converter = 8192
if self.mj in ('mj04', 'mj22', 'mj36'): converter = 4096
l = []
f = './temp/ini/base.ini'
if not os.path.isfile(f): return
i = IniParser()
i.read(f)
for s in i.options('AnalogInput_Signals'):
if 'scale' in s.lower() and i.get('AnalogInput_Signals', s) != 'x': l.append((s, i.getint('AnalogInput_Signals', s)))
#endfor
f = './temp/ini/parameters.ini'
if not os.path.isfile(f): return
i = IniParser()
i.read(f)
for s, a in l:
self.coefficients[s] = i.getfloat(s, 'Range')
#endfor
f = './temp/ini/settings.ini'
if not os.path.isfile(f): return
i = IniParser()
i.read(f)
for s, a in l:
print s, a
print self.coefficients[s]
print i.getfloat('AnalogInput_Coefficients', 'Coefficient%i' % a)
print type(self.coefficients[s])
print type(i.getfloat('AnalogInput_Coefficients', 'Coefficient%i' % a))
capacity = self.coefficients[s]
difference = (i.getfloat('AnalogInput_Coefficients', 'Coefficient%i' % a) / (capacity/converter) - 1) * 100
offset = i.getint('AnalogInput_Offsets', 'Offset%i' % a)
self.coefficients[s] = (capacity, difference, offset)
def get_little_amounts(self):
log.info('little amounts function called')
silo = {}
self.little_amounts_fails = {}
doba = None
vaha = None
navazeno = None
analog = None
adaptivita = None
pozadavek = None
zasobnik = 99 # can delete this line after three month (bacuse of error in info txt from Piestany)
fn = 'control.log'
for l in give_lines(fn, 3):
if 'poddavkovani ' in l.lower():
zasobnik = int(l.split()[-1])
vaha = None
navazeno = None
pozadavek = None
elif 'dobaotevreni ' in l.lower():
pom = float(l.split()[-1])
if pom > 3 or pom < 0: continue # filtrace demetskych hodnot v milisekundach a zapornych hodnot
doba = pom
elif 'mezistav1 ' in l.lower():
vaha = float(l.split()[-1])
elif 'NejblizsiCil2' in l and vaha:
pozadavek = float(l.split()[-1]) - vaha
elif 'PrumernaHodnotaAI[CisloAnalogu] - VychoziVaha' in l:
analog = float(l.split()[-1])
elif 'frakcenamereno2 ' in l.lower() and vaha:
navazeno = float(l.split()[-1]) - vaha
#endif
if doba and navazeno:
log.info('%i, %0.2f, %i' % (navazeno, doba, zasobnik))
if navazeno < 30 or navazeno > 500:
log.info('little amount fail: %i, %0.2f, %i' % (navazeno, doba, zasobnik))
if not zasobnik in self.little_amounts_fails:
self.little_amounts_fails[zasobnik] = []
#endif
self.little_amounts_fails[zasobnik].append('%i, %0.2f' % (navazeno, doba))
doba = None
continue
#endif
if analog: adaptivita = navazeno - analog + vaha
if not zasobnik in silo: silo[zasobnik] = []
silo[zasobnik].append((doba, navazeno, pozadavek, adaptivita))
doba = None
vaha = None
navazeno = None
analog = None
adaptivita = None
#endif
#endfor
#log.info('%s' % silo)
self.silo = {}
self.little_amounts_error = False
f = './temp/ini/adapt.ini'
if not os.path.isfile(f): return
adapt = IniParser()
adapt.read(f)
f = './temp/ini/parameters.ini'
if not os.path.isfile(f): return
parameters = IniParser()
parameters.read(f)
for k in silo:
self.little_amounts_error = adapt.has_option('Aggregate', 'VelocityOfFilling1-%s' % k)
param_velocity = 0
param_offset = 0
if parameters.has_option('Aggregate_Bin1_Silo%s' % k, 'Velocity_Base'):
param_velocity = parameters.getint('Aggregate_Bin1_Silo%s' % k, 'Velocity_Base')
param_offset = parameters.getfloat('Aggregate_Bin1_Silo%s' % k, 'Additional_Time')
else:
self.little_amounts_error = True
#endif
x = 0
sez = {}
while x < 1.5:
if not x in sez: sez[x] = []
for doba, navazeno, pozadavek, adaptivita in silo[k]:
if doba-x == 0: continue #I dont know why
sez[x].append(navazeno/(doba - x))
#endfor
x = x + 0.01
#endwhile
minimum = None
for x in sorted(sez):
sqr_sum = 0
avg = 0
dev = 0
avg = sum(sez[x]) * 1.0 / len(sez[x])
for velocity in sez[x]:
sqr_sum += ((velocity - avg) ** 2) / len(sez[x])
#endfor
dev = sqr_sum ** 0.5
if not minimum or dev < minimum:
#print dev
minimum = dev
self.silo[k] = (x, avg, len(sez[x]), dev, param_velocity, param_offset)
def get_adaptivity(mj):
log.log('get_adaptivity function called')
class Adaptivity():
def __init__(self, silo, scale, coefficient):
self.silo = silo
#log.info('creating new class for %s' % self.silo)
self.scale = scale
self.coefficient = coefficient
self.first = None
#enddef
def check_line(self, line):
if not self.scale in line: return
if not self.first and self.first != 0:
self.first = int(line.split('=')[1])
#log.info('%s first %s' % (self.silo, self.first))
else:
self.adaptivity = (int(line.split('=')[1]) - self.first) * self.coefficient
self.time = get_datetime_from_log(line)
return 'all done'
#endif
#enddef
#enclass
ret = {} # ret[silo] = [(time, adaptivity), ...]
materials = {'AGG': 'Aggregate', 'CEM': 'Cement', 'ADM': 'Admixture', 'WTR': 'Water'}
cache = {} # cache[silo] = {'Scale': scale, 'Coefficient': coefficient}
incompletes = [] # incompletes = [Adativity(), ...]
f = './temp/ini/base.ini'
if not os.path.isfile(f): return ret
base = IniParser()
base.read(f)
f = './temp/ini/settings.ini'
if not os.path.isfile(f): return ret
settings = IniParser()
settings.read(f)
fn = 'recordedsignals.log'
for l in give_lines(fn, 3):
for m in materials:
if m in l and 'O_' in l and 'intoBin' in l and ')=0' in l:
silo = '%s%s' % (m, l.split(m)[1].split('into')[0]) # for example AGG5
if not silo in cache:
bin = 'Bin%s' % l.split('intoBin')[1].split(')=')[0]
scale = '%s_Scale%s' % (materials[m], base.get('%s_%s' % (materials[m], bin), 'CorrespondingScale'))
signal = base.get('AnalogInput_Signals', scale)
coefficient = settings.getfloat('AnalogInput_Coefficients', 'Coefficient%s' % signal)
cache[silo] = {}
cache[silo]['Scale'] = scale
cache[silo]['Coefficient'] = coefficient
#endif
if not silo in ret: ret[silo] = []
# Adding new adaptivities to list
incompletes.append(Adaptivity(silo, cache[silo]['Scale'], cache[silo]['Coefficient']))
break
#endif
#endfor
if 'Scale' in l:
for i in incompletes:
if i.check_line(l):
ret[i.silo].append((i.time, i.adaptivity))
incompletes.remove(i)
#endif
#endfor
#endif
#endfor
# read alternately adaptivity from info.txt
fn = 'control.log'
materials = {'aggregate': 'AGG', 'cement': 'CEM', 'admixture': 'ADM', 'water': 'WTR'}
for l in give_lines(fn, 3):
if not 'adaptivity for' in l: continue
material_type = materials[l.split(' for ')[1].split(' bin ')[0]]
silo_number = l.split(' silo ')[1].split(':')[0]
adaptivity = l.split(': ')[-1].strip()
time = get_datetime_from_log(l)
silo = '%s%s' % (material_type, silo_number)
if not silo in ret: ret[silo] = []
ret[silo].append((time, float(adaptivity)))
#endfor
# create output files
for silo in ret:
log.info('creating file for silo %s' % silo)
a = '../public_html_incoming'
if not os.path.isdir(a):
a = './public_html_incoming'
#endif
f = open('%s/%s%s.csv' % (a, mj, silo), 'w')
for x,y in ret[silo]:
f.write('%s;%i\n' % (x.strftime('%Y-%m-%d %H:%M:%S'), y))
#endfor
f.close()
#endfor
return ret
