Search in sources :

Example 31 with Issue

use of org.vcell.util.Issue in project vcell by virtualcell.

the class ReactionRule method gatherIssues.

// -----------------------------------------------------------------------------------------------------
// 
public void gatherIssues(IssueContext issueContext, List<Issue> issueList) {
    issueContext = issueContext.newChildContext(ContextType.ReactionRule, this);
    if (name == null || name.isEmpty()) {
        issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, "Label is missing", Issue.Severity.ERROR));
    }
    if (name != null && name.startsWith("_reverse_")) {
        String msg = "The prefix '_reverse_' is a BioNetGen reserved keyword. Please rename the " + typeName + ".";
        issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, msg, Issue.Severity.ERROR));
    }
    for (MolecularType mt : model.getRbmModelContainer().getMolecularTypeList()) {
        if (mt.getName().equals(name)) {
            String msg = "Name '" + name + "' name already used for " + mt.getDisplayType() + " '" + mt.getName() + "'.";
            issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, msg, Issue.Severity.WARNING));
            // no need to look further, we won't find another match since Molecule names are unique
            break;
        }
    }
    if (reactantPatterns == null) {
        issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, "Reactant Pattern is null", Issue.Severity.ERROR));
    } else if (reactantPatterns.isEmpty()) {
        String msg = typeName + " " + name + " does not have any Reactants.\n";
        issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, msg, Issue.Severity.ERROR));
    } else {
        checkReactantPatterns(issueContext, issueList);
        for (ReactantPattern rp : reactantPatterns) {
            issueContext = issueContext.newChildContext(ContextType.ReactionRule, this);
            rp.getSpeciesPattern().gatherIssues(issueContext, issueList);
        }
    }
    if (productPatterns == null) {
        issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, "Product Pattern is null", Issue.Severity.ERROR));
    } else if (productPatterns.isEmpty()) {
        String msg = typeName + " " + name + " does not have any Products.\n";
        issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, msg, Issue.Severity.ERROR));
    } else {
        checkProductPatterns(issueContext, issueList);
        for (ProductPattern pp : productPatterns) {
            issueContext = issueContext.newChildContext(ContextType.ReactionRule, this);
            pp.getSpeciesPattern().gatherIssues(issueContext, issueList);
        }
    }
    kineticLaw.gatherIssues(issueContext, issueList);
    if (molecularTypeMappings == null) {
        issueList.add(new Issue(this, issueContext, IssueCategory.KineticsExpressionMissing, MolecularType.typeName + " Mapping is null", Issue.Severity.WARNING));
    }
    // the structure must be in the list of anchors, if anchors are being used -------------------------------
    for (ReactantPattern rp : reactantPatterns) {
        SpeciesPattern sp = rp.getSpeciesPattern();
        for (MolecularTypePattern mtp : sp.getMolecularTypePatterns()) {
            MolecularType mt = mtp.getMolecularType();
            if (mt.isAnchorAll()) {
                // no restrictions for this molecular type
                continue;
            }
            if (!mt.getAnchors().contains(rp.getStructure())) {
                String message = "The Structure " + structure.getName() + " is not allowed for the Molecule " + mt.getDisplayName();
                issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, message, Issue.Severity.ERROR));
            }
        }
    }
    for (ProductPattern pp : productPatterns) {
        SpeciesPattern sp = pp.getSpeciesPattern();
        for (MolecularTypePattern mtp : sp.getMolecularTypePatterns()) {
            MolecularType mt = mtp.getMolecularType();
            if (mt.isAnchorAll()) {
                // no restrictions for this molecular type
                continue;
            }
            if (!mt.getAnchors().contains(pp.getStructure())) {
                String message = "The Structure " + structure.getName() + " is not allowed for the Molecule " + mt.getDisplayName();
                issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, message, Issue.Severity.ERROR));
            }
        }
    }
    for (ReactantPattern rp : reactantPatterns) {
        SpeciesPattern sp = rp.getSpeciesPattern();
        for (MolecularTypePattern mtp : sp.getMolecularTypePatterns()) {
            String name = mtp.getMolecularType().getDisplayName().toLowerCase();
            if (name.equals("trash")) {
                String message = "'Trash' is a reserved NFSim keyword and it cannot be used as a reactant.";
                issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, message, Issue.Severity.WARNING));
            }
        }
    }
    if (bReversible) {
        for (ProductPattern pp : productPatterns) {
            SpeciesPattern sp = pp.getSpeciesPattern();
            for (MolecularTypePattern mtp : sp.getMolecularTypePatterns()) {
                String name = mtp.getMolecularType().getDisplayName().toLowerCase();
                if (name.equals("trash")) {
                    String message = "'Trash' is a reserved NFSim keyword and it cannot be used as a reactant.";
                    issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, message, Issue.Severity.WARNING));
                }
            }
        }
    }
    // match management
    issueContext = issueContext.newChildContext(ContextType.ReactionRule, this);
    Map<String, Integer> matchedReactants = new LinkedHashMap<String, Integer>();
    Map<String, Integer> unmatchedReactants = new LinkedHashMap<String, Integer>();
    Map<String, MolecularTypePattern> rMatches = new LinkedHashMap<String, MolecularTypePattern>();
    for (ReactantPattern rp : reactantPatterns) {
        SpeciesPattern sp = rp.getSpeciesPattern();
        for (MolecularTypePattern mtp : sp.getMolecularTypePatterns()) {
            String thisName = mtp.getMolecularType().getDisplayName();
            if (!mtp.hasExplicitParticipantMatch()) {
                if (unmatchedReactants.containsKey(thisName)) {
                    int newCounter = unmatchedReactants.get(thisName) + 1;
                    unmatchedReactants.put(thisName, newCounter);
                } else {
                    unmatchedReactants.put(thisName, 1);
                }
                continue;
            }
            // the value doesn't matter, it's important to be there
            matchedReactants.put(thisName, 99);
            String key = mtp.getParticipantMatchLabel();
            if (rMatches.containsKey(key)) {
                // no duplicates in reactants allowed
                String message = "Multiple Reactants with the same match id " + key + " are not allowed.";
                issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, message, Issue.Severity.WARNING));
            } else {
                rMatches.put(key, mtp);
            }
        }
    }
    Map<String, Integer> matchedProducts = new LinkedHashMap<String, Integer>();
    Map<String, Integer> unmatchedProducts = new LinkedHashMap<String, Integer>();
    Map<String, MolecularTypePattern> pMatches = new LinkedHashMap<String, MolecularTypePattern>();
    for (ProductPattern pp : productPatterns) {
        SpeciesPattern sp = pp.getSpeciesPattern();
        for (MolecularTypePattern mtp : sp.getMolecularTypePatterns()) {
            String thisName = mtp.getMolecularType().getDisplayName();
            if (!mtp.hasExplicitParticipantMatch()) {
                if (unmatchedProducts.containsKey(thisName)) {
                    int newCounter = unmatchedProducts.get(thisName) + 1;
                    unmatchedProducts.put(thisName, newCounter);
                } else {
                    unmatchedProducts.put(thisName, 1);
                }
                continue;
            }
            matchedProducts.put(thisName, 100);
            String key = mtp.getParticipantMatchLabel();
            if (pMatches.containsKey(key)) {
                // no duplicates in products allowed
                String message = "Multiple Products with the same match id " + key + " are not allowed.";
                issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, message, Issue.Severity.WARNING));
            } else {
                pMatches.put(key, mtp);
            }
        }
    }
    // 2+ unmatched reactants and 1+ unmatched products of same molecule mean error (or vice-versa: 1+ && 2+)
    for (String key : unmatchedReactants.keySet()) {
        if (unmatchedProducts.containsKey(key)) {
            // both maps have the same unmatched molecule, anything larger than 1 and 1 means matching error
            int rCounter = unmatchedReactants.get(key);
            int pCounter = unmatchedProducts.get(key);
            if ((rCounter > 1 && pCounter > 1) || (rCounter > 0 && pCounter > 1)) {
                String message = "Matching missing for Molecule " + key;
                issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, message, Issue.Severity.ERROR));
            }
        }
    }
    // must look also for combination with matches and unmatches for the same molecule like pA rA matched AND also pA and rA unmatched - one of each unmatched
    for (String key : unmatchedReactants.keySet()) {
        if (unmatchedProducts.containsKey(key) && matchedProducts.containsKey(key) && matchedReactants.containsKey(key)) {
            int rCounter = matchedReactants.get(key);
            int pCounter = matchedProducts.get(key);
            if (rCounter == 1 && pCounter == 1) {
                // the cases when either is >1 was addressed above, don't need duplicates
                String message = "Matching missing for Molecule " + key;
                issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, message, Issue.Severity.ERROR));
            }
        }
    }
    // unmatched products must be unambiguous (bond can't be 'possible' and state can't be 'any')
    for (ProductPattern pp : productPatterns) {
        SpeciesPattern sp = pp.getSpeciesPattern();
        for (MolecularTypePattern mtp : sp.getMolecularTypePatterns()) {
            MolecularTypePattern match = getMatchingReactantMolecule(mtp);
            if (match == null) {
                // this product has no reactant to match
                for (MolecularComponentPattern mcp : mtp.getComponentPatternList()) {
                    if (mcp.isAmbiguousBond()) {
                        String message = "Ambiguous bond '" + mcp.getBondType().name() + "' in site '" + mcp.getMolecularComponent().getDisplayName() + "' in unmatched product molecule '" + mtp.getMolecularType().getDisplayName() + "'";
                        issueList.add(new Issue(this, mcp, issueContext, IssueCategory.Identifiers, message, message, Issue.Severity.ERROR));
                    }
                    if (mcp.isAmbiguousState()) {
                        String message = "Ambiguous state 'any' in site '" + mcp.getMolecularComponent().getDisplayName() + "' in unmatched product molecule '" + mtp.getMolecularType().getDisplayName() + "'";
                        issueList.add(new Issue(this, mcp, issueContext, IssueCategory.Identifiers, message, message, Issue.Severity.ERROR));
                    }
                }
            }
        }
    }
    if (bReversible) {
        // same as above for reactants if the rule is reversible
        for (ReactantPattern rp : reactantPatterns) {
            SpeciesPattern sp = rp.getSpeciesPattern();
            for (MolecularTypePattern mtp : sp.getMolecularTypePatterns()) {
                MolecularTypePattern match = getMatchingProductMolecule(mtp);
                if (match == null) {
                    for (MolecularComponentPattern mcp : mtp.getComponentPatternList()) {
                        if (mcp.isAmbiguousBond()) {
                            String message = "Ambiguous bond '" + mcp.getBondType().name() + "' in site '" + mcp.getMolecularComponent().getDisplayName() + "' in unmatched reactant molecule '" + mtp.getMolecularType().getDisplayName() + "'";
                            issueList.add(new Issue(this, mcp, issueContext, IssueCategory.Identifiers, message, message, Issue.Severity.ERROR));
                        }
                        if (mcp.isAmbiguousState()) {
                            String message = "Ambiguous state 'any' in site '" + mcp.getMolecularComponent().getDisplayName() + "' in unmatched reactant molecule '" + mtp.getMolecularType().getDisplayName() + "'";
                            issueList.add(new Issue(this, mcp, issueContext, IssueCategory.Identifiers, message, message, Issue.Severity.ERROR));
                        }
                    }
                }
            }
        }
    }
    // invalid combinations of bonds / states for corresponding reactant/product component pairs
    for (ReactantPattern rp : reactantPatterns) {
        SpeciesPattern sp = rp.getSpeciesPattern();
        for (MolecularTypePattern mtpReactant : sp.getMolecularTypePatterns()) {
            MolecularTypePattern mtpProduct = getMatchingProductMolecule(mtpReactant);
            // mtpProduct may be null (perhaps we don't have yet any explicit or implicit match for this reactant)
            if (mtpProduct == null) {
                continue;
            }
            for (MolecularComponentPattern mcpReactant : mtpReactant.getComponentPatternList()) {
                // search for a match of the component in this product molecule
                MolecularComponentPattern mcpProduct = ReactionRule.getMatchingComponent(mtpProduct, mcpReactant);
                if (mcpProduct == null) {
                    throw new RuntimeException("Reactant Site " + mcpReactant.getDisplayName() + " is missing its matching Product Site.");
                }
                boolean incompatibleBonds = false;
                switch(mcpReactant.getBondType()) {
                    case Exists:
                        switch(mcpProduct.getBondType()) {
                            case Exists:
                                break;
                            case None:
                                incompatibleBonds = true;
                                break;
                            case Possible:
                                incompatibleBonds = true;
                                break;
                            case Specified:
                                incompatibleBonds = true;
                                break;
                        }
                        break;
                    case None:
                        switch(mcpProduct.getBondType()) {
                            case Exists:
                                incompatibleBonds = true;
                                break;
                            case None:
                                break;
                            case Possible:
                                incompatibleBonds = true;
                                break;
                            case Specified:
                                break;
                        }
                        break;
                    case Possible:
                        switch(mcpProduct.getBondType()) {
                            case Exists:
                                incompatibleBonds = true;
                                break;
                            case None:
                                incompatibleBonds = true;
                                break;
                            case Possible:
                                break;
                            case Specified:
                                incompatibleBonds = true;
                                break;
                        }
                        break;
                    case Specified:
                        switch(mcpProduct.getBondType()) {
                            case Exists:
                                incompatibleBonds = true;
                                break;
                            case None:
                                break;
                            case Possible:
                                incompatibleBonds = true;
                                break;
                            case Specified:
                                break;
                        }
                        break;
                }
                // ----- end of bonds switch
                if (incompatibleBonds) {
                    String message = "Reactant " + MolecularComponentPattern.typeName + " " + mcpReactant.getDisplayName();
                    message += " and its matching Product " + MolecularComponentPattern.typeName + " " + mcpProduct.getDisplayName() + " have incompatible Bond Types.";
                    issueList.add(new Issue(this, mcpProduct, issueContext, IssueCategory.Identifiers, message, message, Issue.Severity.ERROR));
                }
                boolean incompatibleStates = true;
                if (mcpReactant.getMolecularComponent().getComponentStateDefinitions().isEmpty()) {
                // nothing to do if there are no States defined
                } else {
                    ComponentStatePattern cspReactant = mcpReactant.getComponentStatePattern();
                    ComponentStatePattern cspProduct = mcpProduct.getComponentStatePattern();
                    if (cspReactant == null) {
                        String message = "Required " + ComponentStatePattern.typeName + " missing for " + MolecularComponentPattern.typeName + " " + mcpReactant.getDisplayName();
                        issueList.add(new Issue(this, mcpReactant, issueContext, IssueCategory.Identifiers, message, message, Issue.Severity.ERROR));
                    }
                    if (cspProduct == null) {
                        String message = "Required " + ComponentStatePattern.typeName + " missing for " + MolecularComponentPattern.typeName + " " + mcpProduct.getDisplayName();
                        issueList.add(new Issue(this, mcpProduct, issueContext, IssueCategory.Identifiers, message, message, Issue.Severity.ERROR));
                    }
                    if (cspReactant.isAny() && cspProduct.isAny()) {
                        incompatibleStates = false;
                    } else if (!cspReactant.isAny() && !cspProduct.isAny()) {
                        incompatibleStates = false;
                    }
                    if (incompatibleStates) {
                        String message = "Reactant " + MolecularComponentPattern.typeName + " " + mcpReactant.getDisplayName();
                        message += " and its matching Product " + MolecularComponentPattern.typeName + " " + mcpProduct.getDisplayName() + " have incompatible States";
                        issueList.add(new Issue(this, mcpProduct, issueContext, IssueCategory.Identifiers, message, message, Issue.Severity.ERROR));
                    }
                }
            }
        }
    }
}
Also used : Issue(org.vcell.util.Issue) MolecularComponentPattern(org.vcell.model.rbm.MolecularComponentPattern) ComponentStatePattern(org.vcell.model.rbm.ComponentStatePattern) SpeciesPattern(org.vcell.model.rbm.SpeciesPattern) LinkedHashMap(java.util.LinkedHashMap) MolecularType(org.vcell.model.rbm.MolecularType) MolecularTypePattern(org.vcell.model.rbm.MolecularTypePattern)

