CEAlign > 0 for exact same structures ?

Jay Vyas <[email protected]>
Newsgroups gmane.comp.java.bio.general
Message-ID <[email protected]>
Hi guys.. Im finding that two structures that are the same give me non-zero
RMSD alignments....   There could be a bug in my code, but this is an
initial notice that I'm pretty sure about.... Any thoughts on the
performance of CEAlign in ideal cases ?


        Atom[] ca1 = GPdbUtils.getAtoms(s1, type);
        Atom[] ca2 = GPdbUtils.getAtoms(s2, type);

        System.out.println("Aligning two sets of atoms, " + ca1.length
+","+ca2.length);
        // get default parameters
        CeParameters params = new CeParameters();

        // set the maximum gap size to unlimited
        // params.setMaxGapSize(-1);
        StructureAlignment algorithm = StructureAlignmentFactory
                .getAlgorithm(CeMain.algorithmName);

        // The results are stored in an AFPChain object
        AFPChain afpChain = algorithm.align(ca1, ca2, params);
        afpChain.setName1("A");
        afpChain.setName2("B");

        return (float) afpChain.getChainRmsd();
_______________________________________________
Biojava-l mailing list  -  [email protected]
http://lists.open-bio.org/mailman/listinfo/biojava-l
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.