refactor gmx_energy

utils/gmx-energy-plugin/ict.yml → utils/molecular-dynamics/gromacs/gmx-energy-tool/ict.yml

@@ -1,12 +1,12 @@
 specVersion: "0.1.0"
 name: gmx_energy
 version: 0.1.0
-container: gmx-energy-plugin
+container: gmx-energy-tool
 entrypoint:
 title: gmx_energy
 description: Wrapper of the GROMACS energy module for extracting energy components from a given GROMACS energy file.
-author: Data Scientist
-contact: [email protected]
+author: Brandon Walker, Nazanin Donyapour
+contact: [email protected], [email protected]
 repository:
 documentation:
 citation:

utils/gmx-energy-plugin/pyproject.toml → utils/molecular-dynamics/gromacs/gmx-energy-tool/pyproject.toml

@@ -8,8 +8,7 @@ readme = "README.md"
 [tool.poetry.dependencies]
 python = ">=3.9,<3.12"
 typer = "^0.7.0"
-cwl-utils = "0.33"
-cwltool = "3.1.20240404144621"
+sophios = "0.1.1"
 
 [tool.poetry.group.dev.dependencies]
 bump2version = "^1.0.1"

utils/molecular-dynamics/gromacs/gmx-energy-tool/tests/test_gmx_energy.py

@@ -0,0 +1,33 @@
+"""Tests for gmx_energy."""
+from pathlib import Path
+
+from sophios.api.pythonapi import Step
+from sophios.api.pythonapi import Workflow
+
+
+def test_gmx_energy() -> None:
+    """Test gmx_energy."""
+    # Define paths
+    cwl_file_str = "gmx_energy_0@[email protected]"
+    cwl_file = Path(__file__).resolve().parent.parent / Path(cwl_file_str)
+    input_energy_path = Path(__file__).resolve().parent / Path("energy.edr")
+
+    # Create the CWL step
+    gmx_energy_step = Step(clt_path=cwl_file)
+    gmx_energy_step.config = '{"terms": ["Potential"]}'
+    gmx_energy_step.input_energy_path = input_energy_path
+    gmx_energy_step.output_xvg_path = "system.xvg"
+
+    # Create the workflow and run it
+    steps = [gmx_energy_step]
+    filename = "gmx_energy_workflow"
+    workflow = Workflow(steps, filename)
+    workflow.run()
+
+    # Define output directory and check for expected output
+    outdir = Path("outdir")
+    files = list(outdir.rglob("system.xvg"))
+
+    assert (
+        files
+    ), f"The file 'system.xvg' does not exist in any subdirectory of '{outdir}'."