Example 32 with Issue

use of org.vcell.util.Issue in project vcell by virtualcell.

the class SpeciesContext method checkComponentStateConsistency.

public void checkComponentStateConsistency(IssueContext issueContext, List<Issue> issueList, MolecularTypePattern mtpThis) {
    issueContext = issueContext.newChildContext(ContextType.SpeciesContext, this);
    MolecularType mtThat = mtpThis.getMolecularType();
    for (MolecularComponentPattern mcpThis : mtpThis.getComponentPatternList()) {
        ComponentStatePattern cspThis = mcpThis.getComponentStatePattern();
        String mcNameThis = mcpThis.getMolecularComponent().getName();
        MolecularComponent[] mcThatList = mtThat.getMolecularComponents(mcNameThis);
        if (mcThatList.length == 0) {
            System.out.println("we already fired an issue about component missing");
            // nothing to do here, we already fired an issue about component missing
            continue;
        } else if (mcThatList.length > 1) {
            String msg = "Multiple " + MolecularComponent.typeName + "s with the same name are not yet supported.";
            issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, msg, Issue.SEVERITY_ERROR));
        } else {
            // found exactly 1 component
            MolecularComponent mcThat = mcThatList[0];
            List<ComponentStateDefinition> csdListThat = mcThat.getComponentStateDefinitions();
            if (csdListThat.size() == 0) {
                // component has no states, we check if mcpThis has any states... it shouldn't
                if (cspThis == null) {
                    // all is well
                    continue;
                }
                if (!cspThis.isAny() || (cspThis.getComponentStateDefinition() != null)) {
                    String msg = MolecularComponentPattern.typeName + " " + mcNameThis + " is in an invalid State.";
                    issueList.add(new Issue(this, mcpThis, issueContext, IssueCategory.Identifiers, msg, null, Issue.SEVERITY_WARNING));
                }
            } else {
                // we check if mcpThis has any of these states... it should!
                if ((cspThis == null) || cspThis.isAny() || (cspThis.getComponentStateDefinition() == null)) {
                    String msg = MolecularComponentPattern.typeName + " " + mcNameThis + " must be in an explicit State.";
                    issueList.add(new Issue(this, mcpThis, issueContext, IssueCategory.Identifiers, msg, null, Issue.Severity.ERROR));
                } else {
                    String csdNameThis = cspThis.getComponentStateDefinition().getName();
                    if (csdNameThis.isEmpty() || (mcThat.getComponentStateDefinition(csdNameThis) == null)) {
                        String msg = "Invalid State " + csdNameThis + " for " + MolecularComponentPattern.typeName + " " + mcNameThis;
                        issueList.add(new Issue(this, mcpThis, issueContext, IssueCategory.Identifiers, msg, null, Issue.SEVERITY_WARNING));
                    }
                }
            }
        }
    }
}
Also used : MolecularType(org.vcell.model.rbm.MolecularType) Issue(org.vcell.util.Issue) MolecularComponentPattern(org.vcell.model.rbm.MolecularComponentPattern) MolecularComponent(org.vcell.model.rbm.MolecularComponent) ComponentStatePattern(org.vcell.model.rbm.ComponentStatePattern) List(java.util.List)

