From 8906ee1e7c76cc694793097abbb0fefe4434bb5a Mon Sep 17 00:00:00 2001
From: jprocter <Jim Procter>
Date: Mon, 21 Aug 2006 12:43:39 +0000
Subject: [PATCH] inverse ShiftList and SequenceI.gapMap->ShiftList methods.

---
 src/jalview/util/ShiftList.java |  185 +++++++++++++++++++++++----------------
 1 file changed, 111 insertions(+), 74 deletions(-)

diff --git a/src/jalview/util/ShiftList.java b/src/jalview/util/ShiftList.java
index 905b098..78b99ef 100644
--- a/src/jalview/util/ShiftList.java
+++ b/src/jalview/util/ShiftList.java
@@ -1,74 +1,111 @@
-package jalview.util;
-
-import java.util.*;
-
-/**
- * ShiftList
- * Simple way of mapping a linear series to a new linear range with new points introduced.
- * Use at your own risk!
- * <p>Title: ShiftList</p>
- *
- * <p>Description: </p>
- *
- * <p>Copyright: Copyright (c) 2004</p>
- *
- * <p>Company: Dundee University</p>
- *
- * @author not attributable
- * @version 1.0
- */
-public class ShiftList
-{
-    Vector shifts;
-  public ShiftList()
-  {
-    shifts = new Vector();
-  }
-
-  /**
-   * addShift
-   * @param pos start position for shift (in original reference frame)
-   * @param shift length of shift
-   */
-  public void addShift(int pos, int shift)
-  {
-    int sidx = 0;
-    int[] rshift=null;
-    while (sidx<shifts.size() && (rshift=(int[]) shifts.elementAt(sidx))[0]<pos)
-      sidx++;
-    if (sidx==shifts.size())
-      shifts.insertElementAt(new int[] { pos, shift}, sidx);
-    else
-      rshift[1]+=shift;
-  }
-
-  /**
-   * shift
-   *
-   * @param pos int
-   * @return int shifted position
-   */
-  public int shift(int pos)
-  {
-    if (shifts.size()==0)
-      return pos;
-    int shifted=pos;
-    int sidx=0;
-    int rshift[];
-    while (sidx<shifts.size()
-           &&
-           (rshift=((int[]) shifts.elementAt(sidx++)))[0]<=pos) {
-      shifted += rshift[1];
-    }
-    return shifted;
-  }
-
-  /**
-   * clear all shifts
-   */
-  public void clear()
-  {
-    shifts.removeAllElements();
-  }
-
-}
+package jalview.util;
+
+import jalview.datamodel.SequenceI;
+
+import java.util.*;
+
+/**
+ * ShiftList
+ * Simple way of mapping a linear series to a new linear range with new points introduced.
+ * Use at your own risk!
+ * <p>Title: ShiftList</p>
+ *
+ * <p>Description: </p>
+ *
+ * <p>Copyright: Copyright (c) 2004</p>
+ *
+ * <p>Company: Dundee University</p>
+ *
+ * @author not attributable
+ * @version 1.0
+ */
+public class ShiftList
+{
+    public Vector shifts;
+  public ShiftList()
+  {
+    shifts = new Vector();
+  }
+
+  /**
+   * addShift
+   * @param pos start position for shift (in original reference frame)
+   * @param shift length of shift
+   */
+  public void addShift(int pos, int shift)
+  {
+    int sidx = 0;
+    int[] rshift=null;
+    while (sidx<shifts.size() && (rshift=(int[]) shifts.elementAt(sidx))[0]<pos)
+      sidx++;
+    if (sidx==shifts.size())
+      shifts.insertElementAt(new int[] { pos, shift}, sidx);
+    else
+      rshift[1]+=shift;
+  }
+
+  /**
+   * shift
+   *
+   * @param pos int
+   * @return int shifted position
+   */
+  public int shift(int pos)
+  {
+    if (shifts.size()==0)
+      return pos;
+    int shifted=pos;
+    int sidx=0;
+    int rshift[];
+    while (sidx<shifts.size()
+           &&
+           (rshift=((int[]) shifts.elementAt(sidx++)))[0]<=pos) {
+      shifted += rshift[1];
+    }
+    return shifted;
+  }
+
+  /**
+   * clear all shifts
+   */
+  public void clear()
+  {
+    shifts.removeAllElements();
+  }
+  /**
+   * invert the shifts
+   * @return ShiftList with inverse shift operations
+   */
+  public ShiftList getInverse() {
+    ShiftList inverse=new ShiftList();
+    if (shifts!=null) {
+      for (int i=0,j=shifts.size(); i<j; i++) {
+        int[] sh=(int[]) shifts.get(i);
+        if (sh!=null)
+          inverse.shifts.add(new int[] {sh[0], -sh[1]});
+      }
+    }
+    return inverse;
+  }
+
+  /**
+   * parse a 1d map of position 1<i<n to L<pos[i]<N
+   * such as that returned from SequenceI.gapMap() 
+   * @param gapMap
+   * @return
+   */
+  public static ShiftList parseMap(int[] gapMap) {
+    ShiftList shiftList = null;
+    if (gapMap!=null && gapMap.length>0) {
+      shiftList=new ShiftList();
+      
+      for (int i=0,p=0; i<gapMap.length; p++,i++) {
+        if (p!=gapMap[i]) {
+          shiftList.addShift(p, gapMap[i]-p);
+          p=gapMap[i];
+        }
+      }
+    }
+    return shiftList;
+  }
+}
-- 
1.7.10.2