Skip to content

Commit

Permalink
improved convergence N figure
Browse files Browse the repository at this point in the history
  • Loading branch information
simongravelle committed Jan 24, 2024
1 parent e4fd582 commit 5b04b75
Show file tree
Hide file tree
Showing 146 changed files with 74 additions and 63 deletions.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
31 changes: 17 additions & 14 deletions examples/illustrations/bulk-water-tip4p-vs-size/analysis.py
Original file line number Diff line number Diff line change
Expand Up @@ -18,8 +18,8 @@
R1_intra = 1e-5
R1_inter = 1e-5
while run:
for N in [25, 39, 62, 99, 158, 251, 398, 631, 1002, 1589, 2521]:
for n in range(10):
for N in [25, 39, 62, 99, 158, 251, 398, 631, 1002, 1589, 2521, 6500, 10000]:
for n in range(20):
# Import MDA universe
datapath = path + "N"+str(N)+"/"
if os.path.exists(datapath+"prod-"+str(n)+".xtc"):
Expand All @@ -28,10 +28,14 @@
name_inter = "N"+str(N)+"_inter_n"+str(n)
n_intra = measure_N(name_intra)
n_inter = measure_N(name_inter)
u = mda.Universe(datapath+"topology.data", datapath+"prod-"+str(n)+".xtc")
hydrogen = u.select_atoms("type 2")
try:
u = mda.Universe(datapath+"topology.data", datapath+"prod-"+str(n)+".xtc")
hydrogen = u.select_atoms("type 2")
except:
u = mda.Universe(datapath+"prod.tpr", datapath+"prod-"+str(n)+".xtc")
hydrogen = u.select_atoms("type HW")
n_hydrogen = hydrogen.n_atoms
if n_intra < n_hydrogen:
if n_intra < n_hydrogen*2:
# Calculated NMR properties
intra = nmrformd.NMR(u, atom_group = hydrogen,
neighbor_group = hydrogen,
Expand All @@ -44,13 +48,12 @@
type_analysis = "inter_molecular",
number_i = 1)
save_result(inter, name = name_inter)
else:
print("N", N, "no calculation")
#run = False
# Print information
R1_intra = measure_R1(name_intra)
R1_inter = measure_R1(name_inter)
print("N =", N, "n", n, "----", n_intra, n_hydrogen,
"R1:", np.round(R1_intra,2), np.round(R1_inter,2),
"T1:", np.round(1/(R1_intra+R1_inter),2))
# Print information
R1_intra = measure_R1(name_intra)
R1_inter = measure_R1(name_inter)
print("N =", N, "n", n, "----", n_intra, n_hydrogen,
"R1:", np.round(R1_intra,2), np.round(R1_inter,2),
"T1:", np.round(1/(R1_intra+R1_inter),2))
#else:
# print("N", N, "no calculation")
print(" ")
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
106 changes: 57 additions & 49 deletions examples/illustrations/bulk-water.ipynb

Large diffs are not rendered by default.

0 comments on commit 5b04b75

Please sign in to comment.