pythia-uq 3.0.0

Toolbox for non-intrusive functional approximation of data via (sparse) general polynomial chaos.

PyThia Uncertainty Quantification Toolbox

The PyThia UQ toolbox uses polynomial chaos surrogates to efficiently generate a surrogate of any (parametric) forward problem. The surrogate is fast to evaluate, allows analytical differentiation and has a built-in global sensitivity analysis via Sobol indices. Assembling the surrogate is done non-intrusive by least-squares regression, hence only training pairs of parameter realizations and evaluations of the forward problem are required to construct the surrogate. No need to compute any nasty interfaces for lagacy code.

Why the Name?

Pythia was the title of the high priestess of the temple of Apollo in Delphi. Hence you could say she used her prophetic abilities to quantify which was uncertain. Moreover, the package is written in python, so...

Installation

The latest stable version of PyThia can be installed using pip

pip install pythia-uq

To install PyThia from source, i.e., if you want to work with the latest (and possibly unstable) changes, simply clone the repository and run the setup script to install PyThia to any environment

cd path/to/pythia-repo/
pip install .

PyThia can then be imported from any location with import pythia.

Documentation

The documentation of the stable package version can be found here: pythia-uq.rtfd.io. For other versions of the documentation, see https://readthedocs.org/projects/pythia-uq/.

How to cite PyThia

There is no official related article to cite PyThia yet. If you make use of PyThia in a publication, please cite it with a BibTeX entry similar to this:

@misc{pythia,
    author = {Farchmin, Nando},
    title = {PyThia Uncertainty Quantification Toolbox},
    howpublished = {vers.: 3.0},
    note = {\url{https://gitlab1.ptb.de/pythia/pythia/}},
    year = {2022},
    month = {12}
}

Want to contribute?

Check out the contribution guidelines on how to create issues or file bug reports and feature requests. Or ever better start developping the PyThia project yourself after reading the development guidelines.

Roadmap and TODOs

After making the project public, there are a few necessary user experience changes left that we should make available:

• create stack overflow hashtag for PyThia

Finally, here is a roadmap of features that we plan to add to pythia in the future.

• integrate tensor train representations of coefficients
• add tensor train regression (VMC)[^vmc-base]
• add exponentiation of tensor trains (expTT)[^expTT]
• add efficient posterior rejection sampling for tensor trains posteriors[^TT-posterior-sampling]

References

Here we list the papers that describe concepts implemented in PyThia for the interested user. In principle PyThia uses a (sparse) polynomial chaos expansion to construct a surrogate of any function via least-squares regression. We first applied the PyThia software package to analyse the sensitivities of a scatterometry experiment [^pythia-scat-A] using global sensitivity analysis via Sobol indices [^sobol-indices]. We also solved the inverse problem for the same experiment [^pythia-scat-B] via Bayesian inversion. To use a minimal but still sufficient amount of random samples for the regression, we integrated weighted least-squares sampling [^wls-sampling] into PyThia.

[^pythia-scat-A]: An efficient approach to global sensitivity analysis and parameter estimation for line gratings [^pythia-scat-B]: Efficient Bayesian inversion for shape reconstruction of lithography masks [^sobol-indices]: Global sensitivity analysis using polynomial chaos expansions [^wls-sampling]: Optimal weighted least-squares methods [^vmc-base]: Variational Monte Carlo - Bridging Concepts of Machine Learning and High Dimensional Partial Differential Equations [^expTT]: Efficient approximation of high-dimensional exponentials by tensornetworks [^TT-posterior-sampling]: Computation of the Response Surface in the Tensor Train data format