ExcitedStates / qfit-3.0

qFit: Automated and unbiased multi-conformer models from X-ray and EM maps.
MIT License
35 stars 11 forks source link

Make BIC routine _actually_ select based on lowest BIC #328

Closed blake-riley closed 1 year ago

blake-riley commented 1 year ago

Pull Request Checklist


Description of the Change

We were previously selecting the number of conformers from the last iteration of the BIC loop, because a "else break" had been commented and removed. As such, every MIQP was performed 5 times (redundantly!), and the last MIQP iteration was the one that made it through (threshold/tdmin = 0.2).

This PR implements the desired behaviour. It calculates the MIQP solutions for all thresholds, and then chooses the solution with lowest BIC.

Discussion points

The number of conformers selected is drastically decreased, and two tests fail: ones which check that the provided examples return an expected number of conformers. These two examples return only one rotamer each now. I do not know why.

Release Notes