annotate icqsol_utils.py @ 3:22bb80483a8d draft

Uploaded
author greg
date Mon, 11 Jan 2016 09:21:40 -0500
parents d7c8a1ad787d
children c3b8705325ea
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
1 import os
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
2 import sys
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
3 import tempfile
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
4
3
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
5 from icqsol.shapes.icqShapeManager import ShapeManager
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
6 from icqsol.bem.icqLaplaceMatrices import LaplaceMatrices
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
7
0
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
8 PLY = 'ply'
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
9 POLYDATA = 'POLYDATA'
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
10 VTK = 'vtk'
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
11
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
12
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
13 def asbool(val):
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
14 return str(val).lower() in ['yes', 'true']
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
15
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
16
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
17 def get_format_and_type(galaxy_ext):
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
18 # Define the output file format and type.
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
19 format = None
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
20 datatype = None
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
21 if galaxy_ext in ['vtkascii', 'vtkbinary']:
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
22 format = VTK
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
23 elif galaxy_ext in ['plyascii', 'plybinary']:
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
24 format = PLY
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
25 if galaxy_ext in ['vtkascii', 'plyascii']:
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
26 datatype = 'ascii'
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
27 elif galaxy_ext in ['vtkbinary', 'plybinary']:
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
28 datatype = 'binary'
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
29 return format, datatype
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
30
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
31
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
32 def get_input_file_path(tmp_dir, input_file, format):
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
33 """
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
34 iCqSol uses file extensions (e.g., .ply, .vtk) when reading and
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
35 writing files, so the Galaxy dataset naming convention of
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
36 setting all file extensions as .dat must be handled.
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
37 """
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
38 file_path = get_temporary_file_path(tmp_dir, format)
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
39 # Remove the file so we can create a symlink.
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
40 os.remove(file_path)
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
41 os.symlink(input_file, file_path)
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
42 return file_path
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
43
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
44
3
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
45 def get_laplace_solver(shape_data, max_edge_length=float('inf')):
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
46 return LaplaceMatrices(shape_data, max_edge_length=max_edge_length)
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
47
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
48
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
49 def get_shape_manager(format, dataset_type):
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
50 # Instantiate a ShapeManager.
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
51 if format == VTK:
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
52 return ShapeManager(file_format=format, vtk_dataset_type=dataset_type)
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
53 else:
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
54 return ShapeManager(file_format=format)
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
55
22bb80483a8d Uploaded
greg
parents: 2
diff changeset
56
0
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
57 def get_temp_dir(prefix='tmp-vtk-', dir=None):
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
58 """
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
59 Return a temporary directory.
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
60 """
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
61 return tempfile.mkdtemp(prefix=prefix, dir=dir)
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
62
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
63
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
64 def get_tempfilename(dir=None, suffix=None):
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
65 """
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
66 Return a temporary file name.
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
67 """
2
d7c8a1ad787d Uploaded
greg
parents: 0
diff changeset
68 if suffix is None:
d7c8a1ad787d Uploaded
greg
parents: 0
diff changeset
69 s = None
d7c8a1ad787d Uploaded
greg
parents: 0
diff changeset
70 elif suffix.startswith('.'):
d7c8a1ad787d Uploaded
greg
parents: 0
diff changeset
71 s = suffix
d7c8a1ad787d Uploaded
greg
parents: 0
diff changeset
72 else:
d7c8a1ad787d Uploaded
greg
parents: 0
diff changeset
73 s = '.%s' % suffix
d7c8a1ad787d Uploaded
greg
parents: 0
diff changeset
74 fd, name = tempfile.mkstemp(suffix=s, dir=dir)
0
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
75 os.close(fd)
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
76 return name
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
77
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
78
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
79 def get_temporary_file_path(tmp_dir, file_extension):
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
80 """
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
81 Return the path to a temporary file with a valid VTK format
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
82 file extension.
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
83 """
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
84 return get_tempfilename(tmp_dir, file_extension)
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
85
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
86
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
87 def stop_err(msg):
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
88 sys.stderr.write("%s\n" % msg)
c8b87f2bb2e8 Uploaded
greg
parents:
diff changeset
89 sys.exit()