# # You may distribute this module under the same terms as perl itself # # POD documentation - main docs before the code =pod =head1 NAME Bio::EnsEMBL::Hive::Queen =head1 SYNOPSIS The Queen of the Hive based job control system =head1 DESCRIPTION The Queen of the Hive based job control system is responsible to 'birthing' the correct number of workers of the right type so that they can find jobs to do. It will also free up jobs of Workers that died unexpectantly so that other workers can claim them to do. Hive based processing is a concept based on a more controlled version of an autonomous agent type system. Each worker is not told what to do (like a centralized control system - like the current pipeline system) but rather queries a central database for jobs (give me jobs). Each worker is linked to an analysis_id, registers its self on creation into the Hive, creates a RunnableDB instance of the Analysis->module, gets $worker->batch_size() jobs from the analysis_job table, does its work, creates the next layer of analysis_job entries by interfacing to the DataflowRuleAdaptor to determine the analyses it needs to pass its output data to and creates jobs on the next analysis database. It repeats this cycle until it has lived its lifetime or until there are no more jobs left. The lifetime limit is just a safety limit to prevent these from 'infecting' a system. The Queens job is to simply birth Workers of the correct analysis_id to get the work down. The only other thing the Queen does is free up jobs that were claimed by Workers that died unexpectantly so that other workers can take over the work. The Beekeeper is in charge of interfacing between the Queen and a compute resource or 'compute farm'. Its job is to query Queens if they need any workers and to send the requested number of workers to open machines via the runWorker.pl script. It is also responsible for interfacing with the Queen to identify worker which died unexpectantly. =head1 CONTACT Contact Jessica Severin on EnsEMBL::Hive implemetation/design detail: jessica@ebi.ac.uk Contact Ewan Birney on EnsEMBL in general: birney@sanger.ac.uk =head1 APPENDIX The rest of the documentation details each of the object methods. Internal methods are usually preceded with a _ =cut package Bio::EnsEMBL::Hive::Queen; use strict; use POSIX; use Bio::EnsEMBL::Utils::Argument; use Bio::EnsEMBL::Utils::Exception; use Bio::EnsEMBL::Hive::Worker; use Bio::EnsEMBL::DBSQL::BaseAdaptor; use Sys::Hostname; use Bio::EnsEMBL::Hive::DBSQL::AnalysisCtrlRuleAdaptor; our @ISA = qw(Bio::EnsEMBL::DBSQL::BaseAdaptor); # # PUBLIC METHODS # ################ =head2 create_new_worker Arg [1] : $analysis_id (optional) Example : Description: If analysis_id is specified it will try to create a worker based on that analysis. If not specified the queen will analyze the hive and pick the analysis that has the most amount of work to be done. It creates an entry in the hive table, and returns a Worker object based on that insert. This guarantees that each worker registered in this queens hive is properly registered. Returntype : Bio::EnsEMBL::Hive::Worker Exceptions : Caller : =cut sub create_new_worker { my ($self, @args) = @_; my ($analysis_id, $beekeeper ,$pid) = rearrange([qw(analysis_id beekeeper process_id) ], @args); my $analStatsDBA = $self->db->get_AnalysisStatsAdaptor; return undef unless($analStatsDBA); return undef if($self->get_hive_current_load() >= 1.5); my $analysisStats; if($analysis_id) { $analysisStats = $analStatsDBA->fetch_by_analysis_id($analysis_id); } else { $analysisStats = $self->_pick_best_analysis_for_new_worker; } return undef unless($analysisStats); $analStatsDBA->decrement_needed_workers($analysisStats->analysis_id); $analysisStats->print_stats; if($analysisStats->status eq 'BLOCKED') { print("Analysis is BLOCKED, can't create workers\n"); return undef; } if($analysisStats->status eq 'DONE') { print("Analysis is DONE, don't need to create workers\n"); return undef; } my $host = hostname; $pid = $$ unless($pid); $beekeeper = '' unless($beekeeper); my $sql = q{INSERT INTO hive (born, last_check_in, process_id, analysis_id, beekeeper, host) VALUES (NOW(), NOW(), ?,?,?,?)}; my $sth = $self->prepare($sql); $sth->execute($pid, $analysisStats->analysis_id, $beekeeper, $host); my $hive_id = $sth->{'mysql_insertid'}; $sth->finish; my $worker = $self->_fetch_by_hive_id($hive_id); $worker=undef unless($worker and $worker->analysis); if($worker and $analysisStats) { $analysisStats->update_status('WORKING'); } return $worker; } sub register_worker_death { my ($self, $worker) = @_; return unless($worker); # if called without a defined cause_of_death, assume catastrophic failure $worker->cause_of_death('FATALITY') unless(defined($worker->cause_of_death)); my $sql = "UPDATE hive SET died=now(), last_check_in=now()"; $sql .= " ,work_done='" . $worker->work_done . "'"; $sql .= " ,cause_of_death='". $worker->cause_of_death ."'"; $sql .= " WHERE hive_id='" . $worker->hive_id ."'"; my $sth = $self->prepare($sql); $sth->execute(); $sth->finish; if($worker->cause_of_death eq "NO_WORK") { $self->db->get_AnalysisStatsAdaptor->update_status($worker->analysis->dbID, "ALL_CLAIMED"); } if($worker->cause_of_death eq "FATALITY") { #print("FATAL DEATH Arrrrgggghhhhhhhh (hive_id=",$worker->hive_id,")\n"); $self->db->get_AnalysisJobAdaptor->reset_dead_jobs_for_worker($worker); } # re-sync the analysis_stats when a worker dies as part of dynamic sync system if($self->synchronize_AnalysisStats($worker->analysis->stats)->status ne 'DONE') { # since I'm dying I should make sure there is someone to take my place after I'm gone ... # above synch still sees me as a 'living worker' so I need to compensate for that $self->db->get_AnalysisStatsAdaptor->increment_needed_workers($worker->analysis->dbID); } } sub worker_check_in { my ($self, $worker) = @_; return unless($worker); my $sql = "UPDATE hive SET last_check_in=now()"; $sql .= " ,work_done='" . $worker->work_done . "'"; $sql .= " WHERE hive_id='" . $worker->hive_id ."'"; my $sth = $self->prepare($sql); $sth->execute(); $sth->finish; $self->synchronize_AnalysisStats($worker->analysis->stats); } sub fetch_overdue_workers { my ($self,$overdue_secs) = @_; $overdue_secs = 3600 unless(defined($overdue_secs)); my $constraint = "h.cause_of_death='' ". "AND (UNIX_TIMESTAMP()-UNIX_TIMESTAMP(h.last_check_in))>$overdue_secs"; return $self->_generic_fetch($constraint); } =head2 synchronize_hive Example : $queen->synchronize_hive(); Description: Runs through all analyses in the system and synchronizes the analysis_stats summary with the states in the analysis_job and hive tables. Then follows by checking all the blocking rules and blocks/unblocks analyses as needed. Exceptions : none Caller : general =cut sub synchronize_hive { my $self = shift; my $start_time = time(); my $allAnalysis = $self->db->get_AnalysisAdaptor->fetch_all; print("analyze ", scalar(@$allAnalysis), "\n"); foreach my $analysis (@$allAnalysis) { my $stats = $analysis->stats; $self->synchronize_AnalysisStats($stats); } foreach my $analysis (@$allAnalysis) { $self->check_blocking_control_rules_for_AnalysisStats($analysis->stats); } print((time() - $start_time), " secs to synchronize_hive\n"); } =head2 synchronize_AnalysisStats Arg [1] : Bio::EnsEMBL::Hive::AnalysisStats object Example : $self->synchronize($analysisStats); Description: Queries the analysis_job and hive tables to get summary counts and rebuilds the AnalysisStats object. Then updates the analysis_stats table with the new summary info Returntype : newly synced Bio::EnsEMBL::Hive::AnalysisStats object Exceptions : none Caller : general =cut sub synchronize_AnalysisStats { my $self = shift; my $analysisStats = shift; return $analysisStats unless($analysisStats); return $analysisStats unless($analysisStats->analysis_id); return $analysisStats if(($analysisStats->status eq 'WORKING') and ($analysisStats->seconds_since_last_update < 3*60)); return $analysisStats if(($analysisStats->status eq 'SYNCHING') and ($analysisStats->seconds_since_last_update < 10*60)); $analysisStats->update_status('SYNCHING'); $analysisStats->total_job_count(0); $analysisStats->unclaimed_job_count(0); $analysisStats->done_job_count(0); $analysisStats->failed_job_count(0); $analysisStats->num_required_workers(0); my $sql = "SELECT status, count(*) FROM analysis_job ". "WHERE analysis_id=? GROUP BY status"; my $sth = $self->prepare($sql); $sth->execute($analysisStats->analysis_id); while (my ($status, $count)=$sth->fetchrow_array()) { my $total = $analysisStats->total_job_count(); $analysisStats->total_job_count($total + $count); if($status eq 'READY') { $analysisStats->unclaimed_job_count($count); my $numWorkers = POSIX::ceil($count * $analysisStats->avg_msec_per_job / 300000); # guess num needed workers by total jobs / (num jobs a worker could do in 5 minutes) $numWorkers=$count if($numWorkers==0); if($analysisStats->hive_capacity>0 and $numWorkers > $analysisStats->hive_capacity) { $numWorkers=$analysisStats->hive_capacity; } $analysisStats->num_required_workers($numWorkers); } if($status eq 'DONE') { $analysisStats->done_job_count($count); } if($status eq 'FAILED') { $analysisStats->failed_job_count($count); } } $sth->finish; if($analysisStats->status ne 'BLOCKED') { $analysisStats->determine_status(); } # # adjust_stats_for_living_workers # if($analysisStats->hive_capacity > 0) { my $sql = "SELECT count(*) FROM hive WHERE cause_of_death='' and analysis_id=?"; $sth = $self->prepare($sql); $sth->execute($analysisStats->analysis_id); my($liveCount)=$sth->fetchrow_array(); $sth->finish; my $numWorkers = $analysisStats->num_required_workers; my $capacityAdjust = ($numWorkers + $liveCount) - $analysisStats->hive_capacity; $numWorkers -= $capacityAdjust if($capacityAdjust > 0); $numWorkers=0 if($numWorkers<0); $analysisStats->num_required_workers($numWorkers); } $analysisStats->update; return $analysisStats; } sub check_blocking_control_rules_for_AnalysisStats { my $self = shift; my $stats = shift; return unless($stats); #print("check ctrl on analysis "); $stats->print_stats; my $ctrlRules = $self->db->get_AnalysisCtrlRuleAdaptor-> fetch_by_ctrled_analysis_id($stats->analysis_id); my $allRulesDone = 1; if(scalar @$ctrlRules > 0) { #print("HAS blocking_ctrl_rules to check\n"); foreach my $ctrlrule (@{$ctrlRules}) { #use this method because the condition_analysis objects can be #network distributed to a different database so use it's adaptor to get #the AnalysisStats object #$ctrlrule->print_rule; my $condAnalysis = $ctrlrule->condition_analysis; my $condStats = $condAnalysis->stats if($condAnalysis); $allRulesDone = 0 unless($condStats and $condStats->status eq 'DONE'); #print(" "); $condStats->print_stats; } if($allRulesDone) { if($stats->status eq 'BLOCKED') { #print(" UNBLOCK analysis : all conditions met\n"); $stats->update_status('LOADING'); #trigger sync } } else { #print(" RE-BLOCK analysis : some conditions failed\n"); $stats->update_status('BLOCKED'); } } } sub get_hive_current_load { my $self = shift; my $sql = "SELECT sum(1/analysis_stats.hive_capacity) FROM hive, analysis_stats ". "WHERE hive.analysis_id=analysis_stats.analysis_id and cause_of_death ='' ". "AND analysis_stats.hive_capacity>0"; my $sth = $self->prepare($sql); $sth->execute(); (my $load)=$sth->fetchrow_array(); $sth->finish; $load=0 unless($load); print("current hive load = $load\n"); return $load; } sub get_num_running_workers { my $self = shift; my $sql = "SELECT count(*) FROM hive WHERE cause_of_death =''"; my $sth = $self->prepare($sql); $sth->execute(); (my $runningCount)=$sth->fetchrow_array(); $sth->finish; $runningCount=0 unless($runningCount); print("current hive num_running_workers = $runningCount\n"); return $runningCount; } =head2 get_num_needed_workers Example : $count = $queen->get_num_needed_workers(); Description: Runs through the analyses in the system which are waiting for workers to be created for them. Calculates the maximum number of workers needed to fill the current needs of the system Exceptions : none Caller : general =cut sub get_num_needed_workers { my $self = shift; my $statsDBA = $self->db->get_AnalysisStatsAdaptor; my $neededAnals = $statsDBA->fetch_by_needed_workers(); my $deeper_stats_list = $statsDBA->fetch_by_status('LOADING', 'BLOCKED'); push @$neededAnals, @$deeper_stats_list; return 0 unless($neededAnals); my $availableLoad = 1.0 - $self->get_hive_current_load(); return 0 if($availableLoad <0.0); my $numWorkers = 0; foreach my $analysis_stats (@{$neededAnals}) { #$analysis_stats->print_stats(); #digging deeper under the surface so need to sync if(($analysis_stats->status eq 'LOADING') or ($analysis_stats->status eq 'BLOCKED')) { $self->synchronize_AnalysisStats($analysis_stats); $self->check_blocking_control_rules_for_AnalysisStats($analysis_stats); } next if($analysis_stats->status eq 'BLOCKED'); next if($analysis_stats->num_required_workers == 0); my $thisLoad = 0.0; if($analysis_stats->hive_capacity>0) { $thisLoad = $analysis_stats->num_required_workers * (1/$analysis_stats->hive_capacity); } if(($analysis_stats->hive_capacity<=0) or ($thisLoad < $availableLoad)) { $numWorkers += $analysis_stats->num_required_workers; $availableLoad -= $thisLoad; $analysis_stats->print_stats(); printf(" %5d workers (%1.3f remaining-hive-load)\n", $numWorkers, $availableLoad); } else { my $workerCount = POSIX::ceil($availableLoad * $analysis_stats->hive_capacity); $numWorkers += $workerCount; $availableLoad -= $workerCount * (1/$analysis_stats->hive_capacity); $analysis_stats->print_stats(); printf(" %5d workers (%1.3f remaining-hive-load) use only %3d workers\n", $numWorkers, $availableLoad, $workerCount); last; } last if($availableLoad <= 0.0); } printf("need $numWorkers workers (availLoad=%1.5f)\n", $availableLoad); return $numWorkers; } sub get_hive_progress { my $self = shift; my $sql = "SELECT sum(done_job_count ), sum(total_job_count) FROM analysis_stats"; my $sth = $self->prepare($sql); $sth->execute(); my ($done, $total) = $sth->fetchrow_array(); $sth->finish; my $completed = 0.0; $completed = ((100.0 * $done)/$total) if($total>0); printf("hive %1.3f%% complete (%d done / %d total)\n", $completed, $done, $total); return $done, $total; } sub print_hive_status { my $self = shift; $self->print_analysis_status; $self->print_running_worker_status; } sub print_analysis_status { my $self = shift; my $allStats = $self->db->get_AnalysisStatsAdaptor->fetch_all(); foreach my $analysis_stats (@{$allStats}) { $analysis_stats->print_stats; } } sub print_running_worker_status { my $self = shift; my $total = 0; print("HIVE LIVE WORKERS====\n"); my $sql = "select logic_name, count(*) from hive, analysis ". "where hive.analysis_id=analysis.analysis_id and hive.cause_of_death='' ". "group by hive.analysis_id"; my $sth = $self->prepare($sql); $sth->execute(); while((my $logic_name, my $count)=$sth->fetchrow_array()) { printf("%20s : %d workers\n", $logic_name, $count); $total += $count; } printf(" %d total workers\n", $total); print("=====================\n"); $sth->finish; } # # INTERNAL METHODS # ################### sub _pick_best_analysis_for_new_worker { my $self = shift; my $statsDBA = $self->db->get_AnalysisStatsAdaptor; return undef unless($statsDBA); my ($stats) = @{$statsDBA->fetch_by_needed_workers(1)}; return $stats if($stats); # ok so no analyses 'need' workers. # see if anything needs an update, in case there are # hidden jobs that haven't made it into the summary stats print("QUEEN: no obvious needed workers, need to dig deeper\n"); my $stats_list = $statsDBA->fetch_by_status('LOADING', 'BLOCKED'); foreach $stats (@$stats_list) { #$stats->print_stats(); $self->synchronize_AnalysisStats($stats); $self->check_blocking_control_rules_for_AnalysisStats($stats); #$stats->print_stats(); return $stats if(($stats->status eq 'READY') and ($stats->num_required_workers > 0)); } ($stats) = @{$statsDBA->fetch_by_needed_workers(1)}; return $stats if($stats); return undef; } =head2 _fetch_by_hive_id Arg [1] : int $id the unique database identifier for the feature to be obtained Example : $feat = $queen->fetch_by_dbID(1234); Description: Returns the feature created from the database defined by the the id $id. Returntype : Bio::EnsEMBL::Hive::Worker Exceptions : thrown if $id is not defined Caller : general =cut sub _fetch_by_hive_id{ my ($self,$id) = @_; unless(defined $id) { throw("fetch_by_dbID must have an id"); } my @tabs = $self->_tables; my $constraint = "h.hive_id = $id"; #return first element of _generic_fetch list my ($obj) = @{$self->_generic_fetch($constraint)}; return $obj; } =head2 _generic_fetch Arg [1] : (optional) string $constraint An SQL query constraint (i.e. part of the WHERE clause) Arg [2] : (optional) string $logic_name the logic_name of the analysis of the features to obtain Example : $fts = $a->_generic_fetch('contig_id in (1234, 1235)', 'Swall'); Description: Performs a database fetch and returns feature objects in contig coordinates. Returntype : listref of Bio::EnsEMBL::SeqFeature in contig coordinates Exceptions : none Caller : BaseFeatureAdaptor, ProxyDnaAlignFeatureAdaptor::_generic_fetch =cut sub _generic_fetch { my ($self, $constraint, $join) = @_; my @tables = $self->_tables; my $columns = join(', ', $self->_columns()); if ($join) { foreach my $single_join (@{$join}) { my ($tablename, $condition, $extra_columns) = @{$single_join}; if ($tablename && $condition) { push @tables, $tablename; if($constraint) { $constraint .= " AND $condition"; } else { $constraint = " $condition"; } } if ($extra_columns) { $columns .= ", " . join(', ', @{$extra_columns}); } } } #construct a nice table string like 'table1 t1, table2 t2' my $tablenames = join(', ', map({ join(' ', @$_) } @tables)); my $sql = "SELECT $columns FROM $tablenames"; my $default_where = $self->_default_where_clause; my $final_clause = $self->_final_clause; #append a where clause if it was defined if($constraint) { $sql .= " WHERE $constraint "; if($default_where) { $sql .= " AND $default_where "; } } elsif($default_where) { $sql .= " WHERE $default_where "; } #append additional clauses which may have been defined $sql .= " $final_clause"; my $sth = $self->prepare($sql); $sth->execute; # print STDERR $sql,"\n"; return $self->_objs_from_sth($sth); } sub _tables { my $self = shift; return (['hive', 'h']); } sub _columns { my $self = shift; return qw (h.hive_id h.analysis_id h.beekeeper h.host h.process_id h.work_done h.born h.last_check_in h.died h.cause_of_death ); } sub _objs_from_sth { my ($self, $sth) = @_; my %column; $sth->bind_columns( \( @column{ @{$sth->{NAME_lc} } } )); my @workers = (); while ($sth->fetch()) { my $worker = new Bio::EnsEMBL::Hive::Worker; $worker->init; $worker->hive_id($column{'hive_id'}); $worker->beekeeper($column{'beekeeper'}); $worker->host($column{'host'}); $worker->process_id($column{'process_id'}); $worker->work_done($column{'work_done'}); $worker->born($column{'born'}); $worker->last_check_in($column{'last_check_in'}); $worker->died($column{'died'}); $worker->cause_of_death($column{'cause_of_death'}); $worker->queen($self); $worker->db($self->db); if($column{'analysis_id'} and $self->db->get_AnalysisAdaptor) { $worker->analysis($self->db->get_AnalysisAdaptor->fetch_by_dbID($column{'analysis_id'})); } push @workers, $worker; } $sth->finish; return \@workers } sub _default_where_clause { my $self = shift; return ''; } sub _final_clause { my $self = shift; return ''; } 1;