X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bp_bin%2Fmax_vals;h=f45b3cac2817941b55d79c9292e50cdde7d41016;hb=5de6112b70b59420b245ce636a8b2e3c90acbe00;hp=4d42fcfc5e89febc4075d03652f3ee411e338413;hpb=5aed4676ad1bf35abbade6215ac90ee591e6d8fe;p=biopieces.git diff --git a/bp_bin/max_vals b/bp_bin/max_vals index 4d42fcf..f45b3ca 100755 --- a/bp_bin/max_vals +++ b/bp_bin/max_vals @@ -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,22 +66,31 @@ 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"; -close $fh; + Maasha::Biopieces::put_record( $new_record, $fh ); +} -Maasha::Biopieces::close_stream( $in ); -Maasha::Biopieces::close_stream( $out ); +close $fh; # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<