diff --git a/docs/.DS_Store b/docs/.DS_Store new file mode 100644 index 0000000..d50adf6 Binary files /dev/null and b/docs/.DS_Store differ diff --git a/docs/build/doctrees/environment.pickle b/docs/build/doctrees/environment.pickle index bf201cd..1545d0f 100644 Binary files a/docs/build/doctrees/environment.pickle and b/docs/build/doctrees/environment.pickle differ diff --git a/docs/build/doctrees/index.doctree b/docs/build/doctrees/index.doctree index 92f56cc..dfb2e3f 100644 Binary files a/docs/build/doctrees/index.doctree and b/docs/build/doctrees/index.doctree differ diff --git a/docs/build/doctrees/modules.doctree b/docs/build/doctrees/modules.doctree index 614127a..6c93083 100644 Binary files a/docs/build/doctrees/modules.doctree and b/docs/build/doctrees/modules.doctree differ diff --git a/docs/build/html/.buildinfo b/docs/build/html/.buildinfo index eae5db8..42691c2 100644 --- a/docs/build/html/.buildinfo +++ b/docs/build/html/.buildinfo @@ -1,4 +1,4 @@ # Sphinx build info version 1 # This file hashes the configuration used when building these files. When it is not found, a full rebuild will be done. -config: af1005e6b5cb943e75ffafa7a50014fd +config: 3e24725cde270cd0e889a75e43e9746b tags: 645f666f9bcd5a90fca523b33c5a78b7 diff --git a/docs/build/html/_modules/index.html b/docs/build/html/_modules/index.html index 5e00664..eba95e4 100644 --- a/docs/build/html/_modules/index.html +++ b/docs/build/html/_modules/index.html @@ -4,13 +4,21 @@
-
system.step += 1
if system.step % sample_frequency == 0:
sample_system.update(system)
- return system
-
system.step += 1
if system.step % sample_frequency == 0:
sample_system.update(system)
- return system
-
velocities[0] += 0.5 * accelerations[0] * timestep_length
velocities[1] += 0.5 * accelerations[1] * timestep_length
return [velocities[0], velocities[1]]
line = self.ax[1, 1].lines[0]
line.set_ydata(iq)
line.set_xdata(x)
- self.ax[1, 1].set_ylim([np.amin(iq) - np.amax(iq) * 0.05, np.amax(iq) + np.amax(iq) * 0.05])
- self.ax[1, 1].set_xlim([np.amin(x), np.amax(x)])
fig, ax = environment(4)
setup_cellview(ax[0, 0], system)
- setup_forceview(ax[0, 1])
+ setup_forceview(ax[1, 1])
setup_pressureview(ax[1, 0])
- setup_tempview(ax[1, 1])
+ setup_tempview(ax[0, 1])
plt.tight_layout()
self.ax = ax
@@ -128,8 +137,8 @@ Source code for pylj.sample
The whole system information.
"""
update_cellview(self.ax[0, 0], system)
- update_forceview(self.ax[0, 1], system)
- update_tempview(self.ax[1, 1], system)
+ update_forceview(self.ax[1, 1], system)
+ update_tempview(self.ax[0, 1], system)
update_pressureview(self.ax[1, 0], system)
self.fig.canvas.draw()
The axes related to each of the panes. For panes=1 this is a single object, for panes=2 it is a 1-D array and
for panes=4 it is a 2-D array.
"""
- fig, ax = plt.subplot()
if panes == 1:
fig, ax = plt.subplots(figsize=(4, 4))
elif panes == 2:
@@ -305,7 +313,7 @@ Source code for pylj.sample
ax.plot([0], color='#34a5daff')
ax.set_yticks([])
ax.set_ylabel('log(I[q])', fontsize=16)
- ax.set_xlabel('log(q)/Å$^{-1}$', fontsize=16)
+ ax.set_xlabel('q/Å$^{-1}$', fontsize=16)
q: array_like
The scattering profile's q for each timestep, to later be averaged.
"""
- hist, bin_edges = np.histogram(system.distances, bins=np.linspace(0, system.box_length/2 + 0.5, 100))
- gr = hist / (system.number_of_particles * (system.number_of_particles / system.box_length ** 2) * np.pi *
- (bin_edges[:-1] + 0.5 / 2.) * 0.5)
- x2 = np.log(np.fft.rfftfreq(len(gr))[5:])
- y2 = np.log(np.fft.rfft(gr)[5:])
+ qw = np.linspace(2 * np.pi /(system.box_length/2)*4, system.box_length/10, 100)
+ i = np.zeros_like(qw)
+ for j in range(0, len(qw)):
+ i[j] = np.sum(3.644 * (np.sin(qw[j] * system.distances))/(qw[j] * system.distances))
+ if i[j] < 0:
+ i[j] = 0
+ x2 = qw
+ y2 = i
average_diff.append(y2)
q.append(x2)
line1 = ax.lines[0]
line1.set_xdata(x2)
line1.set_ydata(y2)
- ax.set_ylim([np.amin(y2) - np.amax(y2) * 0.05, np.amax(y2) + np.amax(y2) * 0.05])
- ax.set_xlim([np.amin(x2), np.amax(x2)])
The whole system information.
"""
line = ax.lines[0]
- line.set_ydata(system.pressure_sample * 1e10)
+ line.set_ydata(system.pressure_sample)
line.set_xdata(np.arange(0, system.step) * system.timestep_length)
ax.set_xlim(0, system.step * system.timestep_length)
- ax.set_ylim(np.amin(system.pressure_sample * 1e10) - np.amax(system.pressure_sample * 1e10) * 0.05,
- np.amax(system.pressure_sample * 1e10) + np.amax(system.pressure_sample * 1e10) * 0.05)