dt_Handle transform, container, molecule1, moleculeX;
/* Read SMIRKS in. */
transform = dt_smirkin(25, "[O:1][H:2]>>[O-:1].[H+:2]");
container = dt_alloc_seq();
molecule1 = dt_smilin(1, "O");
dt_append(container, molecule1);
moleculeX = dt_smilin(12, "Oc1ccc(O)cc1");
dt_append(container, moleculeX);
|