proface-abaqus-pre 0.1.0__tar.gz

proface_abaqus_pre-0.1.0/.github/workflows/python-publish.yml

@@ -0,0 +1,77 @@
+# SPDX-FileCopyrightText: Copyright (c) 2020 GitHub
+#
+# SPDX-License-Identifier: MIT
+
+# Original file downloaded from
+# https://raw.githubusercontent.com/actions/starter-workflows/00795b7feebceed43fe2510ccb976238441d9c07/ci/python-publish.yml
+
+# This workflow will upload a Python Package to PyPI when a release is created
+# For more information see: https://docs.github.com/en/actions/automating-builds-and-tests/building-and-testing-python#publishing-to-package-registries
+
+# This workflow uses actions that are not certified by GitHub.
+# They are provided by a third-party and are governed by
+# separate terms of service, privacy policy, and support
+# documentation.
+
+name: Upload Python Package
+
+on:
+  release:
+    types: [published]
+
+permissions:
+  contents: read
+
+jobs:
+  release-build:
+    runs-on: ubuntu-latest
+
+    steps:
+      - uses: actions/checkout@v4
+
+      - uses: actions/setup-python@v5
+        with:
+          python-version: "3.x"
+
+      - name: Build release distributions
+        run: |
+          # NOTE: put your own distribution build steps here.
+          python -m pip install build
+          python -m build
+
+      - name: Upload distributions
+        uses: actions/upload-artifact@v4
+        with:
+          name: release-dists
+          path: dist/
+
+  pypi-publish:
+    runs-on: ubuntu-latest
+    needs:
+      - release-build
+    permissions:
+      # IMPORTANT: this permission is mandatory for trusted publishing
+      id-token: write
+
+    # Dedicated environments with protections for publishing are strongly recommended.
+    # For more information, see: https://docs.github.com/en/actions/deployment/targeting-different-environments/using-environments-for-deployment#deployment-protection-rules
+    environment:
+      name: pypi
+      # OPTIONAL: uncomment and update to include your PyPI project URL in the deployment status:
+      # url: https://pypi.org/p/YOURPROJECT
+      #
+      # ALTERNATIVE: if your GitHub Release name is the PyPI project version string
+      # ALTERNATIVE: exactly, uncomment the following line instead:
+      url: https://pypi.org/project/proface-abaqus-pre/${{ github.event.release.name }}
+
+    steps:
+      - name: Retrieve release distributions
+        uses: actions/download-artifact@v4
+        with:
+          name: release-dists
+          path: dist/
+
+      - name: Publish release distributions to PyPI
+        uses: pypa/gh-action-pypi-publish@release/v1
+        with:
+          packages-dir: dist/

proface_abaqus_pre-0.1.0/.gitignore

@@ -0,0 +1,17 @@
+# SPDX-FileCopyrightText: 2025 ProFACE developers
+#
+# SPDX-License-Identifier: CC0-1.0
+
+# Python-generated files
+__pycache__/
+*.py[oc]
+build/
+dist/
+wheels/
+*.egg-info
+
+# Virtual environments
+.venv
+
+# hatch-vcs build hook
+src/proface/abaqus/preprocessor/_version.py

proface_abaqus_pre-0.1.0/.pre-commit-config.yaml

@@ -0,0 +1,28 @@
+# SPDX-FileCopyrightText: 2025 ProFACE developers
+#
+# SPDX-License-Identifier: CC0-1.0
+
+# See https://pre-commit.com for more information
+# See https://pre-commit.com/hooks.html for more hooks
+repos:
+  - repo: https://github.com/pre-commit/pre-commit-hooks
+    rev: v5.0.0
+    hooks:
+    # general checks
+    -   id: trailing-whitespace
+    -   id: end-of-file-fixer
+    -   id: mixed-line-ending
+    -   id: check-added-large-files
+    -   id: check-case-conflict
+    # check file formats
+    -   id: check-yaml
+    -   id: check-toml
+  - repo: https://github.com/astral-sh/ruff-pre-commit
+    rev: v0.11.3
+    hooks:
+      - id: ruff
+      - id: ruff-format
+  - repo: https://github.com/fsfe/reuse-tool
+    rev: v5.0.2
+    hooks:
+      - id: reuse

