fimx 0.1.0__tar.gz

fimx-0.1.0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 Bernardo Porto Veronese
+
+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.

fimx-0.1.0/PKG-INFO

@@ -0,0 +1,203 @@
+Metadata-Version: 2.4
+Name: fimx
+Version: 0.1.0
+Summary: Dense Fisher matrix functions on labeled xarray arrays
+Keywords: fisher-matrix,forecasting,cosmology,xarray
+Author: binado
+Author-email: binado <bernardopveronese@gmail.com>
+License-Expression: MIT
+License-File: LICENSE
+Classifier: Development Status :: 3 - Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: Programming Language :: Python :: 3
+Classifier: Topic :: Scientific/Engineering
+Classifier: Typing :: Typed
+Requires-Dist: numpy
+Requires-Dist: xarray
+Requires-Dist: xarray-einstats>=0.10.0
+Requires-Dist: h5netcdf[h5py] ; extra == 'io'
+Requires-Dist: getdist ; extra == 'plotting'
+Requires-Dist: matplotlib ; extra == 'plotting'
+Requires-Python: >=3.12
+Project-URL: Homepage, https://github.com/binado/fimx
+Project-URL: Repository, https://github.com/binado/fimx
+Project-URL: Issues, https://github.com/binado/fimx/issues
+Project-URL: Changelog, https://github.com/binado/fimx/blob/main/CHANGELOG.md
+Provides-Extra: io
+Provides-Extra: plotting
+Description-Content-Type: text/markdown
+
+# fimx
+
+`fimx` provides functions for dense Fisher information matrices stored
+as ordinary `xarray.DataArray` objects. NumPy handles the computation; xarray
+handles labels, metadata containers, and persistence.
+
+```sh
+uv add fimx
+# Optional NetCDF backend:
+uv add 'fimx[io]'
+# Optional Gaussian plotting backend:
+uv add 'fimx[plotting]'
+```
+
+## Quickstart
+
+```python
+from fimx import combine, errors, fix, gaussian_prior, inv, marginalize, matrix
+
+F = matrix([[4, 1], [1, 2]], ["a", "b"])
+fixed = fix(F, "b")  # [[4.0]]: b is held fixed
+reduced = marginalize(F, "b")  # [[3.5]]: uncertainty in b is integrated out
+C = inv(F)  # covariance, with the same row/col labels
+sigma = errors(F)  # sqrt(diag(C)), dimension "parameter"
+
+prior = gaussian_prior({"a": 0.5})  # standard deviation 0.5 -> information 4
+posterior = combine(F, prior)
+```
+
+| Function | Behavior |
+| --- | --- |
+| `matrix(values, parameters)` | Construct and validate a canonical matrix. |
+| `expand(F, parameters)` | Embed in a larger or reordered parameter set, filling with zeros. |
+| `fix(F, parameters)` | Remove names through a principal submatrix. |
+| `marginalize(F, parameters)` | Remove names through the Schur complement. |
+| `inv(F)` | Return the inverse using Cholesky. |
+| `errors(F)` | Return marginalized standard deviations. |
+| `transform(F, jacobian)` | Change variables using `J.T @ F @ J`. |
+| `combine(*matrices)` | Sum independent information over the parameter union. |
+| `gaussian_prior(sigmas)` | Construct diagonal information `1 / sigma**2`. |
+
+## Array contract
+
+- Matrices have dimensions exactly `("row", "col")`, a nonempty square shape,
+  and explicit, identical, ordered coordinates of unique string names.
+  Batches, inferred dimensions, and coordinate repair are unsupported.
+- Values must be real, finite, and symmetric. They are converted to float64.
+- Positive semidefiniteness is not checked, so singular matrices can be built,
+  fixed, transformed, or combined. `inv` and `errors` require positive
+  definite input; `marginalize` only requires the removed block to be.
+- Functions return fresh DataArrays and never mutate their inputs. Attributes,
+  names, and auxiliary coordinates are not preserved.
+- Arrays built directly with xarray work if they meet the contract.
+
+Behavior worth knowing:
+
+- `fix` and `marginalize` accept one name or a sequence and keep the remaining
+  order. Unknown names raise `KeyError`; duplicates or removing every
+  parameter raise `ValueError`.
+- `expand` embeds a matrix in a larger or reordered set of names, filling new
+  entries with zeros. The target must contain every existing name. Matrices
+  expanded to the same names add with plain `+`, whereas `+` on mismatched
+  names silently keeps only the overlap.
+- `combine` keeps the first matrix's order and appends new parameters as they
+  appear. Missing entries contribute zero.
+- `gaussian_prior` takes **standard deviations**, not information values.
+- Failures are explicit: `ValueError` for malformed or non-finite data,
+  `TypeError` for non-DataArray inputs, and `numpy.linalg.LinAlgError` for
+  singular or indefinite matrices. There is no pseudoinverse or regularization.
+
+### Changes of variables
+
+The Jacobian has dimensions `("old", "new")` and orientation
+**`J[i, j] = d theta_i / d phi_j`**, where theta are the old parameters and phi
+the new ones. Every old parameter must appear exactly once (rows are reordered
+to match `F`), and the order of the new parameters sets the output order.
+Rectangular Jacobians are supported.
+
+```python
+import xarray as xr
+from fimx import transform
+
+# theta_a = 2 * phi_x, theta_b = phi_x
+# Rows deliberately appear in reverse order to F.
+J = xr.DataArray(
+    [[1.0], [2.0]],
+    dims=("old", "new"),
+    coords={"old": ["b", "a"], "new": ["x"]},
+)
+G = transform(F, J)  # [[22.0]]
+```
+
+## Plotting
+
+Plotting needs the `plotting` extra (GetDist and Matplotlib, loaded lazily).
+
+### Forecast Datasets
+
+Keep fiducials separate from the Fisher matrix and bundle them with
+`dataset()`. In the mapping of extra variables, 1D arrays use `row` and 2D
+arrays use `(row, col)`; plain arrays follow matrix order and indexed axes are
+aligned by label. Strings, units, and nonfinite metadata are allowed. The names
+`fisher`, `row`, and `col` are reserved.
+
+### Corner plots
+
+```python
+from fimx import dataset, plot
+
+fiducials = xr.DataArray([1.0, 2.0], dims="row", coords={"row": ["a", "b"]})
+survey_a = dataset(F, {"fiducials": fiducials, "units": ["km", "s"]})
+survey_b = dataset(2 * F, {"fiducials": fiducials})
+
+fig = plot({"Survey A": survey_a, "Survey B": survey_b})
+fig.savefig("constraints.pdf")
+
+# Explicit order and GetDist customization:
+fig = plot(
+    {"Survey A": survey_a, "Survey B": survey_b},
+    parameters=["b", "a"],
+    backend="getdist",
+    filled=False,
+    backend_kwargs={"contour_colors": ["C0", "C1"]},
+)
+```
+
+- `plot()` takes a nonempty mapping of labels to Datasets with `fisher` and
+  finite, real `fiducials`. Each forecast is centered on its own fiducials;
+  mapping order sets the overlay order and keys become legend labels.
+- By default, the plot uses the parameters shared by all forecasts, in the
+  first forecast's order. An explicit `parameters` list must be unique and
+  present in every forecast.
+- Each full Fisher matrix is inverted before selecting parameters, so omitted
+  parameters are marginalized and every complete matrix (nuisance blocks
+  included) must be positive definite.
+- The GetDist backend draws analytic Gaussians, with no sampling, and returns
+  a Matplotlib `Figure` without showing or saving it. Filled contours are the
+  default. `backend_kwargs` is forwarded to `triangle_plot`, except `roots`,
+  `params`, `legend_labels`, and `filled`, which `plot()` controls.
+- Additional backends implement `fimx.plotting.base.PlotBackend`.
+
+### Command line
+
+With both extras installed, `fimx-plot` draws a corner plot from NetCDF
+Datasets that contain `fisher` and `fiducials`:
+
+```sh
+uv add 'fimx[io,plotting]'
+fimx-plot --file survey-a.nc --file survey-b.nc \
+    --figure-file constraints.png --figure-dpi 200 \
+    --parameters a b --no-filled --backend getdist \
+    --backend-kwargs '{"contour_colors": ["C0", "C1"]}'
+```
+
+Each file stem becomes a legend label, so stems must be unique. The output
+defaults to `plot.png` at 150 dpi. `--parameters` selects and orders
+parameters, `--no-filled` draws line contours, `--backend` picks a backend, and
+`--backend-kwargs` takes a JSON object. `fimx.io.load_dataset(path)` loads the
+same files from Python.
+
+## Storage
+
+There is no `fimx` file format; use xarray. The `io` extra supplies `h5netcdf`.
+
+```python
+F.to_netcdf("fisher.nc", engine="h5netcdf")
+restored = xr.load_dataarray("fisher.nc", engine="h5netcdf")
+C_restored = inv(restored)
+```
+
+## More
+
+See [CHANGELOG.md](CHANGELOG.md) for release notes. `fimx` is released under
+the MIT license; see [LICENSE](LICENSE).

