Skip to content

Solve issues with snapshot generation in Sprint 12 #175

@jchodera

Description

@jchodera

I must have broken some aspect of snapshot generation in Sprint 12:

ValueError: xyz must be shape (Any, 9388, 3). You supplied  (4, 9465, 3)
Traceback (most recent call last):
  File "/home/server/fah-xchem-prod/fah-xchem/fah_xchem/analysis/structures.py", line 605, in generate_representative_snapshot
    sliced_snapshots, components = self.extract_snapshot(
  File "/home/server/fah-xchem-prod/fah-xchem/fah_xchem/analysis/structures.py", line 428, in extract_snapshot
    trajectory = self.load_trajectory(
  File "/home/server/fah-xchem-prod/fah-xchem/fah_xchem/analysis/structures.py", line 312, in load_trajectory
    trajectory = md.Trajectory(xyz, topology)
  File "/home/server/miniconda/envs/fah-xchem-prod-2022.01.31/lib/python3.9/site-packages/mdtraj/core/trajectory.py", line 1226, in __init__
    self.xyz = xyz
  File "/home/server/miniconda/envs/fah-xchem-prod-2022.01.31/lib/python3.9/site-packages/mdtraj/core/trajectory.py", line 939, in xyz
    value = ensure_type(value, np.float32, 3, 'xyz', shape=shape,
  File "/home/server/miniconda/envs/fah-xchem-prod-2022.01.31/lib/python3.9/site-packages/mdtraj/utils/validation.py", line 148, in ensure_type
    raise error

The good news is that the solute atoms in core22 0.0.20 are now written in sorted (ascending) atom index order, so this should be straightforward to solve in a manner that reduces complexity.

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions