jumbo46cmlfragmentstest.java

化学图形处理软件

    public void testCoordinates2D() throws Exception {
        String cmlString = "<molecule id='m1'><atomArray atomID='a1 a2' x2='0.0 0.1' y2='1.2 1.3'/></molecule>";
        
        IChemFile chemFile = parseCMLString(cmlString);
        org.openscience.cdk.interfaces.IMolecule mol = checkForSingleMoleculeFile(chemFile);
        
        assertEquals(2, mol.getAtomCount());
        assertNotNull(mol.getAtom(0).getPoint2d());
        assertNotNull(mol.getAtom(1).getPoint2d());
        assertNull(mol.getAtom(0).getPoint3d());
        assertNull(mol.getAtom(1).getPoint3d());
    }
  
    public void testCoordinates3D() throws Exception {
        String cmlString = "<molecule id='m1'><atomArray atomID='a1 a2' x3='0.0 0.1' y3='1.2 1.3' z3='2.1 2.5'/></molecule>";
        
        IChemFile chemFile = parseCMLString(cmlString);
        org.openscience.cdk.interfaces.IMolecule mol = checkForSingleMoleculeFile(chemFile);
        
        assertEquals(2, mol.getAtomCount());
        assertNull(mol.getAtom(0).getPoint2d());
        assertNull(mol.getAtom(1).getPoint2d());
        assertNotNull(mol.getAtom(0).getPoint3d());
        assertNotNull(mol.getAtom(1).getPoint3d());
    }
    
    public void testFractional3D() throws Exception {
        String cmlString = "<molecule id='m1'><atomArray atomID='a1 a2' xFract='0.0 0.1' yFract='1.2 1.3' zFract='2.1 2.5'/></molecule>";
        
        IChemFile chemFile = parseCMLString(cmlString);
        org.openscience.cdk.interfaces.IMolecule mol = checkForSingleMoleculeFile(chemFile);
        
        assertEquals(2, mol.getAtomCount());
        assertNull(mol.getAtom(0).getPoint3d());
        assertNull(mol.getAtom(1).getPoint3d());
        assertNotNull(mol.getAtom(0).getFractionalPoint3d());
        assertNotNull(mol.getAtom(1).getFractionalPoint3d());
    }
    
    public void testMissing2DCoordinates() throws Exception {
        String cmlString = 
          "<molecule id='m1'><atomArray><atom id='a1' xy2='0.0 0.1'/><atom id='a2'/><atom id='a3' xy2='0.1 0.0'/></atomArray></molecule>";
          
        IChemFile chemFile = parseCMLString(cmlString);
        org.openscience.cdk.interfaces.IMolecule mol = checkForSingleMoleculeFile(chemFile);
        
        assertEquals(3, mol.getAtomCount());
        IAtom atom1 = mol.getAtom(0);
        IAtom atom2 = mol.getAtom(1);
        IAtom atom3 = mol.getAtom(2);
        
        assertNotNull(atom1.getPoint2d());
        assertNull   (atom2.getPoint2d());
        assertNotNull(atom3.getPoint2d());
    }

    public void testMissing3DCoordinates() throws Exception {
        String cmlString = 
          "<molecule id='m1'><atomArray><atom id='a1' xyz3='0.0 0.1 0.2'/><atom id='a2'/><atom id='a3' xyz3='0.1 0.0 0.2'/></atomArray></molecule>";
          
        IChemFile chemFile = parseCMLString(cmlString);
        org.openscience.cdk.interfaces.IMolecule mol = checkForSingleMoleculeFile(chemFile);
        
        assertEquals(3, mol.getAtomCount());
        IAtom atom1 = mol.getAtom(0);
        IAtom atom2 = mol.getAtom(1);
        IAtom atom3 = mol.getAtom(2);
        
        assertNotNull(atom1.getPoint3d());
        assertNull   (atom2.getPoint3d());
        assertNotNull(atom3.getPoint3d());
    }
    
    public void testCrystal() throws Exception {
        StringBuffer cmlStringB = new StringBuffer("  <molecule id=\"m1\">\n");
        cmlStringB.append("    <crystal z=\"4\">\n");
        cmlStringB.append("      <scalar id=\"sc1\" title=\"a\" errorValue=\"0.001\" units=\"units:angstrom\">4.500</scalar>\n");
        cmlStringB.append("      <scalar id=\"sc2\" title=\"b\" errorValue=\"0.001\" units=\"units:angstrom\">4.500</scalar>\n");
        cmlStringB.append("      <scalar id=\"sc3\" title=\"c\" errorValue=\"0.001\" units=\"units:angstrom\">4.500</scalar>\n");
        cmlStringB.append("      <scalar id=\"sc4\" title=\"alpha\" units=\"units:degrees\">90</scalar>\n");
        cmlStringB.append("      <scalar id=\"sc5\" title=\"beta\" units=\"units:degrees\">90</scalar>\n");
        cmlStringB.append("      <scalar id=\"sc6\" title=\"gamma\" units=\"units:degrees\">90</scalar>\n");
        cmlStringB.append("      <symmetry id=\"s1\" spaceGroup=\"Fm3m\"/>\n");
        cmlStringB.append("    </crystal>\n");
        cmlStringB.append("    <atomArray>\n");
        cmlStringB.append("      <atom id=\"a1\" elementType=\"Na\" formalCharge=\"1\" xyzFract=\"0.0 0.0 0.0\"\n");
        cmlStringB.append("        xy2=\"+23.1 -21.0\"></atom>\n");
        cmlStringB.append("      <atom id=\"a2\" elementType=\"Cl\" formalCharge=\"-1\" xyzFract=\"0.5 0.0 0.0\"></atom>\n");
        cmlStringB.append("    </atomArray>\n");
        cmlStringB.append("  </molecule>\n");
        
        IChemFile chemFile = parseCMLString(cmlStringB.toString());
        org.openscience.cdk.interfaces.ICrystal crystal = checkForCrystalFile(chemFile);
        assertEquals(4, crystal.getZ());
        assertEquals("Fm3m", crystal.getSpaceGroup());
        assertEquals(2, crystal.getAtomCount());
        Vector3d aaxis = crystal.getA();
        assertEquals(4.5, aaxis.x, 0.1);
        assertEquals(0.0, aaxis.y, 0.1);
        assertEquals(0.0, aaxis.z, 0.1);
        Vector3d baxis = crystal.getB();
        assertEquals(0.0, baxis.x, 0.1);
        assertEquals(4.5, baxis.y, 0.1);
        assertEquals(0.0, baxis.z, 0.1);
        Vector3d caxis = crystal.getC();
        assertEquals(0.0, caxis.x, 0.1);
        assertEquals(0.0, caxis.y, 0.1);
        assertEquals(4.5, caxis.z, 0.1);
    }

    public void testMoleculeId() throws Exception {
        String cmlString = "<molecule id='m1'><atomArray><atom id='a1'/></atomArray></molecule>";
        
        IChemFile chemFile = parseCMLString(cmlString);
        org.openscience.cdk.interfaces.IMolecule mol = checkForSingleMoleculeFile(chemFile);

        assertEquals("m1", mol.getID());
    }
    
    private IChemFile parseCMLString(String cmlString) throws Exception {
        IChemFile chemFile = null;
        CMLReader reader = new CMLReader(new ByteArrayInputStream(cmlString.getBytes()));
        chemFile = (IChemFile)reader.read(new org.openscience.cdk.ChemFile());
        return chemFile;
    }

    /**
     * Tests wether the file is indeed a single molecule file
     */
    private org.openscience.cdk.interfaces.IMolecule checkForSingleMoleculeFile(IChemFile chemFile) {
        return checkForXMoleculeFile(chemFile, 1);
    }
    
    private org.openscience.cdk.interfaces.IMolecule checkForXMoleculeFile(IChemFile chemFile, int numberOfMolecules) {
        assertNotNull(chemFile);
        
        assertEquals(chemFile.getChemSequenceCount(), 1);
        org.openscience.cdk.interfaces.IChemSequence seq = chemFile.getChemSequence(0);
        assertNotNull(seq);
        
        assertEquals(seq.getChemModelCount(), 1);
        org.openscience.cdk.interfaces.IChemModel model = seq.getChemModel(0);
        assertNotNull(model);
        
        org.openscience.cdk.interfaces.IMoleculeSet moleculeSet = model.getMoleculeSet();
        assertNotNull(moleculeSet);
        
        assertEquals(moleculeSet.getMoleculeCount(), numberOfMolecules);
        org.openscience.cdk.interfaces.IMolecule mol = null;
        for (int i=0; i<numberOfMolecules; i++) {
            mol = moleculeSet.getMolecule(i);
            assertNotNull(mol);
        }
        return mol;
    }

    private org.openscience.cdk.interfaces.ICrystal checkForCrystalFile(IChemFile chemFile) {
        assertNotNull(chemFile);
        
        assertEquals(chemFile.getChemSequenceCount(), 1);
        org.openscience.cdk.interfaces.IChemSequence seq = chemFile.getChemSequence(0);
        assertNotNull(seq);
        
        assertEquals(seq.getChemModelCount(), 1);
        org.openscience.cdk.interfaces.IChemModel model = seq.getChemModel(0);
        assertNotNull(model);
        
        org.openscience.cdk.interfaces.ICrystal crystal = model.getCrystal();
        assertNotNull(crystal);
        
        return crystal;
    }

}