WSTester updated to work plus hopefully all the other changes that need to go into...

[jabaws.git] / binaries / src / ViennaRNA / Utils / dpzoom.pl

diff --git a/binaries/src/ViennaRNA/Utils/dpzoom.pl b/binaries/src/ViennaRNA/Utils/dpzoom.pl
new file mode 100644 (file)
index 0000000..2fd0d2d
--- /dev/null
+++ b/binaries/src/ViennaRNA/Utils/dpzoom.pl
@@ -0,0 +1,46 @@
+#!/usr/bin/perl -w
+# zoom into a PostScript dot plot produced by RNAfold, RNAalifold or alidot
+# Last changed Time-stamp: <2003-03-24 12:47:26 ivo> 
+sub usage {
+  die "Usage: $0 [-f first] [-l last] [dp_file]\n";
+}
+
+$from=1;
+$to = -1;
+# parse command line
+while ($_ = shift) {
+  /^-[\?h]/ && usage();
+  /^-f$/ && ($from = shift, next );
+  /^-l$/ && ($to   = shift, next );
+  unshift (@ARGV, $_);
+  last;
+}
+
+usage() if eof();
+
+while (<>) {
+  if (/\/sequence \{ \((\S*)[\\\)]/) {
+    $seq = $1;              # empty for new version
+    while (!/\) \} def/) {  # read until end of definition
+      $_ = <>;
+      /(\S*)[\\\)]/;      # ends either with `)' or `\'
+      $seq .= $1;
+    }
+    if ($to>$from) { $seq=substr($seq,$from-1,$to-$from+1) }
+    else { $seq=substr($seq,$from-1) }
+    print  "% Subsequence from $from to $to\n";
+    print "/sequence { (\\\n";
+    for ($p=0; $p<length($seq); $p+=250)  {
+      print substr($seq,$p, 250), "\\\n";
+    }
+    print ") } def\n";
+  }
+  elsif (/(\d+) (\d+) (.*box$)/) {
+    $i=$1; $j=$2;
+    next if (($i<$from) || (($j>$to)&&($to>0)));
+    $i -= ($from-1);
+    $j -= ($from-1);
+    printf "%d %d %s\n", $i, $j, $3;
+  }
+  else {print}
+}