Warning

This document is for an old release of Galaxy. You can alternatively view this page in the latest release if it exists or view the top of the latest release's documentation.

Source code for galaxy.tool_util.cwl.parser

""" This module provides proxy objects around objects from the common
workflow language reference implementation library cwltool. These proxies
adapt cwltool to Galaxy features and abstract the library away from the rest
of the framework.
"""

import base64
import copy
import json
import logging
import os
import pickle
from abc import ABCMeta, abstractmethod
from uuid import uuid4


from galaxy.exceptions import MessageException
from galaxy.util import (
    listify,
    safe_makedirs,
    unicodify,
)
from galaxy.util.bunch import Bunch
from .cwltool_deps import (
    beta_relaxed_fmt_check,
    ensure_cwltool_available,
    getdefault,
    pathmapper,
    process,
    ref_resolver,
    relink_initialworkdir,
    RuntimeContext,
    sourceline,
    StdFsAccess,
)
from .representation import (
    field_to_field_type,
    FIELD_TYPE_REPRESENTATION,
    INPUT_TYPE,
    type_descriptions_for_field_types,
    USE_FIELD_TYPES,
    USE_STEP_PARAMETERS,
)
from .schema import (
    non_strict_non_validating_schema_loader,
    schema_loader,
)
from .util import SECONDARY_FILES_EXTRA_PREFIX

log = logging.getLogger(__name__)

JOB_JSON_FILE = ".cwl_job.json"

DOCKER_REQUIREMENT = "DockerRequirement"
SUPPORTED_TOOL_REQUIREMENTS = [
    "CreateFileRequirement",
    "DockerRequirement",
    "EnvVarRequirement",
    "InitialWorkDirRequirement",
    "InlineJavascriptRequirement",
    "ResourceRequirement",
    "ShellCommandRequirement",
    "ScatterFeatureRequirement",
    "SchemaDefRequirement",
    "SubworkflowFeatureRequirement",
    "StepInputExpressionRequirement",
    "MultipleInputFeatureRequirement",
]


SUPPORTED_WORKFLOW_REQUIREMENTS = SUPPORTED_TOOL_REQUIREMENTS + [
]

PERSISTED_REPRESENTATION = "cwl_tool_object"