fimx-0.1.0/README.md

@@ -0,0 +1,174 @@
+# fimx
+
+`fimx` provides functions for dense Fisher information matrices stored
+as ordinary `xarray.DataArray` objects. NumPy handles the computation; xarray
+handles labels, metadata containers, and persistence.
+
+```sh
+uv add fimx
+# Optional NetCDF backend:
+uv add 'fimx[io]'
+# Optional Gaussian plotting backend:
+uv add 'fimx[plotting]'
+```
+
+## Quickstart
+
+```python
+from fimx import combine, errors, fix, gaussian_prior, inv, marginalize, matrix
+
+F = matrix([[4, 1], [1, 2]], ["a", "b"])
+fixed = fix(F, "b")  # [[4.0]]: b is held fixed
+reduced = marginalize(F, "b")  # [[3.5]]: uncertainty in b is integrated out
+C = inv(F)  # covariance, with the same row/col labels
+sigma = errors(F)  # sqrt(diag(C)), dimension "parameter"
+
+prior = gaussian_prior({"a": 0.5})  # standard deviation 0.5 -> information 4
+posterior = combine(F, prior)
+```
+
+| Function | Behavior |
+| --- | --- |
+| `matrix(values, parameters)` | Construct and validate a canonical matrix. |
+| `expand(F, parameters)` | Embed in a larger or reordered parameter set, filling with zeros. |
+| `fix(F, parameters)` | Remove names through a principal submatrix. |
+| `marginalize(F, parameters)` | Remove names through the Schur complement. |
+| `inv(F)` | Return the inverse using Cholesky. |
+| `errors(F)` | Return marginalized standard deviations. |
+| `transform(F, jacobian)` | Change variables using `J.T @ F @ J`. |
+| `combine(*matrices)` | Sum independent information over the parameter union. |
+| `gaussian_prior(sigmas)` | Construct diagonal information `1 / sigma**2`. |
+
+## Array contract
+
+- Matrices have dimensions exactly `("row", "col")`, a nonempty square shape,
+  and explicit, identical, ordered coordinates of unique string names.
+  Batches, inferred dimensions, and coordinate repair are unsupported.
+- Values must be real, finite, and symmetric. They are converted to float64.
+- Positive semidefiniteness is not checked, so singular matrices can be built,
+  fixed, transformed, or combined. `inv` and `errors` require positive
+  definite input; `marginalize` only requires the removed block to be.
+- Functions return fresh DataArrays and never mutate their inputs. Attributes,
+  names, and auxiliary coordinates are not preserved.
+- Arrays built directly with xarray work if they meet the contract.
+
+Behavior worth knowing:
+
+- `fix` and `marginalize` accept one name or a sequence and keep the remaining
+  order. Unknown names raise `KeyError`; duplicates or removing every
+  parameter raise `ValueError`.
+- `expand` embeds a matrix in a larger or reordered set of names, filling new
+  entries with zeros. The target must contain every existing name. Matrices
+  expanded to the same names add with plain `+`, whereas `+` on mismatched
+  names silently keeps only the overlap.
+- `combine` keeps the first matrix's order and appends new parameters as they
+  appear. Missing entries contribute zero.
+- `gaussian_prior` takes **standard deviations**, not information values.
+- Failures are explicit: `ValueError` for malformed or non-finite data,
+  `TypeError` for non-DataArray inputs, and `numpy.linalg.LinAlgError` for
+  singular or indefinite matrices. There is no pseudoinverse or regularization.
+
+### Changes of variables
+
+The Jacobian has dimensions `("old", "new")` and orientation
+**`J[i, j] = d theta_i / d phi_j`**, where theta are the old parameters and phi
+the new ones. Every old parameter must appear exactly once (rows are reordered
+to match `F`), and the order of the new parameters sets the output order.
+Rectangular Jacobians are supported.
+
+```python
+import xarray as xr
+from fimx import transform
+
+# theta_a = 2 * phi_x, theta_b = phi_x
+# Rows deliberately appear in reverse order to F.
+J = xr.DataArray(
+    [[1.0], [2.0]],
+    dims=("old", "new"),
+    coords={"old": ["b", "a"], "new": ["x"]},
+)
+G = transform(F, J)  # [[22.0]]
+```
+
+## Plotting
+
+Plotting needs the `plotting` extra (GetDist and Matplotlib, loaded lazily).
+
+### Forecast Datasets
+
+Keep fiducials separate from the Fisher matrix and bundle them with
+`dataset()`. In the mapping of extra variables, 1D arrays use `row` and 2D
+arrays use `(row, col)`; plain arrays follow matrix order and indexed axes are
+aligned by label. Strings, units, and nonfinite metadata are allowed. The names
+`fisher`, `row`, and `col` are reserved.
+
+### Corner plots
+
+```python
+from fimx import dataset, plot
+
+fiducials = xr.DataArray([1.0, 2.0], dims="row", coords={"row": ["a", "b"]})
+survey_a = dataset(F, {"fiducials": fiducials, "units": ["km", "s"]})
+survey_b = dataset(2 * F, {"fiducials": fiducials})
+
+fig = plot({"Survey A": survey_a, "Survey B": survey_b})
+fig.savefig("constraints.pdf")
+
+# Explicit order and GetDist customization:
+fig = plot(
+    {"Survey A": survey_a, "Survey B": survey_b},
+    parameters=["b", "a"],
+    backend="getdist",
+    filled=False,
+    backend_kwargs={"contour_colors": ["C0", "C1"]},
+)
+```
+
+- `plot()` takes a nonempty mapping of labels to Datasets with `fisher` and
+  finite, real `fiducials`. Each forecast is centered on its own fiducials;
+  mapping order sets the overlay order and keys become legend labels.
+- By default, the plot uses the parameters shared by all forecasts, in the
+  first forecast's order. An explicit `parameters` list must be unique and
+  present in every forecast.
+- Each full Fisher matrix is inverted before selecting parameters, so omitted
+  parameters are marginalized and every complete matrix (nuisance blocks
+  included) must be positive definite.
+- The GetDist backend draws analytic Gaussians, with no sampling, and returns
+  a Matplotlib `Figure` without showing or saving it. Filled contours are the
+  default. `backend_kwargs` is forwarded to `triangle_plot`, except `roots`,
+  `params`, `legend_labels`, and `filled`, which `plot()` controls.
+- Additional backends implement `fimx.plotting.base.PlotBackend`.
+
+### Command line
+
+With both extras installed, `fimx-plot` draws a corner plot from NetCDF
+Datasets that contain `fisher` and `fiducials`:
+
+```sh
+uv add 'fimx[io,plotting]'
+fimx-plot --file survey-a.nc --file survey-b.nc \
+    --figure-file constraints.png --figure-dpi 200 \
+    --parameters a b --no-filled --backend getdist \
+    --backend-kwargs '{"contour_colors": ["C0", "C1"]}'
+```
+
+Each file stem becomes a legend label, so stems must be unique. The output
+defaults to `plot.png` at 150 dpi. `--parameters` selects and orders
+parameters, `--no-filled` draws line contours, `--backend` picks a backend, and
+`--backend-kwargs` takes a JSON object. `fimx.io.load_dataset(path)` loads the
+same files from Python.
+
+## Storage
+
+There is no `fimx` file format; use xarray. The `io` extra supplies `h5netcdf`.
+
+```python
+F.to_netcdf("fisher.nc", engine="h5netcdf")
+restored = xr.load_dataarray("fisher.nc", engine="h5netcdf")
+C_restored = inv(restored)
+```
+
+## More
+
+See [CHANGELOG.md](CHANGELOG.md) for release notes. `fimx` is released under
+the MIT license; see [LICENSE](LICENSE).

