From: martinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Date: Wed, 3 Sep 2008 05:46:38 +0000 (+0000)
Subject: fixed -S and -L for calc_fixedstep biopiece
X-Git-Url: https://git.donarmstrong.com/?a=commitdiff_plain;h=eb1f4ee65cc0a1002078b0f4db57c2d3de1619c6;p=biopieces.git

fixed -S and -L for calc_fixedstep biopiece

git-svn-id: http://biopieces.googlecode.com/svn/trunk@241 74ccb610-7750-0410-82ae-013aeee3265d
---

diff --git a/code_perl/Maasha/Biopieces.pm b/code_perl/Maasha/Biopieces.pm
index 3344774..b50b3fc 100644
--- a/code_perl/Maasha/Biopieces.pm
+++ b/code_perl/Maasha/Biopieces.pm
@@ -444,6 +444,8 @@ sub get_options
     elsif ( $script eq "calc_fixedstep" )
     {
         @options = qw(
+            score|S
+            log10|L
         );
     }
     elsif ( $script eq "transliterate_seq" )
@@ -2552,6 +2554,8 @@ sub script_calc_fixedstep
             {
                 if ( $beg > $max )
                 {
+                    map { $_ = sprintf( "%.4f", Maasha::Calc::log10( $_ ) ) } @{ $block } if $options->{ "log10" };
+
                     $record->{ "CHR" }     = $chr;
                     $record->{ "CHR_BEG" } = $beg_block;
                     $record->{ "STEP" }    = 1;
@@ -2584,6 +2588,8 @@ sub script_calc_fixedstep
 
         close $fh_in;
 
+        map { $_ = sprintf( "%.4f", Maasha::Calc::log10( $_ ) ) } @{ $block } if $options->{ "log10" };
+
         $record->{ "CHR" }     = $chr;
         $record->{ "CHR_BEG" } = $beg_block;
         $record->{ "STEP" }    = 1;