#!/usr/bin/perl -sl
use strict;
use Benchmark qw|:all|;

my $ch;
my $filename = "ciona_fasta_two";
$/ = ">";


cmpthese( 10, {
              getc => sub {
		open (FH, "<$filename");
		until (eof(FH)) {
		  $ch = getc(FH);
		}
		close FH;
              },
              slurp_substr => sub {
		open (FH, "<$filename");
		my $i = 0;
		while ( <FH>) {
		  while ($ch = substr($_,$i++,1)){ }
		}
		close FH;
              },
              raw_slurp_substr => sub {
		open (FH, '<:raw',"$filename");
		my $i = 0;
		while ( <FH>) {
		  while ($ch = substr($_,$i++,1)){ }
		}
		close FH;
              },
              slurp_regex => sub {
		open (FH, "<$filename");
		while ( <FH>){
		  while ( /(.)/g){ } #the char is in $1
		}
		close FH;
              },
              raw_sysread_onechar => sub {
		open (FH, '<:raw',"$filename");
		while ( sysread( FH, $ch, 1 ) ) {}
		close FH;
              },
              nonraw_sysread_onechar => sub {
		open (FH, '<:raw',"$filename");
		while ( sysread( FH, $ch, 1 ) ) {}
		close FH;
              },
              raw_sysread_buffer => sub {
		my ($read, $buf);
		open (FH, '< :raw', "$filename");
		while ( $read = sysread( FH, $buf, 100*4096) ) {#faster than 1*4096
		  for ( 1 .. $read ) {
		    $ch = substr( $buf, $_, 1 );
		  }
		}
		close FH;
              },
              nonraw_sysread_buffer => sub {
		my ($read, $buf);
		open (FH, "<$filename");
		while ( $read = sysread( FH, $buf, 100*4096) ) {
		  for ( 1 .. $read ) {
		    $ch = substr( $buf, $_, 1 );
		  }
		}
		close FH;
              },
	     }
        );

cmpthese( 10, {
              slurp_substr => sub {
		open (FH, "<$filename");
		my $i = 0;
		while ( <FH>) {
		  while ($ch = substr($_,$i++,1)){
		  }
		}
		close FH;
              },
              slurp_simpleregex => sub {
		my $len=0;
		open (FH, "<$filename");
		while ( <FH>){
		  $_ =~ /(.)$/;
		}
		close FH;
              },
              slurp_length => sub {
		my $len=0;
		open (FH, "<$filename");
		while ( <FH>){
		  $len += length($_);
		}
		close FH;
              },
	     }

);

-----> RESULTS  ----->


                    s/iter raw_sysread_onechar nonraw_sysread_onechar raw_slurp_substr getc nonraw_sysread_buffer slurp_regex raw_sysread_buffer slurp_substr
raw_sysread_onechar    2.97   --  -0%  -4%  -19%  -46%  -51% -53% -70%
nonraw_sysread_onechar 2.95   0%   --  -3%  -19%  -46%  -51% -52% -70%
raw_slurp_substr       2.85   4%   4%   --  -16%  -44%  -49%  -51% -69%
getc                   2.40  24%  23%  19%   --  -33%   -40% -41% -63% 
nonraw_sysread_buffer  1.60  86%  85%  79%  50%    --  -10%  -12% -45%
slurp_regex            1.45 105% 104%  97%  66%   11%    --  -3%  -39%
raw_sysread_buffer     1.41 111% 110% 103%  70%   13%    3%  --   -37%
slurp_substr          0.886 235% 234% 222% 171%   80%   63%  59%   --
           

                  Rate slurp_substr slurp_length  slurp..regex
slurp_substr      1.14/s     --         -96%          -97%
slurp_length      31.2/s   2634%         --           -12%
slurp_simpleregex 35.7/s   3025%         14%            --