Example 33 with Issue

use of org.vcell.util.Issue in project vcell by virtualcell.

the class SpeciesContext method gatherIssues.

public void gatherIssues(IssueContext issueContext, List<Issue> issueList) {
    issueContext = issueContext.newChildContext(ContextType.SpeciesContext, this);
    if (species == null) {
        issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, "Species is null", Issue.SEVERITY_WARNING));
    } else {
        if (speciesPattern != null) {
            checkBondsSufficiency(issueContext, issueList, speciesPattern);
            speciesPattern.checkSpeciesPattern(issueContext, issueList);
            speciesPattern.gatherIssues(issueContext, issueList);
            for (MolecularType mtThat : model.getRbmModelContainer().getMolecularTypeList()) {
                for (MolecularTypePattern mtpThis : speciesPattern.getMolecularTypePatterns()) {
                    if (mtThat.getName().equals(mtpThis.getMolecularType().getName())) {
                        // this should never fire issues!!!
                        checkMolecularTypeConsistency(issueContext, issueList, mtThat, mtpThis);
                        if (!mtThat.compareEqual(mtpThis.getMolecularType())) {
                            String msg = MolecularType.typeName + "s " + mtThat.getName() + " and " + mtpThis.getMolecularType().getName() + " do not compare equal.";
                            issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, msg, Issue.SEVERITY_WARNING));
                        } else {
                        // TODO: components cannot be missing, see if this test is really needed
                        // TODO: the isImplied test is wrong here
                        // String msg = "All components in the " + mtThat.getDisplayType() + " definition must be present. Missing:";
                        // boolean isInvalid = false;
                        // for(int i=0; i< mtpThis.getComponentPatternList().size(); i++) {
                        // MolecularComponentPattern mcp = mtpThis.getComponentPatternList().get(i);
                        // if(mcp.isImplied()) {
                        // if(isInvalid) {
                        // msg += ",";
                        // }
                        // isInvalid = true;
                        // msg += " " + mcp.getMolecularComponent().getName();
                        // }
                        // }
                        // if(isInvalid) {
                        // msg += ".";
                        // issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, msg, Issue.Severity.ERROR));
                        // }
                        }
                    }
                }
            }
            // if the component has states, the component pattern must have ONE of those states
            for (MolecularTypePattern mtpThis : speciesPattern.getMolecularTypePatterns()) {
                checkComponentStateConsistency(issueContext, issueList, mtpThis);
            }
            // the bond must not be ambiguous - seed species must be either unbound or specified (explicit) bond
            for (MolecularTypePattern mtpThis : speciesPattern.getMolecularTypePatterns()) {
                for (MolecularComponentPattern mcpThis : mtpThis.getComponentPatternList()) {
                    if (mcpThis.getBondType() == BondType.Possible || mcpThis.getBondType() == BondType.Exists) {
                        String msg = "The Bond of " + mcpThis.getMolecularComponent().getDisplayName() + " must be 'Unbound' or explicit.";
                        issueList.add(new Issue(this, mcpThis, issueContext, IssueCategory.Identifiers, msg, null, Issue.Severity.ERROR));
                    }
                }
            }
            // the structure must be in the list of anchors, if anchors are being used
            for (MolecularTypePattern mtp : speciesPattern.getMolecularTypePatterns()) {
                MolecularType mt = mtp.getMolecularType();
                if (mt.isAnchorAll()) {
                    // no restrictions for this molecular type
                    continue;
                }
                if (!mt.getAnchors().contains(structure)) {
                    String message = "The Structure " + structure.getName() + " is not allowed for the Molecule " + mt.getDisplayName();
                    issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, message, Issue.Severity.ERROR));
                }
            }
            // check for reserved name 'trash'
            for (MolecularTypePattern mtp : speciesPattern.getMolecularTypePatterns()) {
                String name = mtp.getMolecularType().getDisplayName().toLowerCase();
                if (name.equals("trash")) {
                    String message = "'Trash' is a reserved NFSim keyword and it cannot be used as a seed species.";
                    issueList.add(new Issue(this, issueContext, IssueCategory.Identifiers, message, Issue.Severity.WARNING));
                }
            }
        }
    }
}
Also used : MolecularType(org.vcell.model.rbm.MolecularType) Issue(org.vcell.util.Issue) MolecularComponentPattern(org.vcell.model.rbm.MolecularComponentPattern) MolecularTypePattern(org.vcell.model.rbm.MolecularTypePattern)

