org.xmlcml.cml.tools.AtomMatcher3D Maven / Gradle / Ivy

Go to download
/**
 *    Copyright 2011 Peter Murray-Rust et. al.
 *
 *    Licensed under the Apache License, Version 2.0 (the "License");
 *    you may not use this file except in compliance with the License.
 *    You may obtain a copy of the License at
 *
 *        http://www.apache.org/licenses/LICENSE-2.0
 *
 *    Unless required by applicable law or agreed to in writing, software
 *    distributed under the License is distributed on an "AS IS" BASIS,
 *    WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 *    See the License for the specific language governing permissions and
 *    limitations under the License.
 */

package org.xmlcml.cml.tools;

import java.util.ArrayList;
import java.util.List;

import org.xmlcml.cml.element.CMLAtomSet;
import org.xmlcml.cml.element.CMLMap;
import org.xmlcml.cml.element.CMLMolecule;
import org.xmlcml.euclid.Transform3;

public class AtomMatcher3D extends AtomMatcher {

	public CMLMap match(
			CMLAtomSet atomSet0,CMLAtomSet atomSet1, String title) {
		CMLMap cmlMap = makeMap();
		throw new RuntimeException("NYI");
	}
 
    /** pairwise atom-atom alignment of molecules.
     * 
     * @param identicalMoleculeList list of identical molecules (e.g. by morgan)
     * @return upper triangle of transformations from mol(j) onto mol(i)
     */
	private static List> matchAndAlignMolecules(
			List identicalMoleculeList) {
		AtomMatcher atomMatcher = new MorganAtomMatcher();
//		atomMatcher.setAtomMatchStrategy(Strategy.MATCH_MORGAN);
		CMLMap[][] mapMatrix = atomMatcher.getMoleculeMatch(identicalMoleculeList, identicalMoleculeList);
		List> moleculePairListList = new ArrayList>();
		int i = 0;
		for (CMLMap[] mapx : mapMatrix) {
			List moleculePairList = new ArrayList();
			moleculePairListList.add(moleculePairList);
			int j = 0;
			for (CMLMap map : mapx) {
				if (j > i) {
					MoleculeTool moleculeTooli = MoleculeTool.getOrCreateTool(identicalMoleculeList.get(i));
					MoleculePair moleculePair = moleculeTooli.fitToMoleculeTool(map, identicalMoleculeList.get(j));
					moleculePairList.add(moleculePair);
				}
				j++;
			}
			i++;
		}
		return moleculePairListList;
	}

	public CMLMap match(CMLAtomSet atomSet0, CMLAtomSet atomSet1, Transform3 t3) {
		throw new RuntimeException("NYI");
	}
}