P283370/NuMRI
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

press new

Browse Source
This commit is contained in:
jeremias
2021-03-29 10:45:02 +02:00
parent 842a3fbcc6
commit 74785679cb
5 changed files with 148 additions and 224 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
kalman/graphics/figure2.py
View File

@ -72,7 +72,7 @@ def plot_parameters(dat, input_file, deparameterize=False, ref=None):
if bnd_c['id'] == bnd_set['id']:
ids.append(bnd_c['id'])
current_val.append(bnd_set['value'][0])
labels.append('R_' + str(bnd_c['id']))
labels.append('$R_' + str(bnd_c['id']))
elif 'dirichlet' in bnd_c['type']:
current_val.append(inputfile['boundary_conditions'][1]['parameters']['U'])
@ -105,13 +105,17 @@ def plot_parameters(dat, input_file, deparameterize=False, ref=None):
for i in range(dim):
axes.plot(t, theta[:, i] + 1.5*i, '-', color=col_,label=legends_)
true_level = np.log(true_values[ids[i]]/current_val[i])/np.log(2)
rec_value = np.round(2**theta[-1, i]*current_val[i],1)
cur_key = ids[i]
axes.plot(t, theta[:, i] + 1.5*i, '-', color=col_,label=legends_ + '= ' + str(rec_value) + '/' + str(true_values[cur_key]) + '$')
axes.fill_between(t, theta[:, i] + 1.5*i - np.sqrt(P[:, i, i]),
theta[:, i] + 1.5*i + np.sqrt(P[:, i, i]), alpha=0.3,
color=col_)
true_level = np.log(true_values[ids[i]]/current_val[i])/np.log(2)
axes.plot(t,1.5*i + t*0 + true_level , color=col_,ls='--')
col_ = next(col)
legends_=next(legends)