X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bp_bin%2Fmax_vals;h=52ede7c4fd16f18c5e817d37f0c9e17602dc398e;hb=f2b18307294690cf1376b51357469910b509e695;hp=b8e91b41b482492154d1651f1b0b3d0b0591166d;hpb=74bddf2fe57af278ed3a88de8b7243587024050d;p=biopieces.git diff --git a/bp_bin/max_vals b/bp_bin/max_vals index b8e91b4..52ede7c 100755 --- a/bp_bin/max_vals +++ b/bp_bin/max_vals @@ -1,4 +1,4 @@ -#!/usr/bin/env perl -w +#!/usr/bin/env perl # Copyright (C) 2007-2009 Martin A. Hansen. @@ -26,23 +26,29 @@ # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< +use warnings; use strict; use Maasha::Biopieces; +use Maasha::Calc; # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< -my ( $run_time_beg, $run_time_end, $options, $in, $out, $record, $key, $fh, %max_hash, $max_record ); +my ( $options, $in, $out, $record, $key, $fh, %max_hash, $new_record ); $options = Maasha::Biopieces::parse_options( [ - { long => 'no_stream', short => 'x', type => 'flag', mandatory => 'no', default => undef, allowed => undef, disallowed => undef }, - { long => 'data_out', short => 'o', type => 'file', mandatory => 'no', default => undef, allowed => undef, disallowed => undef }, - { long => 'keys', short => 'k', type => 'list', mandatory => 'yes', default => undef, allowed => undef, disallowed => undef }, + { long => 'no_stream', short => 'x', type => 'flag', mandatory => 'no', default => undef, allowed => undef, disallowed => undef }, + { long => 'data_out', short => 'o', type => 'file', mandatory => 'no', default => undef, allowed => undef, disallowed => undef }, + { long => 'keys', short => 'k', type => 'list', mandatory => 'no', default => undef, allowed => undef, disallowed => undef }, + { long => 'list', short => 'l', type => 'string', mandatory => 'no', default => undef, allowed => undef, disallowed => undef }, ] ); +Maasha::Common::error( qq(both --keys and --list specified) ) if $options->{ "keys" } and $options->{ "list" }; +Maasha::Common::error( qq(no --keys or --list specified) ) if not $options->{ "keys" } and not $options->{ "list" }; + $in = Maasha::Biopieces::read_stream( $options->{ "stream_in" } ); $out = Maasha::Biopieces::write_stream( $options->{ "stream_out" } ); @@ -60,17 +66,29 @@ while ( $record = Maasha::Biopieces::get_record( $in ) ) } } + if ( $options->{ 'list' } and $record->{ $options->{ 'list' } } ) { + $record->{ $options->{ 'list' } . "_MAX" } = sprintf( "%.2f", Maasha::Calc::list_max( [ split ";", $record->{ $options->{ 'list' } } ] ) ); + } + + Maasha::Biopieces::put_record( $record, $out ) if not $options->{ "no_stream" }; } +Maasha::Biopieces::close_stream( $in ); +Maasha::Biopieces::close_stream( $out ); + $fh = Maasha::Biopieces::write_stream( $options->{ "data_out" } ); -foreach $key ( @{ $options->{ "keys" } } ) -{ - $max_record->{ $key . "_MAX" } = $max_hash{ $key }; +foreach $key ( @{ $options->{ "keys" } } ) { + $new_record->{ $key . "_MAX" } = $max_hash{ $key }; } -Maasha::Biopieces::put_record( $max_record, $fh ); +if ( $options->{ "keys" } and $new_record ) +{ + $new_record->{ 'REC_TYPE' } = "MIN"; + + Maasha::Biopieces::put_record( $new_record, $fh ); +} close $fh; @@ -80,20 +98,13 @@ close $fh; BEGIN { - $run_time_beg = Maasha::Biopieces::run_time(); - - Maasha::Biopieces::log_biopiece(); + Maasha::Biopieces::status_set(); } END { - Maasha::Biopieces::close_stream( $in ); - Maasha::Biopieces::close_stream( $out ); - - $run_time_end = Maasha::Biopieces::run_time(); - - Maasha::Biopieces::run_time_print( $run_time_beg, $run_time_end, $options ); + Maasha::Biopieces::status_log(); } @@ -101,4 +112,3 @@ END __END__ -