fimx-0.1.0/pyproject.toml

@@ -0,0 +1,71 @@
+[project]
+name = "fimx"
+version = "0.1.0"
+description = "Dense Fisher matrix functions on labeled xarray arrays"
+readme = "README.md"
+license = "MIT"
+license-files = ["LICENSE*"]
+keywords = [
+    "fisher-matrix",
+    "forecasting",
+    "cosmology",
+    "xarray",
+]
+classifiers = [
+    "Development Status :: 3 - Alpha",
+    "Intended Audience :: Science/Research",
+    "Programming Language :: Python :: 3",
+    "Topic :: Scientific/Engineering",
+    "Typing :: Typed",
+]
+requires-python = ">=3.12"
+dependencies = [
+    "numpy",
+    "xarray",
+    "xarray-einstats>=0.10.0",
+]
+
+[[project.authors]]
+name = "binado"
+email = "bernardopveronese@gmail.com"
+
+[project.optional-dependencies]
+io = ["h5netcdf[h5py]"]
+plotting = [
+    "getdist",
+    "matplotlib",
+]
+
+[project.urls]
+Homepage = "https://github.com/binado/fimx"
+Repository = "https://github.com/binado/fimx"
+Issues = "https://github.com/binado/fimx/issues"
+Changelog = "https://github.com/binado/fimx/blob/main/CHANGELOG.md"
+
+[project.scripts]
+fimx-plot = "fimx.cli:main"
+
+[build-system]
+requires = ["uv_build>=0.12.7,<0.13.0"]
+build-backend = "uv_build"
+
+[dependency-groups]
+dev = [
+    "getdist>=1.7.7",
+    "h5netcdf[h5py]",
+    "ipykernel>=7.1.0",
+    "matplotlib>=3.10.8",
+    "prek>=0.5.3",
+    "pytest>=9.0.2",
+    "ruff>=0.16.0",
+    "ty>=0.0.83",
+]
+
+[tool.ruff]
+extend-exclude = ["examples"]
+
+[tool.ty.src]
+include = [
+    "src",
+    "tests",
+]