Example 34 with Issue

use of org.vcell.util.Issue in project vcell by virtualcell.

the class WorkflowFileTest method main.

public static void main(String[] args) {
    if (args.length != 2) {
        System.out.println("expecting 2 arguments");
        System.out.println("usage: java " + Workflow.class.getSimpleName() + " workingdir workflowInputFile");
        System.out.println("workingdir example: " + "D:\\developer\\eclipse\\workspace\\VCell_5.4_vmicro\\datadir");
        System.out.println("workflowInputFile example: " + "D:\\developer\\eclipse\\workspace\\VCell_5.4_vmicro\\workflow1.txt");
        System.exit(1);
    }
    try {
        File workingDirectory = new File(args[0]);
        LocalWorkspace localWorkspace = new LocalWorkspace(workingDirectory);
        String workflowLanguageText = BeanUtils.readBytesFromFile(new File(args[1]), null);
        Repository repository = new MemoryRepository();
        Workflow workflow = Workflow.parse(repository, localWorkspace, workflowLanguageText);
        TaskContext taskContext = new TaskContext(workflow, repository, localWorkspace);
        WorkflowUtilities.displayWorkflowGraph(taskContext.getWorkflow());
        WorkflowUtilities.displayWorkflowTable(taskContext);
        WorkflowUtilities.displayWorkflowGraphJGraphX(new WorkflowJGraphProxy(taskContext));
        ArrayList<Issue> issues = new ArrayList<Issue>();
        taskContext.getWorkflow().reportIssues(issues, Issue.SEVERITY_INFO, true);
        // 
        // execute the workflow
        // 
        taskContext.getWorkflow().compute(taskContext, new WorkflowUtilities.Progress());
    } catch (Exception e) {
        e.printStackTrace(System.out);
    }
}
Also used : LocalWorkspace(org.vcell.vmicro.workflow.data.LocalWorkspace) TaskContext(org.vcell.workflow.TaskContext) Issue(org.vcell.util.Issue) ArrayList(java.util.ArrayList) Workflow(org.vcell.workflow.Workflow) WorkflowJGraphProxy(org.vcell.vmicro.workflow.jgraphx.WorkflowJGraphProxy) Repository(org.vcell.workflow.Repository) MemoryRepository(org.vcell.workflow.MemoryRepository) MemoryRepository(org.vcell.workflow.MemoryRepository) File(java.io.File)

