IPKM/mdevaluate
6
1
Fork 0
Code Issues Pull Requests Activity

Compare commits

base: IPKM:787882810c4496d16bd76b69b128e77724aa444c
Branches Tags
IPKM:main IPKM:mdeval_dev
...
compare: IPKM:4394f7053017b79b99151cb77112438d50149e81
Branches Tags
IPKM:mdeval_dev IPKM:main

3 Commits
787882810c ... 4394f70530

Author SHA1 Message Date
Sebastian Kloth
4394f70530 Fixed pyproject.toml 2024-01-31 09:02:57 +01:00
Sebastian Kloth
298da3818d Changed parameter normed to density in np.histogram 2024-01-31 09:02:36 +01:00
Sebastian Kloth
d9278eed83 Removed prints from rdf 2024-01-31 08:59:39 +01:00
2 changed files with 2 additions and 5 deletions
Expand all files Collapse all files

2
pyproject.toml
View File

@ -10,6 +10,6 @@ dependencies = [
"pandas", "pandas",
"dask", "dask",
"pathos", "pathos",
"tables" "tables",
"pyedr" "pyedr"
] ]

5
src/mdevaluate/distribution.py
View File

@ -126,9 +126,6 @@ def rdf(
particles_in_volume = int( particles_in_volume = int(
np.max(number_of_neighbors(atoms_a, query_atoms=atoms_b, r_max=bins[-1])) * 1.1 np.max(number_of_neighbors(atoms_a, query_atoms=atoms_b, r_max=bins[-1])) * 1.1
) )
print(atoms_a[:5])
print(atoms_b[:5])
print(" ")
distances, indices = next_neighbors( distances, indices = next_neighbors(
atoms_a, atoms_a,
atoms_b, atoms_b,
@ -309,7 +306,7 @@ def next_neighbor_distribution(
)[1] )[1]
resname_nn = reference.residue_names[nn] resname_nn = reference.residue_names[nn]
count_nn = (resname_nn == atoms.residue_names.reshape(-1, 1)).sum(axis=1) count_nn = (resname_nn == atoms.residue_names.reshape(-1, 1)).sum(axis=1)
return np.histogram(count_nn, bins=bins, normed=normed)[0] return np.histogram(count_nn, bins=bins, density=normed)[0]
def hbonds( def hbonds(