39 |
iters, simuls, fc = get_output(fname, "f=") |
iters, simuls, fc = get_output(fname, "f=") |
40 |
# sort out restarts |
# sort out restarts |
41 |
iters0 = np.asarray(iters) |
iters0 = np.asarray(iters) |
42 |
for k,it in enumerate(iters[0:]): |
for k,it in enumerate(iters[1:]): |
43 |
if iters0[k]<iters0[k-1]: |
print(k,it) |
44 |
iters0[k:] = iters0[k:]+(iters0[k-1]-iters0[k]+1) |
kp1=k+1 |
45 |
|
if iters0[kp1]<iters0[kp1-1]: |
46 |
|
iters0[kp1:] = iters0[kp1:]+(iters0[k]-iters0[kp1]+1) |
47 |
|
|
48 |
fig, ax1 = plt.subplots() |
fig, ax1 = plt.subplots() |
49 |
ax1.semilogy(iters0,(np.asarray(fc)-fc[-1]),'bx-') |
ax1.semilogy(iters0,(np.asarray(fc)-fc[-1]),'bx-') |