fimx-0.1.0/pyproject.toml.orig

@@ -0,0 +1,59 @@
+[project]
+name = "fimx"
+version = "0.1.0"
+description = "Dense Fisher matrix functions on labeled xarray arrays"
+readme = "README.md"
+license = "MIT"
+license-files = ["LICENSE*"]
+authors = [
+    { name = "binado", email = "bernardopveronese@gmail.com" }
+]
+keywords = ["fisher-matrix", "forecasting", "cosmology", "xarray"]
+classifiers = [
+    "Development Status :: 3 - Alpha",
+    "Intended Audience :: Science/Research",
+    "Programming Language :: Python :: 3",
+    "Topic :: Scientific/Engineering",
+    "Typing :: Typed",
+]
+requires-python = ">=3.12"
+dependencies = [
+    "numpy",
+    "xarray",
+    "xarray-einstats>=0.10.0",
+]
+
+[project.optional-dependencies]
+io = ["h5netcdf[h5py]"]
+plotting = ["getdist", "matplotlib"]
+
+[project.urls]
+Homepage = "https://github.com/binado/fimx"
+Repository = "https://github.com/binado/fimx"
+Issues = "https://github.com/binado/fimx/issues"
+Changelog = "https://github.com/binado/fimx/blob/main/CHANGELOG.md"
+
+[project.scripts]
+fimx-plot = "fimx.cli:main"
+
+[build-system]
+requires = ["uv_build>=0.12.7,<0.13.0"]
+build-backend = "uv_build"
+
+[dependency-groups]
+dev = [
+    "getdist>=1.7.7",
+    "h5netcdf[h5py]",
+    "ipykernel>=7.1.0",
+    "matplotlib>=3.10.8",
+    "prek>=0.5.3",
+    "pytest>=9.0.2",
+    "ruff>=0.16.0",
+    "ty>=0.0.83",
+]
+
+[tool.ruff]
+extend-exclude = ["examples"]
+
+[tool.ty.src]
+include = ["src", "tests"]