proface_abaqus_pre-0.1.0/LICENSES/CC0-1.0.txt

@@ -0,0 +1,121 @@
+Creative Commons Legal Code
+
+CC0 1.0 Universal
+
+    CREATIVE COMMONS CORPORATION IS NOT A LAW FIRM AND DOES NOT PROVIDE
+    LEGAL SERVICES. DISTRIBUTION OF THIS DOCUMENT DOES NOT CREATE AN
+    ATTORNEY-CLIENT RELATIONSHIP. CREATIVE COMMONS PROVIDES THIS
+    INFORMATION ON AN "AS-IS" BASIS. CREATIVE COMMONS MAKES NO WARRANTIES
+    REGARDING THE USE OF THIS DOCUMENT OR THE INFORMATION OR WORKS
+    PROVIDED HEREUNDER, AND DISCLAIMS LIABILITY FOR DAMAGES RESULTING FROM
+    THE USE OF THIS DOCUMENT OR THE INFORMATION OR WORKS PROVIDED
+    HEREUNDER.
+
+Statement of Purpose
+
+The laws of most jurisdictions throughout the world automatically confer
+exclusive Copyright and Related Rights (defined below) upon the creator
+and subsequent owner(s) (each and all, an "owner") of an original work of
+authorship and/or a database (each, a "Work").
+
+Certain owners wish to permanently relinquish those rights to a Work for
+the purpose of contributing to a commons of creative, cultural and
+scientific works ("Commons") that the public can reliably and without fear
+of later claims of infringement build upon, modify, incorporate in other
+works, reuse and redistribute as freely as possible in any form whatsoever
+and for any purposes, including without limitation commercial purposes.
+These owners may contribute to the Commons to promote the ideal of a free
+culture and the further production of creative, cultural and scientific
+works, or to gain reputation or greater distribution for their Work in
+part through the use and efforts of others.
+
+For these and/or other purposes and motivations, and without any
+expectation of additional consideration or compensation, the person
+associating CC0 with a Work (the "Affirmer"), to the extent that he or she
+is an owner of Copyright and Related Rights in the Work, voluntarily
+elects to apply CC0 to the Work and publicly distribute the Work under its
+terms, with knowledge of his or her Copyright and Related Rights in the
+Work and the meaning and intended legal effect of CC0 on those rights.
+
+1. Copyright and Related Rights. A Work made available under CC0 may be
+protected by copyright and related or neighboring rights ("Copyright and
+Related Rights"). Copyright and Related Rights include, but are not
+limited to, the following:
+
+  i. the right to reproduce, adapt, distribute, perform, display,
+     communicate, and translate a Work;
+ ii. moral rights retained by the original author(s) and/or performer(s);
+iii. publicity and privacy rights pertaining to a person's image or
+     likeness depicted in a Work;
+ iv. rights protecting against unfair competition in regards to a Work,
+     subject to the limitations in paragraph 4(a), below;
+  v. rights protecting the extraction, dissemination, use and reuse of data
+     in a Work;
+ vi. database rights (such as those arising under Directive 96/9/EC of the
+     European Parliament and of the Council of 11 March 1996 on the legal
+     protection of databases, and under any national implementation
+     thereof, including any amended or successor version of such
+     directive); and
+vii. other similar, equivalent or corresponding rights throughout the
+     world based on applicable law or treaty, and any national
+     implementations thereof.
+
+2. Waiver. To the greatest extent permitted by, but not in contravention
+of, applicable law, Affirmer hereby overtly, fully, permanently,
+irrevocably and unconditionally waives, abandons, and surrenders all of
+Affirmer's Copyright and Related Rights and associated claims and causes
+of action, whether now known or unknown (including existing as well as
+future claims and causes of action), in the Work (i) in all territories
+worldwide, (ii) for the maximum duration provided by applicable law or
+treaty (including future time extensions), (iii) in any current or future
+medium and for any number of copies, and (iv) for any purpose whatsoever,
+including without limitation commercial, advertising or promotional
+purposes (the "Waiver"). Affirmer makes the Waiver for the benefit of each
+member of the public at large and to the detriment of Affirmer's heirs and
+successors, fully intending that such Waiver shall not be subject to
+revocation, rescission, cancellation, termination, or any other legal or
+equitable action to disrupt the quiet enjoyment of the Work by the public
+as contemplated by Affirmer's express Statement of Purpose.
+
+3. Public License Fallback. Should any part of the Waiver for any reason
+be judged legally invalid or ineffective under applicable law, then the
+Waiver shall be preserved to the maximum extent permitted taking into
+account Affirmer's express Statement of Purpose. In addition, to the
+extent the Waiver is so judged Affirmer hereby grants to each affected
+person a royalty-free, non transferable, non sublicensable, non exclusive,
+irrevocable and unconditional license to exercise Affirmer's Copyright and
+Related Rights in the Work (i) in all territories worldwide, (ii) for the
+maximum duration provided by applicable law or treaty (including future
+time extensions), (iii) in any current or future medium and for any number
+of copies, and (iv) for any purpose whatsoever, including without
+limitation commercial, advertising or promotional purposes (the
+"License"). The License shall be deemed effective as of the date CC0 was
+applied by Affirmer to the Work. Should any part of the License for any
+reason be judged legally invalid or ineffective under applicable law, such
+partial invalidity or ineffectiveness shall not invalidate the remainder
+of the License, and in such case Affirmer hereby affirms that he or she
+will not (i) exercise any of his or her remaining Copyright and Related
+Rights in the Work or (ii) assert any associated claims and causes of
+action with respect to the Work, in either case contrary to Affirmer's
+express Statement of Purpose.
+
+4. Limitations and Disclaimers.
+
+ a. No trademark or patent rights held by Affirmer are waived, abandoned,
+    surrendered, licensed or otherwise affected by this document.
+ b. Affirmer offers the Work as-is and makes no representations or
+    warranties of any kind concerning the Work, express, implied,
+    statutory or otherwise, including without limitation warranties of
+    title, merchantability, fitness for a particular purpose, non
+    infringement, or the absence of latent or other defects, accuracy, or
+    the present or absence of errors, whether or not discoverable, all to
+    the greatest extent permissible under applicable law.
+ c. Affirmer disclaims responsibility for clearing rights of other persons
+    that may apply to the Work or any use thereof, including without
+    limitation any person's Copyright and Related Rights in the Work.
+    Further, Affirmer disclaims responsibility for obtaining any necessary
+    consents, permissions or other rights required for any use of the
+    Work.
+ d. Affirmer understands and acknowledges that Creative Commons is not a
+    party to this document and has no duty or obligation with respect to
+    this CC0 or use of the Work.

proface_abaqus_pre-0.1.0/LICENSES/MIT.txt

@@ -0,0 +1,18 @@
+MIT License
+
+Copyright (c) <year> <copyright holders>
+
+Permission is hereby granted, free of charge, to any person obtaining a copy of this software and
+associated documentation files (the "Software"), to deal in the Software without restriction, including
+without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the
+following conditions:
+
+The above copyright notice and this permission notice shall be included in all copies or substantial
+portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT
+LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO
+EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER
+IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE
+USE OR OTHER DEALINGS IN THE SOFTWARE.

proface_abaqus_pre-0.1.0/PKG-INFO

@@ -0,0 +1,78 @@
+Metadata-Version: 2.4
+Name: proface-abaqus-pre
+Version: 0.1.0
+Summary: Abaqus to ProFACE preprocessor
+Project-URL: documentation, https://github.com/ProFACE-dev/proface-abaqus-pre#proface-abaqus-pre
+Project-URL: issues, https://github.com/ProFACE-dev/proface-abaqus-pre/issues
+Project-URL: source, https://github.com/ProFACE-dev/proface-abaqus-pre
+Author-email: Stefano Miccoli <stefano.miccoli@polimi.it>
+License-Expression: MIT
+Classifier: Development Status :: 3 - Alpha
+Classifier: Environment :: Console
+Classifier: Programming Language :: Python
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.13
+Classifier: Programming Language :: Python :: Implementation :: CPython
+Requires-Python: >=3.12
+Requires-Dist: h5py>=3.13.0
+Requires-Dist: numpy>=2.0
+Requires-Dist: proface-pre>=0.2.0
+Requires-Dist: suanpan-abaqus~=0.1.0
+Description-Content-Type: text/markdown
+
+<!--
+SPDX-FileCopyrightText: 2025 ProFACE developers
+
+SPDX-License-Identifier: MIT
+-->
+
+# ProFACE-Abaqus-Pre
+
+ProFACE-Abaqus-Pre is a python package that provides the `proface.preprocessor.abaqus` plugin for use with the [`proface-pre`](https://github.com/ProFACE-dev/proface-pre) command-line interface (CLI).
+
+This plugin enables the conversion of Abaqus binary `.fil` results files to ProFACE `.h5` FEA input files.
+
+## Installation
+
+Install from <https://pypi.org> with
+
+```
+pip install proface-abaqus-pre
+```
+
+## Usage
+
+Running
+
+```
+proface-pre example.toml
+```
+
+will produce `example.h5` from `example.fil`.
+
+`example.toml` has the following format:
+
+```toml
+fea_software = "Abaqus"
+
+[Abaqus.input]
+fil = "example.fil"  # optional if .toml and .fil have the same name
+
+[Abaqus.results.ref_load]
+step = 1  # 0 for last step
+increment = 1  # 0 for last increment
+```
+
+## License
+
+ProFACE-Abaqus-Pre is licensed under the MIT license.
+ProFACE-Pre makes no claims about [ProFACE](https://proface.polimi.it) which is a distinct program with different licensing requirements.
+
+### Disclaimer
+
+This project is not affiliated with, endorsed by, or sponsored by Dassault Systèmes, The HDF Group, or any of their respective products.
+
+- Abaqus is a registered trademark of Dassault Systèmes or its subsidiaries in the United States and/or other countries.
+- HDF5 is a trademark of The HDF Group.
+
+All trademarks and registered trademarks are the property of their respective owners.

proface_abaqus_pre-0.1.0/README.md

@@ -0,0 +1,56 @@
+<!--
+SPDX-FileCopyrightText: 2025 ProFACE developers
+
+SPDX-License-Identifier: MIT
+-->
+
+# ProFACE-Abaqus-Pre
+
+ProFACE-Abaqus-Pre is a python package that provides the `proface.preprocessor.abaqus` plugin for use with the [`proface-pre`](https://github.com/ProFACE-dev/proface-pre) command-line interface (CLI).
+
+This plugin enables the conversion of Abaqus binary `.fil` results files to ProFACE `.h5` FEA input files.
+
+## Installation
+
+Install from <https://pypi.org> with
+
+```
+pip install proface-abaqus-pre
+```
+
+## Usage
+
+Running
+
+```
+proface-pre example.toml
+```
+
+will produce `example.h5` from `example.fil`.
+
+`example.toml` has the following format:
+
+```toml
+fea_software = "Abaqus"
+
+[Abaqus.input]
+fil = "example.fil"  # optional if .toml and .fil have the same name
+
+[Abaqus.results.ref_load]
+step = 1  # 0 for last step
+increment = 1  # 0 for last increment
+```
+
+## License
+
+ProFACE-Abaqus-Pre is licensed under the MIT license.
+ProFACE-Pre makes no claims about [ProFACE](https://proface.polimi.it) which is a distinct program with different licensing requirements.
+
+### Disclaimer
+
+This project is not affiliated with, endorsed by, or sponsored by Dassault Systèmes, The HDF Group, or any of their respective products.
+
+- Abaqus is a registered trademark of Dassault Systèmes or its subsidiaries in the United States and/or other countries.
+- HDF5 is a trademark of The HDF Group.
+
+All trademarks and registered trademarks are the property of their respective owners.

proface_abaqus_pre-0.1.0/pyproject.toml

@@ -0,0 +1,77 @@
+# SPDX-FileCopyrightText: 2025 ProFACE developers
+#
+# SPDX-License-Identifier: MIT
+
+[project]
+dynamic = ["version"]
+name = "proface-abaqus-pre"
+description = "Abaqus to ProFACE preprocessor"
+readme = "README.md"
+authors = [
+    { name = "Stefano Miccoli", email = "stefano.miccoli@polimi.it" }
+]
+mantainers = [
+    { name = "ProFACE developers", email = "proface@polimi.it" }
+]
+classifiers = [
+  "Development Status :: 3 - Alpha",
+  "Environment :: Console",
+  "Programming Language :: Python",
+  "Programming Language :: Python :: 3.12",
+  "Programming Language :: Python :: 3.13",
+  "Programming Language :: Python :: Implementation :: CPython",
+]
+requires-python = ">=3.12"
+dependencies = [
+    "proface-pre>=0.2.0",
+    "h5py>=3.13.0",
+    "numpy>=2.0",
+    "suanpan-abaqus~=0.1.0",
+]
+license = "MIT"
+
+[project.entry-points."proface.preprocessor"]
+abaqus = "proface.abaqus.preprocessor:main"
+
+[project.urls]
+documentation = "https://github.com/ProFACE-dev/proface-abaqus-pre#proface-abaqus-pre"
+issues = "https://github.com/ProFACE-dev/proface-abaqus-pre/issues"
+source = "https://github.com/ProFACE-dev/proface-abaqus-pre"
+
+[build-system]
+requires = ["hatchling", "hatch-vcs"]
+build-backend = "hatchling.build"
+
+#
+# hatch
+#
+[tool.hatch.version]
+source = "vcs"
+[tool.hatch.version.raw-options]
+version_scheme = "release-branch-semver"
+local_scheme = "node-and-date"
+[tool.hatch.build.hooks.vcs]
+version-file = "src/proface/abaqus/preprocessor/_version.py"
+[tool.hatch.build.targets.wheel]
+packages = ["src/proface"]
+
+#
+# ruff
+#
+[tool.ruff]
+line-length = 80
+target-version = "py312"
+
+[tool.ruff.lint]
+select = ["A", "ARG", "B", "BLE", "C4", "C90", "COM818",
+  "E", "W", "EM", "F", "FA", "FBT", "FLY", "G", "I", "ICN",
+  "INP", "INT", "ISC", "N", "NPY", "PERF", "PGH", "PIE",
+  "PLE", "PLR", "PLW", "PT", "PYI", "RET", "RUF", "S",
+  "SIM", "SLF", "SLOT", "T10", "T20", "TCH", "TRY", "TID",
+  "UP", "YTT",
+]
+ignore = ["S101", "PT015", "B011", "TRY400"]
+
+
+[tool.ruff.lint.isort]
+known-first-party = ["proface"]

proface_abaqus_pre-0.1.0/src/proface/abaqus/preprocessor/__init__.py

@@ -0,0 +1,10 @@
+# SPDX-FileCopyrightText: 2025 ProFACE developers
+#
+# SPDX-License-Identifier: MIT
+
+"""ProFACE Preprocessor package for Abaqus FEA"""
+
+from ._version import __version__
+from .translator import AbaqusTranslatorError, main
+
+__all__ = ["AbaqusTranslatorError", "__version__", "main"]

proface_abaqus_pre-0.1.0/src/proface/abaqus/preprocessor/_version.py

@@ -0,0 +1,21 @@
+# file generated by setuptools-scm
+# don't change, don't track in version control
+
+__all__ = ["__version__", "__version_tuple__", "version", "version_tuple"]
+
+TYPE_CHECKING = False
+if TYPE_CHECKING:
+    from typing import Tuple
+    from typing import Union
+
+    VERSION_TUPLE = Tuple[Union[int, str], ...]
+else:
+    VERSION_TUPLE = object
+
+version: str
+__version__: str
+__version_tuple__: VERSION_TUPLE
+version_tuple: VERSION_TUPLE
+
+__version__ = version = '0.1.0'
+__version_tuple__ = version_tuple = (0, 1, 0)

proface_abaqus_pre-0.1.0/src/proface/abaqus/preprocessor/translator.py

@@ -0,0 +1,337 @@
+# SPDX-FileCopyrightText: 2025 ProFACE developers
+#
+# SPDX-License-Identifier: MIT
+
+"""Translator implementation"""
+
+import collections.abc
+import logging
+from pathlib import Path
+
+import h5py
+import numpy as np
+from suanpan.abqfil import AbqFil
+
+from proface.preprocessor import PreprocessorError
+
+from . import __version__
+
+logger = logging.getLogger(__name__)
+
+
+# translation from Abaqus result codes to ProFACE ids
+ABQ_VAR = {"R11": "S", "R401": "SP", "R76": "IVOL", "R8": "COORD"}
+# translation from Abaqus location codes to ProFACE ids
+ABQ_LOC = {0: "integration_point", 4: "nodal_averaged"}
+
+# function for computing results path in h5
+_h5_path = "{var:s}/{loc:s}/{eltype:s}".format
+
+# save results in single precision
+RES_DTYPE = np.dtype("float32")
+
+
+class AbaqusTranslatorError(PreprocessorError):
+    pass
+
+
+def main(
+    *, job: collections.abc.Mapping, job_path: Path, h5: h5py.Group
+) -> None:
+    """main entrypoint for abaqus preprocessor"""
+
+    # runtime type checking
+    if not isinstance(job, collections.abc.Mapping):
+        msg = "'job' must be a mapping"
+        raise TypeError(msg)
+    if not isinstance(job_path, Path):
+        msg = "'job_path' must be a pathlib.Path"
+        raise TypeError(msg)
+    if not isinstance(h5, h5py.Group):
+        msg = "'h5' must be a h5py.Group or h5py.File"
+        raise TypeError(msg)
+
+    logger.info(
+        "\U0001f680 START Abaqus to ProFACE translator, ver. %s",
+        __version__,
+    )  # 🚀
+
+    # compute .fil path
+    if "fil" in job.get("input", {}):
+        filpth = job_path.parent / job["input"]["fil"]
+    else:
+        filpth = job_path.with_suffix(".fil")
+
+    # run translator
+    try:
+        _pre(filpth=filpth, h5=h5, results=job.get("results", {}))
+    except OSError:
+        # caller should treat OSError
+        raise
+    except Exception:
+        msg = "Internal Error"
+        logger.exception(msg)
+        raise AbaqusTranslatorError(msg) from None
+    logger.info("\U0001f3c1 END Abaqus to ProFACE translator")  # 🏁
+
+
+def _pre(filpth, h5: h5py.File, results):
+    """Abaqus translator"""
+
+    logger.info("reading %s", filpth.resolve().as_uri())
+    logger.info("writing %s", Path(h5.filename).resolve().as_uri())
+
+    fil = AbqFil(filpth)
+
+    _write_meta(fil=fil, h5=h5)
+    _write_nodal(fil=fil, h5=h5)
+    _write_element(fil=fil, h5=h5)
+    _write_sets(fil=fil, h5=h5)
+
+    if not results:
+        logger.warning("no results request")
+        return
+    h5_res = h5.create_group("results")
+    for k, v in results.items():
+        if "step" not in v or "increment" not in v:
+            logger.error(
+                "Results request '%s': "
+                "both 'step' and 'increment' must be specified.",
+                k,
+            )
+            continue
+        _write_results(
+            fil=fil,
+            h5=h5,
+            h5_res=h5_res,
+            name=k,
+            step=v["step"],
+            inc=v["increment"],
+        )
+
+
+def _write_meta(fil, h5):
+    #
+    # metadata
+    #
+    logger.info("Abaqus ver. %s", _label(fil.info["ver"]))
+    logger.info(
+        "Analysis run on %s",
+        _label(fil.info["date"]) + " " + _label(fil.info["time"]),
+    )
+    if fil.heading.strip():
+        logger.info("Heading '%s'", _label(fil.heading))
+    logger.info("Number of elements: %9d", fil.info["nelm"])
+    logger.info("Number of nodes   : %9d", fil.info["nnod"])
+
+    h5.attrs["program"] = "Abaqus"
+    h5.attrs["version"] = _label(fil.info["ver"])
+    h5.attrs["run_datetime"] = (
+        _label(fil.info["date"]) + " " + _label(fil.info["time"])
+    )
+    h5.attrs["title"] = _label(fil.heading)
+
+
+def _write_nodal(fil, h5):
+    #
+    # nodal data
+    #
+    h5_nodes = h5.create_group("nodes")
+    h5_nodes.attrs["number"] = fil.info["nnod"]
+    h5_nodes.create_dataset("coordinates", data=fil.coord["coord"])
+    h5_nodes.create_dataset("numbers", data=fil.coord["nnum"])
+
+
+def _write_element(fil, h5):
+    #
+    # element data
+    #
+    h5_elements = h5.create_group("elements")
+    h5_elements.attrs["number"] = fil.info["nelm"]
+    h5_elements.attrs["size"] = fil.info["elsiz"]
+    for elbloc in fil.elm:
+        assert (elbloc["eltyp"] == elbloc["eltyp"][0]).all()
+        eltype = _safe_label(elbloc["eltyp"][0])
+        if eltype in h5_elements:
+            label_m = _mangle_element_label(eltype, h5_elements.keys())
+            logger.warning(
+                "Mapping %s \N{RIGHTWARDS ARROW FROM BAR} %s",
+                eltype,
+                label_m,
+            )
+            eltype = label_m
+        h5_elgroup = h5_elements.create_group(eltype)
+        h5_elgroup.create_dataset("incidences", data=elbloc["ninc"])
+        h5_elgroup.create_dataset("numbers", data=elbloc["elnum"])
+
+        h5_elgroup.create_dataset(
+            "nodes", data=np.unique(h5_elgroup["incidences"])
+        )
+
+
+def _write_sets(fil, h5):
+    #
+    # sets
+    #
+    h5_sets = h5.create_group("sets")
+    h5_sets_node = h5_sets.create_group("node")
+    for k, nset in fil.nset.items():
+        k = _safe_label(fil.label.get(k, k))  # noqa: PLW2901
+        h5_sets_node.create_dataset(k, data=nset)
+    h5_sets_element = h5_sets.create_group("element")
+    for k, elset in fil.elset.items():
+        k = _safe_label(fil.label.get(k, k))  # noqa: PLW2901
+        h5_sets_element.create_dataset(k, data=elset)
+
+
+def _write_results(fil, h5, h5_res, name: str, step: int, inc: int) -> None:  # noqa: PLR0913
+    #
+    # results
+    #
+    try:
+        i = _find_step_inc(fil.step, step, inc)
+    except ValueError as exc:
+        logger.error("Results '%s': %s", name, exc)
+        ## fixme: raise error ?
+        return
+
+    logger.info(
+        "Results '%s': step %d, increment %d, time %#.3g, '%s'",
+        name,
+        *fil.step[i][["step", "incr", "ttime"]],
+        _label(fil.step[i]["subheading"]),
+    )
+    h5_k = h5_res.create_group(name)
+    h5_k.attrs["step"] = fil.step[i]["step"]
+    h5_k.attrs["increment"] = fil.step[i]["incr"]
+    h5_k.attrs["time"] = fil.step[i]["ttime"]
+
+    _write_step_output_blocks(fil, h5, h5_k, i)
+
+
+def _write_step_output_blocks(fil, h5, h5_k, i):  # noqa: C901
+    for flag, elset, eltype, data in fil.get_step(i):
+        # check block flag: 0 is element output
+        if flag != 0:
+            logger.warning("Skipping non element output: flag = %d", flag)
+            continue
+
+        # check block elset
+        if _label(elset) != "":
+            logger.error("Results file with element sets not supported")
+            continue
+
+        # check block element type
+        eltype = _safe_label(eltype)  # noqa: PLW2901
+        if not eltype.startswith("C3D"):
+            logger.warning("Results for element %s ignored", eltype)
+            continue
+
+        # check block location
+        loc = data["loc"][0]
+        assert (data["loc"] == loc).all()
+        match ABQ_LOC.get(loc):
+            case "integration_point":
+                # reshape data to index as [el_num, ip_num]
+                nr_ip = _guess_nr_ip(data)
+                data = data.reshape(-1, nr_ip)  # noqa: PLW2901
+                # data["num"] is elnum across columns
+                if not np.all(
+                    data["num"]
+                    == np.expand_dims(h5["elements"][eltype]["numbers"], -1)
+                ):
+                    msg = f"Inconsistent records for {eltype}: element numbers"
+                    raise ValueError(msg)
+                # data["ipnum"] is 1..nr_ip across rows
+                assert np.all(data["ipnum"] == 1 + np.arange(nr_ip))
+            case "nodal_averaged":
+                if not np.all(data["num"] == h5["elements"][eltype]["nodes"]):
+                    msg = f"Inconsistent records for {eltype}: node numbers"
+                assert (data["ipnum"] == 0).all()
+            case None:
+                logger.warning("Unknown location code %d", loc)
+                continue
+
+        # save block data
+        for name in data.dtype.names:
+            if not name.startswith("R"):
+                continue
+            if name not in ABQ_VAR:
+                logger.warning("Unexpected results code %s", name)
+                continue
+            dset = h5_k.create_dataset(
+                _h5_path(var=ABQ_VAR[name], loc=ABQ_LOC[loc], eltype=eltype),
+                data=data[name].astype(RES_DTYPE),
+            )
+            logger.debug("Wrote %s: %s", dset.name, dset.shape)
+
+
+def _label(lab: bytes) -> str:
+    return lab.decode("ASCII").strip()
+
+
+def _safe_label(lab: bytes) -> str:
+    slabel = _label(lab)
+    if slabel == ".":
+        return "._"
+    return slabel.replace("/", "|")
+
+
+def _mangle_element_label(label: str, groups) -> str:
+    assert label in groups
+    i = 1
+    for s in groups:
+        if s.startswith(f"{label:s}@"):
+            i += 1
+    mangled = f"{label:2}@{i:d}"
+    assert mangled not in groups
+    return mangled
+
+
+def _find_step_inc(stepdata, step, inc):
+    """search for requested step/increment"""
+
+    if len(stepdata) == 0:
+        msg = "no stepdata in file"
+        raise ValueError(msg)
+    logger.debug("Requested: Step %d, increment %d", step, inc)
+    if step > 0 and step not in stepdata["step"]:
+        # explicit request of inesistent step
+        msg = f"step '{step}' not found"
+        raise ValueError(msg)
+    if step == 0:
+        # last step requested
+        step = stepdata["step"][-1]
+
+    c_step = stepdata[stepdata["step"] == step]
+    if inc > 0 and inc not in c_step["incr"]:
+        # explicit request for inesistent increment
+        msg = f"increment '{inc}' not found in step '{step}'"
+        raise ValueError(msg)
+    if inc == 0:
+        # last increment requested
+        inc = c_step["incr"][-1]
+
+    assert stepdata.ndim == 1
+    (i,) = np.nonzero((stepdata["step"] == step) & (stepdata["incr"] == inc))
+    assert np.shape(i) == (1,), f"Multiple step blocks found: {i}"
+    i = i.item()
+
+    logger.debug("Found: Step %d, increment %d at position %s", step, inc, i)
+    assert stepdata[i][["step", "incr"]].item() == (step, inc)
+    return i
+
+
+def _guess_nr_ip(data):
+    """heuristics to find the number of integration points
+    for a loc 0 (integration points) record"""
+
+    # assumption is that number of ip's is smallish
+    num = data["num"]
+    for i, v in enumerate(num):
+        if v != num[0]:
+            return i
+
+    # edge case: record contains a single element
+    assert np.all(num == num[0])
+    return len(data)