[Bio] / FortyEight / rp_peer_sim_status.pl Repository:
ViewVC logotype

View of /FortyEight/rp_peer_sim_status.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.2 - (download) (as text) (annotate)
Fri Nov 7 22:42:29 2008 UTC (11 years, 1 month ago) by olson
Branch: MAIN
CVS Tags: mgrast_dev_08112011, rast_rel_2009_05_18, mgrast_dev_08022011, rast_rel_2014_0912, mgrast_dev_05262011, rast_rel_2008_12_18, mgrast_dev_04082011, rast_rel_2010_0928, mgrast_version_3_2, mgrast_dev_12152011, mgrast_dev_06072011, rast_rel_2009_0925, rast_rel_2010_0526, rast_rel_2014_0729, mgrast_dev_02212011, rast_rel_2010_1206, mgrast_release_3_0, mgrast_dev_03252011, rast_rel_2010_0118, mgrast_rel_2008_1110_v2, rast_rel_2009_02_05, rast_rel_2011_0119, mgrast_release_3_0_4, mgrast_release_3_0_2, mgrast_release_3_0_3, mgrast_release_3_0_1, mgrast_dev_03312011, mgrast_release_3_1_2, mgrast_release_3_1_1, mgrast_release_3_1_0, mgrast_dev_04132011, mgrast_dev_04012011, rast_rel_2009_07_09, rast_rel_2010_0827, mgrast_rel_2008_1110, myrast_33, rast_rel_2011_0928, mgrast_dev_04052011, mgrast_dev_02222011, rast_rel_2009_03_26, mgrast_dev_10262011, rast_rel_2008_11_24, HEAD
Changes since 1.1: +10 -0 lines
First working versoin of the peer sims code.

#
# Determine if the peer sims between the two organisms have been computed.
#

use strict;
use Data::Dumper;
use Carp;
use DB_File;
use strict;
use FIG;
use FIG_Config;
use File::Basename;
use GenomeMeta;
use FileHandle;
use Sim;

my $meta;

@ARGV == 2 or die "Usage: $0 job-id job-id\n";

my @job_ids = @ARGV;

my @jobs;
for my $j (@job_ids)
{
    my $d = "$FIG_Config::rast_jobs/$j";
    my $g = &FIG::file_head("$d/GENOME_ID", 1);
    $g or die "cannot read jobdir $g";
    chomp $g;
    my $simdir = "$d/rp/$g";

    if (! -d "$simdir/sims")
    {
	mkdir("$simdir/sims");
    }
    if (! -f "$simdir/sims/lock")
    {
	open(LF, ">", "$simdir/sims/lock");
	close(LF);
    }

    push(@jobs, [$j, $d, $g, $simdir]);
}

#print Dumper(\@jobs);

my $ok = 0;
my $in_progress = 0;
my $missing = 0;
my $queued = 0;
for my $j1idx (0..@jobs-1)
{
    my($j1_id, $j1_jobdir, $j1_genome, $j1_orgdir) = @{$jobs[$j1idx]};
    for my $j2idx (0..@jobs-1)
    {
	next if $j1idx == $j2idx;
	my($j2_id, $j2_jobdir, $j2_genome, $j2_orgdir) = @{$jobs[$j2idx]};

	# print "check $j1_genome $j2_genome\n";

	if (-f "$j1_orgdir/sims/$j2_genome.queued")
	{
	    $queued++;
	}
	if (-f "$j1_orgdir/sims/$j2_genome.in_progress")
	{
	    $in_progress++;
	}
	elsif (-f "$j1_orgdir/sims/$j2_genome")
	{
	    $ok++;
	}
	else
	{
	    $missing++;
	}
    }
}

if ($missing)
{
    print "missing\n";
}
elsif ($in_progress)
{
    print "in_progress\n";
}
elsif ($queued)
{
    print "queued\n";
}
else
{
    print "complete\n";
}

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3