fimx-0.1.0/src/fimx/__init__.py

@@ -0,0 +1,22 @@
+"""Functions for Fisher information on labeled xarray matrices."""
+
+from .combination import combine, gaussian_prior
+from .construction import matrix
+from .datasets import dataset
+from .inversion import errors, inv
+from .parameters import expand, fix, marginalize, transform
+from .plotting import plot
+
+__all__ = [
+    "combine",
+    "dataset",
+    "errors",
+    "expand",
+    "fix",
+    "gaussian_prior",
+    "inv",
+    "marginalize",
+    "matrix",
+    "plot",
+    "transform",
+]

fimx-0.1.0/src/fimx/cli.py

@@ -0,0 +1,97 @@
+"""Command-line interface for plotting Fisher forecasts."""
+
+from argparse import ArgumentParser, BooleanOptionalAction
+from collections.abc import Sequence
+from json import loads
+from pathlib import Path
+
+from .io import load_dataset
+from .plotting import plot
+
+
+def _parser() -> ArgumentParser:
+    """Build the command-line argument parser."""
+    parser = ArgumentParser(
+        description="Overlay Fisher forecast Datasets in a corner plot."
+    )
+    parser.add_argument(
+        "--file",
+        dest="files",
+        action="append",
+        required=True,
+        type=Path,
+        metavar="PATH",
+        help="NetCDF Dataset to plot; repeat to overlay multiple forecasts.",
+    )
+    parser.add_argument(
+        "--figure-file",
+        type=Path,
+        default=Path("plot.png"),
+        metavar="PATH",
+        help="Output figure path (default: plot.png).",
+    )
+    parser.add_argument(
+        "--figure-dpi",
+        type=int,
+        default=150,
+        metavar="DPI",
+        help="Output resolution in dots per inch (default: 150).",
+    )
+    parser.add_argument(
+        "--parameters",
+        nargs="+",
+        metavar="NAME",
+        help="Parameters to include in the plot, in the requested order.",
+    )
+    parser.add_argument(
+        "--filled",
+        action=BooleanOptionalAction,
+        default=True,
+        help="Fill two-dimensional contours (default; use --no-filled for lines).",
+    )
+    parser.add_argument(
+        "--backend",
+        default="getdist",
+        metavar="NAME",
+        help="Plotting backend (default: getdist).",
+    )
+    parser.add_argument(
+        "--backend-kwargs",
+        type=loads,
+        metavar="JSON",
+        help="JSON object of backend-specific plotting options.",
+    )
+    return parser
+
+
+def main(argv: Sequence[str] | None = None) -> None:
+    """Load Fisher forecast files, plot them, and save the figure.
+
+    Parameters
+    ----------
+    argv : sequence of str, optional
+        Arguments to parse. Defaults to the process command line.
+    """
+    parser = _parser()
+    args = parser.parse_args(argv)
+    datasets = {}
+    for path in args.files:
+        label = path.stem
+        if label in datasets:
+            parser.error(
+                f"Input filenames must have unique stems; duplicate {label!r}."
+            )
+        datasets[label] = load_dataset(path)
+
+    figure = plot(
+        datasets,
+        parameters=args.parameters,
+        filled=args.filled,
+        backend=args.backend,
+        backend_kwargs=args.backend_kwargs,
+    )
+    figure.savefig(args.figure_file, dpi=args.figure_dpi)
+
+
+if __name__ == "__main__":
+    main()