Example 35 with Issue

use of org.vcell.util.Issue in project vcell by virtualcell.

the class WorkflowTest method main.

public static void main(String[] args) {
    if (args.length != 2) {
        System.out.println("expecting 2 arguments");
        System.out.println("usage: java " + Workflow.class.getSimpleName() + " workingdir workflowInputFile");
        System.out.println("workingdir example: " + "D:\\developer\\eclipse\\workspace\\VCell_5.4_vmicro\\datadir");
        System.out.println("workflowInputFile example: " + "D:\\developer\\eclipse\\workspace\\VCell_5.4_vmicro\\workflow1.txt");
        System.exit(1);
    }
    try {
        // PropertyLoader.loadProperties();
        // workflowInputFile "C:\\developer\\eclipse\\workspace\\VCell_5.3_vmicro\\workflow1.txt"
        File workingDirectory = new File(args[0]);
        LocalWorkspace localWorkspace = new LocalWorkspace(workingDirectory);
        Repository repository = new MemoryRepository();
        String workflowLanguageText = BeanUtils.readBytesFromFile(new File(args[1]), null);
        Workflow workflow = Workflow.parse(repository, localWorkspace, workflowLanguageText);
        System.err.println(">>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> using hard-coded example instead <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<");
        // File vfrapFile = new File("D:\\Developer\\eclipse\\workspace_refactor\\VCell_5.4_vmicro\\3D_FRAP_2_ZProjection_Simulation1.vfrap");
        // Workflow workflow = getVFrapSimpleExample(workingDirectory, vfrapFile);
        TaskContext taskContext = new TaskContext(workflow, repository, localWorkspace);
        ArrayList<Issue> issues = new ArrayList<Issue>();
        IssueContext issueContext = new IssueContext();
        workflow.gatherIssues(issueContext, issues);
        // WorkflowJGraphProxy workflowJGraphProxy = new WorkflowJGraphProxy(workflow);
        // displayWorkflowGraphJGraphX(workflowJGraphProxy);
        WorkflowUtilities.displayWorkflowGraph(workflow);
        WorkflowUtilities.displayWorkflowTable(taskContext);
        workflow.reportIssues(issues, Issue.SEVERITY_INFO, true);
        // 
        // execute the workflow
        // 
        workflow.compute(taskContext, new WorkflowUtilities.Progress());
    } catch (Exception e) {
        e.printStackTrace(System.out);
    }
}
Also used : LocalWorkspace(org.vcell.vmicro.workflow.data.LocalWorkspace) TaskContext(org.vcell.workflow.TaskContext) Issue(org.vcell.util.Issue) ArrayList(java.util.ArrayList) Workflow(org.vcell.workflow.Workflow) MemoryRepository(org.vcell.workflow.MemoryRepository) Repository(org.vcell.workflow.Repository) MemoryRepository(org.vcell.workflow.MemoryRepository) IssueContext(org.vcell.util.IssueContext) File(java.io.File)

