use strict;
use warnings;
use 5.012;
use FastBioRanges;
# use PerlIO::gzip; # I just noticed I forgot this 'use' but it worked
+ fine... how?
use Storable qw (retrieve_fd);
use Time::HiRes qw(gettimeofday tv_interval);
my $time = [gettimeofday];
open( my $fastbioranges_fh, "<:gzip",
'fastbioranges.store.gz') or die;
my $fastbioranges = retrieve_fd($fastbioranges_fh) or die;
close($fastbioranges_fh);
say "loaded in ", tv_interval($time), " seconds";
my $n = 1000;
$time = [gettimeofday];
#say "start\tend\tcover";
for ( 1 .. $n ) {
my $start = int rand(877878);
my $size = int rand(7000);
my $end = ( $start + $size ) % 877879 + 1;
my $cover = $fastbioranges->num_ranges_containing( $start, $end );
# say "$start\t$end\t$cover";
}
say "$n queries in ", tv_interval($time), " seconds";
loaded in 0.385292 seconds
1000 queries in 0.005204 seconds
Recall usually my objects are 5-10 time larger and the number of queries is in the millions. The querying is not optimized (all the 'my $...', 'rand' etc.) but still it's very fast.
-
Are you posting in the right place? Check out Where do I post X? to know for sure.
-
Posts may use any of the Perl Monks Approved HTML tags. Currently these include the following:
<code> <a> <b> <big>
<blockquote> <br /> <dd>
<dl> <dt> <em> <font>
<h1> <h2> <h3> <h4>
<h5> <h6> <hr /> <i>
<li> <nbsp> <ol> <p>
<small> <strike> <strong>
<sub> <sup> <table>
<td> <th> <tr> <tt>
<u> <ul>
-
Snippets of code should be wrapped in
<code> tags not
<pre> tags. In fact, <pre>
tags should generally be avoided. If they must
be used, extreme care should be
taken to ensure that their contents do not
have long lines (<70 chars), in order to prevent
horizontal scrolling (and possible janitor
intervention).
-
Want more info? How to link
or How to display code and escape characters
are good places to start.
|