|
330 | 330 | "path_lengths = [len(step.active[0].trajectory) for step in flexible.steps]\n", |
331 | 331 | "plt.hist(path_lengths, bins=40, alpha=0.5);\n", |
332 | 332 | "print(\"Maximum:\", max(path_lengths), \"(\"+str(max(path_lengths)*engine.snapshot_timestep)+\")\")\n", |
333 | | - "print(\"Average:\", \"{0:.2f}\".format(np.mean(path_lengths)), \"(\"+(np.mean(path_lengths)*engine.snapshot_timestep).format(\"%.3f\")+\")\")" |
| 333 | + "print(\"Average:\", \"{0:.2f}\".format(np.mean(path_lengths)), \"(\"+(np.mean(path_lengths)*engine.snapshot_timestep).format(\"%.3f\")+\")\")\n", |
| 334 | + "plt.ylabel(\"Count\")\n", |
| 335 | + "plt.xlabel(\"Path length (Frames)\");" |
334 | 336 | ] |
335 | 337 | }, |
336 | 338 | { |
|
406 | 408 | " xticklabels=tick_labels,\n", |
407 | 409 | " yticklabels=tick_labels, \n", |
408 | 410 | " label_format=\"{:4.2f}\")\n", |
409 | | - "ax = plotter.plot(cmap=\"Blues\")" |
| 411 | + "ax = plotter.plot(cmap=\"Blues\")\n", |
| 412 | + "plt.xlabel(\"$\\phi$\")\n", |
| 413 | + "plt.ylabel(\"$\\psi$\");" |
410 | 414 | ] |
411 | 415 | }, |
412 | 416 | { |
|
442 | 446 | "ax = plotter.plot(xlim=(-np.pi, 0), ylim=(-np.pi/2, np.pi), cmap=\"Blues\")\n", |
443 | 447 | "plotter.plot_trajectory(trajA, '-r', lw=1)\n", |
444 | 448 | "plotter.plot_trajectory(trajB, '-k', lw=1)\n", |
| 449 | + "plt.xlabel(\"$\\phi$\")\n", |
| 450 | + "plt.ylabel(\"$\\psi$\");\n", |
445 | 451 | "#plt.savefig(\"AD_tps_pathdensity.pdf\") # uncomment this line to save as PDF" |
446 | 452 | ] |
447 | 453 | }, |
|
0 commit comments