JAL-2089 patch broken merge to master for Release 2.10.0b1
[jalview.git] / src / jalview / gui / PairwiseAlignPanel.java
index 30effb7..22f1368 100755 (executable)
 /*
- * Jalview - A Sequence Alignment Editor and Viewer
- * Copyright (C) 2006 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle
- *
- * This program is free software; you can redistribute it and/or
- * modify it under the terms of the GNU General Public License
- * as published by the Free Software Foundation; either version 2
+ * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$)
+ * Copyright (C) $$Year-Rel$$ The Jalview Authors
+ * 
+ * This file is part of Jalview.
+ * 
+ * Jalview is free software: you can redistribute it and/or
+ * modify it under the terms of the GNU General Public License 
+ * as published by the Free Software Foundation, either version 3
  * of the License, or (at your option) any later version.
- *
- * This program is distributed in the hope that it will be useful,
- * but WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
- * GNU General Public License for more details.
- *
+ *  
+ * Jalview is distributed in the hope that it will be useful, but 
+ * WITHOUT ANY WARRANTY; without even the implied warranty 
+ * of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
+ * PURPOSE.  See the GNU General Public License for more details.
+ * 
  * You should have received a copy of the GNU General Public License
- * along with this program; if not, write to the Free Software
- * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA
+ * along with Jalview.  If not, see <http://www.gnu.org/licenses/>.
+ * The Jalview Authors are detailed in the 'AUTHORS' file.
  */
 package jalview.gui;
 
-import jalview.analysis.*;
-
-import jalview.datamodel.*;
-
-import jalview.jbgui.*;
-
-import java.awt.event.*;
-
-import java.util.*;
+import jalview.analysis.AlignSeq;
+import jalview.datamodel.Alignment;
+import jalview.datamodel.Sequence;
+import jalview.datamodel.SequenceI;
+import jalview.jbgui.GPairwiseAlignPanel;
+import jalview.util.MessageManager;
+import jalview.viewmodel.AlignmentViewport;
 
