X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bp_bin%2Fmax_vals;h=f45b3cac2817941b55d79c9292e50cdde7d41016;hb=da9fc300ec207fce5d21e883781dab839b378fc7;hp=7c760a5b1a2bffb63050f0abd9f40468d087fa9d;hpb=b607bb9d6e4adcebd9f9e128aacc91a212010ee0;p=biopieces.git diff --git a/bp_bin/max_vals b/bp_bin/max_vals index 7c760a5..f45b3ca 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 ( $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' } = "MAX"; + + 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__ -