[](https://github.com/hippylib/hippylib/actions/workflows/ci.yaml)
[](https://hippylib.readthedocs.io/en/latest/)
[](http://joss.theoj.org/papers/053e0d08a5e9755e7b78898cff6f6208) [](https://doi.org/10.5281/zenodo.596931)
Inverse Problem PYthon library
```
__ ______ _______ _______ __ __ __ __ __
/ | / |/ \ / \ / \ / |/ |/ |/ |
$$ |____ $$$$$$/ $$$$$$$ |$$$$$$$ |$$ \ /$$/ $$ |$$/ $$ |____
$$ \ $$ | $$ |__$$ |$$ |__$$ | $$ \/$$/ $$ |/ |$$ \
$$$$$$$ | $$ | $$ $$/ $$ $$/ $$ $$/ $$ |$$ |$$$$$$$ |
$$ | $$ | $$ | $$$$$$$/ $$$$$$$/ $$$$/ $$ |$$ |$$ | $$ |
$$ | $$ | _$$ |_ $$ | $$ | $$ | $$ |$$ |$$ |__$$ |
$$ | $$ |/ $$ |$$ | $$ | $$ | $$ |$$ |$$ $$/
$$/ $$/ $$$$$$/ $$/ $$/ $$/ $$/ $$/ $$$$$$$/
```
https://hippylib.github.io
`hIPPYlib` implements state-of-the-art scalable algorithms for
deterministic and Bayesian inverse problems governed by partial differential equations (PDEs).
It builds on [FEniCS](http://fenicsproject.org/)
(a parallel finite element element library) for the discretization of the PDE
and on [PETSc](http://www.mcs.anl.gov/petsc/) for scalable and efficient linear
algebra operations and solvers.
For building instructions, see the file `INSTALL.md`. Copyright information
and licensing restrictions can be found in the file `COPYRIGHT`.
The best starting point for new users interested in `hIPPYlib`'s
features are the interactive tutorials in the `tutorial folder`.
Conceptually, `hIPPYlib` can be viewed as a toolbox that provides the
building blocks for experimenting new ideas and developing scalable
algorithms for PDE-constrained deterministic and Bayesian inverse problems.
In `hIPPYlib` the user can express the forward PDE and the likelihood in
weak form using the friendly, compact, near-mathematical notation of
`FEniCS`, which will then automatically generate efficient code for the
discretization. Linear and nonlinear, and stationary and
time-dependent PDEs are supported in `hIPPYlib`.
For stationary problems, gradient and Hessian information can be
automatically generated by `hIPPYlib` using `FEniCS` symbolic differentiation
of the relevant weak forms. For time-dependent problems, instead, symbolic
differentiation can only be used for the spatial terms, and the contribution
to gradients and Hessians arising from the time dynamics needs to be provided
by the user.
Noise and prior covariance operators are modeled as inverses of
elliptic differential operators allowing us to build on existing fast
multigrid solvers for elliptic operators without explicitly
constructing the dense covariance operator.
The key property of the algorithms underlying `hIPPYlib` is that solution
of the deterministic and Bayesian inverse problem is computed
at a cost, measured in forward PDE solves, that is independent of the
parameter dimension.
`hIPPYlib` provides a robust implementation of the inexact
Newton-conjugate gradient algorithm to compute the maximum a posterior
(MAP) point. The gradient and Hessian actions are
computed via their weak form specification in `FEniCS` by
constraining the state and adjoint variables to satisfy the forward
and adjoint problem. The Newton system is solved inexactly by early
termination of CG iterations via Eisenstat-Walker (to prevent
oversolving) and Steihaug (to avoid negative curvature)
criteria. Two globalization techniques are available to the user:
Armijo back-tracking line search and trust region.
In `hIPPYlib`, the posterior covariance is approximated by the
inverse of the Hessian of the negative log posterior evaluated at
the MAP point. This Gaussian approximation is exact when the
parameter-to-observable map is linear; otherwise, its logarithm agrees
to two derivatives with the log posterior at the MAP point, and thus it
can serve as a proposal for Hessian-based Markov chain Monte Carlo (MCMC)
methods. `hIPPYlib` makes the construction of the posterior covariance
tractable by invoking a low-rank approximation of the Hessian of the
log likelihood.
`hIPPYlib` also offers scalable methods for sample generation.
To sample large scale spatially correlated Gaussian random fields from the prior
distribution, `hIPPYlib` implements a new method that strongly relies on the
structure of the covariance operator defined as the inverse of a differential operator:
by exploiting the assembly procedure of finite element matrices `hIPPYlib` constructs a sparse Cholesky-like rectangular decomposition of the precision operator.
To sample from a local Gaussian approximation to the posterior (such as at the MAP point)
`hIPPYlib` exploits the low rank factorization of the Hessian of the
log likelihood to correct samples from the prior distribution.
Finally, to explore the posterior distribution, `hIPPYlib` implements
dimension independent MCMC sampling methods enchanted by Hessian information.
Finally, randomized and probing algorithms are available to compute
the pointwise variance of the prior/posterior distribution and the trace
of the covariance operator.
Raw data
{
"_id": null,
"home_page": "https://hippylib.github.io/",
"name": "hippylib",
"maintainer": "",
"docs_url": null,
"requires_python": ">=3.6,<4",
"maintainer_email": "",
"keywords": "Infinite-dimensional inverse problems,adjoint-based methods,numerical optimization,low-rank approximation,Bayesian inference,uncertainty quantification,sampling",
"author": "Umberto Villa, Noemi Petra, Omar Ghattas",
"author_email": "uvilla@wustl.edu",
"download_url": "",
"platform": null,
"description": "[](https://github.com/hippylib/hippylib/actions/workflows/ci.yaml)\n[](https://hippylib.readthedocs.io/en/latest/)\n[](http://joss.theoj.org/papers/053e0d08a5e9755e7b78898cff6f6208) [](https://doi.org/10.5281/zenodo.596931)\n\n Inverse Problem PYthon library\n\n```\n __ ______ _______ _______ __ __ __ __ __ \n/ | / |/ \\ / \\ / \\ / |/ |/ |/ | \n$$ |____ $$$$$$/ $$$$$$$ |$$$$$$$ |$$ \\ /$$/ $$ |$$/ $$ |____ \n$$ \\ $$ | $$ |__$$ |$$ |__$$ | $$ \\/$$/ $$ |/ |$$ \\ \n$$$$$$$ | $$ | $$ $$/ $$ $$/ $$ $$/ $$ |$$ |$$$$$$$ |\n$$ | $$ | $$ | $$$$$$$/ $$$$$$$/ $$$$/ $$ |$$ |$$ | $$ |\n$$ | $$ | _$$ |_ $$ | $$ | $$ | $$ |$$ |$$ |__$$ |\n$$ | $$ |/ $$ |$$ | $$ | $$ | $$ |$$ |$$ $$/ \n$$/ $$/ $$$$$$/ $$/ $$/ $$/ $$/ $$/ $$$$$$$/ \n``` \n\n\n https://hippylib.github.io\n\n`hIPPYlib` implements state-of-the-art scalable algorithms for\ndeterministic and Bayesian inverse problems governed by partial differential equations (PDEs).\nIt builds on [FEniCS](http://fenicsproject.org/)\n(a parallel finite element element library) for the discretization of the PDE\nand on [PETSc](http://www.mcs.anl.gov/petsc/) for scalable and efficient linear\nalgebra operations and solvers.\n\nFor building instructions, see the file `INSTALL.md`. Copyright information\nand licensing restrictions can be found in the file `COPYRIGHT`.\n\nThe best starting point for new users interested in `hIPPYlib`'s\nfeatures are the interactive tutorials in the `tutorial folder`.\n\nConceptually, `hIPPYlib` can be viewed as a toolbox that provides the\nbuilding blocks for experimenting new ideas and developing scalable\nalgorithms for PDE-constrained deterministic and Bayesian inverse problems.\n\nIn `hIPPYlib` the user can express the forward PDE and the likelihood in\nweak form using the friendly, compact, near-mathematical notation of\n`FEniCS`, which will then automatically generate efficient code for the\ndiscretization. Linear and nonlinear, and stationary and\ntime-dependent PDEs are supported in `hIPPYlib`.\nFor stationary problems, gradient and Hessian information can be\nautomatically generated by `hIPPYlib` using `FEniCS` symbolic differentiation\nof the relevant weak forms. For time-dependent problems, instead, symbolic\ndifferentiation can only be used for the spatial terms, and the contribution\nto gradients and Hessians arising from the time dynamics needs to be provided\nby the user.\n\nNoise and prior covariance operators are modeled as inverses of\nelliptic differential operators allowing us to build on existing fast\nmultigrid solvers for elliptic operators without explicitly\nconstructing the dense covariance operator.\n\nThe key property of the algorithms underlying `hIPPYlib` is that solution\nof the deterministic and Bayesian inverse problem is computed\nat a cost, measured in forward PDE solves, that is independent of the\nparameter dimension.\n\n`hIPPYlib` provides a robust implementation of the inexact\nNewton-conjugate gradient algorithm to compute the maximum a posterior\n(MAP) point. The gradient and Hessian actions are\ncomputed via their weak form specification in `FEniCS` by\nconstraining the state and adjoint variables to satisfy the forward\nand adjoint problem. The Newton system is solved inexactly by early\ntermination of CG iterations via Eisenstat-Walker (to prevent\noversolving) and Steihaug (to avoid negative curvature)\ncriteria. Two globalization techniques are available to the user:\nArmijo back-tracking line search and trust region.\n\nIn `hIPPYlib`, the posterior covariance is approximated by the\ninverse of the Hessian of the negative log posterior evaluated at\nthe MAP point. This Gaussian approximation is exact when the\nparameter-to-observable map is linear; otherwise, its logarithm agrees\nto two derivatives with the log posterior at the MAP point, and thus it\ncan serve as a proposal for Hessian-based Markov chain Monte Carlo (MCMC)\nmethods. `hIPPYlib` makes the construction of the posterior covariance\ntractable by invoking a low-rank approximation of the Hessian of the\nlog likelihood.\n\n`hIPPYlib` also offers scalable methods for sample generation.\nTo sample large scale spatially correlated Gaussian random fields from the prior\ndistribution, `hIPPYlib` implements a new method that strongly relies on the\nstructure of the covariance operator defined as the inverse of a differential operator:\nby exploiting the assembly procedure of finite element matrices `hIPPYlib` constructs a sparse Cholesky-like rectangular decomposition of the precision operator.\nTo sample from a local Gaussian approximation to the posterior (such as at the MAP point)\n`hIPPYlib` exploits the low rank factorization of the Hessian of the\nlog likelihood to correct samples from the prior distribution.\nFinally, to explore the posterior distribution, `hIPPYlib` implements\ndimension independent MCMC sampling methods enchanted by Hessian information.\n\nFinally, randomized and probing algorithms are available to compute\nthe pointwise variance of the prior/posterior distribution and the trace\nof the covariance operator.\n\n\n",
"bugtrack_url": null,
"license": "",
"summary": "an Extensible Software Framework for Large-scale Deterministic and Bayesian Inverse Problems",
"version": "3.1.0",
"split_keywords": [
"infinite-dimensional inverse problems",
"adjoint-based methods",
"numerical optimization",
"low-rank approximation",
"bayesian inference",
"uncertainty quantification",
"sampling"
],
"urls": [
{
"comment_text": "",
"digests": {
"md5": "a70ef0a75ddee1c62416129e46108edb",
"sha256": "1c0fb4894683618b7de12e461b4e186073c0defcb6e7b72834eaa2637e5d65e6"
},
"downloads": -1,
"filename": "hippylib-3.1.0-py3-none-any.whl",
"has_sig": false,
"md5_digest": "a70ef0a75ddee1c62416129e46108edb",
"packagetype": "bdist_wheel",
"python_version": "py3",
"requires_python": ">=3.6,<4",
"size": 122725,
"upload_time": "2022-12-28T18:48:39",
"upload_time_iso_8601": "2022-12-28T18:48:39.696783Z",
"url": "https://files.pythonhosted.org/packages/88/83/d62c2a984c5e30e42e74b0d0bfbce55e363060ddca89510284fe4b0c0b71/hippylib-3.1.0-py3-none-any.whl",
"yanked": false,
"yanked_reason": null
}
],
"upload_time": "2022-12-28 18:48:39",
"github": false,
"gitlab": false,
"bitbucket": false,
"lcname": "hippylib"
}
JSON