Skip to content

Commit

Permalink
add alkane labelling as test
Browse files Browse the repository at this point in the history
  • Loading branch information
lilyminium committed Mar 22, 2024
1 parent 60758a7 commit 50d9b80
Show file tree
Hide file tree
Showing 2 changed files with 56 additions and 2 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -1528,7 +1528,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.11.7"
"version": "3.11.4"
}
},
"nbformat": 4,
Expand Down
56 changes: 55 additions & 1 deletion openff/nagl/tests/label/test_labels.py
Original file line number Diff line number Diff line change
Expand Up @@ -100,6 +100,60 @@ def test_label_with_conformers_on_fly(self, small_dataset):
columns = ["mapped_smiles", "conformers", "n_conformers", "charges"]
assert small_dataset.dataset.schema.names == columns

def test_label_alkane_dataset(self):
# test conformer generation and labelling
# as in examples

training_alkanes = [
'C',
'CC',
'CCC',
'CCCC',
'CC(C)C',
'CCCCC',
'CC(C)CC',
'CCCCCC',
'CC(C)CCC',
'CC(CC)CC',
]

training_dataset = LabelledDataset.from_smiles(
"training_data",
training_alkanes,
mapped=False,
overwrite_existing=True,
)
training_df = training_dataset.to_pandas()
assert training_df.mapped_smiles[0] == "[H:2][C:1]([H:3])([H:4])[H:5]"

label_conformers = LabelConformers(
# create a new 'conformers' with output conformers
conformer_column="conformers",
# create a new 'n_conformers' with number of conformers
n_conformer_column="n_conformers",
n_conformer_pool=500, # initially generate 500 conformers
n_conformers=10, # prune to max 10 conformers
rms_cutoff=0.05,
)

label_am1_charges = LabelCharges(
charge_method="am1-mulliken", # AM1
# use previously generate conformers instead of new ones
use_existing_conformers=True,
# use the 'conformers' column as input for charge assignment
conformer_column="conformers",
# write generated charges to 'target-am1-charges' column
charge_column="target-am1-charges",
)

labellers = [
label_conformers, # generate initial conformers,
label_am1_charges,
]

training_dataset.apply_labellers(labellers)



class TestLabelMultipleDipoles:

Expand Down Expand Up @@ -174,4 +228,4 @@ def test_apply_label(self, dataset_with_conformers_and_charges):

calculated_esps = pydict["esps"]
for esps, lengths in zip(calculated_esps, calculated_esp_lengths):
assert len(esps) == sum(lengths)
assert len(esps) == sum(lengths)

0 comments on commit 50d9b80

Please sign in to comment.