vvoelz / biceps

Bayesian inference of conformational populations
https://github.com/vvoelz/biceps
Other
12 stars 3 forks source link

Cleanup with tests #82

Closed robraddi closed 3 years ago

robraddi commented 3 years ago

The official biceps workflow now passes the ensemble object to PosterSampler without converting to a list. The ensemble object will be needed for future versions.

ensemble = biceps.Ensemble(float(lam), list(energies))
ensemble.initialize_restraints(list(input_data))
sampler = biceps.PosteriorSampler(object(ensemble))

Also added a burn-in time for the sample(). Default is 0.

sampler.sample(int(nsteps), int(burn))
sampler.traj.process_results(str(filename))

Lastly, the sampler object is automatically saved when calling on process_results

Tests were performed for all systems and notebooks were edited to reflect the alterations in the code.