General updates | MT JW Following from lead team discussion - when reasonable let’s add 24h delay after last PR is merged so “latest CI” can run. Request for spiritual guidance: Mysterious CI failures for GBSA energy comparison (for systems created using openff vs. created using OpenMM) , only on macos15, only for a single mol in freesolv, only for a single permutation of settings. Any hints? MT – OS differences for packages are possible. Could try A/B testing openmm versions and platforms. OE/RDK/AT versions also come to mind. If it’s just 1 of many mols, it’s possible that others have differences that are staying just under tolerances. You’ll probably need to diff serialized systems. If you can pass same system through different VMs and get different outputs that’s a smoking gun. I’m seeing some small differences using openmm 8.1.2 and 8.2.0. Recommend looking at nonbonded settings FIRST (especially switching function, cutoffs, 1-4 exceptions). Also compare GBSA parameters. Also wouldn’t ignore possibility that something is being set differently for GBSA. Also I skimmed the OMM 8.2 releasenotes and PRs and didn’t see anything that touches nonbonded values, but I wasn’t looking for GBSA specifically. I have a very weak suspicion that there’s a hard-to-pin-down corner case, maybe something to do with platform selection.
JW – I should check into platform selection on macos 15 - That could explain a lot.
I’ve “picked” OFFTK + OFF Units 0.2.3 main this sprint Issues MT – Start with units main branch, don’t try to use 0.2.3. For mypy issues - We’re a few releases behind pint (and pinned as such) since they move fast and break things. Pint is kinda a relic from pre-typing python - eg if you’re given a class like Quantity, it could be somethat that wraps a scalar or list or numpy array and has fancy stuff with xarray and pandas etc. So it’s highly polymorphic (by design) and getting typing to work is difficult since our usage of it just corresponds to a small subset of pint’s possibilities. So you’ll get a lot of typing errors regarding union-types. One way this could be worked around is by providing stubs (adding type annotations where upstreams don’t provide them). Also, using unit.Quantity to instantiate something confuses unit system, but doing the same by importing Quantity directly is fine. JW – Two things I’m thinking of are:
Gov board told us to go ahead and create PTM parameterization product. A good route forward will be directly assigning charges to some atoms of an already-made protein interchange. Do you estimate that will that hit major roadblocks?
|