[Bio] / Sprout / SHSigGenes.pm Repository:
ViewVC logotype

Diff of /Sprout/SHSigGenes.pm

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.20, Mon Jan 19 21:56:19 2009 UTC revision 1.21, Thu Feb 5 07:24:09 2009 UTC
# Line 203  Line 203 
203          }          }
204          # Only proceed if the filtering parameters are valid.          # Only proceed if the filtering parameters are valid.
205          if ($rhelp->Valid()) {          if ($rhelp->Valid()) {
             # Start the output session.  
             $self->OpenSession($rhelp);  
206              # Now we need to gather and validate the genome sets.              # Now we need to gather and validate the genome sets.
207              $self->PrintLine("Gathering the target genomes.  ");              $self->PrintLine("Gathering the target genomes.  ");
208              my ($givenGenomeID) = $self->GetGenomes('given');              my ($givenGenomeID) = $self->GetGenomes('given');
209                Trace("Given genome is $givenGenomeID.") if T(3);
210              my %targetGenomes = map { $_ => 1 } $self->GetGenomes('target');              my %targetGenomes = map { $_ => 1 } $self->GetGenomes('target');
211                Trace("Target genomes are " . join(", ", sort keys %targetGenomes) . ".") if T(3);
212              $self->PrintLine("Gathering the exclusion genomes.  ");              $self->PrintLine("Gathering the exclusion genomes.  ");
213              my %exclusionGenomes = map { $_ => 1 } $self->GetGenomes('exclusion');              my %exclusionGenomes = map { $_ => 1 } $self->GetGenomes('exclusion');
214                Trace("Exclusion genomes are " . join(", ", sort keys %exclusionGenomes) . ".") if T(3);
215              $self->PrintLine("Validating the genome sets.<br />");              $self->PrintLine("Validating the genome sets.<br />");
216              # Insure the given genome is not in the exclusion set.              # Insure the given genome is not in the exclusion set.
217              if ($exclusionGenomes{$givenGenomeID}) {              if ($exclusionGenomes{$givenGenomeID}) {
218                  $self->SetMessage("The given genome ($givenGenomeID) cannot be in the exclusion set.");                  $self->SetMessage("The given genome ($givenGenomeID) cannot be in the exclusion set.");
219              } else {              } else {
220                    # Start the output session.
221                    $self->OpenSession($rhelp);
222                  # Insure the given genome is in the target set.                  # Insure the given genome is in the target set.
223                  $targetGenomes{$givenGenomeID} = 1;                  $targetGenomes{$givenGenomeID} = 1;
224              }                  Trace("$givenGenomeID added to target set.") if T(3);
225              # Find out if we want to use a statistical analysis.              # Find out if we want to use a statistical analysis.
226              my $statistical = $cgi->param('statistical') || 1;              my $statistical = $cgi->param('statistical') || 1;
227              # Denote we have not yet found any genomes.              # Denote we have not yet found any genomes.
# Line 393  Line 396 
396              $putTimer += time() - $saveTime;              $putTimer += time() - $saveTime;
397          }          }
398      }      }
399        }
400      # Trace the timers.      # Trace the timers.
401      Trace("Time spent: Put = $putTimer, Query = $queryTimer, BBH = $bbhTimer.") if T(3);      Trace("Time spent: Put = $putTimer, Query = $queryTimer, BBH = $bbhTimer.") if T(3);
402      # Return the result count.      # Return the result count.

Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3