+import java.awt.event.ActionEvent;
+import java.util.Vector;
 
 /**
  * DOCUMENT ME!
- *
+ * 
  * @author $author$
  * @version $Revision$
  */
 public class PairwiseAlignPanel extends GPairwiseAlignPanel
 {
 
-    AlignViewport av;
-    Vector sequences;
+  AlignmentViewport av;
+
+  Vector sequences;
+
+  /**
+   * Creates a new PairwiseAlignPanel object.
+   * 
+   * @param av
+   *          DOCUMENT ME!
+   */
+  public PairwiseAlignPanel(AlignmentViewport av)
+  {
+    super();
+    this.av = av;
+
+    sequences = new Vector();
 
-    /**
-     * Creates a new PairwiseAlignPanel object.
-     *
-     * @param av DOCUMENT ME!
-     */
-    public PairwiseAlignPanel(AlignViewport av)
+    SequenceI[] seqs;
+    String[] seqStrings = av.getViewAsString(true);
+
+    if (av.getSelectionGroup() == null)
+    {
+      seqs = av.getAlignment().getSequencesArray();
+    }
+    else
     {
-        super();
-        this.av = av;
+      seqs = av.getSelectionGroup().getSequencesInOrder(av.getAlignment());
+    }
 
-        sequences = new Vector();
+    String type = (av.getAlignment().isNucleotide()) ? AlignSeq.DNA
+            : AlignSeq.PEP;
 
-        SequenceI [] seqs;
-        String []  seqStrings = av.getViewAsString(true);
+    float[][] scores = new float[seqs.length][seqs.length];
+    double totscore = 0;
+    int count = seqs.length;
 
-        if(av.getSelectionGroup()==null)
-        {
-          seqs = av.alignment.getSequencesArray();
-        }
-        else
-        {
-          seqs = av.getSelectionGroup().getSequencesInOrder(av.alignment);
-        }
+    Sequence seq;
 
-        String type=(av.alignment.isNucleotide()) ? AlignSeq.DNA : AlignSeq.PEP;
-        
-        float[][] scores = new float[seqs.length][seqs.length];
-        double totscore = 0;
-        int count = seqs.length;
+    for (int i = 1; i < count; i++)
+    {
+      for (int j = 0; j < i; j++)
+      {
 
-        Sequence seq;
+        AlignSeq as = new AlignSeq(seqs[i], seqStrings[i], seqs[j],
+                seqStrings[j], type);
 
-        for (int i = 1; i < count; i++)
+        if (as.s1str.length() == 0 || as.s2str.length() == 0)
         {
-            for (int j = 0; j < i; j++)
-            {
-
-                AlignSeq as = new AlignSeq(seqs[i], seqStrings[i],
-                        seqs[j], seqStrings[j], type);
-
-                if(as.s1str.length()==0 || as.s2str.length()==0)
-                {
-                  continue;
-                }
-
-                as.calcScoreMatrix();
-                as.traceAlignment();
-
-
-                as.printAlignment(System.out);
-                scores[i][j] = (float) as.getMaxScore() / (float) as.getASeq1().length;
-                totscore = totscore + scores[i][j];
-
-                textarea.append(as.getOutput());
-                seq = new Sequence(as.getS1().getName(),
-                                   as.getAStr1(),
-                                   as.getS1().getStart(),
-                                   as.getS1().getEnd()
-                );
-                sequences.add(seq);
-
-                seq = new Sequence(as.getS2().getName(),
-                                   as.getAStr2(),
-                                   as.getS2().getStart(),
-                                   as.getS2().getEnd() );
-                sequences.add(seq);
-            }
+          continue;
         }
 
-        if (count > 2)
-        {
-            System.out.println(
-                "Pairwise alignment scaled similarity score matrix\n");
-
-            for (int i = 0; i < count; i++)
-            {
-                jalview.util.Format.print(System.out, "%s \n",
-                    ("" + i) + " " +
-                    seqs[i].getName());
-            }
-
-            System.out.println("\n");
-
-            for (int i = 0; i < count; i++)
-            {
-                for (int j = 0; j < i; j++)
-                {
-                    jalview.util.Format.print(System.out, "%7.3f",
-                        scores[i][j] / totscore);
-                }
-            }
-
-            System.out.println("\n");
-        }
+        as.calcScoreMatrix();
+        as.traceAlignment();
+
+        as.printAlignment(System.out);
+        scores[i][j] = (float) as.getMaxScore()
+                / (float) as.getASeq1().length;
+        totscore = totscore + scores[i][j];
+
+        textarea.append(as.getOutput());
+        sequences.add(as.getAlignedSeq1());
+        sequences.add(as.getAlignedSeq2());
+      }
     }
 
-    /**
-     * DOCUMENT ME!
-     *
-     * @param e DOCUMENT ME!
-     */
-    protected void viewInEditorButton_actionPerformed(ActionEvent e)
+    if (count > 2)
     {
-        Sequence[] seq = new Sequence[sequences.size()];
+      System.out
+              .println("Pairwise alignment scaled similarity score matrix\n");
+
+      for (int i = 0; i < count; i++)
+      {
+        jalview.util.Format.print(System.out, "%s \n", ("" + i) + " "
+                + seqs[i].getName());
+      }
 
-        for (int i = 0; i < sequences.size(); i++)
+      System.out.println("\n");
+
+      for (int i = 0; i < count; i++)
+      {
+        for (int j = 0; j < i; j++)
         {
-            seq[i] = (Sequence) sequences.elementAt(i);
+          jalview.util.Format.print(System.out, "%7.3f", scores[i][j]
+                  / totscore);
         }
+      }
 
-        AlignFrame af = new AlignFrame(new Alignment(seq),
-                                           AlignFrame.DEFAULT_WIDTH,
-                                           AlignFrame.DEFAULT_HEIGHT);
+      System.out.println("\n");
+    }
+  }
+
+  /**
+   * DOCUMENT ME!
+   * 
+   * @param e
+   *          DOCUMENT ME!
+   */
+  protected void viewInEditorButton_actionPerformed(ActionEvent e)
+  {
+    Sequence[] seq = new Sequence[sequences.size()];
+
+    for (int i = 0; i < sequences.size(); i++)
+    {
+      seq[i] = (Sequence) sequences.elementAt(i);
+    }
 
-        Desktop.addInternalFrame(af, "Pairwise Aligned Sequences",
+    AlignFrame af = new AlignFrame(new Alignment(seq),
             AlignFrame.DEFAULT_WIDTH, AlignFrame.DEFAULT_HEIGHT);
-    }
+
+    Desktop.addInternalFrame(af,
+            MessageManager.getString("label.pairwise_aligned_sequences"),
+            AlignFrame.DEFAULT_WIDTH, AlignFrame.DEFAULT_HEIGHT);
+  }
 }