GunnerLab / MCCE_Benchmarking

For benchmarking MCCE experiments viz experimental pKas or viz other MCCE experiments (A/B testing)
MIT License
0 stars 0 forks source link

Capture all out of bound messages in pK.out #26

Closed CatChenal closed 7 months ago

CatChenal commented 8 months ago
  1. Check lib/monte.c and lib/mont2.c for all possible "out of bound" info printed in pK.out, often in lieu of actual pKa.

  2. Before collating all pK.out files, save those "oob" lines to a tmp file, e.g.:

    >egrep 'oob msg' /*/pK.out > pk_oob
  3. Get the count:

    >wc -l pk_oob | awk '{print "pk_oob."$1}'
    >> pk_oob.123                                       # => new fname retrieved w/subprocess
  4. Rename oob file: cmd = f"mv pk_oob {fname}; rm pk_oob" # run via subprocess

  5. In code:

for fp in Path(

).glob("pk_oob.*"): n_oob = int(fp.suffix) break


=> Without doing any analysis on oob pKas, their total is known.

5. Determine how `pkout_df` is handling these oob lines.