Mercurial > repos > peterjc > mira4_assembler
annotate tools/mira4/mira4.py @ 0:32f693f6e741 draft
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
author | peterjc |
---|---|
date | Thu, 26 Sep 2013 12:23:42 -0400 |
parents | |
children | df86ed992a1b |
rev | line source |
---|---|
0
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
1 #!/usr/bin/env python |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
2 """A simple wrapper script to call MIRA and collect its output. |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
3 """ |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
4 import os |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
5 import sys |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
6 import subprocess |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
7 import shutil |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
8 import time |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
9 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
10 WRAPPER_VER = "0.0.1" #Keep in sync with the XML file |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
11 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
12 def stop_err(msg, err=1): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
13 sys.stderr.write(msg+"\n") |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
14 sys.exit(err) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
15 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
16 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
17 def get_version(mira_binary): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
18 """Run MIRA to find its version number""" |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
19 # At the commend line I would use: mira -v | head -n 1 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
20 # however there is some pipe error when doing that here. |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
21 cmd = [mira_binary, "-v"] |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
22 try: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
23 child = subprocess.Popen(cmd, |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
24 stdout=subprocess.PIPE, |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
25 stderr=subprocess.STDOUT) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
26 except Exception, err: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
27 sys.stderr.write("Error invoking command:\n%s\n\n%s\n" % (" ".join(cmd), err)) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
28 sys.exit(1) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
29 ver, tmp = child.communicate() |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
30 del child |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
31 return ver.split("\n", 1)[0] |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
32 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
33 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
34 os.environ["PATH"] = "/mnt/galaxy/downloads/mira_4.0rc2_linux-gnu_x86_64_static/bin/:%s" % os.environ["PATH"] |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
35 mira_binary = "mira" |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
36 mira_ver = get_version(mira_binary) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
37 if not mira_ver.strip().startswith("4.0"): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
38 stop_err("This wrapper is for MIRA V4.0, not:\n%s" % mira_ver) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
39 if "-v" in sys.argv: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
40 print "MIRA wrapper version %s," % WRAPPER_VER |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
41 print mira_ver |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
42 sys.exit(0) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
43 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
44 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
45 def log_manifest(manifest): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
46 """Write the manifest file to stderr.""" |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
47 sys.stderr.write("\n%s\nManifest file\n%s\n" % ("="*60, "="*60)) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
48 with open(manifest) as h: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
49 for line in h: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
50 sys.stderr.write(line) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
51 sys.stderr.write("\n%s\nEnd of manifest\n%s\n" % ("="*60, "="*60)) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
52 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
53 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
54 def massage_symlinks(manifest): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
55 """Create FASTQ aliases and edit the manifest to use them. |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
56 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
57 Short term measure for MIRA 4.0RC2 which depends on data file |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
58 extensions to decide the file format, and doesn't like *.dat |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
59 as used in Galaxy. |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
60 """ |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
61 base = os.path.split(manifest)[0] |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
62 with open(manifest) as h: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
63 lines = h.readlines() |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
64 f = 0 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
65 for i, line in enumerate(lines): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
66 if not line.startswith("data ="): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
67 continue |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
68 #Assumes no spaces in filename, would they have to be escaped? |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
69 new_line = "data =" |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
70 for filename in line[6:].strip().split(): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
71 if not filename: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
72 continue |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
73 assert os.path.isfile(filename), filename |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
74 f += 1 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
75 alias = os.path.join(base, "input%i.fastq" % f) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
76 new_line += " " + alias |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
77 cmd = "ln -s %s %s" % (filename, alias) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
78 if os.system(cmd): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
79 stop_err("Problem creating FASTQ alias:\n%s" % cmd) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
80 lines[i] = new_line + "\n" |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
81 with open(manifest, "w") as h: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
82 for line in lines: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
83 #sys.stderr.write(line) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
84 h.write(line) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
85 return True |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
86 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
87 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
88 def collect_output(temp, name): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
89 n3 = (temp, name, name, name) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
90 f = "%s/%s_assembly/%s_d_results" % (temp, name, name) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
91 if not os.path.isdir(f): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
92 log_manifest(manifest) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
93 stop_err("Missing output folder") |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
94 if not os.listdir(f): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
95 log_manifest(manifest) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
96 stop_err("Empty output folder") |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
97 missing = [] |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
98 for old, new in [("%s/%s_out.maf" % (f, name), out_maf), |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
99 ("%s/%s_out.unpadded.fasta" % (f, name), out_fasta)]: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
100 if not os.path.isfile(old): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
101 missing.append(os.path.splitext(old)[-1]) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
102 else: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
103 shutil.move(old, new) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
104 if missing: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
105 log_manifest(manifest) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
106 sys.stderr.write("Contents of %r: %r\n" % (f, os.listdir(f))) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
107 stop_err("Missing output files: %s" % ", ".join(missing)) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
108 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
109 def clean_up(temp, name): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
110 folder = "%s/%s_assembly" % (temp, name) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
111 if os.path.isdir(folder): |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
112 shutil.rmtree(folder) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
113 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
114 #TODO - Run MIRA in /tmp or a configurable directory? |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
115 #Currently Galaxy puts us somewhere safe like: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
116 #/opt/galaxy-dist/database/job_working_directory/846/ |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
117 temp = "." |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
118 #name, out_fasta, out_qual, out_ace, out_caf, out_wig, out_log = sys.argv[1:8] |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
119 name = "MIRA" |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
120 manifest, out_maf, out_fasta, out_log = sys.argv[1:5] |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
121 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
122 #Hack until MIRA v4 lets us specify file format explicitly, |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
123 massage_symlinks(manifest) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
124 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
125 start_time = time.time() |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
126 #cmd_list =sys.argv[8:] |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
127 cmd_list = [mira_binary, manifest] |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
128 cmd = " ".join(cmd_list) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
129 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
130 assert os.path.isdir(temp) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
131 d = "%s_assembly" % name |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
132 assert not os.path.isdir(d), "Path %s already exists" % d |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
133 try: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
134 #Check path access |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
135 os.mkdir(d) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
136 except Exception, err: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
137 log_manifest(manifest) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
138 sys.stderr.write("Error making directory %s\n%s" % (d, err)) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
139 sys.exit(1) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
140 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
141 #print os.path.abspath(".") |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
142 #print cmd |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
143 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
144 handle = open(out_log, "w") |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
145 try: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
146 #Run MIRA |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
147 child = subprocess.Popen(cmd_list, |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
148 stdout=handle, |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
149 stderr=subprocess.STDOUT) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
150 except Exception, err: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
151 log_manifest(manifest) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
152 sys.stderr.write("Error invoking command:\n%s\n\n%s\n" % (cmd, err)) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
153 #TODO - call clean up? |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
154 handle.write("Error invoking command:\n%s\n\n%s\n" % (cmd, err)) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
155 handle.close() |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
156 sys.exit(1) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
157 #Use .communicate as can get deadlocks with .wait(), |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
158 stdout, stderr = child.communicate() |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
159 assert not stdout and not stderr #Should be empty as sent to handle |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
160 run_time = time.time() - start_time |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
161 return_code = child.returncode |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
162 handle.write("\n\nMIRA took %0.2f minutes\n" % (run_time / 60.0)) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
163 print "MIRA took %0.2f minutes" % (run_time / 60.0) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
164 if return_code: |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
165 handle.write("Return error code %i from command:\n" % return_code) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
166 handle.write(cmd + "\n") |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
167 handle.close() |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
168 clean_up(temp, name) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
169 log_manifest(manifest) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
170 stop_err("Return error code %i from command:\n%s" % (return_code, cmd), |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
171 return_code) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
172 handle.close() |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
173 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
174 #print "Collecting output..." |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
175 collect_output(temp, name) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
176 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
177 #print "Cleaning up..." |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
178 clean_up(temp, name) |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
179 |
32f693f6e741
Uploaded v0.0.1 preview0, very much a work in progress, primarily checking mira_datatypes dependency
peterjc
parents:
diff
changeset
|
180 print "Done" |