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.
Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
Read Where should I post X? if you're not absolutely sure you're posting in the right place.
Please read these before you post! —
Posts may use any of the Perl Monks Approved HTML tags:
- a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
Outside of code tags, you may need to use entities for some characters:
| |
For: |
|
Use: |
| & | | & |
| < | | < |
| > | | > |
| [ | | [ |
| ] | | ] |
Link using PerlMonks shortcuts! What shortcuts can I use for linking?
See Writeup Formatting Tips and other pages linked from there for more info.
|
|