Aggregations

Issue (org.vcell.util.Issue)88 ArrayList (java.util.ArrayList)18 IssueContext (org.vcell.util.IssueContext)14 Expression (cbit.vcell.parser.Expression)13 ExpressionException (cbit.vcell.parser.ExpressionException)13 PropertyVetoException (java.beans.PropertyVetoException)9 MolecularType (org.vcell.model.rbm.MolecularType)9 Structure (cbit.vcell.model.Structure)8 VCUnitDefinition (cbit.vcell.units.VCUnitDefinition)8 SymbolTableEntry (cbit.vcell.parser.SymbolTableEntry)7 MolecularComponentPattern (org.vcell.model.rbm.MolecularComponentPattern)7 Model (cbit.vcell.model.Model)6 ModelParameter (cbit.vcell.model.Model.ModelParameter)6 ExpressionBindingException (cbit.vcell.parser.ExpressionBindingException)6 BioModel (cbit.vcell.biomodel.BioModel)5 ModelException (cbit.vcell.model.ModelException)5 Product (cbit.vcell.model.Product)5 Reactant (cbit.vcell.model.Reactant)5 ReactionParticipant (cbit.vcell.model.ReactionParticipant)5 InteriorPoint (org.sbml.jsbml.ext.spatial.InteriorPoint)5