App-gapstat
view release on metacpan or search on metacpan
#$g22 { $gap } = dualvar $prev , $prev ;
$prev = dualvar $. , $_ and $neof = 0 or goto ULOOP if eof && $neof ; # and or xor ã¯ãããéããããã
} continue {
$prev = dualvar $. , $_ ; #
}
& gap1output unless ($o{G}//'') eq 0 ;
exit 0 if ($o{L}//'') eq 0 ;
sub f ($) { my $t = sprintf "%.12f" , $_[0] ; $t =~ s/\.?0*$//r } # 12æ¡ã«ãã¦ãæ«å°¾ããã®0ãåã
# åºå2. longest length 㨠max gap ####
say join "\t", map { UNDERLINE $_ } "|gap|<=" , "maxlen" , "line:content (length)" ;
for ( grep { $_ != "Inf" } sort { $a <=> $b } uniq map { abs $_ } keys %gd ) {
my @chains = defined $lc{$_} ? @{ $lc {$_} } : () ; # ããã ãã® $_ ã¤ã¾ã max abs gap ãæã¤ãã®ã®ããªã¹ããåãåºãã
#my $Length = defined $chains[0] ? $chains [0] +0 : "NA" ; # æ¬å½ã¯(ãã°ã£ã¦ãªããã°)ãã©ã®è¦ç´ ãåã£ã¦ãã¦ãè¯ãããã®æ°å¤é¨å(dualvar) ãåãåºããããã¯æé·é·ãã
my $Length = max map { $_ + 0 } @chains ;
say join "\t" , f $_ , $Length, map{ sprintf "%s (%g)", $_ , $_ } @chains ; # join ã®ä¸ã§ã¯ãããã¹ãã³ã³ãã¯ã¹ãã
}
( run in 1.259 second using v1.01-cache-2.11-cpan-5735350b133 )