[docs]def tool_proxy(tool_path=None, tool_object=None, strict_cwl_validation=True, tool_directory=None, uuid=None): """ Provide a proxy object to cwltool data structures to just grab relevant data. """ ensure_cwltool_available() tool = _to_cwl_tool_object( tool_path=tool_path, tool_object=tool_object, strict_cwl_validation=strict_cwl_validation, tool_directory=tool_directory, uuid=uuid ) return tool
[docs]def tool_proxy_from_persistent_representation(persisted_tool, strict_cwl_validation=True, tool_directory=None): """Load a ToolProxy from a previously persisted representation.""" ensure_cwltool_available() if PERSISTED_REPRESENTATION == "cwl_tool_object": kwds = {"cwl_tool_object": ToolProxy.from_persistent_representation(persisted_tool)} else: raw_process_reference = persisted_tool # ??? kwds = {"raw_process_reference": ToolProxy.from_persistent_representation(raw_process_reference)} uuid = persisted_tool["uuid"] tool = _to_cwl_tool_object(uuid=uuid, strict_cwl_validation=strict_cwl_validation, tool_directory=tool_directory, **kwds) return tool
[docs]def workflow_proxy(workflow_path, strict_cwl_validation=True): ensure_cwltool_available() workflow = _to_cwl_workflow_object(workflow_path, strict_cwl_validation=strict_cwl_validation) return workflow
[docs]def load_job_proxy(job_directory, strict_cwl_validation=True): ensure_cwltool_available() job_objects_path = os.path.join(job_directory, JOB_JSON_FILE) job_objects = json.load(open(job_objects_path)) job_inputs = job_objects["job_inputs"] output_dict = job_objects["output_dict"] # Any reason to retain older tool_path variant of this? Probably not? if "tool_path" in job_objects: tool_path = job_objects["tool_path"] cwl_tool = tool_proxy(tool_path, strict_cwl_validation=strict_cwl_validation) else: persisted_tool = job_objects["tool_representation"] cwl_tool = tool_proxy_from_persistent_representation(persisted_tool=persisted_tool, strict_cwl_validation=strict_cwl_validation) cwl_job = cwl_tool.job_proxy(job_inputs, output_dict, job_directory=job_directory) return cwl_job
def _to_cwl_tool_object(tool_path=None, tool_object=None, cwl_tool_object=None, raw_process_reference=None, strict_cwl_validation=False, tool_directory=None, uuid=None): if uuid is None: uuid = str(uuid4()) schema_loader = _schema_loader(strict_cwl_validation) if raw_process_reference is None and tool_path is not None: assert cwl_tool_object is None assert tool_object is None raw_process_reference = schema_loader.raw_process_reference(tool_path) cwl_tool = schema_loader.tool( raw_process_reference=raw_process_reference, ) elif tool_object is not None: assert raw_process_reference is None assert cwl_tool_object is None # Allow loading tools from YAML... from ruamel import yaml as ryaml as_str = json.dumps(tool_object) tool_object = ryaml.round_trip_load(as_str) path = tool_directory if path is None: path = os.getcwd() uri = ref_resolver.file_uri(path) + "/" sourceline.add_lc_filename(tool_object, uri) raw_process_reference = schema_loader.raw_process_reference_for_object( tool_object, uri=uri ) cwl_tool = schema_loader.tool( raw_process_reference=raw_process_reference, ) else: cwl_tool = cwl_tool_object if isinstance(cwl_tool, int): raise Exception("Failed to load tool.") raw_tool = cwl_tool.tool # Apply Galaxy hacks to CWL tool representation to bridge semantic differences # between Galaxy and cwltool. _hack_cwl_requirements(cwl_tool) check_requirements(raw_tool) return _cwl_tool_object_to_proxy(cwl_tool, uuid, raw_process_reference=raw_process_reference, tool_path=tool_path) def _cwl_tool_object_to_proxy(cwl_tool, uuid, raw_process_reference=None, tool_path=None): raw_tool = cwl_tool.tool if "class" not in raw_tool: raise Exception("File does not declare a class, not a valid Draft 3+ CWL tool.") process_class = raw_tool["class"] if process_class == "CommandLineTool": proxy_class = CommandLineToolProxy elif process_class == "ExpressionTool": proxy_class = ExpressionToolProxy else: raise Exception("File not a CWL CommandLineTool.") top_level_object = tool_path is not None if top_level_object and ("cwlVersion" not in raw_tool): raise Exception("File does not declare a CWL version, pre-draft 3 CWL tools are not supported.") proxy = proxy_class(cwl_tool, uuid, raw_process_reference, tool_path) return proxy def _to_cwl_workflow_object(workflow_path, strict_cwl_validation=None): proxy_class = WorkflowProxy cwl_workflow = _schema_loader(strict_cwl_validation).tool(path=workflow_path) raw_workflow = cwl_workflow.tool check_requirements(raw_workflow, tool=False) proxy = proxy_class(cwl_workflow, workflow_path) return proxy def _schema_loader(strict_cwl_validation): target_schema_loader = schema_loader if strict_cwl_validation else non_strict_non_validating_schema_loader return target_schema_loader def _hack_cwl_requirements(cwl_tool): move_to_hints = [] for i, requirement in enumerate(cwl_tool.requirements): if requirement["class"] == DOCKER_REQUIREMENT: move_to_hints.insert(0, i) for i in move_to_hints: del cwl_tool.requirements[i] cwl_tool.hints.append(requirement) def check_requirements(rec, tool=True): if isinstance(rec, dict): if "requirements" in rec: for r in rec["requirements"]: if tool: possible = SUPPORTED_TOOL_REQUIREMENTS else: possible = SUPPORTED_WORKFLOW_REQUIREMENTS if r["class"] not in possible: raise Exception("Unsupported requirement %s" % r["class"]) for d in rec: check_requirements(rec[d], tool=tool) if isinstance(rec, list): for d in rec: check_requirements(d, tool=tool) class ToolProxy(metaclass=ABCMeta): def __init__(self, tool, uuid, raw_process_reference=None, tool_path=None): self._tool = tool self._uuid = uuid self._tool_path = tool_path self._raw_process_reference = raw_process_reference # remove input parameter formats from CWL files so that cwltool # does not complain they are missing in the input data for input_field in self._tool.inputs_record_schema["fields"]: if 'format' in input_field: del input_field['format'] def job_proxy(self, input_dict, output_dict, job_directory="."): """ Build a cwltool.job.Job describing computation using a input_json Galaxy will generate mapping the Galaxy description of the inputs into a cwltool compatible variant. """ return JobProxy(self, input_dict, output_dict, job_directory=job_directory) @property def id(self): raw_id = self._tool.tool.get("id", None) return raw_id def galaxy_id(self): raw_id = self.id tool_id = None # don't reduce "search.cwl#index" to search if raw_id: tool_id = os.path.basename(raw_id) # tool_id = os.path.splitext(os.path.basename(raw_id))[0] if not tool_id: return self._uuid assert tool_id if tool_id.startswith("#"): tool_id = tool_id[1:] return tool_id @abstractmethod def input_instances(self): """ Return InputInstance objects describing mapping to Galaxy inputs. """ @abstractmethod def output_instances(self): """ Return OutputInstance objects describing mapping to Galaxy inputs. """ @abstractmethod def docker_identifier(self): """ Return docker identifier for embedding in tool description. """ @abstractmethod def description(self): """ Return description to tool. """ @abstractmethod def label(self): """ Return label for tool. """ def to_persistent_representation(self): """Return a JSON representation of this tool. Not for serialization over the wire, but serialization in a database.""" # TODO: Replace this with some more readable serialization, # I really don't like using pickle here. if PERSISTED_REPRESENTATION == "cwl_tool_object": persisted_obj = remove_pickle_problems(self._tool) else: persisted_obj = self._raw_process_reference return { "class": self._class, "pickle": unicodify(base64.b64encode(pickle.dumps(persisted_obj, pickle.HIGHEST_PROTOCOL))), "uuid": self._uuid, } @staticmethod def from_persistent_representation(as_object): """Recover an object serialized with to_persistent_representation.""" if "class" not in as_object: raise Exception("Failed to deserialize tool proxy from JSON object - no class found.") if "pickle" not in as_object: raise Exception("Failed to deserialize tool proxy from JSON object - no pickle representation found.") if "uuid" not in as_object: raise Exception("Failed to deserialize tool proxy from JSON object - no uuid found.") to_unpickle = base64.b64decode(as_object["pickle"]) loaded_object = pickle.loads(to_unpickle) return loaded_object class CommandLineToolProxy(ToolProxy): _class = "CommandLineTool" def description(self): # Don't use description - typically too verbose. return '' def doc(self): # TODO: parse multiple lines and merge - valid in cwl-1.1 doc = self._tool.tool.get('doc') return doc def label(self): label = self._tool.tool.get('label') if label is not None: return label.partition(":")[0] # return substring before ':' else: return '' def input_fields(self): input_records_schema = self._eval_schema(self._tool.inputs_record_schema) if input_records_schema["type"] != "record": raise Exception("Unhandled CWL tool input structure") # TODO: handle this somewhere else? # schemadef_req_tool_param rval = [] for input in input_records_schema["fields"]: input_copy = copy.deepcopy(input) input_type = input.get("type") if isinstance(input_type, list) or isinstance(input_type, dict): rval.append(input_copy) continue if input_type in self._tool.schemaDefs: input_copy["type"] = self._tool.schemaDefs[input_type] rval.append(input_copy) return rval def _eval_schema(self, io_schema): schema_type = io_schema.get("type") if schema_type in self._tool.schemaDefs: io_schema = self._tool.schemaDefs[schema_type] return io_schema def input_instances(self): return [_outer_field_to_input_instance(_) for _ in self.input_fields()] def output_instances(self): outputs_schema = self._eval_schema(self._tool.outputs_record_schema) if outputs_schema["type"] != "record": raise Exception("Unhandled CWL tool output structure") rval = [] for output in outputs_schema["fields"]: rval.append(_simple_field_to_output(output)) return rval def docker_identifier(self): for hint in self.hints_or_requirements_of_class("DockerRequirement"): if "dockerImageId" in hint: return hint["dockerImageId"] else: return hint["dockerPull"] return None def hints_or_requirements_of_class(self, class_name): tool = self._tool.tool reqs_and_hints = tool.get("requirements", []) + tool.get("hints", []) for hint in reqs_and_hints: if hint["class"] == class_name: yield hint def software_requirements(self): # Roughest imaginable pass at parsing requirements, really need to take in specs, handle # multiple versions, etc... tool = self._tool.tool reqs_and_hints = tool.get("requirements", []) + tool.get("hints", []) requirements = [] for hint in reqs_and_hints: if hint["class"] == "SoftwareRequirement": packages = hint.get("packages", []) for package in packages: versions = package.get("version", []) first_version = None if not versions else versions[0] requirements.append((package["package"], first_version)) return requirements @property def requirements(self): return getattr(self._tool, "requirements", []) class ExpressionToolProxy(CommandLineToolProxy): _class = "ExpressionTool" class JobProxy: def __init__(self, tool_proxy, input_dict, output_dict, job_directory): self._tool_proxy = tool_proxy self._input_dict = input_dict self._output_dict = output_dict self._job_directory = job_directory self._final_output = None self._ok = True self._cwl_job = None self._is_command_line_job = None self._normalize_job() def cwl_job(self): self._ensure_cwl_job_initialized() return self._cwl_job @property def is_command_line_job(self): self._ensure_cwl_job_initialized() assert self._is_command_line_job is not None return self._is_command_line_job def _ensure_cwl_job_initialized(self): if self._cwl_job is None: job_args = dict( basedir=self._job_directory, select_resources=self._select_resources, outdir=os.path.join(self._job_directory, "working"), tmpdir=os.path.join(self._job_directory, "cwltmp"), stagedir=os.path.join(self._job_directory, "cwlstagedir"), use_container=False, beta_relaxed_fmt_check=beta_relaxed_fmt_check, ) args = [] kwargs = {} if RuntimeContext is not None: args.append(RuntimeContext(job_args)) else: kwargs = job_args self._cwl_job = next(self._tool_proxy._tool.job( self._input_dict, self._output_callback, *args, **kwargs )) self._is_command_line_job = hasattr(self._cwl_job, "command_line") def _normalize_job(self): # Somehow reuse whatever causes validate in cwltool... maybe? def pathToLoc(p): if "location" not in p and "path" in p: p["location"] = p["path"] del p["path"] runtime_context = RuntimeContext({}) make_fs_access = getdefault(runtime_context.make_fs_access, StdFsAccess) fs_access = make_fs_access(runtime_context.basedir) process.fill_in_defaults(self._tool_proxy._tool.tool["inputs"], self._input_dict, fs_access) process.visit_class(self._input_dict, ("File", "Directory"), pathToLoc) # TODO: Why doesn't fillInDefault fill in locations instead of paths? process.normalizeFilesDirs(self._input_dict) # TODO: validate like cwltool process _init_job. # validate.validate_ex(self.names.get_name("input_record_schema", ""), builder.job, # strict=False, logger=_logger_validation_warnings) def rewrite_inputs_for_staging(self): if hasattr(self._cwl_job, "pathmapper"): pass # DO SOMETHING LIKE THE FOLLOWING? # path_rewrites = {} # for f, p in self._cwl_job.pathmapper.items(): # if not p.staged: # continue # if p.type in ("File", "Directory"): # path_rewrites[p.resolved] = p.target # for key, value in self._input_dict.items(): # if key in path_rewrites: # self._input_dict[key]["location"] = path_rewrites[value] else: stagedir = os.path.join(self._job_directory, "cwlstagedir") safe_makedirs(stagedir) def stage_recursive(value): is_list = isinstance(value, list) is_dict = isinstance(value, dict) log.info(f"handling value {value}, is_list {is_list}, is_dict {is_dict}") if is_list: for val in value: stage_recursive(val) elif is_dict: if "location" in value and "basename" in value: location = value["location"] basename = value["basename"] if not location.endswith(basename): # TODO: sep()[-1] staged_loc = os.path.join(stagedir, basename) if not os.path.exists(staged_loc): os.symlink(location, staged_loc) value["location"] = staged_loc for dict_value in value.values(): stage_recursive(dict_value) else: log.info("skipping simple value...") stage_recursive(self._input_dict) def _select_resources(self, request, runtime_context=None): new_request = request.copy() new_request["cores"] = "$GALAXY_SLOTS" return new_request @property def command_line(self): if self.is_command_line_job: return self.cwl_job().command_line else: return ["true"] @property def stdin(self): if self.is_command_line_job: return self.cwl_job().stdin else: return None @property def stdout(self): if self.is_command_line_job: return self.cwl_job().stdout else: return None @property def stderr(self): if self.is_command_line_job: return self.cwl_job().stderr else: return None @property def environment(self): if self.is_command_line_job: return self.cwl_job().environment else: return {} @property def generate_files(self): if self.is_command_line_job: return self.cwl_job().generatefiles else: return {} def _output_callback(self, out, process_status): self._process_status = process_status if process_status == "success": self._final_output = out else: self._ok = False log.info(f"Output are {out}, status is {process_status}") def collect_outputs(self, tool_working_directory, rcode): if not self.is_command_line_job: cwl_job = self.cwl_job() if RuntimeContext is not None: cwl_job.run( RuntimeContext({}) ) else: cwl_job.run() if not self._ok: raise Exception("Final process state not ok, [%s]" % self._process_status) return self._final_output else: return self.cwl_job().collect_outputs(tool_working_directory, rcode) def save_job(self): job_file = JobProxy._job_file(self._job_directory) job_objects = { # "tool_path": os.path.abspath(self._tool_proxy._tool_path), "tool_representation": self._tool_proxy.to_persistent_representation(), "job_inputs": self._input_dict, "output_dict": self._output_dict, } json.dump(job_objects, open(job_file, "w")) def _output_extra_files_dir(self, output_name): output_id = self.output_id(output_name) return os.path.join(self._job_directory, "outputs", "dataset_%s_files" % output_id) def output_id(self, output_name): output_id = self._output_dict[output_name]["id"] return output_id def output_path(self, output_name): output_id = self._output_dict[output_name]["path"] return output_id def output_directory_contents_dir(self, output_name, create=False): extra_files_dir = self._output_extra_files_dir(output_name) return extra_files_dir def output_secondary_files_dir(self, output_name, create=False): extra_files_dir = self._output_extra_files_dir(output_name) secondary_files_dir = os.path.join(extra_files_dir, SECONDARY_FILES_EXTRA_PREFIX) if create and not os.path.exists(secondary_files_dir): safe_makedirs(secondary_files_dir) return secondary_files_dir def stage_files(self): cwl_job = self.cwl_job() def stageFunc(resolved_path, target_path): log.info(f"resolving {resolved_path} to {target_path}") try: os.symlink(resolved_path, target_path) except OSError: pass if hasattr(cwl_job, "pathmapper"): process.stage_files(cwl_job.pathmapper, stageFunc, ignore_writable=True, symlink=False) if hasattr(cwl_job, "generatefiles"): outdir = os.path.join(self._job_directory, "working") # TODO: Why doesn't cwl_job.generatemapper work? generate_mapper = pathmapper.PathMapper(cwl_job.generatefiles["listing"], outdir, outdir, separateDirs=False) # TODO: figure out what inplace_update should be. inplace_update = cwl_job.inplace_update process.stage_files(generate_mapper, stageFunc, ignore_writable=inplace_update, symlink=False) relink_initialworkdir(generate_mapper, outdir, outdir, inplace_update=inplace_update) # else: expression tools do not have a path mapper. @staticmethod def _job_file(job_directory): return os.path.join(job_directory, JOB_JSON_FILE) class WorkflowProxy: def __init__(self, workflow, workflow_path=None): self._workflow = workflow self._workflow_path = workflow_path self._step_proxies = None @property def cwl_id(self): return self._workflow.tool["id"] def get_outputs_for_label(self, label): outputs = [] for output in self._workflow.tool['outputs']: step, output_name = split_step_references( output["outputSource"], multiple=False, workflow_id=self.cwl_id, ) if step == label: output_id = output["id"] if "#" not in self.cwl_id: _, output_label = output_id.rsplit("#", 1) else: _, output_label = output_id.rsplit("/", 1) outputs.append({ "output_name": output_name, "label": output_label, }) return outputs def tool_reference_proxies(self): """Fetch tool source definitions for all referenced tools.""" references = [] for step in self.step_proxies(): references.extend(step.tool_reference_proxies()) return references def step_proxies(self): if self._step_proxies is None: proxies = [] num_input_steps = len(self._workflow.tool['inputs']) for i, step in enumerate(self._workflow.steps): proxies.append(build_step_proxy(self, step, i + num_input_steps)) self._step_proxies = proxies return self._step_proxies @property def runnables(self): runnables = [] for step in self._workflow.steps: if "run" in step.tool: runnables.append(step.tool["run"]) return runnables def cwl_ids_to_index(self, step_proxies): index = 0 cwl_ids_to_index = {} for input_dict in self._workflow.tool['inputs']: cwl_ids_to_index[input_dict["id"]] = index index += 1 for step_proxy in step_proxies: cwl_ids_to_index[step_proxy.cwl_id] = index index += 1 return cwl_ids_to_index @property def output_labels(self): return [self.jsonld_id_to_label(o['id']) for o in self._workflow.tool['outputs']] def input_connections_by_step(self, step_proxies): cwl_ids_to_index = self.cwl_ids_to_index(step_proxies) input_connections_by_step = [] for step_proxy in step_proxies: input_connections_step = {} for input_proxy in step_proxy.input_proxies: cwl_source_id = input_proxy.cwl_source_id input_name = input_proxy.input_name # Consider only allow multiple if MultipleInputFeatureRequirement is enabled for (output_step_name, output_name) in split_step_references(cwl_source_id, workflow_id=self.cwl_id): if "#" in self.cwl_id: sep_on = "/" else: sep_on = "#" output_step_id = self.cwl_id + sep_on + output_step_name if output_step_id not in cwl_ids_to_index: template = "Output [%s] does not appear in ID-to-index map [%s]." msg = template % (output_step_id, cwl_ids_to_index.keys()) raise AssertionError(msg) if input_name not in input_connections_step: input_connections_step[input_name] = [] input_connections_step[input_name].append({ "id": cwl_ids_to_index[output_step_id], "output_name": output_name, "input_type": "dataset" }) input_connections_by_step.append(input_connections_step) return input_connections_by_step def to_dict(self): name = os.path.basename(self._workflow.tool.get('label') or self._workflow_path or 'TODO - derive a name from ID') steps = {} step_proxies = self.step_proxies() input_connections_by_step = self.input_connections_by_step(step_proxies) index = 0 for i, input_dict in enumerate(self._workflow.tool['inputs']): steps[index] = self.cwl_input_to_galaxy_step(input_dict, i) index += 1 for i, step_proxy in enumerate(step_proxies): input_connections = input_connections_by_step[i] steps[index] = step_proxy.to_dict(input_connections) index += 1 return { 'name': name, 'steps': steps, 'annotation': self.cwl_object_to_annotation(self._workflow.tool), } def find_inputs_step_index(self, label): for i, input in enumerate(self._workflow.tool['inputs']): if self.jsonld_id_to_label(input["id"]) == label: return i raise Exception("Failed to find index for label %s" % label) def jsonld_id_to_label(self, id): if "#" in self.cwl_id: return id.rsplit("/", 1)[-1] else: return id.rsplit("#", 1)[-1] def cwl_input_to_galaxy_step(self, input, i): input_type = input["type"] label = self.jsonld_id_to_label(input["id"]) input_as_dict = { "id": i, "label": label, "position": {"left": 0, "top": 0}, "annotation": self.cwl_object_to_annotation(input), "input_connections": {}, # Should the Galaxy API really require this? - Seems to. "workflow_outputs": self.get_outputs_for_label(label), } if input_type == "File" and "default" not in input: input_as_dict["type"] = "data_input" elif isinstance(input_type, dict) and input_type.get("type") == "array": input_as_dict["type"] = "data_collection_input" input_as_dict["collection_type"] = "list" elif isinstance(input_type, dict) and input_type.get("type") == "record": input_as_dict["type"] = "data_collection_input" input_as_dict["collection_type"] = "record" else: if USE_STEP_PARAMETERS: input_as_dict["type"] = "parameter_input" # TODO: dispatch on actual type so this doesn't always need # to be field - simpler types could be simpler inputs. tool_state = {} tool_state["parameter_type"] = "field" default_set = "default" in input default_value = input.get("default") optional = default_set if isinstance(input_type, list) and "null" in input_type: optional = True if not optional and isinstance(input_type, dict) and "type" in input_type: raise ValueError("'type' detected in non-optional input dictionary.") if default_set: tool_state["default"] = {"src": "json", "value": default_value} tool_state["optional"] = optional input_as_dict["tool_state"] = tool_state else: input_as_dict["type"] = "data_input" # TODO: format = expression.json return input_as_dict def cwl_object_to_annotation(self, cwl_obj): return cwl_obj.get("doc", None) def split_step_references(step_references, workflow_id=None, multiple=True): """Split a CWL step input or output reference into step id and name.""" # Trim off the workflow id part of the reference. step_references = listify(step_references) split_references = [] for step_reference in step_references: if workflow_id is None: # This path works fine for some simple workflows - but not so much # for subworkflows (maybe same for $graph workflows?) assert "#" in step_reference _, step_reference = step_reference.split("#", 1) else: if "#" in workflow_id: sep_on = "/" else: sep_on = "#" expected_prefix = workflow_id + sep_on if not step_reference.startswith(expected_prefix): raise AssertionError(f"step_reference [{step_reference}] doesn't start with {expected_prefix}") step_reference = step_reference[len(expected_prefix):] # Now just grab the step name and input/output name. assert "#" not in step_reference if "/" in step_reference: step_name, io_name = step_reference.split("/", 1) else: # Referencing an input, not a step. # In Galaxy workflows input steps have an implicit output named # "output" for consistency with tools - in cwl land # just the input name is referenced. step_name = step_reference io_name = "output" split_references.append((step_name, io_name)) if multiple: return split_references else: assert len(split_references) == 1 return split_references[0] def build_step_proxy(workflow_proxy, step, index): step_type = step.embedded_tool.tool["class"] if step_type == "Workflow": return SubworkflowStepProxy(workflow_proxy, step, index) else: return ToolStepProxy(workflow_proxy, step, index) class BaseStepProxy: def __init__(self, workflow_proxy, step, index): self._workflow_proxy = workflow_proxy self._step = step self._index = index self._uuid = str(uuid4()) self._input_proxies = None @property def step_class(self): return self.cwl_tool_object.tool["class"] @property def cwl_id(self): return self._step.id @property def cwl_workflow_id(self): return self._workflow_proxy.cwl_id @property def requirements(self): return self._step.requirements @property def hints(self): return self._step.hints @property def label(self): label = self._workflow_proxy.jsonld_id_to_label(self._step.id) return label def galaxy_workflow_outputs_list(self): return self._workflow_proxy.get_outputs_for_label(self.label) @property def cwl_tool_object(self): return self._step.embedded_tool @property def input_proxies(self): if self._input_proxies is None: input_proxies = [] cwl_inputs = self._step.tool["inputs"] for cwl_input in cwl_inputs: input_proxies.append(InputProxy(self, cwl_input)) self._input_proxies = input_proxies return self._input_proxies def inputs_to_dicts(self): inputs_as_dicts = [] for input_proxy in self.input_proxies: inputs_as_dicts.append(input_proxy.to_dict()) return inputs_as_dicts class InputProxy: def __init__(self, step_proxy, cwl_input): self._cwl_input = cwl_input self.step_proxy = step_proxy self.workflow_proxy = step_proxy._workflow_proxy cwl_input_id = cwl_input["id"] cwl_source_id = cwl_input.get("source", None) if cwl_source_id is None: if "valueFrom" not in cwl_input and "default" not in cwl_input: msg = "Workflow step input must define a source, a valueFrom, or a default value. Obtained [%s]." % cwl_input raise MessageException(msg) assert cwl_input_id step_name, input_name = split_step_references( cwl_input_id, multiple=False, workflow_id=step_proxy.cwl_workflow_id ) self.step_name = step_name self.input_name = input_name self.cwl_input_id = cwl_input_id self.cwl_source_id = cwl_source_id scatter_inputs = [split_step_references( i, multiple=False, workflow_id=step_proxy.cwl_workflow_id )[1] for i in listify(step_proxy._step.tool.get("scatter", []))] scatter = self.input_name in scatter_inputs self.scatter = scatter def to_dict(self): as_dict = { "name": self.input_name, } if "linkMerge" in self._cwl_input: as_dict["merge_type"] = self._cwl_input["linkMerge"] if "scatterMethod" in self.step_proxy._step.tool: as_dict["scatter_type"] = self.step_proxy._step.tool.get("scatterMethod", "dotproduct") else: as_dict["scatter_type"] = "dotproduct" if self.scatter else "disabled" if "valueFrom" in self._cwl_input: # TODO: Add a table for expressions - mark the type as CWL 1.0 JavaScript. as_dict["value_from"] = self._cwl_input["valueFrom"] if "default" in self._cwl_input: as_dict["default"] = self._cwl_input["default"] return as_dict class ToolStepProxy(BaseStepProxy): def __init__(self, workflow_proxy, step, index): super().__init__(workflow_proxy, step, index) self._tool_proxy = None @property def tool_proxy(self): # Neeeds to be cached so UUID that is loaded matches UUID generated with to_dict. if self._tool_proxy is None: self._tool_proxy = _cwl_tool_object_to_proxy(self.cwl_tool_object, uuid=str(uuid4())) return self._tool_proxy def tool_reference_proxies(self): # Return a list so we can handle subworkflows recursively. return [self.tool_proxy] def to_dict(self, input_connections): # We need to stub out null entries for things getting replaced by # connections. This doesn't seem ideal - consider just making Galaxy # handle this. tool_state = {} for input_name in input_connections.keys(): tool_state[input_name] = None outputs = self.galaxy_workflow_outputs_list() return { "id": self._index, "tool_uuid": self.tool_proxy._uuid, # TODO: make sure this is respected... "label": self.label, "position": {"left": 0, "top": 0}, "type": "tool", "annotation": self._workflow_proxy.cwl_object_to_annotation(self._step.tool), "input_connections": input_connections, "inputs": self.inputs_to_dicts(), "workflow_outputs": outputs, } class SubworkflowStepProxy(BaseStepProxy): def __init__(self, workflow_proxy, step, index): super().__init__(workflow_proxy, step, index) self._subworkflow_proxy = None def to_dict(self, input_connections): outputs = self.galaxy_workflow_outputs_list() for key, input_connection_list in input_connections.items(): input_subworkflow_step_id = self.subworkflow_proxy.find_inputs_step_index( key ) for input_connection in input_connection_list: input_connection["input_subworkflow_step_id"] = input_subworkflow_step_id return { "id": self._index, "label": self.label, "position": {"left": 0, "top": 0}, "type": "subworkflow", "subworkflow": self.subworkflow_proxy.to_dict(), "annotation": self.subworkflow_proxy.cwl_object_to_annotation(self._step.tool), "input_connections": input_connections, "inputs": self.inputs_to_dicts(), "workflow_outputs": outputs, } def tool_reference_proxies(self): return self.subworkflow_proxy.tool_reference_proxies() @property def subworkflow_proxy(self): if self._subworkflow_proxy is None: self._subworkflow_proxy = WorkflowProxy(self.cwl_tool_object) return self._subworkflow_proxy def remove_pickle_problems(obj): """doc_loader does not pickle correctly""" if hasattr(obj, "doc_loader"): obj.doc_loader = None if hasattr(obj, "embedded_tool"): obj.embedded_tool = remove_pickle_problems(obj.embedded_tool) if hasattr(obj, "steps"): obj.steps = [remove_pickle_problems(s) for s in obj.steps] return obj class WorkflowToolReference(metaclass=ABCMeta): pass class EmbeddedWorkflowToolReference(WorkflowToolReference): pass class ExternalWorkflowToolReference(WorkflowToolReference): pass def _outer_field_to_input_instance(field): field_type = field_to_field_type(field) # Must be a list if in here? if not isinstance(field_type, list): field_type = [field_type] name, label, description = _field_metadata(field) case_name = "_cwl__type_" case_label = "Specify Parameter %s As" % label def value_input(type_description): value_name = "_cwl__value_" value_label = label value_description = description return InputInstance( value_name, value_label, value_description, input_type=type_description.galaxy_param_type, collection_type=type_description.collection_type, ) select_options = [] case_options = [] type_descriptions = type_descriptions_for_field_types(field_type) for type_description in type_descriptions: select_options.append({"value": type_description.name, "label": type_description.label}) input_instances = [] if type_description.uses_param: input_instances.append(value_input(type_description)) case_options.append((type_description.name, input_instances)) # If there is more than one way to represent this parameter - produce a conditional # requesting user to ask for what form they want to submit the data in, else just map # a simple Galaxy parameter. if len(case_options) > 1 and not USE_FIELD_TYPES: case_input = SelectInputInstance( name=case_name, label=case_label, description=False, options=select_options, ) return ConditionalInstance(name, case_input, case_options) else: if len(case_options) > 1: only_type_description = FIELD_TYPE_REPRESENTATION else: only_type_description = type_descriptions[0] return InputInstance( name, label, description, input_type=only_type_description.galaxy_param_type, collection_type=only_type_description.collection_type ) # Older array to repeat handling, now we are just representing arrays as # dataset collections - we should offer a blended approach in the future. # if field_type in simple_map_type_map.keys(): # input_type = simple_map_type_map[field_type] # return {"input_type": input_type, "array": False} # elif field_type == "array": # if isinstance(field["type"], dict): # array_type = field["type"]["items"] # else: # array_type = field["items"] # if array_type in simple_map_type_map.keys(): # input_type = simple_map_type_map[array_type] # return {"input_type": input_type, "array": True} # else: # raise Exception("Unhandled simple field type encountered - [%s]." % field_type) def _field_metadata(field): name = field["name"] label = field.get("label", None) description = field.get("doc", None) return name, label, description def _simple_field_to_output(field): name = field["name"] output_data_class = field["type"] output_instance = OutputInstance( name, output_data_type=output_data_class, output_type=OUTPUT_TYPE.GLOB ) return output_instance class ConditionalInstance: def __init__(self, name, case, whens): self.input_type = INPUT_TYPE.CONDITIONAL self.name = name self.case = case self.whens = whens def to_dict(self): as_dict = dict( name=self.name, type=INPUT_TYPE.CONDITIONAL, test=self.case.to_dict(), when={}, ) for value, block in self.whens: as_dict["when"][value] = [i.to_dict() for i in block] return as_dict class SelectInputInstance: def __init__(self, name, label, description, options): self.input_type = INPUT_TYPE.SELECT self.name = name self.label = label self.description = description self.options = options def to_dict(self): # TODO: serialize options... as_dict = dict( name=self.name, label=self.label or self.name, help=self.description, type=self.input_type, options=self.options, ) return as_dict class InputInstance: def __init__(self, name, label, description, input_type, array=False, area=False, collection_type=None): self.input_type = input_type self.collection_type = collection_type self.name = name self.label = label self.description = description self.required = True self.array = array self.area = area def to_dict(self, itemwise=True): if itemwise and self.array: as_dict = dict( type="repeat", name="%s_repeat" % self.name, title="%s" % self.name, blocks=[ self.to_dict(itemwise=False) ] ) else: as_dict = dict( name=self.name, label=self.label or self.name, help=self.description, type=self.input_type, optional=not self.required, ) if self.area: as_dict["area"] = True if self.input_type == INPUT_TYPE.INTEGER: as_dict["value"] = "0" if self.input_type == INPUT_TYPE.FLOAT: as_dict["value"] = "0.0" elif self.input_type == INPUT_TYPE.DATA_COLLECTON: as_dict["collection_type"] = self.collection_type return as_dict OUTPUT_TYPE = Bunch( GLOB="glob", STDOUT="stdout", ) # TODO: Different subclasses - this is representing different types of things. class OutputInstance: def __init__(self, name, output_data_type, output_type, path=None, fields=None): self.name = name self.output_data_type = output_data_type self.output_type = output_type self.path = path self.fields = fields __all__ = ( 'tool_proxy', 'load_job_proxy', )